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,1971Note 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