Source Code

These source-code files were from the printout at the MIT Museum, with many thanks to Deborah Douglas, the Museum's Curator of Science and Technology, for releasing it to us. The printout was digitally photographed by Paul Fjeld. A team of volunteers transcribed the source code manually from these images or, where appropriate, modified already-transcribed but similar source Luminary 131 (Apollo 13) source code to agree with the images from the printout. Note that the page images presented online are of reduced quality, and that higher-quality images are available. Report any conversion errors or legibility problems in page images to info@sandroid.org. Notations on the program listing read, in part:

	GAP:  ASSEMBLE REVISION 001 OF AGC PROGRAM LMY99 BY NASA 2021112-061
	20'35 OCT. 28,1968
Note that the date is the date of the printout, not the date of the program revision.

048293,000002:                                                                                                  ## Copyright:   Public domain.
048294,000003:                                                                                                  ## Filename:    LATITUDE_LONGITUDE_SUBROUTINES.agc
048295,000004:                                                                                                  ## Purpose:     Part of the source code for Luminary 1A build 099.
048296,000005:                                                                                                  ##              It is part of the source code for the Lunar Module's (LM)
048297,000006:                                                                                                  ##              Apollo Guidance Computer (AGC), for Apollo 11.
048298,000007:                                                                                                  ##
048299,000008:                                                                                                  ## Assembler:   yaYUL
048300,000009:                                                                                                  ## Contact:     Jim Lawton <jim DOT lawton AT gmail DOT com>
048301,000010:                                                                                                  ## Website:     www.ibiblio.org/apollo.
048302,000011:                                                                                                  ## Pages:       1133-1139
048303,000012:                                                                                                  ## Mod history: 2009-05-28 JL   Updated from page images.
048304,000013:                                                                                                  ##              2011-01-06 JL   Fixed interpretive indentation.
048305,000014:                                                                                                  ##              2016-12-17 RSB  Proofed text comments with octopus/ProoferComments
048306,000015:                                                                                                  ##                              and corrected the errors found.
048307,000016:                                                                                                  ##              2017-03-15 RSB  Comment-text fixes identified in 5-way
048308,000017:                                                                                                  ##                              side-by-side diff of Luminary 69/99/116/131/210.
048309,000018: 
048310,000019:                                                                                                  ## This source code has been transcribed or otherwise adapted from digitized
048311,000020:                                                                                                  ## images of a hardcopy from the MIT Museum.  The digitization was performed
048312,000021:                                                                                                  ## by Paul Fjeld, and arranged for by Deborah Douglas of the Museum.  Many
048313,000022:                                                                                                  ## thanks to both.  The images (with suitable reduction in storage size and
048314,000023:                                                                                                  ## consequent reduction in image quality as well) are available online at
048315,000024:                                                                                                  ## www.ibiblio.org/apollo.  If for some reason you find that the images are
048316,000025:                                                                                                  ## illegible, contact me at info@sandroid.org about getting access to the
048317,000026:                                                                                                  ## (much) higher-quality images which Paul actually created.
048318,000027:                                                                                                  ##
048319,000028:                                                                                                  ## Notations on the hardcopy document read, in part:
048320,000029:                                                                                                  ##
048321,000030:                                                                                                  ##    Assemble revision 001 of AGC program LMY99 by NASA 2021112-061
048322,000031:                                                                                                  ##    16:27 JULY 14, 1969
048323,000032: 

Page 1133

