Source Code

These source-code files derive from a printout of Luminary 210 (Apollo 15-17 Lunar Module guidance computer program), from the personal library of original AGC developer Don Eyles, digitally photographed at archive.org, financially sponsored by Jim Lawton, and transcribed to source code by a team of volunteers. This colorized, syntax-highlighted form was created by assembling that transcribed source code. Note that the full page images are available on the Virtual AGC project page at archive.org, while reduced-size images are presented at the VirtualAGC project website. Report or fix any transcription errors at the Virtual AGC project code repository.

Notations on the program listing read, in part:

	GAP: ASSEMBLE REVISION 210 OF AGC PROGRAM LUMINARY 
	BY NASA 2021112-161 17:11 MAR. 19,1971
Note that the date is the date of the printout, not the date of the program revision.

043146,000002:                                                                                                  ## Copyright:   Public domain.
043147,000003:                                                                                                  ## Filename:    LUNAR_AND_SOLAR_EPHEMERIDES_SUBROUTINES.agc
043148,000004:                                                                                                  ## Purpose:     A section of Luminary revision 210.
043149,000005:                                                                                                  ##              It is part of the source code for the Lunar Module's (LM)
043150,000006:                                                                                                  ##              Apollo Guidance Computer (AGC) for Apollo 15-17.
043151,000007:                                                                                                  ##              This file is intended to be a faithful transcription, except
043152,000008:                                                                                                  ##              that the code format has been changed to conform to the
043153,000009:                                                                                                  ##              requirements of the yaYUL assembler rather than the
043154,000010:                                                                                                  ##              original YUL assembler.
043155,000011:                                                                                                  ## Reference:   pp. 982-985
043156,000012:                                                                                                  ## Assembler:   yaYUL
043157,000013:                                                                                                  ## Contact:     Ron Burkey <info@sandroid.org>.
043158,000014:                                                                                                  ## Website:     www.ibiblio.org/apollo/index.html
043159,000015:                                                                                                  ## Mod history: 2016-11-17 JL   Created from Luminary131 version.
043160,000016:                                                                                                  ##              2016-11-29 HG   Transcribed
043161,000017:                                                                                                  ##              2016-12-12 HG   Fix for interpretive operands 8D -> 9D
043162,000018:                                                                                                  ##                                                            5D -> 6D 
043163,000019:                                                                                                  ##              2016-12-26 RSB  Comment-text proofed using ProoferComments
043164,000020:                                                                                                  ##                              and corrected errors found.
043165,000021:                                                                                                  ##              2017-03-13 RSB  Comment-text fixes noted in proofing Luminary 116.
043166,000022: 

Page 982

