GAP: ASSEMBLE REVISION 072 OF AGC PROGRAM ARTEMIS BY NASA 2021114-011 11:40 FEB. 26, 1971Note that the date is the date of the printout, not the date of the program revision.
Source Code
|
|
These source-code files were obtained by digitally photographing
an Artemis 72 (Apollo 15-17 Command Module) program listing from a private
collection and then
modifying pre-existing Comanche 55 (Apollo 11) source files to
incorporate changes. Photography was by Ron Burkey, and transcription
was performed by a team of volunteers. Note that the page images
presented online are of reduced quality, and that additional images
may be 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 072 OF AGC PROGRAM ARTEMIS BY NASA 2021114-011 11:40 FEB. 26, 1971Note that the date is the date of the printout, not the date of the program revision. |
058781,000002: ## Copyright: Public domain.
058782,000003: ## Filename: TIME_OF_FREE_FALL.agc
058783,000004: ## Purpose: Part of the source code for Artemis (i.e., Colossus 3),
058784,000005: ## build 072. This is for the Command Module's (CM)
058785,000006: ## Apollo Guidance Computer (AGC), for
058786,000007: ## Apollo 15-17.
058787,000008: ## Assembler: yaYUL
058788,000009: ## Contact: Jim Lawton <jim DOT lawton AT gmail DOT com>
058789,000010: ## Website: www.ibiblio.org/apollo/index.html
058790,000011: ## Page Scans: www.ibiblio.org/apollo/ScansForConversion/Artemis072/
058791,000012: ## Mod history: 2009-08-20 JL Adapted from corresponding Comanche 055 file.
058792,000013: ## 2009-08-21 JL Added a note about a problem with the page scan on p.1386.
058793,000014: ## 2010-02-20 RSB Un-##'d this header.
058794,000015: ## 2017-02-08 RSB Proofed comment text by diff'ing vs Comanche 55
058795,000016: ## and/or octopus/ProoferComments as most-appropriate.
058796,000017: ## 2017-03-14 RSB Comment-text fixes noted in proofing Luminary 116.
058797,000018:
|
|
Page 1371 |
058799,000020:
058800,000021: # THE TFF SUBROUTINES MAY BE USED IN EITHER EARTH OR MOON CENTERED COORDINATES. THE TFF ROUTINES NEVER
058801,000022: # KNOW WHICH ORIGIN APPLIES. IT IS THE USER WHO KNOWS, AND WHO SUPPLIES RONE, VONE AND 1/SQRT(MU) AT THE
058802,000023: # APPROPRIATE SCALE LEVEL FOR THE PROPER PRIMARY BODY.
058803,000024:
058804,000025: # EARTH ORIGIN POSITION -29 METERS
058805,000026: # VELOCITY -7 METERS/CENTISECOND
058806,000027: # 1/SQRT(MU) +17 SQRT(CS SQ/METERS CUBED)
058807,000028:
058808,000029: # MOON ORIGIN POSITION -27 METERS
058809,000030: # VELOCITY -5 METERS/CENTISECONDS
058810,000031: # 1/SQRT(MU) +14 SQRT(CS SQ/METERS CUBED)
058811,000032:
058812,000033: # ALL DATA PROVIDED TO AND RECEIVED FROM ANY TFF SUBROUTINE WILL BE AT ONE OF THE LEVELS ABOVE. IN ALL CASES,
058813,000034: # THE FREE FALL TIME IS RETURNED IN CENTISECONDS AT (-28). PROGRAM TFF/CONIC WILL GENERATE VONE/RTMU AND
058814,000035: # LEAVE IT IN VONE' AT (+10) IF EARTH ORIGIN AND (+9) IF MOON ORIGIN.
058815,000036:
058816,000037: # THE USER MUST STORE THE STATE VECTOR IN RONE, VONE AND MU IN THE FORM 1/SQRT(MU) IN TFF/RTMU
058817,000038: # AT THE PROPER SCALE BEFORE CALLING TFF/CONIC. SINCE RONE, VONE ARE IN THE EXTENDED VERB STORAGE AREA,
058818,000039: # THE USER MUST ALSO LOCK OUT THE EXTENDED VERBS, AND RELEASE THEM WHEN FINISHED.
058819,000040:
058820,000041: # PROGRAMS CALC/TFF AND CALC/TPER ASSUME THAT THE TERMINAL RADIUS IS LESS THAN THE PRESENT
058821,000042: # RADIUS. THIS RESTRICTION CAN BE REMOVED BY A 15 W CODING CHANGE, BUT AT PRESENT IT IS NOT DEEMED NECESSARY.
058822,000043:
058823,000044: # THE FOLLOWING ERASABLE QUANTITIES ARE USED BY THE TFF ROUTINES, AND ARE LOCATED IN THE PUSH LIST.
058824,000045:
058825,000046: # BELOW E: IS USED FOR EARTH ORIGIN SCALE
058826,000047: # M: IS USED FOR MOON ORIGIN SCALE
058827,000048:
058828,000049: # TFFSW = 119D BIT1 0 = CALCTFF 1 = CALCTPER
058829,000050: 23,3773 TFFDELQ = 10D # Q2-Q1 E: (-16) M: (-15)
058830,000051: 23,3773 RMAG1 = 12D # ABVAL(RN) M E: (-29) M: (-27)
058831,000052: # RPER = 14D PERIGEE RADIUS M E: (-29) M: (-27)
058832,000053: 23,3773 TFFQ1 = 14D # R.V / SQRT(MUE) E: (-16) M: (-15)
058833,000054: # SDELF/2 SIN(THETA) /2
058834,000055: 23,3773 CDELF/2 = 14D # COS(THETA) /2
058835,000056: # RAPO = 16D APOGEE RADIUS M E: (-29) M: (-27)
058836,000057: 23,3773 NRTERM = 16D # TERMINAL RADIUS M E: (-29+NR)
058837,000058: # M: (-27+NR)
058838,000059: 23,3773 RTERM = 18D # TERMINAL RADIUS M E: (-29) M: (-27)
058839,000060: 23,3773 TFFVSQ = 20D # -(V SQUARED/MU) 1/M E: (20) M: (18)
058840,000061: 23,3773 TFF1/ALF = 22D # SEMI MAJ AXIS M E: (-22-2 NA)
058841,000062: # M: (-20-2 NA)
058842,000063: 23,3773 TFFRTALF = 24D # SQRT(ALFA) E: (10+NA) M: (9+NA)
058843,000064: 23,3773 TFFALFA = 26D # ALFA 1/M E: (26-NR) M: (24-NR)
058844,000065: 23,3773 TFFNP = 28D # SEMI LATUS RECTUM M E: (-38+2 NR)
058845,000066: # M: (-36+2 NR)
058846,000067: 23,3773 TFF/RTMU = 30D # 1/SQRT(MU) E: (17) M: (14)
058847,000068: 23,3773 NRMAG = 32D # PRESENT RADIUS M E: (-29+NR)
058848,000069: # M: (-27+NR)
058849,000070: 23,3773 TFFX = 34D
058850,000071: 23,3773 TFFTEM = 36D # TEMPORARY
058851,000072:
|
|
Page 1372 |
058853,000074: # REGISTERS S1, S2 ARE UNTOUCHED BY ANY TFF SUBROUTINE
058854,000075: # INDEX REGISTERS X1, X2 ARE USED BY ALL TFF SUBROUTINES. THEY ARE ESTAB-
058855,000076: # LISHED IN TFF/CONIC AND MUST BE PRESERVED BETWEEN CALLS TO SUBSEQUENT
058856,000077: # SUBROUTINES.
058857,000078: # -NR C(X1) = NORM COUNT OF RMAG
058858,000079: # -NA C(X2) = NORM COUNT OF SQRT(ABS(ALFA))
058859,000080:
|
|
Page 1373 |
058861,000082:
058862,000083: # SUBROUTINE NAME: TFFCONIC DATE: 01.29.67
058863,000084: # MOD NO: 0 LOG SECTION: TIME OF FREE FALL
058864,000085: # MOD BY: RR BAIRNSFATHER
058865,000086: # MOD NO: 1 MOD BY: RR BAIRNSFATHER DATE: 11 APR 67
058866,000087: # MOD NO: 2 MOD BY: RR BAIRNSFATHER DATE: 21 NOV 67 ADD MOON MU.
058867,000088: # MOD NO: 3 MOD BY: RR BAIRNSFATHER DATE: 21 MAR 68 ACCEPT DIFFERENT EARTH/MOON SCALES
058868,000089:
058869,000090: # FUNCTIONAL DESCRIPTION: THIS SUBROUTINE IS CALLED TO COMPUTE THOSE CONIC PARAMETERS REQUIRED BY THE TFF
058870,000091: # SUBROUTINES AND TO ESTABLISH THEM IN THE PUSH LIST AREA. THE PARAMETERS ARE LISTED UNDER OUTPUT.
058871,000092: # THE EQUATIONS ARE
058872,000093: # - - -
058873,000094: # H = RN*VN ANGULAR MOMENTUM
058874,000095: # - -
058875,000096: # LCP = H.H / MU SEMI LATUS RECTUM
058876,000097: # - -
058877,000098: # ALFA = 2/RN - VN.VN / MU RECIPROCAL SEMI MAJ AXIS, SIGNED
058878,000099:
058879,000100: # AND ALFA IS POS FOR ELLIPTIC ORBITS
058880,000101: # 0 FOR PARABOLIC ORBITS
058881,000102: # NEG FOR HYPERBOLIC ORBITS.
058882,000103: # SUBROUTINE ALSO COMPUTES AND SAVES RMAG.
058883,000104:
058884,000105: # CALLING SEQUENCE:
058885,000106: # TFFCONIC EXPECTS CALLER TO ENTER WITH CORRECT GRAVITATIONAL CONSTANT IN MPAC, IN THE FORM
058886,000107: # 1/SQRT(MU). PROGRAM WILL SAVE IN TFF/RTMU . THE SCALE IS DETERMINED BY WHETHER EARTH OR MOON
058887,000108: # ORIGIN IS USED. THE CALLER MUST LOCK OUT THE EXTENDED VERBS BEFORE PROVIDING STATE VECTOR IN RONE,
058888,000109: # VONE AT PROPER SCALE. THE EXTENDED VERBS MUST BE RESTORED WHEN THE CALLER IS FINISHED USING THE
058889,000110: # TFF ROUTINES.
058890,000111:
058891,000112: # ENTRY POINT TFFCONMU EXPECTS THAT TFF/RTMU IS ALREADY LOADED.
058892,000113:
058893,000114: # TO SPECIFY MU: DLOAD CALL IF MU ALREADY STORED: CALL
058894,000115: # YOURMU 1/RTMU E: (17) M: (14) TFFCONMU
058895,000116: # TFFCONIC
058896,000117:
058897,000118: # PUSHLOC = PDL+0, ARBITRARY IF LEQ 18D
058898,000119:
058899,000120: # SUBROUTINES CALLED: NONE
058900,000121:
058901,000122: # NORMAL EXIT MODES: RVQ
058902,000123:
058903,000124: # ALARMS: NONE
058904,000125:
058905,000126: # OUTPUT: THE FOLLOWING ARE STORED IN THE PUSH LIST AREA.
058906,000127: # RMAG1 E:(-29) M:(-27) M RN, PRESENT RADIUS LENGTH.
058907,000128: # NRMAG E: (-29+NR) M RMAG, NORMALIZED
058908,000129: # M: (-27+NR)
058909,000130: # X1 -NR, NORM COUNT
058910,000131: # TFFNP E: (-38+2NR) M LCP, SEMI LATUS RECTUM, WEIGHTED BY NR. FOR VGAMCALC
058911,000132: # M: (-36+2NR)
058912,000133: # TFF/RTMU E:(17) M:(14) 1/SQRT(MU)
058913,000134: # TFFVSQ E:(20) M:(18) 1/M -(V SQ/MU): PRESENT VELOCITY, NORMLIZED. FOR VGAMCALC
058914,000135: # TFFALFA E: (26-NR) 1/M ALFA, WEIGHTED BY NR
058915,000136: # M: (24-NR)
058916,000137: # TFFRTALF E: (10+NA) SQRT(ALFA), NORMALIZED
058917,000138: # M: (9+NA)
|
|
Page 1374 |
058919,000140: # X2 -NA, NORM COUNT
058920,000141: # TFF1/ALF E: (-22-2NA) SIGNED SEMI MAJ AXIS, WEIGHTED BY NA
058921,000142: # M: (-20-2NA)
058922,000143: # PUSHLOC AT PDL+0
058923,000144: # THE FOLLOWING IS STORED IN GENERAL ERASABLE
058924,000145: # VONE' E:(10) M:(9) V/RT(MU), NORMALIZED VELOCITY
058925,000146:
058926,000147: # ERASABLE INITIALIZATION REQUIRED:
058927,000148: # RONE E:(-29) M:(-27) M STATE VECTOR LEFT BY CALLER
058928,000149: # VONE E:(-7) M:(-5) M/CS STATE VECTOR LEFT BY CALLER
058929,000150: # TFF/RTMU E:(17) M:(14) 1/RT(CS SQ/M CUBE) IF ENTER VIA TFFCONMU.
058930,000151:
058931,000152: # DEBRIS: QPRET, PDL+0 ... PDL+3
058932,000153:
058933,000154: 27,2000 SETLOC TOF-FF
058934,000155: 27,2000 BANK
058935,000156:
058936,000157: 27,3151 COUNT* $$/TFF
058937,000158:
058938,000159: 27,3151 00037 TFFCONIC STORE TFF/RTMU # 1/SQRT(MU) E: (17) M: (14)
058939,000160:
058940,000161: 27,3152 53575 TFFCONMU VLOAD UNIT # COME HERE WITH TFFRTMU LOADED.
058941,000162: 27,3153 02320 RONE # SAVED RN. M E: (-29) M: (-27)
058942,000163: 27,3154 77725 PDDL # UR/2 TO PDL+0, +5
058943,000164: 27,3155 00045 36D # MAGNITUDE
058944,000165: 27,3156 00015 STORE RMAG1 # M E: (-29) M: (-27)
058945,000166:
058946,000167: 27,3157 77701 NORM
058947,000168: 27,3160 00047 X1 # -NR
058948,000169: 27,3161 24041 STOVL NRMAG # RMAG M E: (-29+NR) M: (-27+NR)
058949,000170: 27,3162 02326 VONE # SAVED VN. M/CS E: (-7) M: (-5)
058950,000171: 27,3163 77761 VXSC
058951,000172: 27,3164 00037 TFF/RTMU # E: (17) M: (14)
058952,000173: 27,3165 02362 STORE VONE' # VN/SQRT(MU) E:(10) M:(9)
058953,000174:
058954,000175: 27,3166 47361 VXSC VXV
058955,000176: 27,3167 00041 NRMAG # E: (-29+NR) M: (-27+NR)
058956,000177: # UR/2 FROM PDL
058957,000178: 27,3170 47572 VSL1 VSQ # BEFORE: E:(-19+NR) M:(-18+NR)
058958,000179: 27,3171 14035 STODL TFFNP # LC P M E:(-38+2NR) M:(-36+2NR)
058959,000180: # SAVE ALSO FOR VGAMCALC
058960,000181: 27,3172 15326 TFF1/4
058961,000182: 27,3173 63271 DDV PDVL # (2/RMAG) 1/M E:(26-NR) M:(24-NR)
058962,000183: 27,3174 00041 NRMAG # RMAG M E:(-29+NR) M:(-27+NR)
058963,000184: 27,3175 02362 VONE' # SAVED VN. E: (10) M: (9)
058964,000185: 27,3176 57436 VSQ DCOMP # KEEP MPAC+2 HONEST FOR SQRT.
058965,000186: 27,3177 00025 STORE TFFVSQ # -(V SQ/MU) E:(20) M:(18)
058966,000187: # SAVE FOR VGAMCALC
058967,000188: 27,3200 43257 SR* DAD
058968,000189: 27,3201 20573 0 -6,1 # GET -VSQ/MU E:(26-NR) M:(24-NR)
|
|
Page 1375 |
058970,000191: 27,3202 77626 STADR
058971,000192: # 2/RMAG FROM PDL+2
058972,000193: 27,3203 77744 STORE TFFALFA # ALFA 1/M E:(26-NR) M:(24-NR)
058973,000194: 27,3204 41457 SL* PUSH # TEMP SAVE ALFA E:(20) M:(18)
058974,000195: 27,3205 20173 0 -6,1
058975,000196: 27,3206 75446 ABS SQRT # E:(10) M:(9)
058976,000197: 27,3207 77701 NORM
058977,000198: 27,3210 00050 X2 # X2 = -NA
058978,000199: 27,3211 00031 STORE TFFRTALF # SQRT( ABS(ALFA) ) E:(10+NA) M:(9+NA)
058979,000200: 27,3212 75316 DSQ SIGN # NOT SO ACCURATE, BUT OK
058980,000201: # ALFA FROM PDL+2 E:(20) M:(18)
058981,000202: 27,3213 55254 BZE BDDV # SET 1/ALFA =0, TO SHOW SMALL ALFA
058982,000203: 27,3214 57216 +2
058983,000204: 27,3215 15326 TFF1/4
058984,000205: 27,3216 00027 +2 STORE TFF1/ALF # 1/ALFA E:(-22-2 NA) M:(-20-2 NA)
058985,000206: 27,3217 77616 DUMPCNIC RVQ
058986,000207:
058987,000208: # 39 W
058988,000209:
|
|
Page 1376 |
058990,000211:
058991,000212: # SUBROUTINE NAME: TFFRP/RA DATE: 01.17.67
058992,000213: # MOD NO: 0 LOG SECTION: TIME OF FREE FALL
058993,000214: # MOD BY: RR BAIRNSFATHER
058994,000215: # MOD NO: 1 MOD BY: RR BAIRNSFATHER DATE: 11 APR 67
058995,000216: # MOD NO: 2 MOD BY: RR BAIRNSFATHER DATE: 21 MAR 68 ACCEPT DIFFERENT EARTH/MOON SCALES
058996,000217: # ALSO IMPROVE ACCURACY OF RAPO.
058997,000218:
058998,000219: # FUNCTIONAL DESCRIPTION: USED BY CALCTPER AND TFF DISPLAYS TO CALCULATE PERIGEE RADIUS AND ALSO
058999,000220: # APOGEE RADIUS FOR A GENERAL CONIC.
059000,000221: # PROGRAM GIVES PERIGEE RADIUS AS APOGEE RADIUS IS GIVEN BY
059001,000222: # RP = P /(1+E) RA = (1+E) / ALFA
059002,000223: # WHERE 2
059003,000224: # E = 1 - P ALFA
059004,000225: # IF RA IS NEGATIVE OR SHOWS DIVIDE OVERFLOW, THEN RA = POSMAX BECAUSE
059005,000226: # 1. APOGEE RADIUS IS NOT MEANINGFUL FOR HYPERBOLA
059006,000227: # 2. APOGEE RADIUS IS NOT DEFINED FOR PARABOLA
059007,000228: # 3. APOGEE RADIUS EXCEEDS THE SCALING FOR ELLIPSE.
059008,000229:
059009,000230: # THIS SUBROUTINE REQUIRES THE SIGNED RECIPROCAL SEMI MAJ AXIS, ALFA, AND SEMI LATUS RECTUM AS DATA.
059010,000231:
059011,000232: # CALLING SEQUENCE: CALL
059012,000233: # TFFRP/RA
059013,000234:
059014,000235: # PUSHLOC = PDL+0, ARBITRARY IF LEQ 10D
059015,000236: # C(MPAC) UNSPECIFIED
059016,000237:
059017,000238: # SUBROUTINES CALLED: NONE
059018,000239:
059019,000240: # NORMAL EXIT MODE: RVQ
059020,000241: # IF ELLIPSE, WITHIN NORMAL SCALING, RAPO IS CORRECT.
059021,000242: # OTHERWISE, RAPO = POSMAX.
059022,000243:
059023,000244: # ALARMS: NONE
059024,000245:
059025,000246: # OUTPUT: STORED IN PUSH LIST AREA. SCALE OF OUTPUT AGREES WITH DATA SUPPLIED TO TFF/CONIC.
059026,000247: # RPER E:(-29) M:(-27) M PERIGEE RADIUS DESTROYED BY CALCTFF/CALCTPER, TFFTRIG.
059027,000248: # RAPO E:(-29) M:(-27) M APOGEE RADIUS WILL BE DESTROYED BY CALCTFF/CALCTPER
059028,000249: # PUSHLOC AT PDL+0
059029,000250:
059030,000251: # ERASABLE INITIALIZATION REQUIRED:
059031,000252: # TFFALFA E:(26-NR) M 1/SEMI MAJ AXIS LEFT BY TFFCONIC
059032,000253: # M:(24-NR)
059033,000254: # TFFNP E:(-38+2NR) M LC P, SEMI LATUS RECTUM LEFT BY TFFCONIC
059034,000255: # M:(-36+2NR)
059035,000256: # X1 -NR, NORM COUNT OF RMAG LEFT BY TFFCONIC
059036,000257: # X2 -NA, NORM COUNT OF ALFA LEFT BY TFFCONIC
059037,000258:
059038,000259: # DEBRIS: QPRET, PDL+0 ... PDL+1
059039,000260:
|
|
Page 1377 |
059041,000262: 27,3220 RAPO = 16D # APOGEE RADIUS M E:(-29) M:(-27)
059042,000263: 27,3220 RPER = 14D # PERIGEE RADIUS M E:(-29) M:(-27)
059043,000264:
059044,000265: 27,3220 41345 TFFRP/RA DLOAD DMP
059045,000266: 27,3221 00033 TFFALFA # ALFA 1/M E:(26-NR) M:(24-NR)
059046,000267: 27,3222 00035 TFFNP # LC P M E:(-38+2NR) M:(-36+2NR)
059047,000268: 27,3223 57457 SR* DCOMP # ALFA P (-12+NR)
059048,000269: 27,3224 20571 0 -8D,1 # ALFA P (-4)
059049,000270: 27,3225 51415 DAD ABS # (DCOMP GIVES VALID TP RESULT FOR SQRT)
059050,000271: # (ABS PROTECTS SQRT IF E IS VERY NEAR 0)
059051,000272: 27,3226 31771 DP2(-4)
059052,000273: 27,3227 43366 SQRT DAD # E SQ = (1- P ALFA) (-4)
059053,000274: 27,3230 15326 TFF1/4
059054,000275: 27,3231 55206 PUSH BDDV # (1+E) (-2) TO PDL+0
059055,000276: 27,3232 00035 TFFNP # LCP M E:(-38+2NR) M:(-36+2NR)
059056,000277: 27,3233 53657 SR* SR* # (DOES SR THEN SL TO AVOID OVFL)
059057,000278: 27,3234 20601 0,1 # X1=-NR
059058,000279: 27,3235 20572 0 -7,1 # (EFFECTIVE SL)
059059,000280: 27,3236 14017 STODL RPER # PERIGEE RADIUS M E:(-29) M:(-27)
059060,000281: # (1+E) (-2) FROM PDL+0
059061,000282: 27,3237 41005 DMP BOVB
059062,000283: 27,3240 00027 TFF1/ALF # E:(-22-2NA) M:(-20-2NA)
059063,000284: 27,3241 57544 TCDANZIG # CLEAR OVFIND, IF ON.
059064,000285: 27,3242 53654 BZE SL*
059065,000286: 27,3243 57252 MAXRA # SET POSMAX, IF ALFA=0
059066,000287: 27,3244 57603 0 -5,2 # -5+NA
059067,000288: 27,3245 40057 SL* BOV
059068,000289: 27,3246 57576 0,2
059069,000290: 27,3247 57252 MAXRA # SET POSMAX IF OVFL.
059070,000291: 27,3250 77644 BPL # CONTINUE WITH VALID RAPO.
059071,000292: 27,3251 57254 +3
059072,000293: 27,3252 77745 MAXRA DLOAD # RAPO CALC IS NOT VALID. SET RAPO =
059073,000294: 27,3253 15344 NEARONE # POSMAX AS A TAG.
059074,000295: 27,3254 00021 +3 STORE RAPO # APOGEE RADIUS M E:(-29) M:(-27)
059075,000296: 27,3255 77616 DUMPRPRA RVQ
059076,000297:
059077,000298: # 30 W
059078,000299:
|
|
Page 1378 |
059080,000301:
059081,000302: # SUBROUTINE NAME: CALCTPER / CALCTFF DATE: 01.29.67
059082,000303: # MOD NO: 0 LOG SECTION: TIME OF FREE FALL
059083,000304: # MOD BY: RR BAIRNSFATHER
059084,000305: # MOD NO: 1 MOD BY: RR BAIRNSFATHER DATE: 21 MAR 67
059085,000306: # MOD NO: 2 MOD BY: RR BAIRNSFATHER DATE: 14 APR 67
059086,000307: # MOD NO: 3 MOD BY: RR BAIRNSFATHER DATE: 8 JUL 67 NEAR EARTH MUE AND NEG TFF (GONEPAST)
059087,000308: # MOD NO: 4 MOD BY: RR BAIRNSFATHER DATE: 21 NOV 67 ADD VARIABLE MU.
059088,000309: # MOD NO: 5 MOD BY: RR BAIRNSFATHER DATE: 21 MAR 68 ACCEPT DIFFERENT EARTH/MOON SCALES
059089,000310:
059090,000311: # FUNCTIONAL DESCRIPTION: PROGRAM CALCULATES THE FREE-FALL TIME OF FLIGHT FROM PRESENT POSITION RN AND
059091,000312: # VELOCITY VN TO A RADIUS LENGTH SPECIFIED BY RTERM , SUPPLIED BY THE USER. THE POSITION VECTOR
059092,000313: # RN MAY BE ON EITHER SIDE OF THE CONIC, BUT RTERM IS CONSIDERED ON THE INBOUND SIDE.
059093,000314: # THE EQUATIONS ARE
059094,000315:
059095,000316: # Q2 = -SQRT(RTERM (2-RTERM ALFA) - LCP) (INBOUND SIDE)) LEQ +- LCE/SQRT(ALFA)
059096,000317: # - -
059097,000318: # Q1 = RN.VN / SQRT(MU) LEQ +- LCE/SQRT(ALFA)
059098,000319:
059099,000320: # Z = NUM / DEN LEQ +- 1/SQRT(ALFA)
059100,000321:
059101,000322: # WHERE, IF INBOUND
059102,000323: # NUM = RTERM - RN LEQ +- 2 LCE/ALFA
059103,000324: # DEN = Q2+Q1 LEQ +- 2 LCE/SQRT(ALFA)
059104,000325:
059105,000326: # AND, IF OUTBOUND
059106,000327: # NUM = Q2-Q1 LEQ +- 2 LCE/SQRT(ALFA)
059107,000328: # DEN = 2 - ALFA (RTERM + RN) . LEQ +- 2 LCE
059108,000329:
059109,000330: # IF ALFA ZZ < 1.0 (FOR ALL CONICS EXCEPT ELLIPSES HAVING ABS(DEL ECC ANOM) G 90 DEG)
059110,000331:
059111,000332: # THEN X = ALFA Z Z
059112,000333: # AND TFF = (RTERM +RN -2 ZZ T(X) ) Z/SQRT(MU)
059113,000334:
059114,000335: # EXCEPT IF ALFA PNZ, AND IF TFF NEG,
059115,000336: # THEN TFF = 2 PI /(ALFA SQRT(ALFA)) + TFF
059116,000337:
059117,000338: # OR IF ALFA ZZ GEQ 1.0 (FOR ELLIPSES HAVING ABS(DEL ECC ANOM) GEQ 90 DEG)
059118,000339:
059119,000340: # THEN X = 1/ALFA Z Z
059120,000341: # AND TFF = (PI/SQRT(ALFA) -Q2 +Q1 +2(X T(X) -1) /ALFA Z) /ALFA SQRT(MU)
059121,000342:
059122,000343: # WHERE T(X) IS A POLYNOMIAL APPROXIMATION TO THE SERIES
059123,000344: # 2 3 2
059124,000345: # 1/3 - X/5 + X /7 - X /9 ... (X < 1.0)
059125,000346:
059126,000347: # CALLING SEQUENCE: TIME TO RTERM TIME TO PERIGEE
059127,000348: # CALL CALL
059128,000349: # CALCTFF CALCTPER
059129,000350: # C(MPAC) = TERMNL RAD M C(MPAC) = PERIGEE RAD M
059130,000351: # FOR EITHER, E: (-29) M: (-27)
059131,000352: # FOR EITHER, PUSHLOC = PDL+0 , ARBITRARY IF LEQ 8D.
|
|
Page 1379 |
059133,000354:
059134,000355: # SUBROUTINES CALLED: T(X), VIA RTB
059135,000356:
059136,000357: # NORMAL EXIT MODE: RVQ
059137,000358: # HOWEVER, PROGRAM EXITS WITH ONE OF THE FOLLOWING VALUES FOR TFF (-28) CS IN MPAC. USER MUST STORE.
059138,000359: # A. TFF = FLIGHT TIME. NORMAL CASE FOR POSITIVE FLIGHT TIME LESS THAN ONE ORBITAL PERIOD.
059139,000360: # B. (THIS OPTION IS NO LONGER USED.)
059140,000361: # C. TFF = POSMAX. THIS INDICATES THAT THE CONIC FROM THE PRESENT POSITION WILL NOT RETURN TO
059141,000362: # THE SPECIFIED ALTITUDE. ALSO INDICATES OUTBOUND PARABOLA OR HYPERBOLA.
059142,000363:
059143,000364: # OUTPUT: C(MPAC) (-28) CS TIME OF FLIGHT, OR TIME TO PERIGEE
059144,000365: # TFFX (0) X, LEFT FOR ENTRY DISPLAY TFF ROUTINES
059145,000366: # NRTERM E: (-29+NR) M RTERM, WEIGHTED BY NR LEFT FOR ENTRY DISPLAY TFF ROUTINES
059146,000367: # M: (-27+NR)
059147,000368: # TFFTEM E: (-59+2NR) LCP Z Z SGN(SDELF) LEFT FOR ENTRY DISPLAY TFF ROUTINES
059148,000369: # M: (-55+2NR) LCP /ALFA SGN(SDELF) LEFT FOR ENTRY DISPLAY TFF ROUTINES
059149,000370: # NOTE: TFFTEM = PDL 36D AND WILL BE DESTROYED BY .:UNIT:.
059150,000371: # RMAG1 E:(-29) M:(-27) PDL 12 NOT TOUCHED.
059151,000372: # TFFQ1 E:(-16) M:(-15) PDL 14D
059152,000373: # TFFDELQ E:(-16) M:(-15) PDL 10D
059153,000374: # PUSHLOC AT PDL+0
059154,000375:
059155,000376: # ERASABLE INITIALIZATION REQUIRED:
059156,000377: # RONE E:(-29) M:(-27) M STATE VECTOR LEFT BY USER
059157,000378: # VONE' E:(+10) M:(+9) VN/SQRT(MU) LEFT BY TFF/CONIC
059158,000379: # RMAG1 E:(-29) M:(-27) PRESENT RADIUS, M LEFT BY TFFCONIC
059159,000380: # C(MPAC) E:(-29) M:(-27) RTERM, TERMINAL RADIUS LENGTH, M LEFT BY USER
059160,000381:
059161,000382: # THE FOLLOWING ARE STORED IN THE PUSH LIST AREA.
059162,000383: # TFF/RTMU E:(17) M:(14) 1/SQRT(MU) LEFT BY TFFCONIC.
059163,000384: # NRMAG E: (-29+NR) M RMAG, NORMALIZED LEFT BY TFFCONIC
059164,000385: # M: (-27+NR)
059165,000386: # X1 -NR, NORM COUNT LEFT BY TFFCONIC
059166,000387: # TFFNP E: (-38+2NR) M LCP, SEMI LATUS RECTUM, WEIGHT NR LEFT BY TFFCONIC
059167,000388: # M: (-36+2NR)
059168,000389: # TFFALFA E: (26-NR) 1/M ALFA, WEIGHT NR LEFT BY TFFCONIC
059169,000390: # M: (24-NR)
059170,000391: # TFFRTALF E:(10+NA) SQRT(ALFA), NORMALIZED LEFT BY TFFCONIC
059171,000392: # M:(9+NA)
059172,000393: # X2 -NA, NORM COUNT LEFT BY TFFCONIC
059173,000394: # TFF1/ALF E: (-22-2NA) SIGNED SEMIMAJ AXIS, WEIGHTED BY NA LEFT BY TFFCONIC
059174,000395: # M: (-20-2NA)
059175,000396:
059176,000397: # DEBRIS: QPRET, PDL+0 ... PDL+3
059177,000398: # RTERM E:(-29) M:(-27) RTERM, TERMINAL RADIUS LENGTH
059178,000399: # RAPO E:(-29) M:(-27) PDL 16D (=NRTERM)
059179,000400: # RPER E:(-29) M:(-27) PDL 14D (=TFFQ1)
059180,000401:
|
|
Page 1380 |
059182,000403: 27,3256 77614 CALCTPER SETGO # ENTER WITH RPER IN MPAC
059183,000404: 27,3257 03436 TFFSW
059184,000405: 27,3260 57263 +3
059185,000406: 27,3261 77614 CALCTFF CLEAR # ENTER WITH RTERM IN MPAC
059186,000407: 27,3262 03676 TFFSW
059187,000408: 27,3263 00023 +3 STORE RTERM # E: (-29) M: (-27)
059188,000409: 27,3264 77657 SL*
059189,000410: 27,3265 20201 0,1 # X1=-NR
059190,000411: 27,3266 00021 STORE NRTERM # RTERM E: (-29+NR) M: (-27+NR)
059191,000412: 27,3267 44205 DMP BDSU
059192,000413: 27,3270 00033 TFFALFA # ALFA E: (26-NR) M: (24-NR)
059193,000414: 27,3271 15326 TFF1/4
059194,000415: 27,3272 41206 PUSH DMP # (2-ALFA RTERM) (-3) TO PDL+0
059195,000416: 27,3273 00021 NRTERM # E: (-29+NR) M: (-27+NR)
059196,000417: 27,3274 53725 PDDL SR* # RTERM(2-ALFA RTERM) TO PDL+2
059197,000418: # E: (-32+NR) M: (-30+NR)
059198,000419: 27,3275 00035 TFFNP # LC P E:(-38+2NR) M:(-36+2NR)
059199,000420: 27,3276 20573 0 -6,1 # X1 = -NR
059200,000421: 27,3277 43276 DCOMP DAD # DUE TO SHIFTS, KEEP PRECISION FOR SQRT
059201,000422: # RTERM(2-ALFA RTERM) FROM PDL +2
059202,000423: # E: (-32+NR) M: (-30+NR)
059203,000424: 27,3300 77657 SR* # LEAVE E: (-32) M: (-30)
059204,000425: 27,3301 20601 0,1 # X1 = -NR
059205,000426: 27,3302 71214 BOFF DLOAD # CHECK TFF /TPER SWITCH
059206,000427: 27,3303 03756 TFFSW
059207,000428: 27,3304 57306 +2 # IF TFF, CONTINUE
059208,000429: 27,3305 15336 TFFZEROS # IF TPER, SET Q2 = 0
059209,000430: 27,3306 75440 +2 BMN SQRT # E: (-16) M: (-15)
059210,000431:
059211,000432: 27,3307 57441 MAXTFF1 # NO FREE FALL CONIC TO RTERM FROM HERE
059212,000433: # RESET PDL, SET TFF=POSMAX, AND EXIT.
059213,000434:
059214,000435: 27,3310 41076 DCOMP BOVB # RT IS ON INBOUND SIDE. ASSURE OVFIND=0
059215,000436: 27,3311 57544 TCDANZIG # ANY PORT IN A STORM.
059216,000437: 27,3312 24045 STOVL TFFTEM # Q2 E: (-16) M: (-15)
059217,000438: 27,3313 02362 VONE' # VN/SQRT(MU) E: (10) M: (9)
059218,000439: 27,3314 52441 DOT SL3
059219,000440: 27,3315 02320 RONE # SAVED RN. E: (-29) M: (-27)
059220,000441: 27,3316 00017 STORE TFFQ1 # Q1, SAVE FOR GONEPAST TEST.
059221,000442: # E: (-16) M: (-15)
059222,000443: 27,3317 44240 BMN BDSU
059223,000444: 27,3320 57341 INBOUND # USE ALTERNATE Z
059224,000445: 27,3321 00045 TFFTEM # Q2 E: (-16) M: (-15)
059225,000446:
059226,000447: # OUTBOUND Z CALC CONTINUES HERE
059227,000448:
059228,000449: 27,3322 14043 STODL TFFX # NUM=Q2-Q1 E: (-16) M: (-15)
059229,000450: 27,3323 00033 TFFALFA # ALFA E: (26-NR) M: (24-NR)
059230,000451: 27,3324 44205 DMP BDSU
|
|
Page 1381 |
059232,000453: 27,3325 00041 NRMAG # RMAG E: (-29+NR) M: (-27+NR)
059233,000454: # (2-RTERM ALFA) (-3) FROM PDL+0
059234,000455: 27,3326 51406 SAVEDEN PUSH ABS # DEN TO PDL+0 E: (-3) OR (-16)
059235,000456: # M: (-3) OR (-15)
059236,000457: 27,3327 40015 DAD BOV # INDETERMINANCY TEST
059237,000458: 27,3330 31763 LIM(-22) # =1.0-B(-22)
059238,000459: 27,3331 57352 TFFXTEST # GO IF DEN >/= B(-22)
059239,000460: 27,3332 65345 DLOAD PDDL # SET DEN=0 OTHERWISE
059240,000461: 27,3333 15336 TFFZEROS
059241,000462: # XCH ZERO WITH PDL+0
059242,000463: 27,3334 57545 DLOAD DCOMP
059243,000464: 27,3335 00033 TFFALFA # ALFA E: (26-NR) M: (24-NR)
059244,000465: 27,3336 71240 BMN DLOAD # FOR TPER: Z INDET AT DELE/2=0 AND 90.
059245,000466: 27,3337 57446 TFFEL1 # ASSUME 90, AND LEAVE 0 IN PDL: 1/Z=D/N
059246,000467:
059247,000468: # Z INDET. AT PERIGEE FOR PARAB OR HYPERB.
059248,000469: 27,3340 77616 DUMPTFF1 RVQ # RETURN TFF =0
059249,000470:
059250,000471: # INBOUND Z CALC CONTINUES HERE
059251,000472: 27,3341 77745 INBOUND DLOAD # RESET PDL+0
059252,000473: 27,3342 45345 DLOAD DSU # ALTERNATE Z CALC
059253,000474: 27,3343 00023 RTERM # E: (-29) M: (-27)
059254,000475: 27,3344 00015 RMAG1 # E: (-29) M: (-27)
059255,000476: 27,3345 14043 STODL TFFX # NUM=RTERM-RN E: (-29) M: (-27)
059256,000477: 27,3346 00045 TFFTEM # Q2 E: (-16) M: (-15)
059257,000478: 27,3347 52015 DAD GOTO
059258,000479: 27,3350 00017 TFFQ1 # Q1 E: (-16) M:(-15)
059259,000480: 27,3351 57326 SAVEDEN # DEN = Q2+Q1 E: (-16) M: (-15)
059260,000481:
059261,000482: 27,3352 65215 TFFXTEST DAD PDDL # (ABS(DEN) TO PDL+2)) E: (-3) OR (-16)
059262,000483: # M: (-3) OR (-15)
059263,000484: 27,3353 31765 DP(-22) # RESTORE ABS(DEN) TO MPAC
059264,000485: 27,3354 00043 TFFX # NUM E:(-16) OR (-29) M:(-15) OR (-27)
059265,000486: 27,3355 53605 DMP SR*
059266,000487: 27,3356 00031 TFFRTALF # SQRT(ALFA) E: (10+NA) M: (9+NA)
059267,000488: 27,3357 57201 0 -3,2 # X2=-NA
059268,000489: 27,3360 77671 DDV # C(MPAC) =NUM SQRT(ALFA) E:(-3) OR (-16)
059269,000490: # M:(-3) OR (-15)
059270,000491: # ABS(DEN) FROM PDL+2 E:(-3) OR (-16)
059271,000492: # M:(-3) OR (-15)
059272,000493: 27,3361 40145 DLOAD BOV # (THE DLOAD IS SHARED WITH TFFELL)
059273,000494: 27,3362 00043 TFFX # NUM E: (-16) OR (-29) M:(-15) OR (-27)
059274,000495: 27,3363 57444 TFFELL # USE EQN FOR DELE GEQ 90, LEQ -90
059275,000496:
059276,000497: # OTHERWISE, CONTINUE FOR GENERAL CONIC FOR TFF EQN
059277,000498:
059278,000499: 27,3364 45471 DDV STADR
059279,000500: # DEN FROM PDL+0 E: (-3) OR (-16)
059280,000501: # M: (-3) OR (-15)
059281,000502: 27,3365 77732 STORE TFFTEM # Z SAVE FOR SIGN OF SDELF.
|
|
Page 1382 |
059283,000504: # E: (-13) M: (-12)
059284,000505: 27,3366 63406 PUSH DSQ # Z TO PDL+0
059285,000506: 27,3367 41206 PUSH DMP # Z SQ TO PDL+2 E: (-26) M: (-24)
059286,000507: 27,3370 00035 TFFNP # LC P E: (-38+2NR) M: (-36+2NR)
059287,000508: 27,3371 75261 SL SIGN
059288,000509: 27,3372 20206 5
059289,000510: 27,3373 00045 TFFTEM # AFFIX SIGN FOR SDELF (ENTRY DISPLAY)
059290,000511: 27,3374 14045 STODL TFFTEM # P ZSQ E: (-59+2NR) M: (-55+2NR)
059291,000512: # (ARG IS USED IN TFF/TRIG)
059292,000513: # ZSQ FROM PDL+2 E: (-26) M: (-24)
059293,000514: 27,3375 41206 PUSH DMP # RESTORE PUSH LOC
059294,000515: 27,3376 00033 TFFALFA # ALFA E: (26-NR) M: (24-NR)
059295,000516: 27,3377 77657 SL*
059296,000517: 27,3400 20201 0,1 # X1=-NR
059297,000518: 27,3401 00043 STORE TFFX # X
059298,000519: 27,3402 41234 RTB DMP
059299,000520: 27,3403 57526 T(X) # POLY
059300,000521: # ZSQ FROM PDL+2 E: (-26) M: (-24)
059301,000522: 27,3404 44302 SR2 BDSU # 2 ZSQ T(X) E: (-29) M: (-27)
059302,000523: 27,3405 00023 RTERM # RTERM E: (-29) M: (-27)
059303,000524: 27,3406 41215 DAD DMP
059304,000525: 27,3407 00015 RMAG1 # E: (-29) M: (-27)
059305,000526: # Z FROM PDL+0 E: (-13) M: (-12)
059306,000527: 27,3410 51042 SR3 BPL # TFF SQRT(MU) E: (-45) M: (-42)
059307,000528: 27,3411 57432 ENDTFF # (NO PUSH UP)
059308,000529: 27,3412 75206 PUSH SIGN # TFF SQRT(MU) TO PDL+0
059309,000530: 27,3413 00017 TFFQ1 # Q1 FOR GONEPAST TEST
059310,000531: 27,3414 71244 BPL DLOAD # GONE PAST ?
059311,000532: 27,3415 57436 NEGTFF # YES. TFF < 0.
059312,000533: 27,3416 00027 TFF1/ALF # 1/ALFA E: (-22-2NA) M: (-20-2NA)
059313,000534: 27,3417 51076 DCOMP BPL # ALFA > 0 ?
059314,000535: 27,3420 57436 NEGTFF # NO. TFF IS NEGATIVE.
059315,000536:
059316,000537: # CORRECT FOR ORBITAL PERIOD.
059317,000538:
059318,000539: 27,3421 77676 DCOMP # YES. CORRECT FOR ORB PERIOD.
059319,000540: 27,3422 56205 DMP DDV
059320,000541: 27,3423 31761 PI/16 # 2 PI (-5)
059321,000542: 27,3424 00031 TFFRTALF # SQRT(ALFA) E: (10+NA) M: (9+NA)
059322,000543: 27,3425 53657 SL* SL*
059323,000544: 27,3426 57602 0 -4,2 # X2=-NA
059324,000545: 27,3427 57602 0 -4,2
059325,000546: 27,3430 43257 SL* DAD
059326,000547: 27,3431 57576 0,2
059327,000548: # TFF SQRT(MU) FROM PDL+0 E:(-45) M:(-42)
059328,000549: 27,3432 40005 ENDTFF DMP BOV # TFF SQRT(MU) IN MPAC E:(-45) M:(-42)
059329,000550: 27,3433 00037 TFF/RTMU # E: (17) M: (14)
059330,000551: 27,3434 57442 MAXTFF # SET POSMAX IF OVFL.
059331,000552:
059332,000553: 27,3435 77616 DUMPTFF2 RVQ # RETURN TFF (-28) CS IN MPAC.
059333,000554:
|
|
Page 1383 |
059335,000556: 27,3436 77745 NEGTFF DLOAD
059336,000557: # TFF SQRT(MU) FROM PDL+0, NEGATIVE.
059337,000558: 27,3437 77650 GOTO
059338,000559: 27,3440 57432 ENDTFF
059339,000560:
059340,000561: 27,3441 77745 MAXTFF1 DLOAD # RESET PDL
059341,000562: 27,3442 43545 MAXTFF DLOAD RVQ
059342,000563: 27,3443 15344 NEARONE
059343,000564:
059344,000565: # TIME OF FLIGHT ELLIPSE WHEN DEL (ECCENTRIC ANOM) GEQ 90 AND LEQ -90.
059345,000566:
059346,000567: # NUM FROM TFFX. E: (-16) OR (-29)
059347,000568: # M: (-15) OR (-27)
059348,000569: 27,3444 77712 TFFELL SL2 # NUM E:(-14) OR (-27) M:(-13) OR (-25)
059349,000570: 27,3445 41465 BDDV PUSH # TEMP SAVE D/N IN PDL+0
059350,000571: # DEN FROM PDL+0 E:(-3)/( 16) M:(-3)/(-15)
059351,000572: # N/D TO PDL+0 E: (11) M: (10)
059352,000573: 27,3446 45345 TFFEL1 DLOAD DSU # (ENTER WITH D/N=0 IN PDL+0)
059353,000574: 27,3447 00045 TFFTEM # Q2 E: (-16) M: (-15)
059354,000575: 27,3450 00017 TFFQ1 # Q1 E: (-16) M: (-15)
059355,000576: 27,3451 14013 STODL TFFDELQ # Q2-Q1 E: (-16) M: (-15)
059356,000577: # D/N FROM PDL+0
059357,000578: 27,3452 77626 STADR
059358,000579: 27,3453 77732 STORE TFFTEM # D/N E: (11) M: (10)
059359,000580: 27,3454 53605 DMP SL*
059360,000581: 27,3455 00027 TFF1/ALF # 1/ALFA E: (-22-2NA) M: (-20-2NA)
059361,000582: 27,3456 57576 0,2 # 1/ALFA Z E: (-11-NA) M: (-10-NA)
059362,000583: 27,3457 41206 PUSH DMP # TO PDL+0
059363,000584: 27,3460 00045 TFFTEM # 1/Z E: (11) M: (10)
059364,000585: 27,3461 41057 SL* BOVB
059365,000586: 27,3462 57576 0,2 # X2= -NA
059366,000587: 27,3463 45761 SIGNMPAC # IN CASE X= 1.0, CONTINUE
059367,000588: 27,3464 00043 STORE TFFX # X=1/ALFA ZSQ
059368,000589: 27,3465 41234 RTB DMP
059369,000590: 27,3466 57526 T(X) # POLY
059370,000591: 27,3467 00043 TFFX
059371,000592: 27,3470 45242 SR3 DSU
059372,000593: 27,3471 31767 DP2(-3)
059373,000594: 27,3472 41405 DMP PUSH # 2(X T(X)-1) /Z ALFA E: (-15-NA)
059374,000595: # M: (-14-NA)
059375,000596: # 1/ALFA Z FROM PDL+0 E: (-11-NA)
059376,000597: # M: (-10-NA)
059377,000598: 27,3473 41345 DLOAD DMP # GET SIGN FOR SDELF
059378,000599: 27,3474 00045 TFFTEM # 1/Z E: (11) M: (10)
059379,000600: 27,3475 00015 RMAG1 # E: (-29) M: (-27)
059380,000601: 27,3476 43312 SL2 DAD
059381,000602: 27,3477 00017 TFFQ1 # Q1 E: (-16) M: (-15)
059382,000603: 27,3500 14045 STODL TFFTEM # (Q1+R 1/Z) =SGN OF SDELF E:(-16) M:(-15
059383,000604: 27,3501 00035 TFFNP # LC P E: (-38+2NR) M: (-36+2NR)
059384,000605: 27,3502 53605 DMP SL* # CALC FOR ARG FOR TFF/TRIG.
|
|
Page 1384 |
059386,000607: 27,3503 00027 TFF1/ALF # 1/ALFA E:(-22-2NA) M:(-20-2NA)
059387,000608: 27,3504 57575 1,2 # X2=-NA
059388,000609: 27,3505 53765 SIGN SL*
059389,000610: 27,3506 00045 TFFTEM # AFFIX SIGN FOR SDELF
059390,000611: 27,3507 57576 0,2
059391,000612: 27,3510 14045 STODL TFFTEM # P/ALFA E:(-59+2NR) M:(-55+2NR)
059392,000613: # (ARG FOR USE IN TFF/TRIG)
059393,000614: 27,3511 00027 TFF1/ALF # 1/ALFA E:(-22-2NA) M:(-20-2NA)
059394,000615: 27,3512 41366 SQRT DMP
059395,000616: 27,3513 31761 PI/16 # PI (-4)
059396,000617: 27,3514 77615 DAD
059397,000618: # 2(XT(X)-1)/Z ALFA FROM PDL E:(-15-NA)
059398,000619: # M:(-14-NA)
059399,000620: 27,3515 45257 SL* DSU
059400,000621: 27,3516 57577 0 -1,2
059401,000622: 27,3517 00013 TFFDELQ # Q2-Q1 E: (-16) M: (-15)
059402,000623: 27,3520 53605 DMP SL*
059403,000624: 27,3521 00027 TFF1/ALF # 1/ALFA E:(-22-2NA) M:(-20-2NA)
059404,000625: 27,3522 57601 0 -3,2
059405,000626: 27,3523 52057 SL* GOTO
059406,000627: 27,3524 57602 0 -4,2
059407,000628: 27,3525 57432 ENDTFF # TFF SQRT(MU) IN MPAC E:(-45) M:(-42)
059408,000629:
|
|
Page 1385 |
059410,000631:
059411,000632: # PROGRAM NAME: T(X) DATE: 01.17.67
059412,000633: # MOD NO: 0 LOG SECTION: TIME OF FREE FALL
059413,000634: # MOD BY: RR BAIRNSFATHER
059414,000635:
059415,000636: # FUNCTIONAL DESCRIPTION: THE POLYNOMIAL T(X) IS USED BY TIME OF FLIGHT SUBROUTINES CALCTFF AND
059416,000637: # CALCTPER TO APPROXIMATE THE SERIES
059417,000638: # 2 3
059418,000639: # 1/3 -X/5 +X /7 -X /9 ...
059419,000640:
059420,000641: # WHERE X = ALFA Z Z IF ALFA Z Z LEQ 1
059421,000642: # X = 1/(ALFA Z Z) IF ALFA Z Z G 1
059422,000643:
059423,000644: # ALSO X IS NEG FOR HYPERBOLIC ORBITS
059424,000645: # X = 0 FOR PARABOLIC ORBITS
059425,000646: # X IS POSITIVE FOR ELLIPTIC ORBITS
059426,000647: # FOR FLIGHT 278, THE POLYNOMIAL T(X) IS FITTED OVER THE RANGE (0,+1) AND HAS A MAXIMUM
059427,000648: # DEVIATION FROM THE SERIES OF 2 E-5 (T(X) IS A CHEBYCHEV TYPE FIT AND WAS OBTAINED USING
059428,000649: # MAC PROGRAM AUTCURFIT294RRB AND IS VALID TO THE SAME TOLERANCE OVER THE RANGE (-.08,+1). )
059429,000650:
059430,000651: # CALLING SEQUENCE: RTB
059431,000652: # T(X)
059432,000653: # C(MPAC) = X
059433,000654:
059434,000655: # SUBROUTINES CALLED: NONE
059435,000656:
059436,000657: # NORMAL EXIT MODE: TC DANZIG
059437,000658:
059438,000659: # ALARMS: NONE
059439,000660:
059440,000661: # OUTPUT: C(MPAC) = T(X)
059441,000662:
059442,000663: # ERASABLE INITIALIZATION REQUIRED:
059443,000664: # C(MPAC) = X
059444,000665:
059445,000666: # DEBRIS: NONE
059446,000667:
059447,000668: 27,3526 07171 T(X) TC POLY
059448,000669: 27,3527 00004 DEC 4 B-14 # N-1
059449,000670: 27,3530 12525 12525 2DEC 3.333333333 E-1
059450,000671: 27,3532 71463 57703 2DEC* -1.999819135 E-1*
059451,000672: 27,3534 04423 17645 2DEC* 1.418148467 E-1*
059452,000673: 27,3536 74604 43667 2DEC* -1.01310997 E-1*
059453,000674: 27,3540 01626 37256 2DEC* 5.609004986 E-2*
059454,000675: 27,3542 77404 52071 2DEC* -1.536156925 E-2*
059455,000676:
059456,000677: 27,3544 06030 ENDT(X) TC DANZIG
059457,000678:
059458,000679: 27,3545 27,3544 TCDANZIG = ENDT(X)
059459,000680:
|
|
Page 1386 |
059461,000682: # TFF CONSTANTS
059462,000683:
059463,000684: 34,2000 SETLOC TOF-FF1
059464,000685: 34,2000 BANK
059465,000686:
059466,000687: 34,3756 COUNT* $$/TFF
059467,000688: # NOTE _ NOTE _ ADJUSTED MUE FOR NEAR EARTH TRAJ.
059468,000689:
|
|
(JL,2009-08-21) The Artemis072 page scan for the following three lines is garbled (lines 0660, 0661, 0662 in the original printout), due to what looks like overprinting in the original printout. I checked the generated octal (which is not overprinted) against the Comanche055 listing and they are both identical. |
059473,000694:
059474,000695: # MUE = 3.990815471 E10 M CUBE/CS SQ
059475,000696: # RTMUE = 1.997702549 E5 B-18* MODIFIED EARTH MU
059476,000697:
059477,000698: 34,3756 24775 30424 1/RTMU 2DEC* .5005750271 E-5 B17* # MODIFIED EARTH MU
059478,000699:
|
|
(JL,2009-08-21) End of garbled section. |
059480,000701:
059481,000702: # NOTE _ NOTE _ ADJUSTED MUE FOR NEAR EARTH TRAJ.
059482,000703: # MUM = 4.902778 E8 M CUBE/CS SQ
059483,000704: # RTMUM 2DEC* 2.21422176 E4 B-18*
059484,000705: 34,3760 06220 37553 PI/16 2DEC 3.141592653 B-4
059485,000706: 34,3762 37777 37700 LIM(-22) 2OCT 3777737700 # 1.0 -B(-22)
059486,000707: 34,3764 00000 00100 DP(-22) 2OCT 0000000100 # B(-22)
059487,000708: 34,3766 04000 00000 DP2(-3) 2DEC 1 B-3
059488,000709: 34,3770 02000 00000 DP2(-4) 2DEC 1 B-4 # 1/16
059489,000710:
059490,000711: # RPAD1 2DEC 6373338 B-29 M (-29) = 20 909 901.57 FT
059491,000712: 34,3772 22,3172 RPAD1 = RPAD
059492,000713: 34,3772 26,3343 NEARONE EQUALS DPPOSMAX
059493,000714: 34,3772 26,3335 TFFZEROS EQUALS HI6ZEROS
059494,000715: 34,3772 26,3325 TFF1/4 EQUALS HIDP1/4
End of include-file TIME_OF_FREE_FALL.agc. Parent file is MAIN.agc