048325,000034:                                                                                                  #  SUBROUTINE TO CONVERT RAD VECTOR AT GIVEN TIME TO LAT,LONG AND ALT
048326,000035: 
048327,000036:                                                                                                  #  CALLING SEQUENCE
048328,000037:                                                                                                  #        L-1     CALL
048329,000038:                                                                                                  #        L               LAT-LONG
048330,000039: 
048331,000040:                                                                                                  #  SUBROUTINES USED
048332,000041:                                                                                                  #        R-TO-RP, ARCTAN, SETGAMMA, SETRE
048333,000042: 
048334,000043:                                                                                                  #  ERASABLE INIT. REQ.
048335,000044:                                                                                                  #        AXO, -AYO, AZO, TEPHEM (SET AT LAUNCH TIME)
048336,000045:                                                                                                  #        ALPHAV = POSITION VECTOR METERS B-29
048337,000046:                                                                                                  #        MPAC -- TIME (CSECS B-28)
048338,000047:                                                                                                  #        ERADFLAG =1, TO COMPUTE EARTH RADIUS, =0 FOR FIXED EARTH RADIUS
048339,000048:                                                                                                  #        LUNAFLAG=0 FOR EARTH, 1 FOR MOON
048340,000049: 
048341,000050:                                                                                                  #  OUTPUT
048342,000051:                                                                                                  #        LATITUDE IN LAT         (REVS. B-0)
048343,000052:                                                                                                  #        LONGITUDE IN LONG       (REVS. B-0)
048344,000053:                                                                                                  #        ALTITUDE IN ALT         METERS B-29
048345,000054: 
048346,000055: 30,3775                                           BANK     30                                    
048347,000056: 13,2000                                           SETLOC   LATLONG                               
048348,000057: 13,2000                                           BANK                                           
048349,000058: 
048350,000059: 13,2351                                           COUNT*   $$/LT-LG                              
048351,000060: 13,2351  E4,1431                                  EBANK=   ALPHAV                                
048352,000061: 13,2351           40220        LAT-LONG           STQ      SETPD                                 
048353,000062: 13,2352           03674                                    INCORPEX                              
048354,000063: 13,2353           00001                                    0D                                    
048355,000064: 13,2354           24007                           STOVL    6D                                    #  SAVE TIME IN 6-7D FOR R-TO-RP
048356,000065: 13,2355           02032                                    ALPHAV                                
048357,000066: 13,2356           51406                           PUSH     ABVAL                                 #  0-5D= R FOR R-TO-RP
048358,000067: 13,2357           16070                           STODL    ALPHAM                                #  ABS. VALUE OF R FOR ALT FORMULA BELOW
048359,000068: 13,2360           24007                                    ZEROVEC                               #  SET MPAC=0 FOR EARTH, NON-ZERO FOR MOON
048360,000069: 13,2361           71414                           BOFF     COS                                   #  USE COS(0) TO GET NON-ZERO IN MPAC
048361,000070: 13,2362           01743                                    LUNAFLAG                              #  0=EARTH, 1=MOON
048362,000071: 13,2363           26364                                    CALLRTRP                              
048363,000072: 13,2364           77624        CALLRTRP           CALL                                           
048364,000073: 13,2365           51670                                    R-TO-RP                               #  RP VECTOR CONVERTED FROM R B-29
048365,000074: 13,2366           77656                           UNIT                                           #  UNIT RP B-1
048366,000075: 13,2367           36032                           STCALL   ALPHAV                                #  U2= 1/2 SINL FOR SETRE SUBR BELOW
048367,000076: 13,2370           26550                                    SETGAMMA                              #        SET GAMMA=B2/A2 FOR EARTH, =1 FOR MOON
048368,000077: 13,2371           77624                           CALL                                           #        SCALED B-1
048369,000078: 13,2372           26560                                    SETRE                                 #  CALC RE METERS B-29
048370,000079: 13,2373           63545                           DLOAD    DSQ                                   
048371,000080: 13,2374           02032                                    ALPHAV                                
048372,000081: 13,2375           63525                           PDDL     DSQ                                   
048373,000082: 13,2376           02034                                    ALPHAV     +2                         
048374,000083: 13,2377           75415                           DAD      SQRT                                  

Page 1134