043168,000024:                                                                                                  #  NAME - LSPOS  - LOCATE SUN AND MOON                   DATE - 25 OCT 67
043169,000025:                                                                                                  #  MOD NO.2  (DATE OCT 70)
043170,000026:                                                                                                  #  MOD BY NEVILLE                                        ASSEMBLY SUNDANCE
043171,000027: 
043172,000028:                                                                                                  #  FUNCTIONAL DESCRIPTION
043173,000029: 
043174,000030:                                                                                                  #         COMPUTES UNIT POSITION VECTOR OF THE SUN AND MOON IN THE BASIC REFERENCE SYSTEM. THE SUN VECTOR S IS
043175,000031:                                                                                                  #  LOCATED VIA TWO ANGLES. THE FIRST ANGLE(OBLIQUITY) IS THE ANGLE BETWEEN THE EARTH EQUATOR AND THE ECLIPTIC.  THE
043176,000032:                                                                                                  #  SECOND ANGLE IS THE LONGITUDE OF THE SUN MEASURED IN THE ECLIPTIC.
043177,000033:                                                                                                  #  THE POSITION VECTOR OF THE SUN IS
043178,000034:                                                                                                  #         -
043179,000035:                                                                                                  #         S=(COS(LOS), COS(OBL)*SIN(LOS), SIN(OBL)*SIN(LOS)), WHERE
043180,000036: 
043181,000037:                                                                                                  #       LOS=LOS +LOS *T-(C *SIN(2PI*T)/365.24 +C *COS(2PI*T)/365.24)
043182,000038:                                                                                                  #              0    R     0                     1
043183,000039:                                                                                                  #       LOS  (RAD) IS THE LONGITUDE OF THE SUN FOR MIDNIGHT JUNE 30TH OF THE PARTICULAR YEAR.
043184,000040:                                                                                                  #          0
043185,000041:                                                                                                  #       LOS  (RAD/DAY) IS THE MEAN RATE FOR THE PARTICULAR YEAR.
043186,000042:                                                                                                  #          R
043187,000043:                                                                                                  #  LOS  AND LOS  ARE STORED AS LOSC AND LOSR IN RATESP.
043188,000044:                                                                                                  #     0        R
043189,000045:                                                                                                  #  COS(OBL) AND SIN(OBL) ARE STORED IN THE MATRIX KONMAT.
043190,000046:                                                                                                  #  T, TIME MEASURED IN DAYS (24 HOURS), IS STORED IN TIMEP.
043191,000047:                                                                                                  #  C  AND C  ARE FUDGE FACTORS TO MINIMIZE THE DEVIATION.  THEY ARE STORED AS ONE CONSTANT(CMOD), SINCE
043192,000048:                                                                                                  #   0      1                               2  2 1/2
043193,000049:                                                                                                  #  C *SIN(X)+C *COS(X) CAN BE WRITTEN AS (C +C )   *SIN(X+PHI), WHERE PHI=ARCTAN(C /C ).
043194,000050:                                                                                                  #   0         1                            0  1                                   1  0
043195,000051: 
043196,000052:                                                                                                  #  THE MOON IS LOCATED VIA FOUR ANGLES, THE FIRST IS THE OBLIQUITY. THE SECOND IS THE MEAN LONGITUDE OF THE MOON,
043197,000053:                                                                                                  #  MEASURED IN THE ECLIPTIC FROM THE MEAN EQUINOX TO THE MEAN ASCENDING NODE OF THE LUNAR ORBIT, AND THEN ALONG THE
043198,000054:                                                                                                  #  ORBIT. THE THIRD ANGLE IS THE ANGLE BETWEEN THE ECLIPTIC AND THE LUNAR ORBIT. THE FOURTH ANGLE IS THE LONGITUDE
043199,000055:                                                                                                  #  OF THE NODE OF THE MOON, MEASURED IN THE LUNAR ORBIT. LET THESE ANGLES BE OBL,LOM,IM, AND LON RESPECTIVELY.
043200,000056: 
043201,000057:                                                                                                  #  THE SIMPLIFIED POSITION VECTOR OF THE MOON IS
043202,000058:                                                                                                  #    -
043203,000059:                                                                                                  #    M=(COS(LOM), COS(OBL)*SIN(LOM)-SIN(OBL)*SIN(IM)*SIN(LOM-LON), SIN(OBL)*SIN(LOM)+COS(OBL)*SIN(IM)*SIN(LOM-LON))
043204,000060: 
043205,000061:                                                                                                  #    WHERE
043206,000062:                                                                                                  #       LOM=LOM +LOM *T-(A *SIN 2PI*T/27.5545+A *COS(2PI*T/27.5545)+B *SIN 2PI*T/32+B *COS(2PI*T/32)), AND
043207,000063:                                                                                                  #              0    R     0                    1                     0               1
043208,000064:                                                                                                  #        LON=LON +LON
043209,000065:                                                                                                  #               0    R
043210,000066:                                                                                                  #  A , A , B  AND B  ARE STORED AS AMOD AND BMOD (SEE DESCRIPTION OF CMOD, ABOVE).  COS(OBL), SIN(OBL)*SIN(IM),
043211,000067:                                                                                                  #   0   1   0      1
043212,000068:                                                                                                  #  SIN(OBL), AND COS(OBL)*SIN(IM) ARE STORED IN KONMAT AS K1, K2, K3 AND K4, RESPECTIVELY. LOM , LOM , LON , LON
043213,000069:                                                                                                  #  ARE STORED AS LOM0, LOMR, LON0, AND LONR IN RATESP.                                        0     R     0     R
043214,000070:                                                                                                  #  THE THREE PHIS ARE STORED AS AARG, BARG, AND CARG(SUN).  ALL CONSTANTS ARE UPDATED BY YEAR.
043215,000071: 
043216,000072:                                                                                                  #  CALLING SEQUENCE
043217,000073: 

Page 983

