Source Code

These source-code files were obtained by digitally photographing an Artemis 72 program listing from a private collection and then modifying pre-existing Comanche 55 (Apollo 11) source files to incorporate changes. 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, 1971

	THIS AGC PROGRAM SHALL ALSO BE REFERRED TO AS
				COLOSSUS 3
	THIS PROGRAM IS INTENDED FOR USE IN THE CM AS SPECIFIED
	IN REPORT R-577.  THIS PROGRAM WAS PREPARED UNDER OSR
	PROJECT 55-23890, SPONSORED BY THE MANNED SPACECRAFT
	CENTER OF THE NATIONAL AERONAUTICS AND SPACE
	ADMINISTRATION THROUGH CONTRACT NAS 9-4065 WITH THE
	CHARLES STARK DRAPER LABORATORY, MASSACHUSETTS INSTITUTE OF
	TECHNOLOGY, CAMBRIDGE, MASS.

058521,000002:                                                                                                 #  Copyright:     Public domain.
058522,000003:                                                                                                 #  Filename:      TIME_OF_FREE_FALL.agc
058523,000004:                                                                                                 #  Purpose:       Part of the source code for Artemis (i.e., Colossus 3),
058524,000005:                                                                                                 #                 build 072.  This is for the Command Module's (CM)
058525,000006:                                                                                                 #                 Apollo Guidance Computer (AGC), we believe for
058526,000007:                                                                                                 #                 Apollo 15-17.
058527,000008:                                                                                                 #  Assembler:     yaYUL
058528,000009:                                                                                                 #  Contact:       Jim Lawton <jim DOT lawton AT gmail DOT com>
058529,000010:                                                                                                 #  Website:       www.ibiblio.org/apollo/index.html
058530,000011:                                                                                                 #  Page scans:    www.ibiblio.org/apollo/ScansForConversion/Artemis072/
058531,000012:                                                                                                 #  Mod history:   2009-08-20 JL   Adapted from corresponding Comanche 055 file.
058532,000013:                                                                                                 #                 2009-08-21 JL   Added a note about a problem with the page scan on p.1386.
058533,000014:                                                                                                 #                 2010-02-20 RSB  Un-##'d this header.
058534,000015: 
058535,000016:                                                                                                 # # Page 1371
058536,000017: 
058537,000018:                                                                                                 #  THE TFF SUBROUTINES MAY BE USED IN EITHER EARTH OR MOON CENTERED COORDINATES. THE TFF ROUTINES NEVER
058538,000019:                                                                                                 #  KNOW WHICH ORIGIN APPLIES. IT IS THE USER WHO KNOWS, AND WHO SUPPLIES  RONE, VONE, AND 1/SQRT(MU)  AT THE
058539,000020:                                                                                                 #  APPROPRiATE SCALE LEVEL FOR THE PROPER PRIMARY BODY.
058540,000021: 
058541,000022:                                                                                                 #         EARTH ORIGIN    POSITION        -29     METERS
058542,000023:                                                                                                 #                         VELOCITY        -7      METERS/CENTISECOND
058543,000024:                                                                                                 #                         1/SQRT(MU)      +17     SQRT(CS SQ/METERS CUBED)
058544,000025: 
058545,000026:                                                                                                 #         MOON ORIGIN     POSITION        -27     METERS
058546,000027:                                                                                                 #                         VELOCITY        -5      METERS/CENTISECONDS
058547,000028:                                                                                                 #                         1/SQRT(MU)      +14     SQRT(CS SQ/METERS CUBED)
058548,000029: 
058549,000030:                                                                                                 #  ALL DATA PROVIDED TO AND RECEIVED FROM ANY TFF SUBROUTINE WILL BE AT ONE OF THE LEVELS ABOVE. IN ALL CASES,
058550,000031:                                                                                                 #  THE FREE FALL TIME IS RETURNED IN CENTISECONDS AT (-28). PROGRAM TFF/CONIC WILL GENERATE VONE/RTMU AND
058551,000032:                                                                                                 #  LEAVE IT IN VONE' AT (+10) IF EARTH ORIGIN AND (+9) IF MOON ORIGIN.
058552,000033: 
058553,000034:                                                                                                 #  THE USER MUST STORE THE STATE VECTOR IN  RONE, VONE, AND MU  IN THE FORM 1/SQRT(MU) IN TFF/RTMU
058554,000035:                                                                                                 #  AT THE PROPER SCALE BEFORE CALLING TFF/CONIC. SINCE  RONE, VONE  ARE IN THE EXTENDED VERB STORAGE AREA,
058555,000036:                                                                                                 #  THE USER MUST ALSO LOCK OUT THE EXTENDED VERBS, AND RELEASE THEM WHEN FINISHED.
058556,000037: 
058557,000038:                                                                                                 #  PROGRAMS CALC/TFF AND CALC/TPER ASSUME THAT THE TERMINAL RADIUS IS LESS THAN THE PRESENT
058558,000039:                                                                                                 #  RADIUS. THIS RESTRICTION CAN BE REMOVED BY A 15 W CODING CHANGE, BUT AT PRESENT IT IS NOT DEEMED NECESSARY.
058559,000040: 
058560,000041:                                                                                                 #  THE FOLLOWING ERASABLE QUANTITIES ARE USED BY THE TFF ROUTINES, AND ARE LOCATED IN THE PUSH LIST.
058561,000042: 
058562,000043:                                                                                                 #                                         BELOW   E:  IS USED FOR EARTH ORIGIN SCALE
058563,000044:                                                                                                 #                                                 M:  IS USED FOR MOON  ORIGIN SCALE
058564,000045: 
058565,000046:                                                                                                 # TFFSW           =       119D    # BIT1  0 = CALCTFF             1 = CALCTPER
058566,000047: ???????                        TFFDELQ           =        10D                                   #         Q2-Q1                   E: (-16)  M: (-15)
058567,000048: ???????                        RMAG1             =        12D                                   #         ABVAL(RN)  M            E: (-29)  M: (-27)
058568,000049:                                                                                                 # RPER            =       14D     #       PERIGEE RADIUS  M       E: (-29)  M: (-27)
058569,000050: ???????                        TFFQ1             =        14D                                   #         R.V / SQRT(MUE)         E: (-16)  M: (-15)
058570,000051:                                                                                                 # SDELF/2                 #       SIN(THETA) /2
058571,000052: ???????                        CDELF/2           =        14D                                   #         COS(THETA) /2
058572,000053:                                                                                                 # RAPO            =       16D     #       APOGEE RADIUS  M        E: (-29)  M: (-27)
058573,000054: ???????                        NRTERM            =        16D                                   #         TERMINAL RADIUS M       E: (-29+NR)
058574,000055:                                                                                                 #                                 M: (-27+NR)
058575,000056: ???????                        RTERM             =        18D                                   #         TERMINAL RADIUS  M      E: (-29)  M: (-27)
058576,000057: ???????                        TFFVSQ            =        20D                                   #         -(V SQUARED/MU)  1/M    E: (20)   M: (18)
058577,000058: ???????                        TFF1/ALF          =        22D                                   #         SEMI MAJ AXIS  M        E: (-22-2 NA)
058578,000059:                                                                                                 #                                 M: (-20-2 NA)
058579,000060: ???????                        TFFRTALF          =        24D                                   #         SQRT(ALFA)              E: (10+NA) M: (9+NA)
058580,000061: ???????                        TFFALFA           =        26D                                   #         ALFA  1/M               E: (26-NR) M: (24-NR)
058581,000062: ???????                        TFFNP             =        28D                                   #         SEMI LATUS RECTUM  M    E: (-38+2 NR)
058582,000063:                                                                                                 #                                 M: (-36+2 NR)
058583,000064: ???????                        TFF/RTMU          =        30D                                   #         1/SQRT(MU)              E: (17)   M: (14)
058584,000065: ???????                        NRMAG             =        32D                                   #         PRESENT RADIUS  M       E: (-29+NR)
058585,000066:                                                                                                 #                                 M: (-27+NR)
058586,000067: ???????                        TFFX              =        34D                                   
058587,000068: ???????                        TFFTEM            =        36D                                   #         TEMPORARY
058588,000069: 
058589,000070:                                                                                                 # # Page 1372
058590,000071:                                                                                                 #                 REGISTERS S1, S2 ARE UNTOUCED BY ANY TFF SUBROUTINE
058591,000072:                                                                                                 #                 INDEX REGISTERS X1, X2 ARE USED BY ALL TFF SUBROUTINES. THEY ARE ESTAB-
058592,000073:                                                                                                 #                 LISHED IN TFF/CONIC AND MUST BE PRESERVED BETWEEN CALLS TO SUBSEQUENT
058593,000074:                                                                                                 #                 SUBROUTINES.
058594,000075:                                                                                                 #                 -NR                             C(X1) = NORM COUNT OF RMAG
058595,000076:                                                                                                 #                 -NA                             C(X2) = NORM COUNT OF SQRT(ABS(ALFA))
058596,000077: 
058597,000078:                                                                                                 # # Page 1373
058598,000079: 
058599,000080:                                                                                                 #  SUBROUTINE NAME:       TFFCONIC                                        DATE:  01.29.67
058600,000081:                                                                                                 #  MOD NO:  0                                                             LOG SECTION:  TIME OF FREE FALL
058601,000082:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
058602,000083:                                                                                                 #  MOD NO:  1     MOD BY:  RR BAIRNSFATHER        DATE: 11 APR 67
058603,000084:                                                                                                 #  MOD NO:  2     MOD BY:  RR BAIRNSFATHER        DATE: 21 NOV 67         ADD MOON MU.
058604,000085:                                                                                                 #  MOD NO:  3     MOD BY:  RR BAIRNSFATHER        DATE: 21 MAR 68         ACCEPT DIFFERENT EARTH/MOON SCALES
058605,000086: 
058606,000087:                                                                                                 #  FUNCTIONAL DESCRIPTION:        THIS SUBROUTINE IS CALLED TO COMPUTE THOSE CONIC PARAMETERS REQUIRED BY THE TFF
058607,000088:                                                                                                 #         SUBROUTINES AND TO ESTABLISH THEM IN THE PUSH LIST AREA. THE PARAMETERS ARE LISTED UNDER OUTPUT.
058608,000089:                                                                                                 #         THE EQUATIONS ARE:
058609,000090:                                                                                                 #                 _   __ __
058610,000091:                                                                                                 #                 H = RN*VN                       ANGULAR MOMENTUM
058611,000092:                                                                                                 #                       _ _
058612,000093:                                                                                                 #                 LCP = H.H / MU                  SEMI LATUS RECTUM
058613,000094:                                                                                                 #                               __ __
058614,000095:                                                                                                 #                 ALFA = 2/RN - VN.VN / MU        RECIPROCAL SEMI MAJ AXIS, SIGNED
058615,000096: 
058616,000097:                                                                                                 #         AND ALFA IS POS FOR ELLIPTIC ORBITS
058617,000098:                                                                                                 #                      0  FOR PARABOLIC ORBITS
058618,000099:                                                                                                 #                     NEG FOR HYPERBOLIC ORBITS.
058619,000100:                                                                                                 #         SUBROUTINE ALSO COMPUTES AND SAVES RMAG.
058620,000101: 
058621,000102:                                                                                                 #  CALLING SEQUENCE:
058622,000103:                                                                                                 #         TFFCONIC  EXPECTS CALLER TO ENTER WITH CORRECT GRAVITATIONAL CONSTANT IN MPAC, IN THE FORM
058623,000104:                                                                                                 #         1/SQRT(MU). THE PROGRAM WILL SAVE IN  TFF/RTMU  . THE SCALE IS DETERMINED BY WHETHER EARTH OR MOON
058624,000105:                                                                                                 #         ORIGIN IS USED. THE CALLER MUST LOCK OUT THE EXTENDED VERBS BEFORE PROVIDING STATE VECTOR IN RONE,
058625,000106:                                                                                                 #         VONE  AT PROPER SCALE. THE EXTENDED VERBS MUST BE RESTORED WHEN THE CALLER IS FINISHED USING THE
058626,000107:                                                                                                 #         TFF ROUTINES.
058627,000108: 
058628,000109:                                                                                                 #         ENTRY POINT  TFFCONMU  EXPECTS THAT  TFF/RTMU  IS ALREADY LOADED.
058629,000110: 
058630,000111:                                                                                                 #         TO SPECIFY MU:  DLOAD   CALL                    IF MU ALREADY STORED:   CALL
058631,000112:                                                                                                 #                                 YOURMU          1/RTMU  E: (17)  M: (14)                TFFCONMU
058632,000113:                                                                                                 #                                 TFFCONIC
058633,000114: 
058634,000115:                                                                                                 #         PUSHLOC = PDL+0, ARBITRARY IF LEQ 18D
058635,000116: 
058636,000117:                                                                                                 #  SUBROUTINES CALLED:  NONE
058637,000118: 
058638,000119:                                                                                                 #  NORMAL EXIT MODES:  RVQ
058639,000120: 
058640,000121:                                                                                                 #  ALARMS:  NONE
058641,000122: 
058642,000123:                                                                                                 #  OUTPUT:        THE FOLLOWING ARE STORED IN THE PUSH LIST AREA.
058643,000124:                                                                                                 #                 RMAG1           E:(-29) M:(-27) M  RN, PRESENT RADIUS LENGTH.
058644,000125:                                                                                                 #                 NRMAG           E: (-29+NR)     M  RMAG, NORMALIZED
058645,000126:                                                                                                 #                                 M: (-27+NR)
058646,000127:                                                                                                 #                 X1                              -NR, NORM COUNT
058647,000128:                                                                                                 #                 TFFNP           E: (-38+2NR)    M  LCP, SEMI LATUS RECTUM, WEIGHTED BY NR.      FOR VGAMCALC.
058648,000129:                                                                                                 #                                 M: (-36+2NR)
058649,000130:                                                                                                 #                 TFF/RTMU        E:(17)  M:(14)  1/SQRT(MU)
058650,000131:                                                                                                 #                 TFFVSQ          E:(20)  M:(18)  1/M  -(V SQ/MU):  PRESENT VELOCITY, NORMLIZED.  FOR VGAMCALC
058651,000132:                                                                                                 #                 TFFALFA         E: (26-NR)      1/M  ALFA, WEIGHTED BY NR
058652,000133:                                                                                                 #                                 M: (24-NR)
058653,000134:                                                                                                 #                 TFFRTALF        E: (10+NA)      SQRT(ALFA), NORMALIZED
058654,000135:                                                                                                 #                                 M: (9+NA)
058655,000136:                                                                                                 # # Page 1374
058656,000137:                                                                                                 #                 X2                              -NA, NORM COUNT
058657,000138:                                                                                                 #                 TFF1/ALF        E: (-22-2NA)    SIGNED SEMI MAJ AXIS, WEIGHTED BY NA
058658,000139:                                                                                                 #                                 M: (-20-2NA)
058659,000140:                                                                                                 #                 PUSHLOC AT PDL+0
058660,000141:                                                                                                 #                 THE FOLLOWING IS STORED IN GENERAL ERASABLE
058661,000142:                                                                                                 #                 VONE'           E:(10)  M:(9)   V/RT(MU), NORMALIZED VELOCITY
058662,000143: 
058663,000144:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
058664,000145:                                                                                                 #                 RONE            E:(-29) M:(-27) M       STATE VECTOR    LEFT BY CALLER
058665,000146:                                                                                                 #                 VONE            E:(-7)  M:(-5)  M/CS    STATE VECTOR    LEFT BY CALLER
058666,000147:                                                                                                 #                 TFF/RTMU        E:(17)  M:(14)  1/RT(CS SQ/M CUBE)      IF ENTER VIA  TFFCONMU.
058667,000148: 
058668,000149:                                                                                                 #  DEBRIS:        QPRET           PDL+0 ... PDL+3
058669,000150: 
058670,000151: 27,2000                                          SETLOC   TOF-FF                                
058671,000152: 27,2000                                          BANK                                           
058672,000153: 
058673,000154: 27,3151                                          COUNT*   $$/TFF                                
058674,000155: 
058675,000156: 27,3151           00037        TFFCONIC          STORE    TFF/RTMU                              #  1/SQRT(MU)             E: (17)  M: (14)
058676,000157: 
058677,000158: 27,3152           53575        TFFCONMU          VLOAD    UNIT                                  #  COME HERE WITH TFFRTMU LOADED.
058678,000159: 27,3153           02320                                   RONE                                  #  SAVED RN.  M           E: (-29)  M: (-27)
058679,000160: 27,3154           77725                          PDDL                                           #  UR/2 TO PDL+0, +5
058680,000161: 27,3155           00045                                   36D                                   #  MAGNITUDE
058681,000162: 27,3156           00015                          STORE    RMAG1                                 #  M      E: (-29)  M: (-27)
058682,000163: 
058683,000164: 27,3157           77701                          NORM                                           
058684,000165: 27,3160           00047                                   X1                                    #  -NR
058685,000166: 27,3161           24041                          STOVL    NRMAG                                 #  RMAG  M        E: (-29+NR)  M: (-27+NR)
058686,000167: 27,3162           02326                                   VONE                                  #  SAVED VN.  M/CS        E: (-7)  M: (-5)
058687,000168: 27,3163           77761                          VXSC                                           
058688,000169: 27,3164           00037                                   TFF/RTMU                              #  E: (17)  M: (14)
058689,000170: 27,3165           02362                          STORE    VONE'                                 #  VN/SQRT(MU)  E:(10)  M:(9)
058690,000171: 
058691,000172: 27,3166           47361                          VXSC     VXV                                   
058692,000173: 27,3167           00041                                   NRMAG                                 #  E: (-29+NR)  M: (-27+NR)
058693,000174:                                                                                                 #  UR/2 FROM PDL
058694,000175: 27,3170           47572                          VSL1     VSQ                                   #  BEFORE:  E:(-19+NR)  M:(-18+NR)
058695,000176: 27,3171           14035                          STODL    TFFNP                                 #  LC P  M  E:(-38+2NR)  M:(-36+2NR)
058696,000177:                                                                                                 #  SAVE ALSO FOR VGAMCALC
058697,000178: 27,3172           15326                                   TFF1/4                                
058698,000179: 27,3173           63271                          DDV      PDVL                                  #  (2/RMAG)  1/M  E:(26-NR)  M:(24-NR)
058699,000180: 27,3174           00041                                   NRMAG                                 #  RMAG  M   E:(-29+NR)  M:(-27+NR)
058700,000181: 27,3175           02362                                   VONE'                                 #  SAVED VN.   E: (10)   M: (9)
058701,000182: 27,3176           57436                          VSQ      DCOMP                                 #  KEEP MPAC+2 HONEST FOR SQRT.
058702,000183: 27,3177           00025                          STORE    TFFVSQ                                #  -(V SQ/MU)  E:(20)  M:(18)
058703,000184:                                                                                                 #  SAVE FOR VGAMCALC
058704,000185: 27,3200           43257                          SR*      DAD                                   
058705,000186: 27,3201           20573                                   0          -6,1                       #  GET -VSQ/MU  E:(26-NR)  M:(24-NR)
058706,000187:                                                                                                 # # Page 1375
058707,000188: 27,3202           77626                          STADR                                          
058708,000189:                                                                                                 #  2/RMAG  FROM PDL+2
058709,000190: 27,3203           77744                          STORE    TFFALFA                               #  ALFA  1/M  E:(26-NR)  M:(24-NR)
058710,000191: 27,3204           41457                          SL*      PUSH                                  #  TEMP SAVE ALFA  E:(20)  M:(18)
058711,000192: 27,3205           20173                                   0          -6,1                       
058712,000193: 27,3206           75446                          ABS      SQRT                                  #  E:(10)  M:(9)
058713,000194: 27,3207           77701                          NORM                                           
058714,000195: 27,3210           00050                                   X2                                    #  X2 = -NA
058715,000196: 27,3211           00031                          STORE    TFFRTALF                              #  SQRT( ABS(ALFA) )  E:(10+NA)  M:(9+NA)
058716,000197: 27,3212           75316                          DSQ      SIGN                                  #  NOT SO ACCURATE, BUT OK
058717,000198:                                                                                                 #  ALFA FROM PDL+2  E:(20)  M:(18)
058718,000199: 27,3213           55254                          BZE      BDDV                                  #  SET 1/ALFA =0, TO SHOW SMALL ALFA
058719,000200: 27,3214           57216                                   +2                                    
058720,000201: 27,3215           15326                                   TFF1/4                                
058721,000202: 27,3216           00027                 +2       STORE    TFF1/ALF                              #  1/ALFA  E:(-22-2 NA)  M:(-20-2 NA)
058722,000203: 27,3217           77616        DUMPCNIC          RVQ                                            
058723,000204: 
058724,000205:                                                                                                 #                         39 W
058725,000206: 
058726,000207:                                                                                                 # # Page 1376
058727,000208: 
058728,000209:                                                                                                 #  SUBROUTINE NAME:       TFFRP/RA                                        DATE: 01.17.67
058729,000210:                                                                                                 #  MOD NO:  0                                                             LOG SECTION:  TIME OF FREE FALL
058730,000211:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
058731,000212:                                                                                                 #  MOD NO:  1     MOD BY:  RR BAIRNSFATHER        DATE: 11 APR 67
058732,000213:                                                                                                 #  MOD NO:  2     MOD BY:  RR BAIRNSFATHER        DATE: 21 MAR 68         ACCEPT DIFFERENT EARTH/MOON SCALES
058733,000214:                                                                                                 #                                                                         ALSO IMPROVE ACCURACY OF RAPO.
058734,000215: 
058735,000216:                                                                                                 #  FUNCTIONAL DESCRIPTION:        USED BY CALCTPER AND TFF DISPLAYS TO CALCULATE PERIGEE RADIUS AND ALSO
058736,000217:                                                                                                 #         APOGEE RADIUS FOR A GENERAL CONIC.
058737,000218:                                                                                                 #         PROGRAM GIVES PERIGEE RADIUS AS         APOGEE RADIUS IS GIVEN BY
058738,000219:                                                                                                 #                 RP = P /(1+E)                           RA = (1+E) / ALFA
058739,000220:                                                                                                 #         WHERE    2
058740,000221:                                                                                                 #                 E  = 1 - P ALFA
058741,000222:                                                                                                 #         IF RA IS NEGATIVE OR SHOWS DIVIDE OVERFLOW, THEN RA = POSMAX BECAUSE
058742,000223:                                                                                                 #                 1. APOGEE RADIUS IS NOT MEANINGFUL FOR HYPERBOLA
058743,000224:                                                                                                 #                 2. APOGEE RADIUS IS NOT DEFINED FOR PARABOLA
058744,000225:                                                                                                 #                 3. APOGEE RADIUS EXCEEDS THE SCALING FOR ELLIPSE.
058745,000226: 
058746,000227:                                                                                                 #         THIS SUBROUTINE REQUIRED THE SIGNED RECIPROCAL SEMI MAJ AXIS, ALFA, AND SEMI LATUS RECTUM AS DATA.
058747,000228: 
058748,000229:                                                                                                 #  CALLING SEQUENCE:      CALL
058749,000230:                                                                                                 #                                 TFFRP/RA
058750,000231: 
058751,000232:                                                                                                 #         PUSHLOC = PDL+0, ARBITRARY IF LEQ 10D
058752,000233:                                                                                                 #         C(MPAC) UNSPECIFIED
058753,000234: 
058754,000235:                                                                                                 #  SUBROUTINES CALLED:    NONE
058755,000236: 
058756,000237:                                                                                                 #  NORMAL EXIT MODE:      RVQ
058757,000238:                                                                                                 #         IF ELLIPSE, WITHIN NORMAL SCALING, RAPO IS CORRECT.
058758,000239:                                                                                                 #         OTHERWISE, RAPO = POSMAX.
058759,000240: 
058760,000241:                                                                                                 #  ALARMS:        NONE
058761,000242: 
058762,000243:                                                                                                 #  OUTPUT:        STORED IN PUSH LIST AREA. SCALE OF OUTPUT AGREES WITH DATA SUPPLIED TO  TFF/CONIC.
058763,000244:                                                                                                 #         RPER    E:(-29)  M:(-27)        M       PERIGEE RADIUS                  DESTROYED BY CALCTFF/CALCTPER, TFFTRIG.
058764,000245:                                                                                                 #         RAPO    E:(-29)  M:(-27)        M       APOGEE RADIUS                   WILL BE DESTROYED BY CALCTFF/CALCTPER
058765,000246:                                                                                                 #         PUSHLOC AT PDL+0
058766,000247: 
058767,000248:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
058768,000249:                                                                                                 #         TFFALFA E:(26-NR)               M       1/SEMI MAJ AXIS                 LEFT BY TFFCONIC
058769,000250:                                                                                                 #                 M:(24-NR)
058770,000251:                                                                                                 #         TFFNP   E:(-38+2NR)             M       LC P, SEMI LATUS RECTUM         LEFT BY TFFCONIC
058771,000252:                                                                                                 #                 M:(-36+2NR)
058772,000253:                                                                                                 #         X1                              -NR, NORM COUNT OF RMAG                 LEFT BY TFFCONIC
058773,000254:                                                                                                 #         X2                              -NA, NORM COUNT OF ALFA                 LEFT BY TFFCONIC
058774,000255: 
058775,000256:                                                                                                 #  DEBRIS:        QPRET, PDL+0 ... PDL+1
058776,000257: 
058777,000258:                                                                                                 # # Page 1377
058778,000259: 27,3220                        RAPO              =        16D                                   #  APOGEE RADIUS  M  E:(-29)  M:(-27)
058779,000260: 27,3220                        RPER              =        14D                                   #  PERIGEE RADIUS  M  E:(-29)  M:(-27)
058780,000261: 
058781,000262: 27,3220           41345        TFFRP/RA          DLOAD    DMP                                   
058782,000263: 27,3221           00033                                   TFFALFA                               #  ALFA  1/M      E:(26-NR)  M:(24-NR)
058783,000264: 27,3222           00035                                   TFFNP                                 #  LC P  M        E:(-38+2NR)  M:(-36+2NR)
058784,000265: 27,3223           57457                          SR*      DCOMP                                 #  ALFA P (-12+NR)
058785,000266: 27,3224           20571                                   0          -8D,1                      #  ALFA P (-4)
058786,000267: 27,3225           51415                          DAD      ABS                                   #  (DCOMP GIVES VALID TP RESULT FOR SQRT)
058787,000268:                                                                                                 #  (ABS PROTECTS SQRT IF E IS VERY NEAR 0)
058788,000269: 27,3226           31771                                   DP2(-4)                               
058789,000270: 27,3227           43366                          SQRT     DAD                                   #  E SQ = (1- P ALFA)     (-4)
058790,000271: 27,3230           15326                                   TFF1/4                                
058791,000272: 27,3231           55206                          PUSH     BDDV                                  #  (1+E)  (-2)  TO PDL+0
058792,000273: 27,3232           00035                                   TFFNP                                 #  LCP  M         E:(-38+2NR)  M:(-36+2NR)
058793,000274: 27,3233           53657                          SR*      SR*                                   #  (DOES SR THEN SL TO AVOID OVFL)
058794,000275: 27,3234           20601                                   0,1                                   #  X1=-NR
058795,000276: 27,3235           20572                                   0          -7,1                       #  (EFFECTIVE SL)
058796,000277: 27,3236           14017                          STODL    RPER                                  #  PERIGEE RADIUS  M      E:(-29)  M:(-27)
058797,000278:                                                                                                 #  (1+E)  (-2)  FROM PDL+0
058798,000279: 27,3237           41005                          DMP      BOVB                                  
058799,000280: 27,3240           00027                                   TFF1/ALF                              #  E:(-22-2NA)  M:(-20-2NA)
058800,000281: 27,3241           57544                                   TCDANZIG                              #  CLEAR OVFIND, IF ON.
058801,000282: 27,3242           53654                          BZE      SL*                                   
058802,000283: 27,3243           57252                                   MAXRA                                 #  SET POSMAX IF ALFA=0
058803,000284: 27,3244           57603                                   0          -5,2                       #  -5+NA
058804,000285: 27,3245           40057                          SL*      BOV                                   
058805,000286: 27,3246           57576                                   0,2                                   
058806,000287: 27,3247           57252                                   MAXRA                                 #  SET POSMAX IF OVFL.
058807,000288: 27,3250           77644                          BPL                                            #  CONTINUE WITH VALID RAPO.
058808,000289: 27,3251           57254                                   +3                                    
058809,000290: 27,3252           77745        MAXRA             DLOAD                                          #  RAPO CALC IS NOT VALID. SET RAPO =
058810,000291: 27,3253           15344                                   NEARONE                               #  POSMAX AS A TAG.
058811,000292: 27,3254           00021                 +3       STORE    RAPO                                  #  APOGEE RADIUS  M  E:(-29)  M:(-27)
058812,000293: 27,3255           77616        DUMPRPRA          RVQ                                            
058813,000294: 
058814,000295:                                                                                                 #                         30 W
058815,000296: 
058816,000297:                                                                                                 # # Page 1378
058817,000298: 
058818,000299:                                                                                                 #  SUBROUTINE NAME:       CALCTPER / CALCTFF                              DATE:  01.29.67
058819,000300:                                                                                                 #  MOD NO:  0                                                             LOG SECTION:  TIME OF FREE FALL
058820,000301:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
058821,000302:                                                                                                 #  MOD NO:  1     MOD BY:  RR BAIRNSFATHER        DATE: 21 MAR 67
058822,000303:                                                                                                 #  MOD NO:  2     MOD BY:  RR BAIRNSFATHER        DATE: 14 APR 67
058823,000304:                                                                                                 #  MOD BY:  3     MOD BY:  RR BAIRNSFATHER        DATE: 8 JUL 67          NEAR EARTH MUE AND NEG TFF (GONEPAST)
058824,000305:                                                                                                 #  MOD BY:  4     MOD BY:  RR BAIRNSFATHER        DATE: 21 NOV 67         ADD VARIABLE MU.
058825,000306:                                                                                                 #  MOD BY:  5     MOD BY:  RR BAIRNSFATHER        DATE: 21 MAR 68         ACCEPT DIFFERENT EARTH/MOON SCALES
058826,000307: 
058827,000308:                                                                                                 #  FUNCTIONAL DESCRIPTION:        PROGRAM CALCULATES THE FREE-FALL TIME OF FLIGHT FROM PRESENT POSITION  RN  AND
058828,000309:                                                                                                 #         VELOCITY  VN  TO A RADIUS LENGTH SPECIFIED BY  RTERM  , SUPPLIED BY THE USER. THE POSITION VECTOR
058829,000310:                                                                                                 #         RN  MAY BE ON EITHER SIDE OF THE CONIC, BUT  RTERM  IS CONSIDERED ON THE INBOUND SIDE.
058830,000311:                                                                                                 #         THE EQUATIONS ARE:
058831,000312: 
058832,000313:                                                                                                 #                 Q2 = -SQRT(RTERM (2-RTERM ALFA) - LCP)  (INBOUND SIDE)  LEQ +- LCE/SQRT(ALFA)
058833,000314:                                                                                                 #                      _  _
058834,000315:                                                                                                 #                 Q1 = RN.VN / SQRT(MU)                                   LEQ +- LCE/SQRT(ALFA)
058835,000316: 
058836,000317:                                                                                                 #                 Z = NUM / DEN                                           LEQ +- 1/SQRT(ALFA)
058837,000318: 
058838,000319:                                                                                                 #         WHERE, IF INBOUND
058839,000320:                                                                                                 #                 NUM = RTERM - RN                                        LEQ +- 2 LCE/ALFA
058840,000321:                                                                                                 #                 DEN = Q2+Q1                                             LEQ +- 2 LCE/SQRT(ALFA)
058841,000322: 
058842,000323:                                                                                                 #         AND, IF OUTBOUND
058843,000324:                                                                                                 #                 NUM = Q2-Q1                                             LEQ +- 2 LCE/SQRT(ALFA)
058844,000325:                                                                                                 #                 DEN = 2 - ALFA (RTERM + RN) .                           LEQ +- 2 LCE
058845,000326: 
058846,000327:                                                                                                 #         IF      ALFA ZZ < 1.0           (FOR ALL CONICS EXCEPT ELLIPSES HAVING ABS(DEL ECC ANOM) G 90 DEG)
058847,000328: 
058848,000329:                                                                                                 #         THEN    X = ALFA Z Z
058849,000330:                                                                                                 #         AND     TFF = (RTERM +RN -2 ZZ T(X) ) Z/SQRT(MU)
058850,000331: 
058851,000332:                                                                                                 #                 EXCEPT  IF ALFA PNZ, AND IF TFF NEG,
058852,000333:                                                                                                 #                 THEN    TFF = 2 PI /(ALFA SQRT(ALFA)) + TFF
058853,000334: 
058854,000335:                                                                                                 #         OR      IF ALFA ZZ GEQ 1.0      (FOR ELLIPSES HAVING ABS(DEL ECC ANOM) GEQ 90 DEG)
058855,000336: 
058856,000337:                                                                                                 #         THEN    X = 1/ALFA Z Z
058857,000338:                                                                                                 #         AND     TFF = (PI/SQRT(ALFA) -Q2 +Q1 +2(X T(X) -1) /ALFA Z) /ALFA SQRT(MU)
058858,000339: 
058859,000340:                                                                                                 #         WHERE   T(X) IS A POLYNOMIAL APPROXIMATION TO THE SERIES
058860,000341:                                                                                                 #                              2      3             2
058861,000342:                                                                                                 #                 1/3 - X/5 + X /7 - X /9 ...     (X  < 1.0)
058862,000343: 
058863,000344:                                                                                                 #  CALLING SEQUENCE:      TIME TO RTERM                   TIME TO PERIGEE
058864,000345:                                                                                                 #                         CALL                            CALL
058865,000346:                                                                                                 #                                 CALCTFF                         CALCTPER
058866,000347:                                                                                                 #                         C(MPAC) = TERMNL RAD M          C(MPAC) = PERIGEE RAD M
058867,000348:                                                                                                 #         FOR EITHER, E: (-29)    M: (-27)
058868,000349:                                                                                                 #         FOR EITHER, PUSHLOC = PDL+0 , ARBITRARY IF LEQ 8D.
058869,000350:                                                                                                 # # Page 1379
058870,000351: 
058871,000352:                                                                                                 #  SUBROUTINES CALLED:    T(X), VIA RTB
058872,000353: 
058873,000354:                                                                                                 #  NORMAL EXIT MODE:      RVQ
058874,000355:                                                                                                 #         HOWEVER, PROGRAM EXITS WITH ONE OF THE FOLLOWING VALUES FOR TFF (-28) CS IN MPAC. USER MUST STORE.
058875,000356:                                                                                                 #                 A. TFF = FLIGHT TIME. NORMAL CASE FOR POSITIVE FLIGHT TIME LESS THAN ONE ORBITAL PERIOD.
058876,000357:                                                                                                 #                 B. (THIS OPTION IS NO LONGER USED.)
058877,000358:                                                                                                 #                 C. TFF = POSMAX.  THIS INDICATES THAT THE CONIC FROM THE PRESENT POSITION WILL NOT RETURN TO
058878,000359:                                                                                                 #                    THE SPECIFIED ALTITUDE.  ALSO INDICATES OUTBOUND PARABOLA OR HYPERBOLA.
058879,000360: 
058880,000361:                                                                                                 #  OUTPUT:        C(MPAC)         (-28) CS        TIME OF FLIGHT, OR TIME TO PERIGEE
058881,000362:                                                                                                 #                 TFFX            (0)             X,                                      LEFT FOR ENTRY DISPLAY TFF ROUTINES
058882,000363:                                                                                                 #                 NRTERM          E: (-29+NR)     M  RTERM, WEIGHTED BY NR                LEFT FOR ENTRY DISPLAY TFF ROUTINES
058883,000364:                                                                                                 #                                 M: (-27+NR)
058884,000365:                                                                                                 #                 TFFTEM          E: (-59+2NR)    LCP Z Z SGN(SDELF)                      LEFT FOR ENTRY DISPLAY TFF ROUTINES
058885,000366:                                                                                                 #                                 M: (-55+2NR)    LCP /ALFA SGN(SDELF)                    LEFT FOR ENTRY DISPLAY TFF ROUTINES
058886,000367:                                                                                                 #                 NOTE:   TFFTEM = PDL 36D AND WILL BE DESTROYED BY .'UNIT'.
058887,000368:                                                                                                 #                 RMAG1           E:(-29) M:(-27) PDL 12 NOT TOUCHED.
058888,000369:                                                                                                 #                 TFFQ1           E:(-16) M:(-15) PDL 14D
058889,000370:                                                                                                 #                 TFFDELQ         E:(-16) M:(-15) PDL 10D
058890,000371:                                                                                                 #                 PUSHLOC AT PDL+0
058891,000372: 
058892,000373:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
058893,000374:                                                                                                 #                 RONE            E:(-29) M:(-27) M  STATE VECTOR                         LEFT BY USER
058894,000375:                                                                                                 #                 VONE'           E:(+10) M:(+9)  VN/SQRT(NU)                             LEFT BY TFF/CONIC
058895,000376:                                                                                                 #                 RMAG1           E:(-29) M:(-27) PRESENT RADIUS, M                       LEFT BY TFFCONIC
058896,000377:                                                                                                 #                 C(MPAC)         E:(-29) M:(-27) RTERM, TERMINAL RADIUS LENGTH, M        LEFT BY USER
058897,000378: 
058898,000379:                                                                                                 #                 THE FOLLOWING ARE STORED IN THE PUSH LIST AREA.
058899,000380:                                                                                                 #                 TFF/RTMU        E:(17) M:(14)   1/SQRT(MU)                              LEFT BY TFFCONIC.
058900,000381:                                                                                                 #                 NRMAG           E: (-29+NR)     M  RMAG, NORMALIZED                     LEFT BY TFFCONIC
058901,000382:                                                                                                 #                                 M: (-27+NR)
058902,000383:                                                                                                 #                 X1                              -NR, NORM COUNT                         LEFT BY TFFCONIC
058903,000384:                                                                                                 #                 TFFNP           E: (-38+2NR)    M  LCP, SEMI LATUS RECTUM, WEIGHT NR    LEFT BY TFFCONIC
058904,000385:                                                                                                 #                                 M: (-36+2N4)
058905,000386:                                                                                                 #                 TFFALFA         E: (26-NR)      1/M  ALFA, WEIGHT NR                    LEFT BY TFFCONIC
058906,000387:                                                                                                 #                                 M: (24-NR)
058907,000388:                                                                                                 #                 TFFRTALF        E:(10+NA)       SQRT(ALFA), NORMALIZED                  LEFT BY TFFCONIC
058908,000389:                                                                                                 #                                 M:(9+NA)
058909,000390:                                                                                                 #                 X2                              -NA, NORM COUNT                         LEFT BY TFFCONIC
058910,000391:                                                                                                 #                 TFF1/ALF        E: (-22-2NA)    SIGNED SEMIMAJ AXIS, WEIGHTED BY NA     LEFT BY TFFCONIC
058911,000392:                                                                                                 #                                 M: (-20-2NA)
058912,000393: 
058913,000394:                                                                                                 #  DEBRIS:        QPRET, PDL+0 ... PDL+3
058914,000395:                                                                                                 #                 RTERM           E:(-29) M(-27)  RTERM, TERMINAL RADIUS LENGTH
058915,000396:                                                                                                 #                 RAPO            E:(-29) M(-27)  PDL 16D (=NRTERM)
058916,000397:                                                                                                 #                 RPER            E:(-29) M(-27)  PDL 14D (=TFFQ1)
058917,000398: 
058918,000399:                                                                                                 # # Page 1380
058919,000400: 27,3256           77614        CALCTPER          SETGO                                          #  ENTER WITH RPER IN MPAC
058920,000401: 27,3257           03436                                   TFFSW                                 
058921,000402: 27,3260           57263                                   +3                                    
058922,000403: 27,3261           77614        CALCTFF           CLEAR                                          #  ENTER WITH RTERM IN MPAC
058923,000404: 27,3262           03676                                   TFFSW                                 
058924,000405: 27,3263           00023                 +3       STORE    RTERM                                 #  E: (-29)  M: (-27)
058925,000406: 27,3264           77657                          SL*                                            
058926,000407: 27,3265           20201                                   0,1                                   #  X1=-NR
058927,000408: 27,3266           00021                          STORE    NRTERM                                #  RTERM  E: (-29+NR)  M: (-27+NR)
058928,000409: 27,3267           44205                          DMP      BDSU                                  
058929,000410: 27,3270           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
058930,000411: 27,3271           15326                                   TFF1/4                                
058931,000412: 27,3272           41206                          PUSH     DMP                                   #  (2-ALFA RTERM)  (-3)  TO PDL+0
058932,000413: 27,3273           00021                                   NRTERM                                #  E: (-29+NR)  M: (-27+NR)
058933,000414: 27,3274           53725                          PDDL     SR*                                   #  RTERM(2-ALFA RTERM) TO PDL+2
058934,000415:                                                                                                 #  E: (-32+NR)  M: (-30+NR)
058935,000416: 27,3275           00035                                   TFFNP                                 #  LC P  E:(-38+2NR)  M:(-36+2NR)
058936,000417: 27,3276           20573                                   0          -6,1                       #  X1 = -NR
058937,000418: 27,3277           43276                          DCOMP    DAD                                   #  DUE TO SHIFTS, KEEP PRECISION FOR SQRT
058938,000419:                                                                                                 #  RTERM(2-ALFA RTERM) FROM PDL +2
058939,000420:                                                                                                 #  E: (-32+NR)  M: (-30+NR)
058940,000421: 27,3300           77657                          SR*                                            #  LEAVE  E: (-32)  M: (-30)
058941,000422: 27,3301           20601                                   0,1                                   #  X1 = -NR
058942,000423: 27,3302           71214                          BOFF     DLOAD                                 #  CHECK TFF /TPER SWITCH
058943,000424: 27,3303           03756                                   TFFSW                                 
058944,000425: 27,3304           57306                                   +2                                    #  IF TFF, CONTINUE
058945,000426: 27,3305           15336                                   TFFZEROS                              #  IF TPER, SET Q2 = 0
058946,000427: 27,3306           75440                 +2       BMN      SQRT                                  #  E: (-16)  M: (-15)
058947,000428: 
058948,000429: 27,3307           57441                                   MAXTFF1                               #  NO FREE FALL CONIC TO RTERM FROM HERE
058949,000430:                                                                                                 #  RESET PDL, SET TFF=POSMAX, AND EXIT.
058950,000431: 
058951,000432: 27,3310           41076                          DCOMP    BOVB                                  #  RT IS ON INBOUND SIDE.  ASSURE OVFIND=0
058952,000433: 27,3311           57544                                   TCDANZIG                              #  ANY PORT IN A STORM.
058953,000434: 27,3312           24045                          STOVL    TFFTEM                                #  Q2  E: (-16)  M: (-15)
058954,000435: 27,3313           02362                                   VONE'                                 #  VN/SQRT(MU)  E: (10)  M: (9)
058955,000436: 27,3314           52441                          DOT      SL3                                   
058956,000437: 27,3315           02320                                   RONE                                  #  SAVED RN.  E: (-29)  M: (-27)
058957,000438: 27,3316           00017                          STORE    TFFQ1                                 #  Q1, SAVE FOR GONEPAST TEST.
058958,000439:                                                                                                 #  E: (-16)  M: (-15)
058959,000440: 27,3317           44240                          BMN      BDSU                                  
058960,000441: 27,3320           57341                                   INBOUND                               #  USE ALTERNATE Z
058961,000442: 27,3321           00045                                   TFFTEM                                #  Q2  E: (-16)  M: (-15)
058962,000443: 
058963,000444:                                                                                                 #  OUTBOUND Z CALC CONTINUES HERE
058964,000445: 
058965,000446: 27,3322           14043                          STODL    TFFX                                  #  NUM=Q2-Q1  E: (-16)  M: (-15)
058966,000447: 27,3323           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
058967,000448: 27,3324           44205                          DMP      BDSU                                  
058968,000449:                                                                                                 # # Page 1381
058969,000450: 27,3325           00041                                   NRMAG                                 #  RMAG  E: (-29+NR) M: (-27+NR)
058970,000451:                                                                                                 #  (2-RTERM ALFA)  (-3) FROM PDL+0
058971,000452: 27,3326           51406        SAVEDEN           PUSH     ABS                                   #  DEN TO PDL+0   E: (-3) OR (-16)
058972,000453:                                                                                                 #                M: (-3) OR (-15)
058973,000454: 27,3327           40015                          DAD      BOV                                   #  INDETERMINANCY TEST
058974,000455: 27,3330           31763                                   LIM(-22)                              #  =1.0-B(-22)
058975,000456: 27,3331           57352                                   TFFXTEST                              #  GO IF DEN >/= B(-22)
058976,000457: 27,3332           65345                          DLOAD    PDDL                                  #  SET DEN=0 OTHERWISE
058977,000458: 27,3333           15336                                   TFFZEROS                              
058978,000459:                                                                                                 #  XCH ZERO WITH PDL+0
058979,000460: 27,3334           57545                          DLOAD    DCOMP                                 
058980,000461: 27,3335           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
058981,000462: 27,3336           71240                          BMN      DLOAD                                 #  FOR TPER:  Z INDET AT DELE/2=0 AND 90.
058982,000463: 27,3337           57446                                   TFFEL1                                #  ASSUME 90, AND LEAVE 0 IN PDL: 1/Z=D/N
058983,000464: 
058984,000465:                                                                                                 #  Z INDET. AT PERIGEE FOR PARAB OR HYPERB.
058985,000466: 27,3340           77616        DUMPTFF1          RVQ                                            #  RETURN TFF =0
058986,000467: 
058987,000468:                                                                                                 #  INBOUND Z CALC CONTINUES HERE
058988,000469: 27,3341           77745        INBOUND           DLOAD                                          #  RESET PDL+0
058989,000470: 27,3342           45345                          DLOAD    DSU                                   #  ALTERNATE Z CALC
058990,000471: 27,3343           00023                                   RTERM                                 #  E: (-29)  M: (-27)
058991,000472: 27,3344           00015                                   RMAG1                                 #  E: (-29)  M: (-27)
058992,000473: 27,3345           14043                          STODL    TFFX                                  #  NUM=RTERM-RN  E: (-29)  M: (-27)
058993,000474: 27,3346           00045                                   TFFTEM                                #  Q2  E: (-16)  M: (-15)
058994,000475: 27,3347           52015                          DAD      GOTO                                  
058995,000476: 27,3350           00017                                   TFFQ1                                 #  Q1  E: (-16)  M:(-15)
058996,000477: 27,3351           57326                                   SAVEDEN                               #  DEN = Q2+Q1  E: (-16)  M: (-15)
058997,000478: 
058998,000479: 27,3352           65215        TFFXTEST          DAD      PDDL                                  #   (ABS(DEN) TO PDL+2))  E: (-3) OR (-16)
058999,000480:                                                                                                 #                         M: (-3) OR (-15)
059000,000481: 27,3353           31765                                   DP(-22)                               #  RESTORE ABS(DEN) TO MPAC
059001,000482: 27,3354           00043                                   TFFX                                  #  NUM  E:(-16) OR (-29)  M:(-15) OR (-27)
059002,000483: 27,3355           53605                          DMP      SR*                                   
059003,000484: 27,3356           00031                                   TFFRTALF                              #  SQRT(ALFA)  E: (10+NA)  M: (9+NA)
059004,000485: 27,3357           57201                                   0          -3,2                       #  X2=-NA
059005,000486: 27,3360           77671                          DDV                                            #  C(MPAC) =NUM SQRT(ALFA)        E:(-3) OR (-16)
059006,000487:                                                                                                 #                                 M:(-3) OR (-15)
059007,000488:                                                                                                 #  ABS(DEN) FROM PDL+2    E:(-3) OR (-16)
059008,000489:                                                                                                 #                         M:(-3) OR (-15)
059009,000490: 27,3361           40145                          DLOAD    BOV                                   #  (THE DLOAD IS SHARED WITH TFFELL)
059010,000491: 27,3362           00043                                   TFFX                                  #  NUM  E: (-16) OR (-29)  M:(-15) OR (-27)
059011,000492: 27,3363           57444                                   TFFELL                                #  USE EQN FOR DELE GEQ 90, LEQ -90
059012,000493: 
059013,000494:                                                                                                 #  OTHERWISE, CONTINUE FOR GENERAL CONIC FOR TFF EQN
059014,000495: 
059015,000496: 27,3364           45471                          DDV      STADR                                 
059016,000497:                                                                                                 #  DEN FROM PDL+0 E: (-3) OR (-16)
059017,000498:                                                                                                 #                         M: (-3) OR (-15)
059018,000499: 27,3365           77732                          STORE    TFFTEM                                #  Z  SAVE FOR SIGN OF SDELF.
059019,000500:                                                                                                 # # Page 1382
059020,000501:                                                                                                 #  E: (-13)  M: (-12)
059021,000502: 27,3366           63406                          PUSH     DSQ                                   #  Z TO PDL+0
059022,000503: 27,3367           41206                          PUSH     DMP                                   #  Z SQ TO PDL+2  E: (-26)  M: (-24)
059023,000504: 27,3370           00035                                   TFFNP                                 #  LC P  E: (-38+2NR)  M: (-36+NR)
059024,000505: 27,3371           75261                          SL       SIGN                                  
059025,000506: 27,3372           20206                                   5                                     
059026,000507: 27,3373           00045                                   TFFTEM                                #  AFFIX SIGN FOR SDELF (ENTRY DISPLAY)
059027,000508: 27,3374           14045                          STODL    TFFTEM                                #  P ZSQ  E: (-59+2NR)  M: (-55+2NR)
059028,000509:                                                                                                 #  (ARG IS USED IN TFF/TRIG)
059029,000510:                                                                                                 #  ZSQ FROM PDL+2  E: (-26)  M: (-24)
059030,000511: 27,3375           41206                          PUSH     DMP                                   #  RESTORE PUSH LOC
059031,000512: 27,3376           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
059032,000513: 27,3377           77657                          SL*                                            
059033,000514: 27,3400           20201                                   0,1                                   #  X1=-NR
059034,000515: 27,3401           00043                          STORE    TFFX                                  #  X
059035,000516: 27,3402           41234                          RTB      DMP                                   
059036,000517: 27,3403           57526                                   T(X)                                  #   POLY
059037,000518:                                                                                                 #  ZSQ FROM PDL+2  E: (-26)  M: (-24)
059038,000519: 27,3404           44302                          SR2      BDSU                                  #  2 ZSQ T(X)  E: (-29)  M: (-27)
059039,000520: 27,3405           00023                                   RTERM                                 #  RTERM  E: (-29)  M: (-27)
059040,000521: 27,3406           41215                          DAD      DMP                                   
059041,000522: 27,3407           00015                                   RMAG1                                 #  E: (-29)  M: (-27)
059042,000523:                                                                                                 #  Z FROM PDL+0  E: (-13)  M: (-12)
059043,000524: 27,3410           51042                          SR3      BPL                                   #  TFF SQRT(MU)  E: (-45)  M: (-42)
059044,000525: 27,3411           57432                                   ENDTFF                                #  (NO PUSH UP)
059045,000526: 27,3412           75206                          PUSH     SIGN                                  #  TFF SQRT(MU)  TO PDL+0
059046,000527: 27,3413           00017                                   TFFQ1                                 #  Q1 FOR GONEPAST TEST
059047,000528: 27,3414           71244                          BPL      DLOAD                                 #  GONE PAST ?
059048,000529: 27,3415           57436                                   NEGTFF                                #  YES. TFF < 0.
059049,000530: 27,3416           00027                                   TFF1/ALF                              #  1/ALFA  E: (-22-2NA)  M: (-20-2NA)
059050,000531: 27,3417           51076                          DCOMP    BPL                                   #  ALFA > 0 ?
059051,000532: 27,3420           57436                                   NEGTFF                                #  NO. TFF IS NEGATIVE.
059052,000533: 
059053,000534:                                                                                                 #  CORRECT FOR ORBITAL PERIOD.
059054,000535: 
059055,000536: 27,3421           77676                          DCOMP                                          #  YES.  CORRECT FOR ORB PERIOD.
059056,000537: 27,3422           56205                          DMP      DDV                                   
059057,000538: 27,3423           31761                                   PI/16                                 #  2 PI (-5)
059058,000539: 27,3424           00031                                   TFFRTALF                              #  SQRT(ALFA)  E: (10+NA)  M: (9+NA)
059059,000540: 27,3425           53657                          SL*      SL*                                   
059060,000541: 27,3426           57602                                   0          -4,2                       #  X2=-NA
059061,000542: 27,3427           57602                                   0          -4,2                       
059062,000543: 27,3430           43257                          SL*      DAD                                   
059063,000544: 27,3431           57576                                   0,2                                   
059064,000545:                                                                                                 #  TFF SQRT(MU) FROM PDL+0        E:(-45) M:(-42)
059065,000546: 27,3432           40005        ENDTFF            DMP      BOV                                   #  TFF SQRT(MU) IN MPAC           E:(-45) M:(-42)
059066,000547: 27,3433           00037                                   TFF/RTMU                              #  E: (17)  M: (14)
059067,000548: 27,3434           57442                                   MAXTFF                                #  SET POSMAX IN OVFL.
059068,000549: 
059069,000550: 27,3435           77616        DUMPTFF2          RVQ                                            #  RETURN TFF     (-28) CS IN MPAC.
059070,000551: 
059071,000552:                                                                                                 # # Page 1383
059072,000553: 27,3436           77745        NEGTFF            DLOAD                                          
059073,000554:                                                                                                 #  TFF SQRT(MU) FROM PDL+0, NEGATIVE.
059074,000555: 27,3437           77650                          GOTO                                           
059075,000556: 27,3440           57432                                   ENDTFF                                
059076,000557: 
059077,000558: 27,3441           77745        MAXTFF1           DLOAD                                          #  RESET PDL
059078,000559: 27,3442           43545        MAXTFF            DLOAD    RVQ                                   
059079,000560: 27,3443           15344                                   NEARONE                               
059080,000561: 
059081,000562:                                                                                                 #  TIME OF FLIGHT ELLIPSE WHEN DEL (ECCENTRIC ANOM) GEQ 90 AND LEQ -90.
059082,000563: 
059083,000564:                                                                                                 #  NUM FROM TFFX. E: (-16) OR (-29)
059084,000565:                                                                                                 #                         M: (-15) OR (-27)
059085,000566: 27,3444           77712        TFFELL            SL2                                            #  NUM  E:(-14) OR (-27)  M:(-13) OR (-25)
059086,000567: 27,3445           41465                          BDDV     PUSH                                  #  TEMP SAVE D/N IN PDL+0
059087,000568:                                                                                                 #  DEN FROM PDL+0  E:(-3)/(-16)  M:(-3)/(-15)
059088,000569:                                                                                                 #  N/D TO PDL+0  E: (11)  M: (10)
059089,000570: 27,3446           45345        TFFEL1            DLOAD    DSU                                   #  (ENTER WITH D/N=0 IN PDL+0)
059090,000571: 27,3447           00045                                   TFFTEM                                #  Q2  E: (-16)  M: (-15)
059091,000572: 27,3450           00017                                   TFFQ1                                 #  Q1  E: (-16)  M: (-15)
059092,000573: 27,3451           14013                          STODL    TFFDELQ                               #  Q2-Q1  E: (-16)  M: (-15)
059093,000574:                                                                                                 #  D/N FROM PDL+0
059094,000575: 27,3452           77626                          STADR                                          
059095,000576: 27,3453           77732                          STORE    TFFTEM                                #  D/N  E: (11)  M: (10)
059096,000577: 27,3454           53605                          DMP      SL*                                   
059097,000578: 27,3455           00027                                   TFF1/ALF                              #  1/ALFA  E: (-22-2NA)  M: (-20-2NA)
059098,000579: 27,3456           57576                                   0,2                                   #  1/ALFA Z  E: (-11-NA)  M: (-10-NA)
059099,000580: 27,3457           41206                          PUSH     DMP                                   #  TO PDL+0
059100,000581: 27,3460           00045                                   TFFTEM                                #  1/Z  E: (11)  M: (10)
059101,000582: 27,3461           41057                          SL*      BOVB                                  
059102,000583: 27,3462           57576                                   0,2                                   #  X2= -NA
059103,000584: 27,3463           45761                                   SIGNMPAC                              #  IN CASE X= 1.0, CONTINUE
059104,000585: 27,3464           00043                          STORE    TFFX                                  #  X=1/ALFA ZSQ
059105,000586: 27,3465           41234                          RTB      DMP                                   
059106,000587: 27,3466           57526                                   T(X)                                  #  POLY
059107,000588: 27,3467           00043                                   TFFX                                  
059108,000589: 27,3470           45242                          SR3      DSU                                   
059109,000590: 27,3471           31767                                   DP2(-3)                               
059110,000591: 27,3472           41405                          DMP      PUSH                                  #  2(X T(X)-1) /Z ALFA    E: (-15-NA)
059111,000592:                                                                                                 #                         M: (-14-NA)
059112,000593:                                                                                                 #  1/ALFA Z FROM PDL+0    E: (-11-NA)
059113,000594:                                                                                                 #                         M: (-10-NA)
059114,000595: 27,3473           41345                          DLOAD    DMP                                   #  GET SIGN FOR SDELF
059115,000596: 27,3474           00045                                   TFFTEM                                #  1/Z  E: (11)  M: (10)
059116,000597: 27,3475           00015                                   RMAG1                                 #  E: (-29)  M: (-27)
059117,000598: 27,3476           43312                          SL2      DAD                                   
059118,000599: 27,3477           00017                                   TFFQ1                                 #  Q1  E: (-16)  M: (-15)
059119,000600: 27,3500           14045                          STODL    TFFTEM                                #  (Q1+R 1/Z) =SGN OF SDELF  E:(-16) M:(-15)
059120,000601: 27,3501           00035                                   TFFNP                                 #  LC P  E: (-38+2NR)  M: (-36+2NR)
059121,000602: 27,3502           53605                          DMP      SL*                                   #  CALC FOR ARG FOR TFF/TRIG.
059122,000603:                                                                                                 # # Page 1384
059123,000604: 27,3503           00027                                   TFF1/ALF                              #  1/ALFA  E:(-22-2NA)  M:(-20-2NA)
059124,000605: 27,3504           57575                                   1,2                                   #  X2=-NA
059125,000606: 27,3505           53765                          SIGN     SL*                                   
059126,000607: 27,3506           00045                                   TFFTEM                                #  AFFIX SIGN FOR SDELF
059127,000608: 27,3507           57576                                   0,2                                   
059128,000609: 27,3510           14045                          STODL    TFFTEM                                #  P/ALFA  E:(-59+2NR)  M:(-55+2NR)
059129,000610:                                                                                                 #  (ARG FOR USE IN TFF/TRIG)
059130,000611: 27,3511           00027                                   TFF1/ALF                              #  1/ALFA  E:(-22-2NA)  M:(-20-2NA)
059131,000612: 27,3512           41366                          SQRT     DMP                                   
059132,000613: 27,3513           31761                                   PI/16                                 #  PI (-4)
059133,000614: 27,3514           77615                          DAD                                            
059134,000615:                                                                                                 #  2(XT(X)-1)/Z ALFA FROM PDL     E:(-15-NA)
059135,000616:                                                                                                 #                                 M:(-14-NA)
059136,000617: 27,3515           45257                          SL*      DSU                                   
059137,000618: 27,3516           57577                                   0          -1,2                       
059138,000619: 27,3517           00013                                   TFFDELQ                               #  Q2-Q1  E: (-16)  M: (-15)
059139,000620: 27,3520           53605                          DMP      SL*                                   
059140,000621: 27,3521           00027                                   TFF1/ALF                              #  1/ALFA  E:(-22-2NA)  M:(-20-2NA)
059141,000622: 27,3522           57601                                   0          -3,2                       
059142,000623: 27,3523           52057                          SL*      GOTO                                  
059143,000624: 27,3524           57602                                   0          -4,2                       
059144,000625: 27,3525           57432                                   ENDTFF                                #  TFF SQRT(MU) IN MPAC E:(-45) M:(-42)
059145,000626: 
059146,000627:                                                                                                 # # Page 1385
059147,000628: 
059148,000629:                                                                                                 #  PROGRAM NAME:          T(X)                                            DATE:  01.17.67
059149,000630:                                                                                                 #  MOD NO:  0                                                             LOG SECTION:  TIME OF FREE FALL
059150,000631:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
059151,000632: 
059152,000633:                                                                                                 #  FUNCTIONAL DESCRIPTION:        THE POLYNOMIAL T(X) IS USED BY TIME OF FLIGHT SUBROUTINES  CALCTFF AND
059153,000634:                                                                                                 #         CALCTPER TO APPROXIMATE THE SERIES
059154,000635:                                                                                                 #                            2     3
059155,000636:                                                                                                 #                 1/3 -X/5 +X /7 -X /9 ...
059156,000637: 
059157,000638:                                                                                                 #         WHERE   X = ALFA Z Z            IF ALFA Z Z LEQ 1
059158,000639:                                                                                                 #                 X = 1/(ALFA Z Z)        IF ALFA Z Z  G  1
059159,000640: 
059160,000641:                                                                                                 #         ALSO    X IS NEG FOR HYPERBOLIC ORBITS
059161,000642:                                                                                                 #                 X = 0 FOR PARABOLIC ORBITS
059162,000643:                                                                                                 #                 X IS POSITIVE FOR ELLIPTIC ORBITS
059163,000644:                                                                                                 #         FOR FLIGHT 278, THE POLYNOMIAL T(X) IS FITTED OVER THE RANGE  (0,+1)  AND HAS A MAXIMUM
059164,000645:                                                                                                 #         DEVIATION FROM THE SERIES OF  2 E-5.  (T(X) IS A CHEBYCHEV TYPE FIT AND WAS OBTAINED USING
059165,000646:                                                                                                 #         MAC PROGRAM AUTCURFIT294RRB AND IS VALID TO THE SAME TOLERANCE OVER THE RANGE  (-.08,+1). )
059166,000647: 
059167,000648:                                                                                                 #  CALLING SEQUENCE:      RTB
059168,000649:                                                                                                 #                                 T(X)
059169,000650:                                                                                                 #         C(MPAC) = X
059170,000651: 
059171,000652:                                                                                                 #  SUBROUTINES CALLED:  NONE
059172,000653: 
059173,000654:                                                                                                 #  NORMAL EXIT MODE:    TC DANZIG
059174,000655: 
059175,000656:                                                                                                 #  ALARMS:  NONE
059176,000657: 
059177,000658:                                                                                                 #  OUTPUT:  C(MPAC) = T(X)
059178,000659: 
059179,000660:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
059180,000661:                                                                                                 #         C(MPAC) = X
059181,000662: 
059182,000663:                                                                                                 #  DEBRIS:  NONE
059183,000664: 
059184,000665: 27,3526           07171        T(X)              TC       POLY                                  
059185,000666: 27,3527           00004                          DEC      4          B-14                       #  N-1
059186,000667: 27,3530           12525 12525                    2DEC     3.333333333 E-1                        
059187,000668: 27,3532           71463 57703                    2DEC*    -1.999819135 E-1*                       
059188,000669: 27,3534           04423 17645                    2DEC*    1.418148467 E-1*                       
059189,000670: 27,3536           74604 43667                    2DEC*    -1.01310997 E-1*                       
059190,000671: 27,3540           01626 37256                    2DEC*    5.609004986 E-2*                       
059191,000672: 27,3542           77404 52071                    2DEC*    -1.536156925 E-2*                       
059192,000673: 
059193,000674: 27,3544           06030        ENDT(X)           TC       DANZIG                                
059194,000675: 
059195,000676: 27,3545  27,3544               TCDANZIG          =        ENDT(X)                               
059196,000677: 
059197,000678:                                                                                                 # # Page 1386
059198,000679:                                                                                                 #  TFF CONSTANTS
059199,000680: 
059200,000681: 34,2000                                          SETLOC   TOF-FF1                               
059201,000682: 34,2000                                          BANK                                           
059202,000683: 
059203,000684: 34,3756                                          COUNT*   $$/TFF                                
059204,000685:                                                                                                 #                                                 # NOTE _  NOTE _ ADJUSTED MUE FOR NEAR EARTH TRAJ.
059205,000686: 