048376,000085: 13,2400           76405                           DMP      SL1R                                  
048377,000086: 13,2401           00011                                    GAMRP                                 
048378,000087: 13,2402           14021                           STODL    COSTH                                 #  COS(LAT) B-1
048379,000088: 13,2403           02036                                    ALPHAV     +4                         
048380,000089: 13,2404           34023                           STCALL   SINTH                                 #  SIN(LAT) B-1
048381,000090: 13,2405           26510                                    ARCTAN                                
048382,000091: 13,2406           15121                           STODL    LAT                                   #  LAT B0
048383,000092: 13,2407           02032                                    ALPHAV                                
048384,000093: 13,2410           14021                           STODL    COSTH                                 #  COS(LONG) B-1
048385,000094: 13,2411           02034                                    ALPHAV     +2                         
048386,000095: 13,2412           34023                           STCALL   SINTH                                 #  SIN(LONG) B-1
048387,000096: 13,2413           26510                                    ARCTAN                                
048388,000097: 13,2414           15123                           STODL    LONG                                  #  LONG. REVS B-0 IN RANGE -1/2 TO 1/2
048389,000098: 13,2415           02070                                    ALPHAM                                
048390,000099: 13,2416           77625                           DSU                                            #  ALT= R-RE METERS B-29
048391,000100: 13,2417           03673                                    ERADM                                 
048392,000101: 13,2420           35125                           STCALL   ALT                                   #  EXIT WITH ALT METERS B-29
048393,000102: 13,2421           03674                                    INCORPEX                              

Page 1135

048395,000104:                                                                                                  #  SUBROUTINE TO CONVERT LAT,LONG,ALT AT GIVEN TIME TO RADIUS VECTOR
048396,000105: 
048397,000106:                                                                                                  #  CALLING SEQUENCE
048398,000107:                                                                                                  #        L-1     CALL
048399,000108:                                                                                                  #        L               LALOTORV
048400,000109: 
048401,000110:                                                                                                  #  SUBROUTINES USED
048402,000111:                                                                                                  #        SETGAMMA, SETRE, RP-TO-R
048403,000112: 
048404,000113:                                                                                                  #  ERASABLE INIT. REQ.
048405,000114:                                                                                                  #        AXO, AYO, AZO, TEPHEM SET AT LAUNCH TIME
048406,000115:                                                                                                  #        LAT -- LATITUDE         (REVS B0)
048407,000116:                                                                                                  #        LONG -- LONGITUDE       (REVS B0)
048408,000117:                                                                                                  #        ALT -- ALTITUDE         (METERS) B-29
048409,000118:                                                                                                  #        MPAC -- TIME            (CSECS B-28)
048410,000119:                                                                                                  #        ERADFLAG =1 TO COMPUTE EARTH RADIUS, =0 FOR FIXED EARTH RADIUS
048411,000120:                                                                                                  #        LUNAFLAG=0 FOR EARTH, 1 FOR MOON
048412,000121: 
048413,000122:                                                                                                  #  OUTPUT
048414,000123:                                                                                                  #        R-VECTOR IN ALPHAV      (METERS B-29)
048415,000124: 
048416,000125: 13,2422           40220        LALOTORV           STQ      SETPD                                 #  LAT,LONG,ALT TO R VECTOR
048417,000126: 13,2423           03674                                    INCORPEX                              
048418,000127: 13,2424           00001                                    0D                                    
048419,000128: 13,2425           34007                           STCALL   6D                                    #  6-7D= TIME FOR RP-TO-R
048420,000129: 13,2426           26550                                    SETGAMMA                              #  GAMMA=B2/A2 FOR EARTH, 1 FOR MOON B-1
048421,000130: 13,2427           73545                           DLOAD    SIN                                   #                COS(LONG)COS(LAT) IN MPAC
048422,000131: 13,2430           01121                                    LAT                                   #      UNIT RP = SIN(LONG)COS(LAT)    2-3D
048423,000132: 13,2431           65275                           DMPR     PDDL                                  #  PD 2          GAMMA*SIN(LAT)       0-1D
048424,000133: 13,2432           00011                                    GAMRP                                 
048425,000134: 13,2433           01121                                    LAT                                   #         0-1D = GAMMA*SIN(LAT) B-2
048426,000135: 13,2434           65346                           COS      PDDL                                  #  PD4    2-3D = COS(LAT) B-1 TEMPORARILY
048427,000136: 13,2435           01123                                    LONG                                  
048428,000137: 13,2436           57356                           SIN      DMPR                                  #  PD 2
048429,000138: 13,2437           71525                           PDDL     COS                                   #  PD 4   2-3D = SIN(LONG)COS(LAT) B-2
048430,000139: 13,2440           01121                                    LAT                                   
048431,000140: 13,2441           71525                           PDDL     COS                                   #  PD 6   4-5D = COS(LAT) B-1 TEMPORARILY
048432,000141: 13,2442           01123                                    LONG                                  
048433,000142: 13,2443           55475                           DMPR     VDEF                                  #  PD 4   MPAC = COS(LONG)COS(LAT) B-2
048434,000143: 13,2444           41456                           UNIT     PUSH                                  #  0-5D= UNIT RP FOR RP-TO-R SUBR.
048435,000144: 13,2445           36032                           STCALL   ALPHAV                                #  ALPHAV +4= SINL FOR SETRE SUBR.
048436,000145: 13,2446           26560                                    SETRE                                 #  RE METERS B-29
048437,000146: 13,2447           43145                           DLOAD    BOFF                                  #  SET MPAC=0 FOR EARTH, NON-ZERO FOR MOON
048438,000147: 13,2450           24007                                    ZEROVEC                               
048439,000148: 13,2451           01743                                    LUNAFLAG                              
048440,000149: 13,2452           26454                                    CALLRPRT                              
048441,000150: 13,2453           77746                           COS                                            #  USE COS(0) TO GET NON-ZERO IN MPAC
048442,000151: 13,2454           77624        CALLRPRT           CALL                                           
048443,000152: 13,2455           55716                                    RP-TO-R                               #  EXIT WITH UNIT R VECTOR IN MPAC
048444,000153: 13,2456           16032                           STODL    ALPHAV                                
048445,000154: 13,2457           03673                                    ERADM                                 