043219,000075:                                                                                                  #    CALL LSPOS.  RETURN IS VIA CPRET.
043220,000076: 
043221,000077:                                                                                                  #  ALARMS OR ABORTS
043222,000078: 
043223,000079:                                                                                                  #    NONE
043224,000080: 
043225,000081:                                                                                                  #  ERASABLE INITIALIZATION REQUIRED
043226,000082: 
043227,000083:                                                                                                  #    TEPHEM - TIME FROM MIDNIGHT 1 JULY PRECEDING THE LAUNCH TO THE TIME OF THE LAUNCH (WHEN THE AGC CLOCK WENT
043228,000084:                                                                                                  #  TO ZERO). TEPHEM IS TP WITH UNITS OF CENTI-SECONDS.
043229,000085:                                                                                                  #    TIME2 AND TIME1 ARE IN MPAC AND MPAC +1 WHEN PROGRAM IS CALLED.
043230,000086: 
043231,000087:                                                                                                  #  OUTPUT
043232,000088: 
043233,000089:                                                                                                  #    UNIT POSITIONAL VECTOR OF SUN IN VSUN.   (SCALED B-1)
043234,000090:                                                                                                  #    UNIT POSITIONAL VECTOR OF MOON IN VMOON. (SCALED B-1)
043235,000091: 
043236,000092:                                                                                                  #  SUBROUTINES USED
043237,000093: 
043238,000094:                                                                                                  #    NONE
043239,000095: 
043240,000096:                                                                                                  #  DEBRIS
043241,000097: 
043242,000098:                                                                                                  #    CURRENT CORE SET, WORK AREA AND FREEFLAG
043243,000099: 04,3030                                           BANK     04                                    
043244,000100: 15,2000                                           SETLOC   EPHEM                                 
043245,000101: 15,2000                                           BANK                                           
043246,000102: 
043247,000103: 15,3646  E5,1714                                  EBANK=   VSUN                                  
043248,000104: 15,3646                                           COUNT*   $$/EPHEM                              
043249,000105: 15,3646  15,3646               LUNPOS             EQUALS   LSPOS                                 
043250,000106: 
043251,000107: 15,3646           54201        LSPOS              SETPD    SR                                    
043252,000108: 15,3647           00001                                    0                                     
043253,000109: 15,3650           20617                                    14D                                   #  TP
043254,000110: 15,3651           56371                           TAD      DDV                                   
043255,000111: 15,3652           01707                                    TEPHEM                                #  TIME OF LAUNCH
043256,000112: 15,3653           32001                                    CSTODAY                               #  24 HOURS-8640000 CENTI-SECS/DAY B-32
043257,000113: 15,3654           00031                           STORE    TIMEP                                 #  T IN DAYS
043258,000114: 15,3655           77170                           AXT,1    AXT,2                                 
043259,000115: 15,3656           00000                                    0                                     
043260,000116: 15,3657           00000                                    0                                     
043261,000117: 15,3660           77614                           CLEAR                                          
043262,000118: 15,3661           00274                                    FREEFLAG                              #  SWITCH BIT
043263,000119: 15,3662           77745        POSITA             DLOAD                                          
043264,000120: 15,3663           12004                                    KONMAT     +2                         #  ZEROS
043265,000121: 15,3664           00027                           STORE    GTMP                                  
043266,000122: 15,3665           40745        POSITB             DLOAD    DMP*                                  
043267,000123: 15,3666           00031                                    TIMEP                                 #  T
043268,000124: 15,3667           12030                                    VAL67      +4,1                       #  1/27 OR 1/32 OR 1/365

Page 984