(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.

059210,000691: 
059211,000692:                                                                                                 # MUE             =       3.990815471 E10         # M CUBE/CS SQ
059212,000693:                                                                                                 # RTMUE           =       1.997702549 E5 B-18*    # MODIFIED EARTH MU
059213,000694: 
059214,000695: 34,3756           24775 30424  1/RTMU            2DEC*    .5005750271 E-5        B17*            #  MODIFIED EARTH MU
059215,000696: 

(JL,2009-08-21) End of garbled section.

059217,000698: 
059218,000699:                                                                                                 #                                                 # NOTE _  NOTE _ ADJUSTED MUE FOR NEAR EARTH TRAJ.
059219,000700:                                                                                                 # MUM             =       4.902778 E8             # M CUBE/CS SQ
059220,000701:                                                                                                 # RTMUM           2DEC*   2.21422176 E4 B-18*
059221,000702: 34,3760           06220 37553  PI/16             2DEC     3.141592653 B-4                        
059222,000703: 34,3762           37777 37700  LIM(-22)          2OCT     3777737700                            #  1.0 -B(-22)
059223,000704: 34,3764           00000 00100  DP(-22)           2OCT     0000000100                            #  B(-22)
059224,000705: 34,3766           04000 00000  DP2(-3)           2DEC     1          B-3                        
059225,000706: 34,3770           02000 00000  DP2(-4)           2DEC     1          B-4                        #  1/16
059226,000707: 
059227,000708:                                                                                                 #  RPAD1          2DEC    6373338 B-29            # M (-29) = 20 909 901.57 FT
059228,000709: 34,3772  22,3172               RPAD1             =        RPAD                                  
059229,000710: 34,3772  26,3343               NEARONE           EQUALS   DPPOSMAX                              
059230,000711: 34,3772  26,3335               TFFZEROS          EQUALS   HI6ZEROS                              
059231,000712: 34,3772  26,3325               TFF1/4            EQUALS   HIDP1/4                               

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