Page 1136

048447,000156: 13,2460           74215                           DAD      VXSC                                  #  (RE + ALT)(UNIT R) METERS B-30
048448,000157: 13,2461           01125                                    ALT                                   
048449,000158: 13,2462           02032                                    ALPHAV                                
048450,000159: 13,2463           77772                           VSL1                                           #  R METERS B-29
048451,000160: 13,2464           36032                           STCALL   ALPHAV                                #  EXIT WITH R IN METERS B-29
048452,000161: 13,2465           03674                                    INCORPEX                              
048453,000162: 
048454,000163:                                                                                                  #  SUBROUTINE TO COMPUTE EARTH RADIUS
048455,000164: 
048456,000165:                                                                                                  #  INPUT
048457,000166:                                                                                                  #        1/2 SIN LAT IN ALPHAV +4
048458,000167: 
048459,000168:                                                                                                  #  OUTPUT
048460,000169:                                                                                                  #        EARTH RADIUS IN ERADM AND MPAC (METERS B-29)
048461,000170: 
048462,000171: 13,2466           63545        GETERAD            DLOAD    DSQ                                   
048463,000172: 13,2467           02036                                    ALPHAV     +4                         #  SIN**2(L)
048464,000173: 13,2470           44352                           SL1      BDSU                                  
048465,000174: 13,2471           24005                                    DP1/2                                 #  COS**2(L)
048466,000175: 13,2472           44275                           DMPR     BDSU                                  
048467,000176: 13,2473           26507                                    EE                                    
048468,000177: 13,2474           24005                                    DP1/2                                 
048469,000178: 13,2475           75465                           BDDV     SQRT                                  
048470,000179: 13,2476           26503                                    B2XSC                                 
048471,000180: 13,2477           77622                           SR4R                                           
048472,000181: 13,2500           03673                           STORE    ERADM                                 
048473,000182: 13,2501           77616                           RVQ                                            
048474,000183: 
048475,000184:                                                                                                  #  THE FOLLOWING CONSTANTS WERE COMPUTED WITH A=6378166, B=6356784 METERS
048476,000185:                                                                                                  #  B2XSC = B**2 SCALED B-51
048477,000186:                                                                                                  #  B2/A2 = B**2/A**2 SCALED B-1
048478,000187:                                                                                                  #  EE = (1-B**2/A**2) SCALED B-0
048479,000188: 
048480,000189: 13,2502           00446 00305  B2XSC              2DEC     .0179450689                            #  B**2 SCALED B-51
048481,000190: 13,2504  12,2004               DP1/2              =        XUNIT                                 
048482,000191: 13,2504           17711 05254  B2/A2              2DEC     .9933064884 B-1                        #  GAMMA= B**2/A**2 B-1
048483,000192: 13,2506           00155 25250  EE                 2DEC     6.6935116  E-3                        #  (1-B**2/A**2) B-0
048484,000193: 