043270,000126: 15,3670           42661                           SL       DAD*                                  
043271,000127: 15,3671           20212                                    9D                                    
043272,000128: 15,3672           12026                                    VAL67      +2,1                       #  AARG
043273,000129: 15,3673           40756                           SIN      DMP*                                  #  SIN(T/27+PHI) OR T/32 OR T/365
043274,000130: 15,3674           12024                                    VAL67,1                               #  (A0**2+A1**2)**1/2 SIN(X+PHIA)
043275,000131: 15,3675           62015                           DAD      INCR,1                                #  PLUS
043276,000132: 15,3676           00027                                    GTMP                                  #  (B0**2+B1**2)**1/2 SIN(X+PHIB)
043277,000133: 15,3677           77771                           DEC      -6         B-14                       
043278,000134: 15,3700           00027                           STORE    GTMP                                  #  OR (C0**2+C1**2)**1/2 SIN(X+PHIC)
043279,000135: 15,3701           77614                           BOFSET                                         
043280,000136: 15,3702           00054                                    FREEFLAG                              
043281,000137: 15,3703           33665                                    POSITB                                
043282,000138: 15,3704           40745        POSITD             DLOAD    DMP*                                  
043283,000139: 15,3705           00031                                    TIMEP                                 #  T
043284,000140: 15,3706           67752                                    RATESP,2                              #  LOMR,LOSR,LONR
043285,000141: 15,3707           42661                           SL       DAD*                                  
043286,000142: 15,3710           20207                                    6D                                    
043287,000143: 15,3711           67744                                    RATESP     +6,2                       #  LOMO,LOSO,LONO
043288,000144: 15,3712           77625                           DSU                                            
043289,000145: 15,3713           00027                                    GTMP                                  
043290,000146: 15,3714           10021                           STORE    STMP,2                                #  LOM,LOS,LON
043291,000147: 15,3715           63135                           SLOAD    INCR,2                                
043292,000148: 15,3716           00050                                    X2                                    
043293,000149: 15,3717           77775                           DEC      -2         B-14                       
043294,000150: 15,3720           53015                           DAD      BZE                                   
043295,000151: 15,3721           32003                                    RCB-13                                #  PLUS 2
043296,000152: 15,3722           33752                                    POSITE                                #  2ND
043297,000153: 15,3723           77644                           BPL                                            
043298,000154: 15,3724           33662                                    POSITA                                #  1ST
043299,000155: 15,3725           45345        POSITF             DLOAD    DSU                                   #  3RD
043300,000156: 15,3726           00021                                    STMP                                  #  LOM
043301,000157: 15,3727           00025                                    STMP       +4                         #  LON
043302,000158: 15,3730           65356                           SIN      PDDL                                  #  SIN(LOM-LON)
043303,000159: 15,3731           00021                                    STMP                                  
043304,000160: 15,3732           65356                           SIN      PDDL                                  #  SIN LOM
043305,000161: 15,3733           00021                                    STMP                                  
043306,000162: 15,3734           55546                           COS      VDEF                                  #  COS LOM
043307,000163: 15,3735           53521                           MXV      UNIT                                  
043308,000164: 15,3736           12002                                    KONMAT                                #  K1,K2,K3,K4,
043309,000165: 15,3737           02723                           STORE    VMOON                                 
043310,000166: 15,3740           65345                           DLOAD    PDDL                                  
043311,000167: 15,3741           12004                                    KONMAT     +2                         #  ZERO
043312,000168: 15,3742           00023                                    STMP       +2                         
043313,000169: 15,3743           65356                           SIN      PDDL                                  #  SIN LOS
043314,000170: 15,3744           00023                                    STMP       +2                         
043315,000171: 15,3745           55546                           COS      VDEF                                  #  COS LOS
043316,000172: 15,3746           53521                           MXV      UNIT                                  
043317,000173: 15,3747           12002                                    KONMAT                                
043318,000174: 15,3750           02715                           STORE    VSUN                                  
043319,000175: 15,3751           77616                           RVQ                                            
043320,000176: 

Page 985

043322,000178: 15,3752           77745        POSITE             DLOAD                                          
043323,000179: 15,3753           12004                                    KONMAT     +2                         #  ZEROS
043324,000180: 15,3754           00027                           STORE    GTMP                                  
043325,000181: 15,3755           77650                           GOTO                                           
043326,000182: 15,3756           33704                                    POSITD                                
043327,000183: 05,2000                                           SETLOC   EPHEM1                                
043328,000184: 05,2000                                           BANK                                           
043329,000185: 
043330,000186: 05,3507                                           COUNT*   $$/EPHEM                              
043331,000187: 
043332,000188: 05,3507                        STMP               EQUALS   16D                                   
043333,000189: 
043334,000190: 05,3507                        GTMP               EQUALS   22D                                   
043335,000191: 
043336,000192: 05,3507                        TIMEP              EQUALS   24D                                   

End of include-file LUNAR_AND_SOLAR_EPHEMERIDES_SUBROUTINES.agc.  Parent file is MAIN.agc