Source Code

These source-code files were obtained by digitally photographing the Colossus 237 program listing from the private collection of original AGC developer Fred Martin, and then modifying pre-existing Colossus 249 (Apollo 9) source files to incorporate changes. Note that the page images presented online are of reduced quality, and that additional images are available. Report any conversion errors or legibility problems in page images to info@sandroid.org. Notations on the program listing read, in part:

	GAP: ASSEMBLE REVISION 237 OF AGC PROGRAM COLOSSUS 
	BY NASA 2021111-031  1:17 AUG. 24, 1968

056523,000002:                                                                                                 #  Copyright:    Public domain.
056524,000003:                                                                                                 #  Filename:      TIME_OF_FREE_FALL.agc
056525,000004:                                                                                                 #  Purpose:      Part of the source code for Colossus build 237.
056526,000005:                                                                                                 #                This is for the Command Module's (CM) Apollo Guidance
056527,000006:                                                                                                 #                Computer (AGC), we believe for Apollo 8.
056528,000007:                                                                                                 #  Assembler:    yaYUL
056529,000008:                                                                                                 #  Contact:      Jim Lawton <jim DOT lawton AT gmail DOT com>
056530,000009:                                                                                                 #  Website:      www.ibiblio.org/apollo/index.html
056531,000010:                                                                                                 #  Page scans:   www.ibiblio.org/apollo/ScansForConversion/Colossus237/
056532,000011:                                                                                                 #  Mod history:  2011-04-17 JL    Adapted from corresponding Colossus 249 file.
056533,000012:                                                                                                 #                 2011-04-17 JL   Removed temporary line.
056534,000013: 
056535,000014:                                                                                                 # # Page 1337
056536,000015:                                                                                                 #           THE TFF SUBROUTINES MAY BE USED IN EITHER EARTH OR MOON CENTERED COORDINATES. THE TFF ROUTINES NEVER
056537,000016:                                                                                                 #  KNOW WHICH ORIGIN APPLIES. IT IS THE USER WHO KNOWS, AND WHO SUPPLIES  RONE, VONE AND 1/SQRT(MU)  AT THE
056538,000017:                                                                                                 #  APPROPRIATE SCALE LEVEL FOR THE PROPER PRIMARY BODY.
056539,000018: 
056540,000019:                                                                                                 #           EARTH ORIGIN           POSITION      -29       METERS
056541,000020:                                                                                                 #                                  VELOCITY      -7        METERS/CENTISECOND
056542,000021:                                                                                                 #                                  1/SQRT(MU)    +17       SQRT(CS SQ/METERS CUBED)
056543,000022: 
056544,000023:                                                                                                 #           MOON ORIGIN            POSITION      -27       METERS
056545,000024:                                                                                                 #                                  VELOCITY      -5        METERS/CENTISECONDS
056546,000025:                                                                                                 #                                  1/SQRT(MU)    +14       SQRT(CS SQ/METERS CUBED)
056547,000026: 
056548,000027:                                                                                                 #  ALL DATA PROVIDED TO AND RECEIVED FROM ANY TFF SUBROUTINE WILL BE AT ONE OF THE LEVELS ABOVE. IN ALL CASES,
056549,000028:                                                                                                 #  THE FREE FALL TIME IS RETURNED IN CENTISECONDS AT (-28). PROGRAM  TFF/CONIC  WILL GENERATE VONE/RTMU  AND
056550,000029:                                                                                                 #  LEAVE IT IN  VONE' AT (+10) IF EARTH ORIGIN AND (+9) IF MOON ORIGIN.
056551,000030:                                                                                                 #           THE USER MUST STORE THE STATE VECTOR IN  RONE, VONE  AND MU IN  THE FORM 1/SQRT(MU) IN  TFF/RTMU
056552,000031:                                                                                                 #  AT THE PROPER SCALE BEFORE CALLING  TFF/CONIC. SINCE  RONE, VONE  ARE IN THE EXTENDED VERB STORAGE AREA,
056553,000032:                                                                                                 #  THE USER MUST ALSO LOCK OUT THE EXTENDED VERBS, AND RELEASE THEM WHEN FINISHED.
056554,000033:                                                                                                 #           PROGRAMS     CALC/TFF  AND  CALC/TPER  ASSUME THAT THE TERMINAL RADIUS IS LESS THAN THE PRESENT
056555,000034:                                                                                                 #  RADIUS. THIS RESTRICTION CAN BE REMOVED BY A  15 W CODING CHANGE, BUT AT PRESENT IT IS NOT DEEMED NECESSARY.
056556,000035: 
056557,000036:                                                                                                 #           THE FOLLOWING ERASABLE QUANTITIES ARE USED BY THE TFF ROUTINES, AND ARE LOCATED IN THE PUSH LIST.
056558,000037: 
056559,000038: 
056560,000039:                                                                                                 #                                                               BELOW      E:   IS USED FOR EARTH ORIGIN SCALE
056561,000040:                                                                                                 #                                                                          M:   IS USED FOR MOON  ORIGIN SCALE
056562,000041: 
056563,000042:                                                                                                 #  TFFSW          =       119D  BIT1      # 0 = CALCTFF        1 = CALCTPER
056564,000043: ???????                        TFFDELQ           =        10D                                   #  Q2-Q1              E: (-16)  M: (-15)
056565,000044: ???????                        RMAG1             =        12D                                   #  ABVAL(RN)  M       E: (-29)  M: (-27)
056566,000045:                                                                                                 #  RPER           =       14D             # PERIGEE RADIUS M   E: (-29)  M: (-27)
056567,000046: ???????                        TFFQ1             =        14D                                   #  R.V / SQRT(MUE)    E: (-16)  M: (-15)
056568,000047:                                                                                                 #  SDELF/2                                # SIN(THETA) /2
056569,000048: ???????                        CDELF/2           =        14D                                   #  COS(THETA) /2
056570,000049:                                                                                                 #  RAPO           =       16D             # APOGEE RADIUS  M   E: (-29)  M: (-27)
056571,000050: ???????                        NRTERM            =        16D                                   #  TERMINAL RADIUS M  E: (-29+NR)
056572,000051:                                                                                                 #                     M: (-27+NR)
056573,000052: ???????                        RTERM             =        18D                                   #  TERMINAL RADIUS M  E: (-29)  M: (-27)
056574,000053: ???????                        TFFVSQ            =        20D                                   #  -(V SQUARED/MU)  1/M  E: (20)  M: (18)
056575,000054: ???????                        TFF1/ALF          =        22D                                   #  SEMI MAJ AXIS  M  E: (-22-2 NA)
056576,000055:                                                                                                 #                    M: (-20-2 NA)
056577,000056: ???????                        TFFRTALF          =        24D                                   #  SQRT(ALFA)  E: (10+NA)  M: (9+NA)
056578,000057: ???????                        TFFALFA           =        26D                                   #  ALFA  1/M   E: (26-NR)  M: (24-NR)
056579,000058: ???????                        TFFNP             =        28D                                   #  SEMI LATUS RECTUM  M  E: (-38+2 NR)
056580,000059:                                                                                                 #                        M: (-36+2 NR)
056581,000060: ???????                        TFF/RTMU          =        30D                                   #  1/SQRT(MU)         E: (17)   M: (14)
056582,000061: ???????                        NRMAG             =        32D                                   #  PRESENT RADIUS  M  E: (-29+NR)
056583,000062:                                                                                                 #                     M: (-27+NR)
056584,000063: ???????                        TFFX              =        34D                                   
056585,000064: ???????                        TFFTEM            =        36D                                   #  TEMPORARY
056586,000065:                                                                                                 # # Page 1338
056587,000066:                                                                                                 #                                          REGISTERS S1, S2 ARE UNTOUCED BY ANY TFF SUBROUTINE
056588,000067:                                                                                                 #                                          INDEX REGISTERS X1, X2 ARE USED BY ALL TFF SUBROUTINES. THEY ARE ESTAB-
056589,000068:                                                                                                 #                                          LISHED IN TFF/CONIC AND MUST BE PRESERVED BETWEEN CALLS TO SUBSEQUENT
056590,000069:                                                                                                 #                                          SUBROUTINES.
056591,000070:                                                                                                 #                                          -NR                             C(X1) = NORM COUNT OF RMAG
056592,000071:                                                                                                 #                                          -NA                             C(X2)= NORM COUNT OF SQRT(ABS(ALFA))
056593,000072: 
056594,000073:                                                                                                 # # Page 1339
056595,000074: 
056596,000075: 
056597,000076:                                                                                                 #  SUBROUTINE NAME:     TFFCONIC                                           DATE:    01.29.67
056598,000077:                                                                                                 #  MOD NO:  0                                                              LOG SECTION:    TIME OF FREE FALL
056599,000078:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
056600,000079:                                                                                                 #  MOD NO:  1           MOD BY:  RR BAIRNSFATHER       DATE:  11 APR 67
056601,000080:                                                                                                 #  MOD NO:  2           MOD BY:  RR BAIRNSFATHER       DATE:  21 NOV 67     ADD MOON MU.
056602,000081:                                                                                                 #  MOD NO:  3           MOD BY:  RR BAIRNSFATHER       DATE:  21 MAR 68     ACCEPT DIFFERENT EARTH/MOON SCALES
056603,000082:                                                                                                 #  FUNCTIONAL DESCRIPTION:      THIS SUBROUTINE IS CALLED TO COMPUTE THOSE CONIC PARAMETERS REQUIRED BY THE TFF
056604,000083:                                                                                                 #           SUBROUTINES AND TO ESTABLISH THEM IN THE PUSH LIST AREA. THE PARAMETERS ARE LISTED UNDER OUTPUT.
056605,000084:                                                                                                 #           THE EQUATIONS ARE
056606,000085:                                                                                                 #                  -   -  -
056607,000086:                                                                                                 #                  H = RN*VN                                               ANGULAR MOMENTUM
056608,000087:                                                                                                 #                        - -
056609,000088:                                                                                                 #                  LCP = H.H / MU                                          SEMI LATUS RECTUM
056610,000089:                                                                                                 #                                -  -
056611,000090:                                                                                                 #                  ALFA = 2/RN - VN.VN / MU                                RECIPROCAL SEMI MAJ AXIS, SIGNED
056612,000091: 
056613,000092:                                                                                                 #           AND ALFA IS POS FOR ELLIPTIC  ORBITS
056614,000093:                                                                                                 #                        0  FOR PARABOLIC ORBITS
056615,000094:                                                                                                 #                       NEG FOR HYPERBOLIC ORBITS.
056616,000095:                                                                                                 #           SUBROUTINE ALSO COMPUTES AND SAVES  RMAG.
056617,000096:                                                                                                 #  CALLING SEQUENCE:
056618,000097:                                                                                                 #           TFFCONIC  EXPECTS CALLER TO ENTER WITH CORRECT GRAVITATIONAL CONSTANT IN MPAC, IN THE FORM
056619,000098:                                                                                                 #           1/SQRT(MU). PROGRAM WILL SAVE IN  TFF/RTMU  . THE SCALE IS DETERMINED BY WHETHER EARTH OR MOON
056620,000099:                                                                                                 #           ORIGIN IS USED.  THE CALLER MUST LOCK OUT THE EXTENDED VERBS BEFORE PROVIDING STATE VECTOR IN RONE,
056621,000100:                                                                                                 #           VONE  AT PROPER SCALE. THE EXTENDED VERBS MUST BE RESTORED WHEN THE CALLER IS FINISHED USING THE
056622,000101:                                                                                                 #           TFF ROUTINES.
056623,000102:                                                                                                 #           ENTRY POINT  TFFCONMU  EXPECTS THAT  TFF/RTMU  IS ALREADY LOADED.
056624,000103:                                                                                                 #           TO SPECIFY MU:      DLOAD  CALL                IF MU ALREADY STORED:    CALL
056625,000104:                                                                                                 #                                      YOURMU       1/RTMU   E: (17)  M: (14)              TFFCONMU
056626,000105:                                                                                                 #                                      TFFCONIC
056627,000106:                                                                                                 #           PUSHLOC = PDL+0, ARBITRARY IF LEQ 18D
056628,000107: 
056629,000108:                                                                                                 #  SUBROUTINES CALLED:  NONE
056630,000109:                                                                                                 #  NORMAL EXIT MODES:   RVQ
056631,000110:                                                                                                 #  ALARMS:  NONE
056632,000111:                                                                                                 #  OUTPUT:  THE FOLLOWING ARE STORED IN THE PUSH LIST AREA.
056633,000112:                                                                                                 #           RMAG1  E:(-29) M:(-27) M  RN, PRESENT RADIUS LENGTH.
056634,000113:                                                                                                 #           NRMAG  E: (-29+NR)     M  RMAG, NORMALIZED
056635,000114:                                                                                                 #                  M: (-27+NR)
056636,000115:                                                                                                 #           X1                        -NR, NORM COUNT
056637,000116:                                                                                                 #           TFFNP  E: (-38+2NR)    M  LCP, SEMI LATUS RECTUM, WEIGHTED BY NR.       FOR VGAMCALC
056638,000117:                                                                                                 #                  M: (-36+2NR)
056639,000118:                                                                                                 #           TFF/RTMU  E:(17)  M:(14)       1/SQRT(MU)
056640,000119:                                                                                                 #           TFFVSQ E:(20)  M:(18)  1/M   -(V SQ/MU): PRESENT VELOCITY,NORMLIZED.    FOR VGAMCALC
056641,000120:                                                                                                 #           TFFALFA  E: (26-NR)    1/M   ALFA, WEIGHTED BY NR
056642,000121:                                                                                                 #                    M: (24-NR)
056643,000122:                                                                                                 #           TFFRTALF  E: (10+NA)   SQRT(ALFA), NORMALIZED
056644,000123:                                                                                                 #                     M: (9+NA)
056645,000124:                                                                                                 # # Page 1340
056646,000125:                                                                                                 #           X2                     -NA, NORMCOUNT
056647,000126:                                                                                                 #           TFF1/ALF  E: (-22-2NA) SIGNED SEMI MAJ AXIS, WEIGHTED BY NA
056648,000127:                                                                                                 #                     M: (-20-2NA)
056649,000128:                                                                                                 #           PUSHLOC AT PDL+0
056650,000129:                                                                                                 #           THE FOLLOWING  IS STORED IN GENERAL ERASABLE
056651,000130:                                                                                                 #           VONE'  E:(10)  M:(9)   V/RT(MU), NORMALIZED VELOCITY
056652,000131:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
056653,000132:                                                                                                 #           RONE   E:(-29)  M:(-27)  M      STATE VECTOR                   LEFT BY CALLER
056654,000133:                                                                                                 #           VONE   E:(-7)   M:(-5)   M/CS   STATE VECTOR                   LEFT BY CALLER
056655,000134:                                                                                                 #           TFF/RTMU  E:(17)  M:(14)      1/RT(CS SQ/M CUBE)               IF ENTER VIA TFFCONMU.
056656,000135:                                                                                                 #  DEBRIS:  QPRET,    PDL+0 ... PDL+3
056657,000136: 
056658,000137: 
056659,000138: 
056660,000139: 33,3773                                          BANK     33                                    
056661,000140: 27,2000                                          SETLOC   TOF-FF                                
056662,000141: 27,2000                                          BANK                                           
056663,000142: 
056664,000143: 27,2721                                          COUNT*   $$/TFF                                
056665,000144: 
056666,000145: 27,2721           00037        TFFCONIC          STORE    TFF/RTMU                              #  1/SQRT(MU)       E: (17)   M: (14)
056667,000146: 
056668,000147: 27,2722           53575        TFFCONMU          VLOAD    UNIT                                  #  COME HERE WITH TFFRTMU LOADED.
056669,000148: 27,2723           02327                                   RONE                                  #  SAVED RN.  M     E: (-29)  M: (-27)
056670,000149: 27,2724           77725                          PDDL                                           #  UR/2 TO PDL+0, +5
056671,000150: 27,2725           00045                                   36D                                   #  MAGNITUDE
056672,000151: 27,2726           00015                          STORE    RMAG1                                 #  M    E: (-29)   M: (-27)
056673,000152: 
056674,000153: 27,2727           77701                          NORM                                           
056675,000154: 27,2730           00047                                   X1                                    #  -NR
056676,000155: 27,2731           24041                          STOVL    NRMAG                                 #  RMAG  M   E: (-29+NR)  M: (-27+NR)
056677,000156: 27,2732           02335                                   VONE                                  #  SAVED VN.  M/CS   E: (-7)   M: (-5)
056678,000157: 27,2733           77761                          VXSC                                           
056679,000158: 27,2734           00037                                   TFF/RTMU                              #  E: (17)   M: (14)
056680,000159: 27,2735           02372                          STORE    VONE'                                 #  VN/SQRT(MU)   E: (10)   M: (9)
056681,000160: 
056682,000161: 27,2736           47361                          VXSC     VXV                                   
056683,000162: 27,2737           00041                                   NRMAG                                 #  E: (-29+NR)  M: (-27+NR)
056684,000163:                                                                                                 #                                           UR/2 FROM PDL
056685,000164: 27,2740           47572                          VSL1     VSQ                                   #  BEFORE:  E:(-19+NR)  M:(-18+NR)
056686,000165: 27,2741           14035                          STODL    TFFNP                                 #  LC P  M   E:(-38+2NR)  M:(-36+2NR)
056687,000166:                                                                                                 #                                           SAVE ALSO FOR VGAMCALC
056688,000167: 27,2742           15322                                   TFF1/4                                
056689,000168: 27,2743           63271                          DDV      PDVL                                  #  (2/RMAG)  1/M   E:(26-NR)   M:(24-NR)
056690,000169: 27,2744           00041                                   NRMAG                                 #  RMAG  M   E:(-29+NR)  M:(-27+NR)
056691,000170: 27,2745           02372                                   VONE'                                 #  SAVED VN.   E: (10)   M: (9)
056692,000171: 27,2746           57436                          VSQ      DCOMP                                 #  KEEP MPAC+2 HONEST FOR SQRT.
056693,000172: 27,2747           00025                          STORE    TFFVSQ                                #  -(V SQ/MU)   E:(20)   M:(18)
056694,000173:                                                                                                 #                                           SAVE FOR VGAMCALC
056695,000174: 27,2750           43257                          SR*      DAD                                   
056696,000175:                                                                                                 # # Page 1341
056697,000176: 27,2751           20573                                   0          -6,1                       #  GET -VSQ/MU   E:(26-NR)   M:(24-NR)
056698,000177: 27,2752           77626                          STADR                                          
056699,000178:                                                                                                 #                                           2/RMAG   FROM PDL+2
056700,000179: 27,2753           77744                          STORE    TFFALFA                               #  ALFA  1/M   E:(26-NR)   M:(24-NR)
056701,000180: 27,2754           41457                          SL*      PUSH                                  #  TEMP SAVE ALFA    E:(20)   M:(18)
056702,000181: 27,2755           20173                                   0          -6,1                       
056703,000182: 27,2756           75446                          ABS      SQRT                                  #  E:(10)   M:(9)
056704,000183: 27,2757           77701                          NORM                                           
056705,000184: 27,2760           00050                                   X2                                    #  X2 = -NA
056706,000185: 27,2761           00031                          STORE    TFFRTALF                              #  SQRT( ABS(ALFA) )  E:(10+NA)   M:(9+NA)
056707,000186: 27,2762           75316                          DSQ      SIGN                                  #  NOT SO ACCURATE, BUT OK
056708,000187:                                                                                                 #                                           ALFA FROM PDL+2   E:(20)   M:(18)
056709,000188: 27,2763           55254                          BZE      BDDV                                  #  SET 1/ALFA =0, TO SHOW SMALL ALFA
056710,000189: 27,2764           56766                                   +2                                    
056711,000190: 27,2765           15322                                   TFF1/4                                
056712,000191: 27,2766           00027        +2                STORE    TFF1/ALF                              #  1/ALFA   E:(-22-2NA)   M:(-20-2NA)
056713,000192: 27,2767           77616        DUMPCNIC          RVQ                                            
056714,000193:                                                                                                 #                                                                 39 W
056715,000194: 
056716,000195:                                                                                                 # # Page 1342
056717,000196:                                                                                                 #  SUBROUTINE NAME:     TFFRP/RA                                           DATE:  01.17.67
056718,000197:                                                                                                 #  MOD NO:  0                                                              LOG SECTION:    TIME OF FREE FALL
056719,000198:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
056720,000199:                                                                                                 #  MOD NO:  1           MOD BY:  RR BAIRNSFATHER      DATE:  11 APR 67
056721,000200:                                                                                                 #  MOD NO:  2           MOD BY:  RR BAIRNSFATHER      DATE:  21 MAR 68     ACCEPT DIFFERENT EARTH/MOON SCALES
056722,000201:                                                                                                 #                                                                          ALSO IMPROVE ACCURACY OF RAPO.
056723,000202:                                                                                                 #  FUNCTIONAL DESCRIPTION:         USED BY CALCTPER AND TFF DISPLAYS TO CALCULATE PERIGEE RADIUS AND ALSO
056724,000203:                                                                                                 #           APOGEE RADIUS FOR A GENERAL CONIC.
056725,000204:                                                                                                 #           PROGRAM GIVES PERIGEE RADIUS AS               APOGEE RADIUS IS GIVEN BY
056726,000205:                                                                                                 #                  RP = P /(1+E)                                 RA = (1+E) / ALFA
056727,000206:                                                                                                 #           WHERE   2
056728,000207:                                                                                                 #                  E  = 1 - P ALFA
056729,000208:                                                                                                 #           IF RA IS NEGATIVE OR SHOWS DIVIDE OVERFLOW, THEN RA = POSMAX BECAUSE
056730,000209:                                                                                                 #                  1. APOGEE RADIUS IS NOT MEANINGFUL FOR HYPERBOLA
056731,000210:                                                                                                 #                  2. APOGEE RADIUS IS NOT DEFINED FOR PARABOLA
056732,000211:                                                                                                 #                  3. APOGEE RADIUS EXCEEDS THE SCALING FOR ELLIPSE.
056733,000212:                                                                                                 #           THIS SUBROUTINE REQUIRES THE SIGNED  RECIPROCAL SEMI MAJ AXIS,  ALFA,  AND SEMI LATUS RECTUM AS DATA.
056734,000213:                                                                                                 #  CALLING SEQUENCE:  CALL
056735,000214:                                                                                                 #                            TFFRP/RA
056736,000215:                                                                                                 #           PUSHLOC = PDL+0, ARBITRARY IF LEQ 10D
056737,000216:                                                                                                 #           C(MPAC) UNSPECIFIED
056738,000217: 
056739,000218:                                                                                                 #  SUBROUTINES CALLED:  NONE
056740,000219:                                                                                                 #  NORMAL EXIT MODE:    RVQ
056741,000220:                                                                                                 #           IF ELLIPSE, WITHIN NORMAL SCALING, RAPO IS CORRECT.
056742,000221:                                                                                                 #           OTHERWISE, RAPO = POSMAX.
056743,000222:                                                                                                 #  ALARMS:  NONE
056744,000223:                                                                                                 #  OUTPUT:  STORED IN PUSH LIST AREA. SCALE OF OUTPUT AGREES WITH DATA SUPPLIED TO  TFF/CONIC.
056745,000224:                                                                                                 #           RPER   E:(-29)  M:(-27)  M     PERIGEE RADIUS                  DESTROYED BY CALCTFF/CALCTPER, TFFTRIG.
056746,000225:                                                                                                 #           RAPO   E:(-29)  M:(-27)  M     APOGEE  RADIUS                  WILL BE DESTROYED BY CALCTFF/CALCTPER
056747,000226:                                                                                                 #           PUSHLOC AT PDL+0
056748,000227:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
056749,000228:                                                                                                 #           TFFALFA E:(26-NR)      M   1/SEMI MAJ AXIS                     LEFT BY TFFCONIC
056750,000229:                                                                                                 #                   M:(24-NR)
056751,000230:                                                                                                 #           TFFNP  E: (-38+2NR)    M  LC P, SEMI LATUS RECTUM              LEFT BY TFFCONIC
056752,000231:                                                                                                 #                  M: (-36+2NR)
056753,000232:                                                                                                 #           X1                     -NR, NORM COUNT OF RMAG                 LEFT BY TFFCONIC
056754,000233:                                                                                                 #           X2                     -NA, NORM COUNT OF ALFA                 LEFT BY TFFCONIC
056755,000234:                                                                                                 #  DEBRIS:  QPRET,    PDL+0 ... PDL+1
056756,000235: 
056757,000236:                                                                                                 # # Page 1343
056758,000237: 
056759,000238: 27,2770                        RAPO              =        16D                                   #  APOGEE RADIUS  M  E:(-29)  M:(-27)
056760,000239: 27,2770                        RPER              =        14D                                   #  PERIGEE RADIUS  M  E:(-29)  M:(-27)
056761,000240: 
056762,000241: 27,2770           41345        TFFRP/RA          DLOAD    DMP                                   
056763,000242: 27,2771           00033                                   TFFALFA                               #  ALFA  1/M  E:(26-NR)   M:(24-NR)
056764,000243: 27,2772           00035                                   TFFNP                                 #  LC P  M   E:(-38+2NR)   M:(-36+2NR)
056765,000244: 27,2773           57457                          SR*      DCOMP                                 #  ALFA P (-12+NR)
056766,000245: 27,2774           20571                                   0          -8D,1                      #  ALFA P  (-4)
056767,000246: 27,2775           51415                          DAD      ABS                                   #  (DCOMP GIVES VALID TP RESULT FOR SQRT)
056768,000247:                                                                                                 #                                           (ABS PROTECTS SQRT IF E IS VERY NEAR 0)
056769,000248: 27,2776           17330                                   DP2(-4)                               
056770,000249: 27,2777           43366                          SQRT     DAD                                   #  E SQ = (1- P ALFA)   (-4)
056771,000250: 27,3000           15322                                   TFF1/4                                
056772,000251: 27,3001           55206                          PUSH     BDDV                                  #  (1+E)  (-2)  TO PDL+0
056773,000252: 27,3002           00035                                   TFFNP                                 #  LCP  M  E:(-38+2NR)   M:(-36+2NR)
056774,000253: 27,3003           53657                          SR*      SR*                                   #  (DOES SR THEN SL TO AVOID OVFL)
056775,000254: 27,3004           20601                                   0,1                                   #  X1=-NR
056776,000255: 27,3005           20572                                   0          -7,1                       #  (EFFECTIVE SL)
056777,000256: 27,3006           14017                          STODL    RPER                                  #  PERIGEE RADIUS  M   E:(-29)   M:(-27)
056778,000257:                                                                                                 #                                           (1+E)  (-2)   FROM PDL+0
056779,000258: 27,3007           41005                          DMP      BOVB                                  
056780,000259: 27,3010           00027                                   TFF1/ALF                              #  E:(-22-2NA)  M:(-20-2NA)
056781,000260: 27,3011           57314                                   TCDANZIG                              #  CLEAR OVFIND, IF ON.
056782,000261: 27,3012           53654                          BZE      SL*                                   
056783,000262: 27,3013           57022                                   MAXRA                                 #  SET POSMAX IF ALFA=0
056784,000263: 27,3014           57603                                   0          -5,2                       #  -5+NA
056785,000264: 27,3015           40057                          SL*      BOV                                   
056786,000265: 27,3016           57576                                   0,2                                   
056787,000266: 27,3017           57022                                   MAXRA                                 #  SET POSMAX IF OVFL.
056788,000267: 27,3020           77644                          BPL                                            #  CONTINUE WITH VALID RAPO.
056789,000268: 27,3021           57024                                   +3                                    
056790,000269: 27,3022           77745        MAXRA             DLOAD                                          #  RAPO CALC IS NOT VALID. SET RAPO =
056791,000270: 27,3023           17334                                   NEARONE                               #  POSMAX AS A TAG.
056792,000271: 27,3024           00021        +3                STORE    RAPO                                  #  APOGEE RADIUS  M   E:(-29)   M:(-27)
056793,000272: 27,3025           77616        DUMPRPRA          RVQ                                            
056794,000273: 
056795,000274:                                                                                                 #                                                                  30 W
056796,000275: 
056797,000276:                                                                                                 # # Page 1344
056798,000277:                                                                                                 #  SUBROUTINE NAME:     CALCTPER  /  CALCTFF                               DATE:    01.29.67
056799,000278:                                                                                                 #  MOD NO:  0                                                              LOG SECTION:    TIME OF FREE FALL
056800,000279:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
056801,000280:                                                                                                 #  MOD NO:  1           MOD BY:  RR BAIRNSFATHER       DATE:  21 MAR 67
056802,000281:                                                                                                 #  MOD NO:  2           MOD BY:  RR BAIRNSFATHER       DATE:  14 APR 67
056803,000282:                                                                                                 #  MOD BY:  3           MOD BY:  RR BAIRNSFATHER       DATE:  8JUL 67       NEAR EARTH MUE AND NEG TFF (GONEPAST)
056804,000283:                                                                                                 #  MOD BY:  4           MOD BY:  RR BAIRNSFATHER       DATE:  21 NOV 67     ADD VARIABLE MU.
056805,000284:                                                                                                 #  MOD BY:  5           MOD BY:  RR BAIRNSFATHER       DATE:  21 MAR 68     ACCEPT DIFFERENT EARTH/MOON SCALES
056806,000285:                                                                                                 #  FUNCTIONAL DESCRIPTION:      PROGRAM CALCULATES THE FREE-FALL TIME OF FLIGHT FROM PRESENT POSITION  RN  AND
056807,000286:                                                                                                 #           VELOCITY  VN  TO A RADIUS LENGTH SPECIFIED BY  RTERM  , SUPPLIED BY THE USER. THE POSITION VECTOR
056808,000287:                                                                                                 #           RN  MAY BE ON EITHER SIDE OF THE CONIC, BUT  RTERM  IS CONSIDERED ON THE INBOUND SIDE.
056809,000288:                                                                                                 #           THE EQUATIONS ARE
056810,000289: 
056811,000290:                                                                                                 #                  Q2 = -SQRT(RTERM (2-RTERM ALFA) - LCP)    (INBOUND SIDE))        LEQ +- LCE/SQRT(ALFA)
056812,000291:                                                                                                 #                        -  -
056813,000292:                                                                                                 #                  Q1 = RN.VN / SQRT(MU)                                            LEQ +- LCE/SQRT(ALFA)
056814,000293: 
056815,000294:                                                                                                 #                  Z = NUM / DEN                                                    LEQ +- 1/SQRT(ALFA)
056816,000295: 
056817,000296:                                                                                                 #           WHERE, IF INBOUND
056818,000297:                                                                                                 #                  NUM = RTERM -RN                                                  LEQ +- 2 LCE/ALFA
056819,000298:                                                                                                 #                  DEN = Q2+Q1                                                      LEQ +- 2 LCE/SQRT(ALFA)
056820,000299: 
056821,000300:                                                                                                 #           AND, IF OUTBOUND
056822,000301:                                                                                                 #                  NUM = Q2-Q1                                                      LEQ +- 2 LCE/SQRT(ALFA)
056823,000302:                                                                                                 #                  DEN = 2 - ALFA (RTERM + RN) .                                    LEQ +- 2 LCE
056824,000303: 
056825,000304:                                                                                                 #           IF      ALFA ZZ  <  1.0               (FOR ALL CONICS EXCEPT ELLIPSES HAVING ABS(DEL ECC ANOM) G 90 DEG)
056826,000305: 
056827,000306:                                                                                                 #           THEN   X = ALFA Z Z
056828,000307:                                                                                                 #           AND    TFF = (RTERM +RN -2 ZZ T(X) ) Z/SQRT(MU)
056829,000308: 
056830,000309:                                                                                                 #                  EXCEPT IF ALFA PNZ, AND IF TFF NEG,
056831,000310:                                                                                                 #                  THEN   TFF = 2 PI /(ALFA SQRT(ALFA)) + TFF
056832,000311: 
056833,000312:                                                                                                 #           OR IF  ALFA ZZ  GEQ  1.0             (FOR ELLIPSES HAVING ABS(DEL ECC ANOM) GEQ 90 DEG)
056834,000313: 
056835,000314:                                                                                                 #           THEN   X = 1/ALFA Z Z
056836,000315:                                                                                                 #           AND    TFF = (PI/SQRT(ALFA) -Q2 +Q1 +2(X T(X) -1) /ALFA Z) /ALFA SQRT(MU)
056837,000316: 
056838,000317:                                                                                                 #           WHERE  T(X) IS A POLYNOMIAL APPROXIMATION TO THE SERIES
056839,000318:                                                                                                 #                               2      3                     2
056840,000319:                                                                                                 #                  1/3 - X/5 + X /7 - X /9 ...             (X   <  1.0)
056841,000320: 
056842,000321:                                                                                                 #  CALLING SEQUENCE:  TIME TO RTERM                           TIME TO PERIGEE
056843,000322:                                                                                                 #                     CALL                                    CALL
056844,000323:                                                                                                 #                            CALCTFF                                 CALCTPER
056845,000324:                                                                                                 #                     C(MPAC) = TERMNL RAD  M                 C(MPAC) = PERIGEE RAD  M
056846,000325:                                                                                                 #           FOR EITHER,   E: (-29)    M: (-27)
056847,000326:                                                                                                 #           FOR EITHER, PUSHLOC = PDL+0 , ARBITRARY IF LEQ 8D.
056848,000327:                                                                                                 # # Page 1345
056849,000328:                                                                                                 #  SUBROUTINES CALLED:  T(X),  VIA RTB
056850,000329:                                                                                                 #  NORMAL EXIT MODE:  RVQ
056851,000330:                                                                                                 #           HOWEVER, PROGRAM EXITS WITH ONE OF THE FOLLOWING VALUES FOR TFF (-28) CS IN MPAC. USER MUST STORE.
056852,000331:                                                                                                 #                  A.  TFF= FLIGHT TIME. NORMAL CASE FOR POSITIVE FLIGHT TIME LESS THAN ONE ORBITAL PERIOD.
056853,000332:                                                                                                 #                  B.  (THIS OPTION IS NO LONGER USED.)
056854,000333:                                                                                                 #                  C.  TFF = POSMAX.  THIS INDICATES THAT THE CONIC FROM THE PRESENT POSITION WILL NOT RETURN TO
056855,000334:                                                                                                 #                      THE SPECIFIED ALTITUDE.  ALSO INDICATES OUTBOUND PARABOLA OR HYPERBOLA.
056856,000335:                                                                                                 #  OUTPUT:  C(MPAC)  (-28)  CS     TIME OF FLIGHT, OR TIME TO PERIGEE
056857,000336:                                                                                                 #           TFFX    (0)            X,                                      LEFT FOR ENTRY DISPLAY TFF ROUTINES
056858,000337:                                                                                                 #           NRTERM  E: (-29+NR)    M  RTERM, WEIGHTED BY NR                LEFT FOR ENTRY DISPLAY TFF ROUTINES
056859,000338:                                                                                                 #                   M: (-27+NR)
056860,000339:                                                                                                 #           TFFTEM  E: (-59+2NR)   LCP Z Z SGN(SDELF)                      LEFT FOR ENTRY DISPLAY TFF ROUTINES
056861,000340:                                                                                                 #                   M: (-55+2NR)   LCP /ALFA SGN(SDELF)                    LEFT FOR ENTRY DISPLAY TFF ROUTINES
056862,000341:                                                                                                 #           NOTE:  TFFTEM = PDL 36D AND WILL BE DESTROYED BY .:UNIT:.
056863,000342:                                                                                                 #           RMAG1  E:(-29)  M:(-27)  PDL 12 NOT TOUCHED.
056864,000343:                                                                                                 #           TFFQ1  E:(-16)  M:(-15)  PDL 14D
056865,000344:                                                                                                 #           TFFDELQ E:(-16) M:(-15)  PDL 10D
056866,000345:                                                                                                 #           PUSHLOC AT PDL+0
056867,000346:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
056868,000347:                                                                                                 #           RONE   E:(-29)  M:(-27)  M    STATE VECTOR                     LEFT BY USER
056869,000348:                                                                                                 #           VONE'  E:(+10)  M:(+9)   VN/SQRT(NU)                           LEFT BY TFF/CONIC
056870,000349:                                                                                                 #           RMAG1  E:(-29)  M:(-27)  PRESENT RADIUS,  M                    LEFT BY TFFCONIC
056871,000350:                                                                                                 #           C(MPAC)E:(-29)  M:(-27)  RTERM, TERMINAL RADIUS LENGTH, M      LEFT BY USER
056872,000351: 
056873,000352:                                                                                                 #           THE FOLLOWING ARE STORED IN THE PUSH LIST AREA.
056874,000353:                                                                                                 #           TFF/RTMU  E:(17)  M:(14)  1/SQRT(MU)                           LEFT BY TFFCONIC.
056875,000354:                                                                                                 #           NRMAG  E: (-29+NR)  M  RMAG, NORMALIZED                        LEFT BY TFFCONIC
056876,000355:                                                                                                 #                  M: (-27+NR)
056877,000356:                                                                                                 #           X1                     -NR, NORM COUNT                         LEFT BY TFFCONIC
056878,000357:                                                                                                 #           TFFNP  E: (-38+2NR)   M   LCP, SEMI LATUS RECTUM, WEIGHT  NR    LEFT BY TFFCONIC
056879,000358:                                                                                                 #                  M: (-36+2N4)
056880,000359:                                                                                                 #           TFFALFA  E: (26-NR)    1/M  ALFA, WEIGHT  NR                   LEFT BY TFFCONIC
056881,000360:                                                                                                 #                    M: (24-NR)
056882,000361:                                                                                                 #           TFFRTALF  E:(10+NA)    SQRT(ALFA), NORMALIZED                  LEFT BY TFFCONIC
056883,000362:                                                                                                 #                     M:(9+NA)
056884,000363:                                                                                                 #           X2                     -NA, NORM COUNT                         LEFT BY TFFCONIC
056885,000364:                                                                                                 #           TFF1/ALF  E: (-22-2NA)  SIGNED SEMIMAJ AXIS, WEIGHTED BY NA    LEFT BY TFFCONIC
056886,000365:                                                                                                 #                     M: (-20-2NA)
056887,000366:                                                                                                 #  DEBRIS:  QPRET,   PDL+0 ... PDL+3
056888,000367:                                                                                                 #           RTERM  E:(-29)  M(-27)  RTERM, TERMINAL RADIUS LENGTH
056889,000368:                                                                                                 #           RAPO   E:(-29)  M(-27)  PDL 16D  (=NRTERM)
056890,000369:                                                                                                 #           RPER   E:(-29)  M(-27)  PDL 14D   (=TFFQ1)
056891,000370: 
056892,000371:                                                                                                 # # Page 1346
056893,000372: 27,3026           77614        CALCTPER          SETGO                                          #  ENTER WITH RPER  IN MPAC
056894,000373: 27,3027           03436                                   TFFSW                                 
056895,000374: 27,3030           57033                                   +3                                    
056896,000375: 27,3031           77614        CALCTFF           CLEAR                                          #  ENTER WITH RTERM IN MPAC
056897,000376: 27,3032           03676                                   TFFSW                                 
056898,000377: 27,3033           00023        +3                STORE    RTERM                                 #  E: (-29)  M: (-27)
056899,000378: 27,3034           77657                          SL*                                            
056900,000379: 27,3035           20201                                   0,1                                   #  X1=-NR
056901,000380: 27,3036           00021                          STORE    NRTERM                                #  RTERM  E: (-29+NR)  M: (-27+NR)
056902,000381: 27,3037           44205                          DMP      BDSU                                  
056903,000382: 27,3040           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
056904,000383: 27,3041           15322                                   TFF1/4                                
056905,000384: 27,3042           41206                          PUSH     DMP                                   #  (2-ALFA RTERM)  (-3)  TO PDL+0
056906,000385: 27,3043           00021                                   NRTERM                                #  E: (-29+NR)  M: (-27+NR)
056907,000386: 27,3044           53725                          PDDL     SR*                                   #  RTERM(2-ALFA RTERM) TO PDL+2
056908,000387:                                                                                                 #                                           E: (-32+NR)   M: (-30+NR)
056909,000388: 27,3045           00035                                   TFFNP                                 #  LC P  E:(-38+2NR)  M:(-36+2NR)
056910,000389: 27,3046           20573                                   0          -6,1                       #  X1 = -NR
056911,000390: 27,3047           43276                          DCOMP    DAD                                   #  DUE TO SHIFTS, KEEP PRECISION FOR SQRT
056912,000391:                                                                                                 #                                           RTERM(2-ALFA RTERM) FROM PDL +2
056913,000392:                                                                                                 #                                           E: (-32+NR)  M: (-30+NR)
056914,000393: 27,3050           77657                          SR*                                            #  LEAVE  E: (-32)  M: (-30)
056915,000394: 27,3051           20601                                   0,1                                   #  X1 = -NR
056916,000395: 27,3052           71214                          BOFF     DLOAD                                 #  CHECK TFF / TPER SWITCH
056917,000396: 27,3053           03756                                   TFFSW                                 
056918,000397: 27,3054           57056                                   +2                                    #  IF TFF, CONTINUE
056919,000398: 27,3055           15332                                   TFFZEROS                              #  IF TPER, SET Q2 = 0
056920,000399: 27,3056           75440        +2                BMN      SQRT                                  #  E: (-16)  M: (-15)
056921,000400: 
056922,000401: 27,3057           57211                                   MAXTFF1                               #  NO FREE FALL CONIC TO RTERM FROM HERE
056923,000402:                                                                                                 #                                           RESET PDL, SET TFF=POSMAX, AND EXIT.
056924,000403: 
056925,000404: 27,3060           41076                          DCOMP    BOVB                                  #  RT IS ON INBOUND SIDE.  ASSURE OVFIND=0
056926,000405: 27,3061           57314                                   TCDANZIG                              #  ANY PORT IN A STORM.
056927,000406: 27,3062           24045                          STOVL    TFFTEM                                #  Q2   E: (-16)  M: (-15)
056928,000407: 27,3063           02372                                   VONE'                                 #  VN/SQRT(MU) E: (10)  M: (9)
056929,000408: 27,3064           52441                          DOT      SL3                                   
056930,000409: 27,3065           02327                                   RONE                                  #  SAVED RN.  E: (-29)  M: (-27)
056931,000410: 27,3066           00017                          STORE    TFFQ1                                 #  Q1, SAVE FOR GONEPAST TEST.
056932,000411:                                                                                                 #                                           E: (-16)  M: (-15)
056933,000412: 27,3067           44240                          BMN      BDSU                                  
056934,000413: 27,3070           57111                                   INBOUND                               #  USE ALTERNATE Z
056935,000414: 27,3071           00045                                   TFFTEM                                #  Q2  E: (-16)  M: (-15)
056936,000415: 
056937,000416:                                                                                                 #                                          OUTBOUND  Z CALC CONTINUES HERE
056938,000417: 
056939,000418: 27,3072           14043                          STODL    TFFX                                  #  NUM=Q2-Q1  E: (-16)  M: (-15)
056940,000419: 27,3073           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
056941,000420: 27,3074           44205                          DMP      BDSU                                  
056942,000421:                                                                                                 # # Page 1347
056943,000422: 27,3075           00041                                   NRMAG                                 #  RMAG  E: (-29+NR)  M: (-27+NR)
056944,000423:                                                                                                 #                                           (2-RTERM ALFA)  (-3) FROM PDL+0
056945,000424: 27,3076           51406        SAVEDEN           PUSH     ABS                                   #  DEN TO PDL+0  E:(-3) OR (-16)
056946,000425:                                                                                                 #                                                         M:(-3) OR (-15)
056947,000426: 27,3077           40015                          DAD      BOV                                   #  INDETERMINANCY TEST
056948,000427: 27,3100           17322                                   LIM(-22)                              #  =1.0-B(-22)
056949,000428: 27,3101           57122                                   TFFXTEST                              #  GO IF DEN >/= B(-22)
056950,000429: 27,3102           65345                          DLOAD    PDDL                                  #  SET DEN=0 OTHERWISE
056951,000430: 27,3103           15332                                   TFFZEROS                              
056952,000431:                                                                                                 #                                           XCH ZERO WITH PDL+0
056953,000432: 27,3104           57545                          DLOAD    DCOMP                                 
056954,000433: 27,3105           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
056955,000434: 27,3106           71240                          BMN      DLOAD                                 #  FOR TPER: Z INDET AT DELE/2=0 AND 90.
056956,000435: 27,3107           57216                                   TFFEL1                                #  ASSUME 90, AND LEAVE 0 IN PDL: 1/Z=D/N
056957,000436: 
056958,000437:                                                                                                 #                                           Z INDET. AT PERIGEE FOR PARAB OR HYPERB.
056959,000438: 27,3110           77616        DUMPTFF1          RVQ                                            #  RETURN  TFF =0
056960,000439: 
056961,000440:                                                                                                 #                                          INBOUND  Z  CALC CONTINUES HERE
056962,000441: 27,3111           77745        INBOUND           DLOAD                                          #  RESET PDL+0
056963,000442: 27,3112           45345                          DLOAD    DSU                                   #  ALTERNATE Z CALC
056964,000443: 27,3113           00023                                   RTERM                                 #  E: (-29)  M: (-27)
056965,000444: 27,3114           00015                                   RMAG1                                 #  E: (-29)  M: (-27)
056966,000445: 27,3115           14043                          STODL    TFFX                                  #  NUM=RTERM-RN  E: (-29)  M: (-27)
056967,000446: 27,3116           00045                                   TFFTEM                                #  Q2  E: (-16)  M: (-15)
056968,000447: 27,3117           52015                          DAD      GOTO                                  
056969,000448: 27,3120           00017                                   TFFQ1                                 #  Q1  E: (-16)  M: (-15)
056970,000449: 27,3121           57076                                   SAVEDEN                               #  DEN = Q2+Q1   E: (-16)  M: (-15)
056971,000450: 
056972,000451: 27,3122           65215        TFFXTEST          DAD      PDDL                                  #   (ABS(DEN) TO PDL+2))  E: (-3) OR (-16)
056973,000452:                                                                                                 #                                                                  M: (-3) OR (-15)
056974,000453: 27,3123           17324                                   DP(-22)                               #  RESTORE ABS(DEN) TO MPAC
056975,000454: 27,3124           00043                                   TFFX                                  #  NUM E:(-16) OR (-29)  M:(-15) OR (-27)
056976,000455: 27,3125           53605                          DMP      SR*                                   
056977,000456: 27,3126           00031                                   TFFRTALF                              #  SQRT(ALFA)  E: (10+NA)  M: (9+NA)
056978,000457: 27,3127           57201                                   0          -3,2                       #  X2=-NA
056979,000458: 27,3130           77671                          DDV                                            #  C(MPAC) =NUM SQRT(ALFA)  E:(-3) OR (-16)
056980,000459:                                                                                                 #                                                                    M:(-3) OR (-15)
056981,000460:                                                                                                 #                                           ABS(DEN) FROM PDL+2  E:(-3) OR (-16)
056982,000461:                                                                                                 #                                                                M:(-3) OR (-15)
056983,000462: 27,3131           40145                          DLOAD    BOV                                   #  (THE DLOAD IS SHARED WITH TFFELL)
056984,000463: 27,3132           00043                                   TFFX                                  #  NUM  E: (-16) OR (-29)  M:(-15) OR (-27)
056985,000464: 27,3133           57214                                   TFFELL                                #  USE EQN FOR DELE GEQ 90, LEQ -90
056986,000465: 
056987,000466:                                                                                                 #                                          OTHERWISE, CONTINUE FOR GENERAL CONIC FOR TFF EQN
056988,000467: 
056989,000468: 27,3134           45471                          DDV      STADR                                 
056990,000469:                                                                                                 #                                           DEN FROM PDL+0   E: (-3) OR (-16)
056991,000470:                                                                                                 #                                                            M: (-3) OR (-15)
056992,000471: 27,3135           77732                          STORE    TFFTEM                                #  Z  SAVE FOR SIGN OF SDELF.
056993,000472:                                                                                                 # # Page 1348
056994,000473:                                                                                                 #                                           E: (-13)  M: (-12)
056995,000474: 27,3136           63406                          PUSH     DSQ                                   #  Z TO PDL+0
056996,000475: 27,3137           41206                          PUSH     DMP                                   #  Z SQ TO PDL+2  E: (-26)  M: (-24)
056997,000476: 27,3140           00035                                   TFFNP                                 #  LC P  E: (-38+2NR)  M: (-36+NR)
056998,000477: 27,3141           75261                          SL       SIGN                                  
056999,000478: 27,3142           20206                                   5                                     
057000,000479: 27,3143           00045                                   TFFTEM                                #  AFFIX SIGN FOR SDELF (ENTRY DISPLAY)
057001,000480: 27,3144           14045                          STODL    TFFTEM                                #  P ZSQ  E: (-59+2NR)  M: (-55+2NR)
057002,000481:                                                                                                 #                                           (ARG IS USED IN TFF/TRIG)
057003,000482:                                                                                                 #                                           ZSQ FROM PDL+2   E: (-26)  M: (-24)
057004,000483: 27,3145           41206                          PUSH     DMP                                   #  RESTORE PUSH LOC
057005,000484: 27,3146           00033                                   TFFALFA                               #  ALFA  E: (26-NR)  M: (24-NR)
057006,000485: 27,3147           77657                          SL*                                            
057007,000486: 27,3150           20201                                   0,1                                   #  X1=-NR
057008,000487: 27,3151           00043                          STORE    TFFX                                  #  X
057009,000488: 27,3152           41234                          RTB      DMP                                   
057010,000489: 27,3153           57276                                   T(X)                                  #   POLY
057011,000490:                                                                                                 #                                           ZSQ FROM PDL+2  E: (-26)  M: (-24)
057012,000491: 27,3154           44302                          SR2      BDSU                                  #  2 ZSQ T(X)  E: (-29)  M: (-27)
057013,000492: 27,3155           00023                                   RTERM                                 #  RTERM  E: (-29)  M: (-27)
057014,000493: 27,3156           41215                          DAD      DMP                                   
057015,000494: 27,3157           00015                                   RMAG1                                 #  E: (-29)  M: (-27)
057016,000495:                                                                                                 #                                           Z FROM PDL+0  E: (-13)  M: (-12)
057017,000496: 27,3160           51042                          SR3      BPL                                   #  TFF SQRT(MU)   E: (-45)  M: (-42)
057018,000497: 27,3161           57202                                   ENDTFF                                #  (NO PUSH UP)
057019,000498: 27,3162           75206                          PUSH     SIGN                                  #  TFF SQRT(MU) TO PDL+0
057020,000499: 27,3163           00017                                   TFFQ1                                 #  Q1 FOR GONEPAST TEST
057021,000500: 27,3164           71244                          BPL      DLOAD                                 #  GONE PAST ?
057022,000501: 27,3165           57206                                   NEGTFF                                #  YES. TFF < 0 .
057023,000502: 27,3166           00027                                   TFF1/ALF                              #  1/ALFA  E: (-22-2NA)  M: (-20-2NA)
057024,000503: 27,3167           51076                          DCOMP    BPL                                   #  ALFA > 0 ?
057025,000504: 27,3170           57206                                   NEGTFF                                #  NO. TFF IS NEGATIVE.
057026,000505: 
057027,000506:                                                                                                 #                                          CORRECT FOR ORBITAL PERIOD.
057028,000507: 
057029,000508: 27,3171           77676                          DCOMP                                          #  YES. CORRECT FOR ORB PERIOD.
057030,000509: 27,3172           56205                          DMP      DDV                                   
057031,000510: 27,3173           17320                                   PI/16                                 #  2 PI (-5)
057032,000511: 27,3174           00031                                   TFFRTALF                              #  SQRT(ALFA)  E: (10+NA)  M: (9+NA)
057033,000512: 27,3175           53657                          SL*      SL*                                   
057034,000513: 27,3176           57602                                   0          -4,2                       #  X2=-NA
057035,000514: 27,3177           57602                                   0          -4,2                       
057036,000515: 27,3200           43257                          SL*      DAD                                   
057037,000516: 27,3201           57576                                   0,2                                   
057038,000517:                                                                                                 #                                           TFF SQRT(MU) FROM PDL+0  E:(-45) M:(-42)
057039,000518: 27,3202           40005        ENDTFF            DMP      BOV                                   #  TFF SQRT(MU) IN MPAC     E:(-45) M:(-42)
057040,000519: 27,3203           00037                                   TFF/RTMU                              #  E: (17)  M: (14)
057041,000520: 27,3204           57212                                   MAXTFF                                #  SET POSMAX IF OVFL.
057042,000521: 
057043,000522: 27,3205           77616        DUMPTFF2          RVQ                                            #  RETURN  TFF   (-28) CS IN MPAC.
057044,000523:                                                                                                 # # Page 1349
057045,000524: 27,3206           77745        NEGTFF            DLOAD                                          
057046,000525:                                                                                                 #                                           TFF SQRT(MU)  FROM PDL+0, NEGATIVE.
057047,000526: 27,3207           77650                          GOTO                                           
057048,000527: 27,3210           57202                                   ENDTFF                                
057049,000528: 
057050,000529: 27,3211           77745        MAXTFF1           DLOAD                                          #  RESET PDL
057051,000530: 27,3212           43545        MAXTFF            DLOAD    RVQ                                   
057052,000531: 27,3213           17334                                   NEARONE                               
057053,000532: 
057054,000533:                                                                                                 #  TIME OF FLIGHT ELLIPSE WHEN DEL (ECCENTRIC ANOM) GEQ 90 AND LEQ -90.
057055,000534: 
057056,000535:                                                                                                 #                                           NUM FROM TFFX.  E: (-16) OR (-29)
057057,000536:                                                                                                 #                                                           M: (-15) OR (-27)
057058,000537: 27,3214           77712        TFFELL            SL2                                            #  NUM  E:(-14) OR (-27)  M:(-13) OR (-25)
057059,000538: 27,3215           41465                          BDDV     PUSH                                  #  TEMP SAVE D/N IN PDL+0
057060,000539:                                                                                                 #                                           DEN FROM PDL+0 E:(-3)/( 16) M:(-3)/(-15)
057061,000540:                                                                                                 #                                           N/D TO PDL+0  E: (11)  M: (10)
057062,000541: 27,3216           45345        TFFEL1            DLOAD    DSU                                   #  (ENTER WITH D/N=0 IN PDL+0)
057063,000542: 27,3217           00045                                   TFFTEM                                #  Q2  E: (-16)  M: (-15)
057064,000543: 27,3220           00017                                   TFFQ1                                 #  Q1  E: (-16)  M: (-15)
057065,000544: 27,3221           14013                          STODL    TFFDELQ                               #  Q2-Q1  E: (-16)  M: (-15)
057066,000545:                                                                                                 #                                           D/N FROM PDL+0
057067,000546: 27,3222           77626                          STADR                                          
057068,000547: 27,3223           77732                          STORE    TFFTEM                                #  D/N  E: (11)  M: (10)
057069,000548: 27,3224           53605                          DMP      SL*                                   
057070,000549: 27,3225           00027                                   TFF1/ALF                              #  1/ALFA E: (-22-2NA)  M: (-20-2NA)
057071,000550: 27,3226           57576                                   0,2                                   #  1/ALFA Z  E: (-11-NA)  M: (-10-NA)
057072,000551: 27,3227           41206                          PUSH     DMP                                   #  TO PDL+0
057073,000552: 27,3230           00045                                   TFFTEM                                #  1/Z  E: (11)  M: (10
057074,000553: 27,3231           41057                          SL*      BOVB                                  
057075,000554: 27,3232           57576                                   0,2                                   #  X2= -NA
057076,000555: 27,3233           45707                                   SIGNMPAC                              #  IN CASE X= 1.0, CONTINUE
057077,000556: 27,3234           00043                          STORE    TFFX                                  #  X=1/ALFA ZSQ
057078,000557: 27,3235           41234                          RTB      DMP                                   
057079,000558: 27,3236           57276                                   T(X)                                  #  POLY
057080,000559: 27,3237           00043                                   TFFX                                  
057081,000560: 27,3240           45242                          SR3      DSU                                   
057082,000561: 27,3241           17326                                   DP2(-3)                               
057083,000562: 27,3242           41405                          DMP      PUSH                                  #  2(X T(X)-1) /Z ALFA  E: (-15-NA)
057084,000563:                                                                                                 #                                                                M: (-14-NA)
057085,000564:                                                                                                 #                                           1/ALFA Z FROM PDL+0  E: (-11-NA)
057086,000565:                                                                                                 #                                                                M: (-10-NA)
057087,000566: 27,3243           41345                          DLOAD    DMP                                   #  GET SIGN FOR SDELF
057088,000567: 27,3244           00045                                   TFFTEM                                #  1/Z  E: (11)  M: (10)
057089,000568: 27,3245           00015                                   RMAG1                                 #  E: (-29)  M: (-27)
057090,000569: 27,3246           43312                          SL2      DAD                                   
057091,000570: 27,3247           00017                                   TFFQ1                                 #  Q1  E: (-16)  M: (-15)
057092,000571: 27,3250           14045                          STODL    TFFTEM                                #  (Q1+R 1/Z)  =SGN OF SDELF E:(-16) M:(-15
057093,000572: 27,3251           00035                                   TFFNP                                 #  LC P E: (-38+2NR) M: (-36+2NR)
057094,000573: 27,3252           53605                          DMP      SL*                                   #  CALC FOR ARG FOR TFF/TRIG.
057095,000574:                                                                                                 # # Page 1350
057096,000575: 27,3253           00027                                   TFF1/ALF                              #  1/ALFA   E:(-22-2NA)  M:(-20-2NA)
057097,000576: 27,3254           57575                                   1,2                                   #  X2=-NA
057098,000577: 27,3255           53765                          SIGN     SL*                                   
057099,000578: 27,3256           00045                                   TFFTEM                                #  AFFIX SIGN FOR SDELF
057100,000579: 27,3257           57576                                   0,2                                   
057101,000580: 27,3260           14045                          STODL    TFFTEM                                #  P/ALFA  E:(-59+2NR)  M:(-55+2NR)
057102,000581:                                                                                                 #                                           (ARG FOR USE IN TFF/TRIG)
057103,000582: 27,3261           00027                                   TFF1/ALF                              #  1/ALFA E:(-22-2NA)  M:(-20-2NA)
057104,000583: 27,3262           41366                          SQRT     DMP                                   
057105,000584: 27,3263           17320                                   PI/16                                 #  PI (-4)
057106,000585: 27,3264           77615                          DAD                                            
057107,000586:                                                                                                 #                                           2(XT(X)-1)/Z ALFA FROM PDL E:(-15-NA)
057108,000587:                                                                                                 #                                                                      M:(-14-NA)
057109,000588: 27,3265           45257                          SL*      DSU                                   
057110,000589: 27,3266           57577                                   0          -1,2                       
057111,000590: 27,3267           00013                                   TFFDELQ                               #  Q2-Q1  E: (-16)  M: (-15)
057112,000591: 27,3270           53605                          DMP      SL*                                   
057113,000592: 27,3271           00027                                   TFF1/ALF                              #  1/ALFA  E:(-22-2NA)  M:(-20-2NA)
057114,000593: 27,3272           57601                                   0          -3,2                       
057115,000594: 27,3273           52057                          SL*      GOTO                                  
057116,000595: 27,3274           57602                                   0          -4,2                       
057117,000596: 27,3275           57202                                   ENDTFF                                #  TFF SQRT(MU) IN MPAC E:(-145) M:(-42)
057118,000597: 
057119,000598:                                                                                                 # # Page 1351
057120,000599:                                                                                                 #  PROGRAM NAME:      T(X)                                                 DATE:    01.17.67
057121,000600:                                                                                                 #  MOD NO:  0                                                              LOG SECTION:    TIME OF FREE FALL
057122,000601:                                                                                                 #  MOD BY:  RR BAIRNSFATHER
057123,000602:                                                                                                 #  FUNCTIONAL DESCRIPTION:      THE POLYNOMIAL T(X) IS USED BY TIME OF FLIGHT SUBROUTINES  CALCTFF AND
057124,000603:                                                                                                 #           CALCTPER TO APPROXIMATE THE SERIES
057125,000604:                                                                                                 #                             2     3
057126,000605:                                                                                                 #                  1/3 -X/5 +X /7 -X /9 ...
057127,000606: 
057128,000607:                                                                                                 #           WHERE  X = ALFA Z Z       IF ALFA Z Z  LEQ 1
057129,000608:                                                                                                 #                  X = 1/(ALFA Z Z )  IF ALFA Z Z   G  1
057130,000609: 
057131,000610:                                                                                                 #           ALSO   X IS NEG FOR HYPERBOLIC ORBITS
057132,000611:                                                                                                 #                  X = 0 FOR PARABOLIC ORBITS
057133,000612:                                                                                                 #                  X IS POSITIVE FOR ELLIPTIC ORBITS
057134,000613:                                                                                                 #           FOR FLIGHT 278, THE POLYNOMIAL  T(X) IS FITTED OVER THE RANGE  (0,+1)  AND HAS A MAXIMUM
057135,000614:                                                                                                 #           DEVIATION FROM THE SERIES OF  2 E-5       (T(X) IS A CHEBYCHEV  TYPE FIT AND WAS OBTAINED USING
057136,000615:                                                                                                 #           MAC PROGRAM AUTCURFIT294RRB  AND IS VALID TO THE SAME TOLERANCE OVER THE RANGE  (-.08,+1). )
057137,000616:                                                                                                 #  CALLING SEQUENCE:  RTB
057138,000617:                                                                                                 #                            T(X)
057139,000618:                                                                                                 #           C(MPAC) = X
057140,000619: 
057141,000620:                                                                                                 #  SUBROUTINES CALLED:  NONE
057142,000621:                                                                                                 #  NORMAL EXIT MODE:    TC   DANZIG
057143,000622:                                                                                                 #  ALARMS:  NONE
057144,000623:                                                                                                 #  OUTPUT:  C(MPAC) = T(X)
057145,000624:                                                                                                 #  ERASABLE INITIALIZATION REQUIRED:
057146,000625:                                                                                                 #           C(MPAC)  =  X
057147,000626:                                                                                                 #  DEBRIS:  NONE
057148,000627: 
057149,000628: 27,3276           07171        T(X)              TC       POLY                                  
057150,000629: 27,3277           00004                          DEC      4          B-14                       #  N-1
057151,000630: 27,3300           12525 12525                    2DEC     3.333333333 E-1                        
057152,000631: 27,3302           71463 57703                    2DEC*    -1.999819135 E-1*                       
057153,000632: 27,3304           04423 17645                    2DEC*    1.418148467 E-1*                       
057154,000633: 27,3306           74604 43667                    2DEC*    -1.01310997 E-1*                       
057155,000634: 27,3310           01626 37256                    2DEC*    5.609004986 E-2*                       
057156,000635: 27,3312           77404 52071                    2DEC*    -1.536156925 E-2*                       
057157,000636: 27,3314           06030        ENDT(X)           TC       DANZIG                                
057158,000637: 
057159,000638: 27,3315  27,3314               TCDANZIG          =        ENDT(X)                               
057160,000639: 
057161,000640:                                                                                                 # # Page 1352
057162,000641:                                                                                                 #  TFF CONSTANTS
057163,000642: 32,3776                                          BANK     32                                    
057164,000643: 
057165,000644: 27,2000                                          SETLOC   TOF-FF1                               
057166,000645: 27,2000                                          BANK                                           
057167,000646: 
057168,000647:                                                                                                 #                                                               NOTE _  NOTE _ ADJUSTED MUE FOR NEAR EARTH TRAJ.
057169,000648: 
057170,000649:                                                                                                 #  MUE            =       3.990 815 471 E10  M CUBE/CS SQ
057171,000650:                                                                                                 #  RTMUE          =       1.997702549 E5 B-18*     MODIFIED EARTH MU
057172,000651: 
057173,000652: 27,3315           24775 30424  1/RTMU            2DEC*    .5005750271 E-5        B17*            #  MODIFIED EARTH MU
057174,000653:                                                                                                 #                                                               NOTE _  NOTE _ ADJUSTED MUE FOR NEAR EARTH TRAJ.
057175,000654: 
057176,000655:                                                                                                 #  MUM            =       4.902 778  E8     M CUBE /CS SQ
057177,000656: 
057178,000657:                                                                                                 #  RTMUM          2DEC*   2.21422176 E4 B-18*
057179,000658: 27,3317           06220 37553  PI/16             2DEC     3.141592653 B-4                        
057180,000659: 27,3321           37777 37700  LIM(-22)          2OCT     37777      37700                      #  1.0 -B(-22)
057181,000660: 27,3323           00000 00100  DP(-22)           2OCT     00000      00100                      #  B(-22)
057182,000661: 27,3325           04000 00000  DP2(-3)           2DEC     1          B-3                        
057183,000662: 27,3327           02000 00000  DP2(-4)           2DEC     1          B-4                        #  1/16
057184,000663:                                                                                                 #  RPAD1    2DEC   6373338 B-29    M  (-29)  =20 909 901.57 FT
057185,000664: 
057186,000665: 27,3331  22,3310               RPAD1             =        RPAD                                  
057187,000666: 27,3331           00305 11205  R300K             2DEC     6464778    B-29                       #   (-29) M
057188,000667: 27,3333           37777 37777  NEARONE           2DEC     .999999999                            
057189,000668: 27,3335  26,3331               TFFZEROS          EQUALS   HI6ZEROS                              
057190,000669: 27,3335  26,3321               TFF1/4            EQUALS   HIDP1/4                               

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