Page 1137

048486,000195:                                                                                                  #  ARCTAN SUBROUTINE
048487,000196: 
048488,000197:                                                                                                  #  CALLING SEQUENCE
048489,000198:                                                                                                  #        SIN THETA IN SINTH B-1
048490,000199:                                                                                                  #        COS THETA IN COSTH B-1
048491,000200:                                                                                                  #        CALL ARCTAN
048492,000201: 
048493,000202:                                                                                                  #  OUTPUT
048494,000203:                                                                                                  #        ARCTAN THETA IN MPAC AND THETA B-0 IN RANGE -1/2 TO +1/2
048495,000204: 
048496,000205: 13,2510           77600        ARCTAN             BOV                                            
048497,000206: 13,2511           26512                                    CLROVFLW                              
048498,000207: 13,2512           63545        CLROVFLW           DLOAD    DSQ                                   
048499,000208: 13,2513           00023                                    SINTH                                 
048500,000209: 13,2514           63525                           PDDL     DSQ                                   
048501,000210: 13,2515           00021                                    COSTH                                 
048502,000211: 13,2516           77615                           DAD                                            
048503,000212: 13,2517           75454                           BZE      SQRT                                  
048504,000213: 13,2520           26536                                    ARCTANXX                              #  ATAN=0/0  SET THETA=0
048505,000214: 13,2521           40065                           BDDV     BOV                                   
048506,000215: 13,2522           00023                                    SINTH                                 
048507,000216: 13,2523           26543                                    ATAN=90                               
048508,000217: 13,2524           67542                           SR1      ASIN                                  
048509,000218: 13,2525           00025                           STORE    THETA                                 
048510,000219: 13,2526           50125                           PDDL     BMN                                   
048511,000220: 13,2527           00021                                    COSTH                                 
048512,000221: 13,2530           26532                                    NEGCOS                                
048513,000222: 13,2531           43545                           DLOAD    RVQ                                   
048514,000223: 13,2532           57545        NEGCOS             DLOAD    DCOMP                                 
048515,000224: 13,2533           43244                           BPL      DAD                                   
048516,000225: 13,2534           26540                                    NEGOUT                                
048517,000226: 13,2535           24005                                    DP1/2                                 
048518,000227: 13,2536           00025        ARCTANXX           STORE    THETA                                 
048519,000228: 13,2537           77616                           RVQ                                            
048520,000229: 
048521,000230: 13,2540           52025        NEGOUT             DSU      GOTO                                  
048522,000231: 13,2541           24005                                    DP1/2                                 
048523,000232: 13,2542           26536                                    ARCTANXX                              
048524,000233: 13,2543           75345        ATAN=90            DLOAD    SIGN                                  
048525,000234: 13,2544           11037                                    LODP1/4                               
048526,000235: 13,2545           00023                                    SINTH                                 
048527,000236: 13,2546           00025                           STORE    THETA                                 
048528,000237: 13,2547           77616                           RVQ                                            
048529,000238: 
048530,000239: 13,2550  12,2006               2DZERO             =        DPZERO                                
048531,000240: 

Page 1138

048533,000242:                                                                                                  #  ..... SETGAMMA SUBROUTINE .....
048534,000243:                                                                                                  #  SUBROUTINE TO SET GAMMA FOR THE LAT-LONG AND LALOTORV SUBROUTINES
048535,000244: 
048536,000245:                                                                                                  #  GAMMA = B**2/A**2 FOR EARTH (B-1)
048537,000246:                                                                                                  #  GAMMA = 1 FOR MOON (B-1)
048538,000247: 
048539,000248:                                                                                                  #  CALLING SEQUENCE
048540,000249:                                                                                                  #        L       CALL
048541,000250:                                                                                                  #        L+1             SETGAMMA
048542,000251: 
048543,000252:                                                                                                  #  INPUT
048544,000253:                                                                                                  #        LUNAFLAG=0 FOR EARTH, =1 FOR MOON
048545,000254: 
048546,000255:                                                                                                  #  OUTPUT
048547,000256:                                                                                                  #        GAMMA IN GAMRP (B-1)
048548,000257: 
048549,000258: 13,2550           43145        SETGAMMA           DLOAD    BOFF                                  #  BRANCH FOR EARTH
048550,000259: 13,2551           26505                                    B2/A2                                 #  EARTH GAMMA
048551,000260: 13,2552           01743                                    LUNAFLAG                              
048552,000261: 13,2553           26556                                    SETGMEX                               
048553,000262: 13,2554           77735                           SLOAD                                          
048554,000263: 13,2555           24005                                    1B1                                   #  MOON GAMMA
048555,000264: 13,2556           00011        SETGMEX            STORE    GAMRP                                 
048556,000265: 13,2557           77616                           RVQ                                            
048557,000266: 13,2560                        GAMRP              =        8D                                    
048558,000267: 

