Source Code
|
These source-code files are part of a reconstructed copy of Artemis 071, the
first release of the Apollo 15 through 17 Command Module (CM) Apollo Guidance
Computer (AGC) software.
The reconstruction began with source code of Artemis 072 previously transcribed
from a digitized copy of that program. The code was then updated by undoing
changes between the two versions. The reconstruction was verified by matching
memory-bank checksums to those listed in drawing 2021154-. Note that page
numbers in the reconstructed code match those on the Artemis 072 printout,
although the changed code would likely have changed page numbers for a real
Artemis 071 listing.
Comments from the original source code are prefixed with a single '#' symbol,
whereas comments added later are prefixed by "##" or "###". Report any errors
noted by creating an
issue report at the Virtual AGC
project's GitHub repository.
|
052750,000002: ## Copyright: Public domain.
052751,000003: ## Filename: LATITUDE_LONGITUDE_SUBROUTINES.agc
052752,000004: ## Purpose: A section of Artemis revision 071.
052753,000005: ## It is part of the reconstructed source code for the first
052754,000006: ## release of the flight software for the Command Module's
052755,000007: ## (CM) Apollo Guidance Computer (AGC) for Apollo 15 through
052756,000008: ## 17. The code has been recreated from a copy of Artemis 072.
052757,000009: ## It has been adapted such that the resulting bugger words
052758,000010: ## exactly match those specified for Artemis 071 in NASA
052759,000011: ## drawing 2021154-, which gives relatively high confidence
052760,000012: ## that the reconstruction is correct.
052761,000013: ## Reference: 1236
052762,000014: ## Assembler: yaYUL
052763,000015: ## Contact: Ron Burkey <info@sandroid.org>.
052764,000016: ## Website: www.ibiblio.org/apollo/index.html
052765,000017: ## Mod history: 2019-08-14 MAS Created from Artemis 072.
052766,000018:
|
|
Page 1236
|
052768,000020: # SUBROUTINE TO CONVERT RAD VECTOR AT GIVEN TIME TO LAT,LONG AND ALT
052769,000021:
052770,000022: # CALLING SEQUENCE
052771,000023:
052772,000024: # L-1 CALL
052773,000025: # L LAT-LONG
052774,000026:
052775,000027: # SUBROUTINES USED
052776,000028:
052777,000029: # R-TO-RP, ARCTAN, SETGAMMA, SETRE
052778,000030:
052779,000031: # ERASABLE INIT. REQ.
052780,000032:
052781,000033: # AXO, -AYO, AZO, TEPHEM (SET AT LAUNCH TIME)
052782,000034: # ALPHAV = POSITION VECTOR METERS B-29
052783,000035: # MPAC -- TIME (CSECS B-28)
052784,000036: # ERADFLAG =1, TO COMPUTE EARTH RADIUS, =0 FOR FIXED EARTH RADIUS
052785,000037: # LUNAFLAG=0 FOR EARTH, 1 FOR MOON
052786,000038:
052787,000039: # OUTPUT
052788,000040:
052789,000041: # LATITUDE IN LAT (REVS. B-0)
052790,000042: # LONGITUDE IN LONG (REVS. B-0)
052791,000043: # ALTITUDE IN ALT METERS B-29
052792,000044:
052793,000045: 13,2000 SETLOC LATLONG
052794,000046: 13,2000 BANK
052795,000047:
052796,000048: 13,2453 COUNT* $$/LT-LG
052797,000049: 13,2453 E4,1551 EBANK= ALPHAV
052798,000050: 13,2453 40220 LAT-LONG STQ SETPD
052799,000051: 13,2454 02241 INCORPEX
052800,000052: 13,2455 00001 0D
052801,000053: 13,2456 24007 STOVL 6D # SAVE TIME IN 6-7D FOR R-TO-RP
052802,000054: 13,2457 02152 ALPHAV
052803,000055: 13,2460 51406 PUSH ABVAL # 0-5D= R FOR R-TO-RP
052804,000056: 13,2461 16301 STODL ALPHAM # ABS. VALUE OF R FOR ALT FORMULA BELOW
052805,000057: 13,2462 11762 ZEROVEC # SET MPAC=0 FOR EARTH, NON-ZERO FOR MOON
052806,000058: 13,2463 71414 BOFF COS # USE COS(0) TO GET NON-ZERO IN MPAC
052807,000059: 13,2464 01743 LUNAFLAG # 0=EARTH,1=MOON
052808,000060: 13,2465 26466 CALLRTRP
052809,000061: 13,2466 77624 CALLRTRP CALL
052810,000062: 13,2467 55370 R-TO-RP # RP VECTOR CONVERTED FROM R B-29
052811,000063: 13,2470 77656 UNIT # UNIT RP B-1
052812,000064: 13,2471 36152 STCALL ALPHAV # U2= 1/2 SINL FOR SETRE SUBR BELOW
052813,000065: 13,2472 35766 SETGAMMA # SET GAMMA=B2/A2 FOR EARTH, =1 FOR MOON
052814,000066: 13,2473 77624 CALL # SCALED B-1
052815,000067: 13,2474 26654 SETRE # CALC RE METERS B-29
052816,000068: 13,2475 63545 DLOAD DSQ
052817,000069: 13,2476 02152 ALPHAV
052818,000070: 13,2477 63525 PDDL DSQ
052819,000071: 13,2500 02154 ALPHAV +2
052820,000072: 13,2501 75415 DAD SQRT
052821,000073: 13,2502 76405 DMP SL1R
|
|
Page 1237
|
052823,000075: 13,2503 00011 GAMRP
052824,000076: 13,2504 14021 STODL COSTH # COS(LAT) B-1
052825,000077: 13,2505 02156 ALPHAV +4
052826,000078: 13,2506 34023 STCALL SINTH # SIN(LAT) B-1
052827,000079: 13,2507 26614 ARCTAN
052828,000080: 13,2510 14747 STODL LAT # LAT B0
052829,000081: 13,2511 02152 ALPHAV
052830,000082: 13,2512 14021 STODL COSTH # COS(LONG) B-1
052831,000083: 13,2513 02154 ALPHAV +2
052832,000084: 13,2514 34023 STCALL SINTH # SIN(LONG) B-1
052833,000085: 13,2515 26614 ARCTAN
052834,000086: 13,2516 14751 STODL LONG # LONG. REVS B-0 IN RANGE -1/2 TO 1/2
052835,000087: 13,2517 02301 ALPHAM
052836,000088: 13,2520 77625 DSU # ALT= R-RE METERS B-29
052837,000089: 13,2521 02240 ERADM
052838,000090: 13,2522 34753 STCALL ALT # EXIT WITH ALT METERS B-29
052839,000091: 13,2523 02241 INCORPEX
052840,000092:
|
|
Page 1238
|
052842,000094: # SUBROUTINE TO CONVERT LAT,LONG,ALT AT GIVEN TIME TO RADIUS VECTOR
052843,000095:
052844,000096: # CALLING SEQUENCE
052845,000097:
052846,000098: # L-1 CALL
052847,000099: # L LALOTORV
052848,000100:
052849,000101: # SUBROUTINES USED
052850,000102:
052851,000103: # SETGAMMA, SETRE, RP-TO-R
052852,000104:
052853,000105: # ERASABLE INIT. REQ.
052854,000106:
052855,000107: # AXO, AYO, AZO, TEPHEM SET AT LAUNCH TIME
052856,000108: # LAT -- LATITUDE (REVS B0)
052857,000109: # LONG -- LONGITUDE (REVS B0)
052858,000110: # ALT -- ALTITUDE (METERS) B-29
052859,000111: # MPAC -- TIME (CSECS B-28)
052860,000112: # ERADFLAG =1 TO COMPUTE EARTH RADIUS, =0 FOR FIXED EARTH RADIUS
052861,000113: # LUNAFLAG=0 FOR EARTH, 1 FOR MOON
052862,000114:
052863,000115: # OUTPUT
052864,000116:
052865,000117: # R-VECTOR IN ALPHAV (METERS B-29)
052866,000118:
052867,000119: 13,2524 40220 LALOTORV STQ SETPD # LAT,LONG,ALT TO R VECTOR
052868,000120: 13,2525 02241 INCORPEX
052869,000121: 13,2526 00001 0D
052870,000122: 13,2527 34007 STCALL 6D # 6-7D = TIME FOR RP-TO-R
052871,000123: 13,2530 35766 SETGAMMA # GAMMA = B2/A2 FOR EARTH, 1 FOR MOON B-1
052872,000124: 13,2531 73545 DLOAD SIN # COS(LONG)COS(LAT) IN MPAC
052873,000125: 13,2532 00747 LAT # UNIT RP = SIN(LONG)COS(LAT) 2-3D
052874,000126: 13,2533 65275 DMPR PDDL # PD 2 GAMMA*SIN(LAT) 0-1D
052875,000127: 13,2534 00011 GAMRP
052876,000128: 13,2535 00747 LAT # 0-1D = GAMMA*SIN(LAT) B-2
052877,000129: 13,2536 65346 COS PDDL # PD4 2-3D = COS(LAT) B-1 TEMPORARILY
052878,000130: 13,2537 00751 LONG
052879,000131: 13,2540 57356 SIN DMPR # PD 2
052880,000132: 13,2541 71525 PDDL COS # PD 4 2-3D = SIN(LONG)COS(LAT) B-2
052881,000133: 13,2542 00747 LAT
052882,000134: 13,2543 71525 PDDL COS # PD 6 4-5D = COS(LAT) B-1 TEMPORARILY
052883,000135: 13,2544 00751 LONG
052884,000136: 13,2545 55475 DMPR VDEF # PD4 MPAC = COS(LONG)COS(LAT) B-2
052885,000137: 13,2546 41456 UNIT PUSH # 0-5D = UNIT RP FOR RP-TO-R SUBR.
052886,000138: 13,2547 36152 STCALL ALPHAV # ALPHAV +4= SINL FOR SETRE SUBR.
052887,000139: 13,2550 26654 SETRE # RE METERS B-29
052888,000140: 13,2551 43145 DLOAD BOFF # SET MPAC = 0 FOR EARTH, NON-ZERO FOR MOON
052889,000141: 13,2552 11762 ZEROVEC
052890,000142: 13,2553 01743 LUNAFLAG
052891,000143: 13,2554 26556 CALLRPRT
052892,000144: 13,2555 77746 COS # USE COS(0) TO GET NON-ZERO IN MPAC
052893,000145: 13,2556 77624 CALLRPRT CALL
052894,000146: 13,2557 55345 RP-TO-R # EXIT WITH UNIT R VECTOR IN MPAC
052895,000147: 13,2560 16152 STODL ALPHAV
052896,000148: 13,2561 02240 ERADM
|
|
Page 1239
|
052898,000150: 13,2562 74215 DAD VXSC # (RE + ALT)(UNIT R) METERS B-30
052899,000151: 13,2563 00753 ALT
052900,000152: 13,2564 02152 ALPHAV
052901,000153: 13,2565 77772 VSL1 # R METERS B-29
052902,000154: 13,2566 36152 STCALL ALPHAV # EXIT WITH R IN METERS B-29
052903,000155: 13,2567 02241 INCORPEX
052904,000156:
052905,000157: # SUBROUTINE TO COMPUTE EARTH RADIUS
052906,000158:
052907,000159: # INPUT
052908,000160:
052909,000161: # 1/2 SIN LAT IN ALPHAV +4
052910,000162:
052911,000163: # OUTPUT
052912,000164:
052913,000165: # EARTH RADIUS IN ERADM AND MPAC (METERS B-29)
052914,000166:
052915,000167: 13,2570 63545 GETERAD DLOAD DSQ
052916,000168: 13,2571 02156 ALPHAV +4 # SIN**2(L)
052917,000169: 13,2572 44352 SL1 BDSU
052918,000170: 13,2573 11760 DP1/2 # COS**2(L)
052919,000171: 13,2574 44275 DMPR BDSU
052920,000172: 13,2575 26611 EE
052921,000173: 13,2576 11760 DP1/2
052922,000174: 13,2577 75465 BDDV SQRT
052923,000175: 13,2600 26605 B2XSC
052924,000176: 13,2601 77622 SR4R
052925,000177: 13,2602 02240 STORE ERADM
052926,000178: 13,2603 77616 RVQ
052927,000179:
052928,000180: # THE FOLLOWING CONSTANTS WERE COMPUTED WITH A=6378166, B=6356784 METERS
052929,000181: # B2XSC = B**2 SCALED B-51
052930,000182: # B2/A2 = B**2/A**2 SCALED B-1
052931,000183: # EE = (1-B**2/A**2) SCALED B-0
052932,000184:
052933,000185: 13,2604 00446 00305 B2XSC 2DEC .0179450689 # B**2 SCALED B-51
052934,000186: 13,2606 04,3757 DP1/2 = XUNIT
052935,000187: 13,2606 17711 05254 B2/A2 2DEC .9933064884 B-1 # GAMMA= B**2/A**2 B-1
052936,000188: 13,2610 00155 25250 EE 2DEC 6.6935116 E-3 # (1-B**2/A**2) B-0
052937,000189: 13,2612 00302 17755 ERAD 2DEC 6373338 B-29 # PAD RADIUS
052938,000190:
|
|
Page 1240
|
052940,000192: # ARCTAN SUBROUTINE
052941,000193:
052942,000194: # CALLING SEQUENCE
052943,000195:
052944,000196: # SIN THETA IN SINTH B-1
052945,000197: # COS THETA IN COSTH B-1
052946,000198: # CALL ARCTAN
052947,000199:
052948,000200: # OUTPUT
052949,000201: # ARCTAN THETA IN MPAC AND THETA B-0 IN RANGE -1/2 TO +1/2
052950,000202:
052951,000203: 13,2614 77600 ARCTAN BOV
052952,000204: 13,2615 26616 CLROVFLW
052953,000205: 13,2616 63545 CLROVFLW DLOAD DSQ
052954,000206: 13,2617 00023 SINTH
052955,000207: 13,2620 63525 PDDL DSQ
052956,000208: 13,2621 00021 COSTH
052957,000209: 13,2622 77615 DAD
052958,000210: 13,2623 75454 BZE SQRT
052959,000211: 13,2624 26642 ARCTANXX # ATAN=0/0 SET THETA=0
052960,000212: 13,2625 40065 BDDV BOV
052961,000213: 13,2626 00023 SINTH
052962,000214: 13,2627 26647 ATAN=90
052963,000215: 13,2630 67542 SR1 ASIN
052964,000216: 13,2631 00025 STORE THETA
052965,000217: 13,2632 50125 PDDL BMN
052966,000218: 13,2633 00021 COSTH
052967,000219: 13,2634 26636 NEGCOS
052968,000220: 13,2635 43545 DLOAD RVQ
052969,000221: 13,2636 57545 NEGCOS DLOAD DCOMP
052970,000222: 13,2637 43244 BPL DAD
052971,000223: 13,2640 26644 NEGOUT
052972,000224: 13,2641 11760 DP1/2
052973,000225: 13,2642 00025 ARCTANXX STORE THETA
052974,000226: 13,2643 77616 RVQ
052975,000227:
052976,000228: 13,2644 52025 NEGOUT DSU GOTO
052977,000229: 13,2645 11760 DP1/2
052978,000230: 13,2646 26642 ARCTANXX
052979,000231: 13,2647 75345 ATAN=90 DLOAD SIGN
052980,000232: 13,2650 33412 LODP1/4
052981,000233: 13,2651 00023 SINTH
052982,000234: 13,2652 00025 STORE THETA
052983,000235: 13,2653 77616 RVQ
052984,000236:
052985,000237: 13,2654 04,3761 2DZERO = DPZERO
052986,000238:
|
|
Page 1241
|
052988,000240: # ..... SETGAMMA SUBROUTINE .....
052989,000241: # SUBROUTINE TO SET GAMMA FOR THE LAT-LONG AND LALOTORV SUBROUTINES
052990,000242:
052991,000243: # GAMMA = B**2/A**2 FOR EARTH (B-1)
052992,000244: # GAMMA = 1 FOR MOON (B-1)
052993,000245:
052994,000246: # CALLING SEQUENCE
052995,000247: # L CALL
052996,000248: # L+1 SETGAMMA
052997,000249:
052998,000250: # INPUT
052999,000251: # LUNAFLAG=0 FOR EARTH, =1 FOR MOON
053000,000252:
053001,000253: # OUTPUT
053002,000254: # GAMMA IN GAMRP (B-1)
053003,000255:
053004,000256: 16,2000 SETLOC LATLONG1
053005,000257: 16,2000 BANK
053006,000258: 16,3766 COUNT* $$/LT-LG
053007,000259:
053008,000260: 16,3766 43145 SETGAMMA DLOAD BOFF # BRANCH FOR EARTH
053009,000261: 16,3767 26607 B2/A2 # EARTH GAMMA
053010,000262: 16,3770 01743 LUNAFLAG
053011,000263: 16,3771 35774 SETGMEX
053012,000264: 16,3772 77735 SLOAD
053013,000265: 16,3773 11760 1B1 # MOON GAMMA
053014,000266: 16,3774 00011 SETGMEX STORE GAMRP
053015,000267: 16,3775 77616 RVQ
053016,000268: 16,3776 GAMRP = 8D
053017,000269: 13,2000 SETLOC LATLONG
053018,000270: 13,2000 BANK
053019,000271:
|
|
Page 1242
|
053021,000273: # ..... SETRE SUBROUTINE .....
053022,000274: # SUBROUTINE TO SET RE (EARTH OR MOON RADIUS)
053023,000275:
053024,000276: # RE = RM FOR MOON
053025,000277: # RE = RREF FOR FIXED EARTH RADIUS OR COMPUTED RF FOR FISCHER ELLIPSOID
053026,000278:
053027,000279: # CALLING SEQUENCE
053028,000280: # L CALL
053029,000281: # L+1 SETRE
053030,000282:
053031,000283: # SUBROUTINES USED
053032,000284: # GETERAD
053033,000285:
053034,000286: # INPUT
053035,000287: # ERADFLAG = 0 FOR FIXED RE, 1 FOR COMPUTED RE
053036,000288: # ALPHAV +4 = 1/2 SINL IF GETERAD IS CALLED
053037,000289: # LUNAFLAG = 0 FOR EARTH, =1 FOR MOON
053038,000290:
053039,000291: # OUTPUT
053040,000292: # ERADM = 504RM FOR MOON (METERS B-29)
053041,000293: # ERADM = ERAD OR COMPUTED RF FOR EARTH (METERS B-29)
053042,000294:
053043,000295: 13,2654 71220 SETRE STQ DLOAD
053044,000296: 13,2655 00051 SETREX
053045,000297: 13,2656 26701 504RM
053046,000298: 13,2657 71214 BON DLOAD # BRANCH FOR MOON
053047,000299: 13,2660 01703 LUNAFLAG
053048,000300: 13,2661 26671 TSTRLSRM
053049,000301: 13,2662 26613 ERAD
053050,000302: 13,2663 45014 BOFF CALL # ERADFLAG=0 FOR FIXED RE, 1 FOR COMPUTED
053051,000303: 13,2664 00742 ERADFLAG
053052,000304: 13,2665 26667 SETRXX
053053,000305: 13,2666 26570 GETERAD
053054,000306: 13,2667 36240 SETRXX STCALL ERADM # EXIT WITH RE OR RM METERS B-29
053055,000307: 13,2670 00051 SETREX
053056,000308: 13,2671 77214 TSTRLSRM BON VLOAD # ERADFLAG=0, SET R0=RLS
053057,000309: 13,2672 00702 ERADFLAG # =1 R0=RM
053058,000310: 13,2673 26667 SETRXX
053059,000311: 13,2674 02026 RLS
053060,000312: 13,2675 64446 ABVAL SR2R # SCALE FROM B-27 TO B-29
053061,000313: 13,2676 77650 GOTO
053062,000314: 13,2677 26667 SETRXX
053063,000315: 13,2700 0000051 SETREX = S2
053064,000316: 13,2700 00065 01265 504RM 2DEC 1738090 B-29 # METERS B-29 (MOON RADIUS)
End of include-file LATITUDE_LONGITUDE_SUBROUTINES.agc. Parent file is MAIN.agc