Page 1139

048560,000269:                                                                                                  #  ..... SETRE SUBROUTINE .....
048561,000270:                                                                                                  #  SUBROUTINE TO SET RE (EARTH OR MOON RADIUS)
048562,000271: 
048563,000272:                                                                                                  #        RE = RM FOR MOON
048564,000273:                                                                                                  #        RE = RREF FOR FIXED EARTH RADIUS OR COMPUTED RF FOR FISCHER ELLIPSOID
048565,000274: 
048566,000275:                                                                                                  #  CALLING SEQUENCE
048567,000276:                                                                                                  #        L       CALL
048568,000277:                                                                                                  #        L+1             SETRE
048569,000278: 
048570,000279:                                                                                                  #  SUBROUTINES USED
048571,000280:                                                                                                  #        GETERAD
048572,000281: 
048573,000282:                                                                                                  #  INPUT
048574,000283:                                                                                                  #        ERADFLAG = 0 FOR FIXED RE, 1 FOR COMPUTED RE
048575,000284:                                                                                                  #        ALPHAV +4 = 1/2 SINL IF GETERAD IS CALLED
048576,000285:                                                                                                  #        LUNAFLAG = 0 FOR EARTH, =1 FOR MOON
048577,000286: 
048578,000287:                                                                                                  #  OUTPUT
048579,000288:                                                                                                  #        ERADM = 504RM FOR MOON (METERS B-29)
048580,000289:                                                                                                  #        ERADM = ERAD OR COMPUTED RF FOR EARTH (METERS B-29)
048581,000290: 
048582,000291: 13,2560           71220        SETRE              STQ      DLOAD                                 
048583,000292: 13,2561           00051                                    SETREX                                
048584,000293: 13,2562           10003                                    504RM                                 
048585,000294: 13,2563           71214                           BON      DLOAD                                 #  BRANCH FOR MOON
048586,000295: 13,2564           01703                                    LUNAFLAG                              
048587,000296: 13,2565           26575                                    TSTRLSRM                              
048588,000297: 13,2566           10001                                    ERAD                                  
048589,000298: 13,2567           45014                           BOFF     CALL                                  #  ERADFLAG=0 FOR FIXED RE, 1 FOR COMPUTED
048590,000299: 13,2570           00742                                    ERADFLAG                              
048591,000300: 13,2571           26573                                    SETRXX                                
048592,000301: 13,2572           26466                                    GETERAD                               
048593,000302: 13,2573           37673        SETRXX             STCALL   ERADM                                 #  EXIT WITH RE OR RM METERS B-29
048594,000303: 13,2574           00051                                    SETREX                                
048595,000304: 13,2575           77214        TSTRLSRM           BON      VLOAD                                 #  ERADFLAG=0, SET R0=RLS
048596,000305: 13,2576           00702                                    ERADFLAG                              #          =1      R0=RM
048597,000306: 13,2577           26573                                    SETRXX                                
048598,000307: 13,2600           02023                                    RLS                                   
048599,000308: 13,2601           64446                           ABVAL    SR2R                                  #  SCALE FROM B-27 TO B-29
048600,000309: 13,2602           77650                           GOTO                                           
048601,000310: 13,2603           26573                                    SETRXX                                
048602,000311: 13,2604  0000051               SETREX             =        S2                                    
048603,000312: 
048604,000313: 

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