Source Code
|
These source-code files were transcribed from a printout in Don Eyles's personal
collection, scanned by archive.org, and financially sponsored by Peter McDermott.
A team of volunteers performed the transcription and proof-reading. The scanned
page images are available at
the Virtual AGC Project website, as well as higher-quality (but much larger)
images at
the Virtual AGC Project's collection in the Internet Archive. Report any problems by creating
"issues" at
the Virtual AGC Project's GitHub Repository. Notations on the program listing read, in part:
YUL SYSTEM FOR AGC: NEW PROGRAM SHEPATIN BY EYLES
DEC 13, 1966
THIS PROGRAM WAS ASSEMBLED AS A VERSION OF REVISION 37
OF PROGRAM SUNBURST BY LEM GROUP
Note that the date is the date of the printout, not the date of the program revision.
|
032761,000002: ## Copyright: Public domain.
032762,000003: ## Filename: ORBITAL_INTEGRATION_PROGRAM.agc
032763,000004: ## Purpose: A section of Sunburst revision 37, or Shepatin revision 0.
032764,000005: ## It is part of an early development version of the software
032765,000006: ## for Apollo Guidance Computer (AGC) on the unmanned Lunar
032766,000007: ## Module (LM) flight Apollo 5. Sunburst 37 was the program
032767,000008: ## upon which Don Eyles's offline development program Shepatin
032768,000009: ## was based; the listing herein transcribed was actually for
032769,000010: ## the equivalent revision 0 of Shepatin.
032770,000011: ## This file is intended to be a faithful transcription, except
032771,000012: ## that the code format has been changed to conform to the
032772,000013: ## requirements of the yaYUL assembler rather than the
032773,000014: ## original YUL assembler.
032774,000015: ## Reference: pp. 723-745
032775,000016: ## Assembler: yaYUL
032776,000017: ## Contact: Ron Burkey <info@sandroid.org>.
032777,000018: ## Website: www.ibiblio.org/apollo/index.html
032778,000019: ## Mod history: 2017-05-24 MAS Created from Sunburst 120.
032779,000020: ## 2017-07-12 MAS Updated for Sunburst 37.
032780,000021: ## 2017-06-23 RSB Proofed comment text with
032781,000022: ## octopus/ProoferComments.
032782,000023:
|
|
Page 723 |
032784,000025: # FBR3 SETS UP A TIMESTEP CALL TO KEPLER.
032785,000026:
032786,000027:
032787,000028:
032788,000029: 07,2477 50545 FBR3 DLOAD SR3
032789,000030: 07,2500 02165 H
032790,000031: 07,2501 43262 SR3R DAD
032791,000032: 07,2502 02045 TC
032792,000033: 07,2503 16161 STODL TAU
032793,000034:
032794,000035: 07,2504 17703 EARTHTAB +9D
032795,000036: 07,2505 54205 DMP SRR
032796,000037: 07,2506 02163 DT/2
032797,000038: 07,2507 21615 12D
032798,000039: 07,2510 77615 DAD
032799,000040: 07,2511 02047 TET
032800,000041: 07,2512 02047 STORE TET
032801,000042:
|
|
Page 724 |
032803,000044: # THIS ORBITAL KEPLER SUBROUTINE FINDS THE POSITION AND VELOCITY OF THE VEHICLE AFTER TIME FOUND IN GIVENT
032804,000045: # SINCE RECTIFICATION TO POSITION RRECT AND VELOCITY VRECT. THE RESULTING POSITION AND VELOCITY ARE LEFT IN
032805,000046:
032806,000047: # FOUNDR AND FOUNDV, RESPECTIVELY.
032807,000048:
032808,000049:
032809,000050:
032810,000051: 07,2513 40375 KEPLER VLOAD SETPD # UNIT OF RECTIFICATION POSITION TO 0
032811,000052: 07,2514 02001 RRECT
032812,000053: 07,2515 00001 0
032813,000054: 07,2516 65256 UNIT PDDL # AND LENGTH OF ORIGINAL IN 6
032814,000055: 07,2517 00045 36D
032815,000056: 07,2520 77715 PDVL # LENGTH OF POSITION AT RECTIFICATION.
032816,000057: 07,2521 02007 VRECT
032817,000058: 07,2522 61436 VSQ ROUND
032818,000059: 07,2523 45205 DMP DSU
032819,000060: 07,2524 00007 6 # LENGTH OF POSITION AT RECTIFICATION
032820,000061: 07,2525 17630 DP1/4
032821,000062: 07,2526 41532 SL2R PUSH # A4 TO REGISTER 8
032822,000063: 07,2527 57502 SR2 DCOMP
032823,000064: 07,2530 56215 DAD DDV
032824,000065: 07,2531 17630 DP1/4
032825,000066: 07,2532 00007 6
032826,000067: 07,2533 50315 PDVL DOT # ALPHA TO REGISTER 10
032827,000068: 07,2534 02001 RRECT
032828,000069: 07,2535 02007 VRECT
032829,000070:
032830,000071: 07,2536 63306 ROUND PDVL # A1 TO REGISTER 12
032831,000072: 07,2537 02031 RCV
032832,000073: 07,2540 50256 UNIT DOT
032833,000074: 07,2541 02037 VCV
032834,000075: 07,2542 77725 PDDL # IR/2.VC IN 14D.
032835,000076: 07,2543 00045 36D
032836,000077: 07,2544 16155 STODL ALPHAM # RC IN ALPHAM
032837,000078: 07,2545 02163 DT/2
032838,000079: 07,2546 54442 SR3 SR3R
032839,000080: 07,2547 77071 DDV AXT,2
032840,000081: 07,2550 02155 ALPHAM
032841,000082: 07,2551 00012 10D # MAX ITERATION COUNT IS 10
032842,000083:
|
|
Page 725 |
032844,000085: 07,2552 56325 PDDL DDV # Q IN 16
032845,000086: 07,2553 17630 DP1/4
032846,000087: 07,2554 02155 ALPHAM
032847,000088: 07,2555 57225 DSU DMPR
032848,000089: 07,2556 00013 10D # 1/4RC : ALPHA
032849,000090: 07,2557 00021 16D # Q( )
032850,000091: 07,2560 57275 DMPR DMPR
032851,000092: 07,2561 00021 16D # QQ( )
032852,000093:
032853,000094: 07,2562 17632 DP1/3
032854,000095: 07,2563 52452 SL3 SL3
032855,000096: 07,2564 57325 PDDL DMPR
032856,000097: 07,2565 00017 14D
032857,000098: 07,2566 00021 16D
032858,000099: 07,2567 41412 SL4 PUSH # 16(UR/2.VC)Q IN 20
032859,000100: 07,2570 57221 BDSU DMPR
032860,000101: 07,2571 17656 DP1/2
032861,000102: 07,2572 45221 BDSU DSU
032862,000103: 07,2573 17656 DP1/2
032863,000104: 07,2574 76405 DMP SL1R
032864,000105: 07,2575 66215 DAD SSP
032865,000106: 07,2576 02051 XKEP
032866,000107: 07,2577 00052 S2
032867,000108: 07,2600 00001 1
032868,000109: 07,2601 02051 STORE XKEP
032869,000110:
|
|
Page 726 |
032871,000112: # ITERATING EQUATIONS - GIVEN X IN MPAC, FIND TIME OF FLIGHT.
032872,000113:
032873,000114:
032874,000115:
032875,000116: 07,2602 61516 KTIMEN+1 DSQ ROUND # FORM ALPHA X-SQUARED AND CALL S AND C
032876,000117: 07,2603 66405 DMP SL2R
032877,000118: 07,2604 00013 10D
032878,000119: 07,2605 45001 SETPD CALL # SET PD INDICATOR TO 16
032879,000120: 07,2606 00021 16D
032880,000121: 07,2607 16677 S(X)C(X)
032881,000122: 07,2610 42405 DMP SL4
032882,000123: 07,2611 02051 XKEP
032883,000124: 07,2612 72405 DMP SL1
032884,000125: 07,2613 02051 XKEP
032885,000126:
032886,000127: 07,2614 76405 DMP SL1R
032887,000128: 07,2615 02051 XKEP
032888,000129: 07,2616 00030 STORE 23D # A3
032889,000130: 07,2617 65275 DMPR PDDL
032890,000131: 07,2620 00011 8D
032891,000132: 07,2621 02051 XKEP
032892,000133: 07,2622 42405 DMP SL4
032893,000134: 07,2623 00021 16D # VALUE OF C
032894,000135: 07,2624 41352 SL1 DMP
032895,000136: 07,2625 02051 XKEP
032896,000137: 07,2626 77732 SL2R
032897,000138: 07,2627 00026 STORE 21D # A2
032898,000139: 07,2630 74405 DMP SR1R
032899,000140: 07,2631 00015 12D # A1
032900,000141: 07,2632 77615 DAD
032901,000142: 07,2633 57325 PDDL DMPR
032902,000143: 07,2634 00007 6
032903,000144: 07,2635 02051 XKEP
032904,000145: 07,2636 77615 DAD
032905,000146: 07,2637 44206 PUSH BDSU # COMPARE WITH GIVEN TIME OF FLIGHT
032906,000147: 07,2640 02161 GIVENT
032907,000148: 07,2641 00021 STORE 16D # DIFFERENCE TO REGISTER 16
032908,000149: 07,2642 77776 EXIT
032909,000150:
|
|
Page 727 |
032911,000152: 07,2643 06103 DUMPDUMP TC INTPRET # FOR DUMP ONLY *******
032912,000153: 07,2644 45246 ABS DSU
032913,000154: 07,2645 17622 KEPSILON # SEE IF WITHIN EPSILON OF GIVEN TIME.
032914,000155: 07,2646 61040 BMN TIX,2 # IF SO, GET R AND V AND EXIT.
032915,000156: 07,2647 16764 GETRANDV
032916,000157: 07,2650 16653 GETNEWX
032917,000158: 07,2651 77650 GOTO
032918,000159: 07,2652 16764 GETRANDV
032919,000160: 07,2653 41345 GETNEWX DLOAD DMP
032920,000161: 07,2654 00013 10D # ALPHA
032921,000162: 07,2655 00030 23D # A3
032922,000163: 07,2656 44332 SL2R BDSU
032923,000164: 07,2657 02051 XKEP
032924,000165: 07,2660 46405 DMP SL4R
032925,000166: 07,2661 00015 12D # A1
032926,000167: 07,2662 14023 STODL 18D
032927,000168: 07,2663 00026 21D # A2
032928,000169: 07,2664 76405 DMP SL1R
032929,000170: 07,2665 00011 8D # A4
032930,000171:
032931,000172: 07,2666 77615 DAD
032932,000173: 07,2667 65215 DAD PDDL
032933,000174: 07,2670 00007 6 # R0
032934,000175: 07,2671 00021 16D
032935,000176: 07,2672 43271 DDV DAD
032936,000177: 07,2673 00023 18D
032937,000178: 07,2674 02051 XKEP
032938,000179: 07,2675 36051 STCALL XKEP
032939,000180: 07,2676 16602 KTIMEN+1
032940,000181:
|
|
Page 728 |
032942,000183: # SUBROUTINE FOR COMPUTING THE UNIVERSAL CONIC FUNCTIONS S(X) AND C(X). THE ACTUAL OUTPUT OF THIS ROUTINE
032943,000184: # CONSISTS OF SCALED VERSIONS DEFINED AS FOLLOWS -
032944,000185:
032945,000186: # S (X) = S(64X) C (X) = C(64X)/4
032946,000187: # S S
032947,000188:
032948,000189: # IT IS ASSUMED THAT THE INPUT ARRIVES IN MPAC,MPAC+1 AND THAT IT LIES BETWEEN -30/64 AND 40/64. UPON EXIT,
032949,000190: # S(X) WILL BE LEFT IN MPAC,MPAC+1 AND C(X) ON TOP OF THE PUSHDOWN LIST.
032950,000191:
032951,000192:
032952,000193:
032953,000194: 07,2677 00043 S(X)C(X) STORE 34D # X TO 34D
032954,000195: 07,2700 63434 RTB DSQ
032955,000196: 07,2701 16726 A(X)
032956,000197: 07,2702 41506 ROUND PUSH
032957,000198: 07,2703 00045 STORE 36D # A SQUARED TO 36D
032958,000199: 07,2704 76405 DMP SL1R
032959,000200: 07,2705 00043 34D # 2 2
032960,000201: 07,2706 57221 BDSU DMPR # C (X) = A (.25 - 2XA ) TO PD LIST
032961,000202: 07,2707 17630 POS1/4 # S
032962,000203: 07,2710 00045 36D
032963,000204: 07,2711 64525 PDDL SR2R # 2
032964,000205: 07,2712 47125 PDDL RTB # A /4 TO PD LIST
032965,000206: 07,2713 00043 34D
032966,000207: 07,2714 16745 B(X)
032967,000208: 07,2715 61516 DSQ ROUND
032968,000209:
032969,000210: 07,2716 57325 PDDL DMPR # B SQUARED TO PD LIST
032970,000211: 07,2717 00043 34D
032971,000212: 07,2720 00045 36D
032972,000213: 07,2721 57221 BDSU DMPR # 2 2 2
032973,000214: 07,2722 17634 POS1/16 # LEAVES S (X)=B (.0625-A X)+A /4 IN MPAC
032974,000215: 07,2723 43415 DAD ITCQ # S
032975,000216: 07,2724 74631 63145 -1/12 2DEC -.1 # DONT MOVE.
032976,000217:
|
|
Page 729 |
032978,000219: # A AND B POLYNOMIALS WHOSE COEFFICIENTS WERE OBTAINED WITH THE *AUTOCURVEFIT* PROGRAM.
032979,000220: 07,2726 07265 A(X) TC POLY
032980,000221: 07,2727 00004 DEC 4 B-14
032981,000222: 07,2730 26501 07463 2DEC 7.071067810 E-1
032982,000223:
032983,000224: 07,2732 60724 60210 2DEC -4.714045180 E-1
032984,000225:
032985,000226: 07,2734 03010 26256 2DEC 9.42808914 E-2
032986,000227:
032987,000228: 07,2736 77554 74242 2DEC -8.9791893 E-3
032988,000229:
032989,000230: 07,2740 00010 05475 2DEC 4.989987 E-4
032990,000231:
032991,000232: 07,2742 77777 66460 2DEC -1.79357 E-5
032992,000233:
032993,000234: 07,2744 06115 TC DANZIG
032994,000235: 07,2745 07265 B(X) TC POLY
032995,000236: 07,2746 00004 DEC 4 B-14
032996,000237: 07,2747 32101 17270 2DEC 8.164965793 E-1
032997,000238:
032998,000239: 07,2751 65431 40174 2DEC -3.265986572 E-1
032999,000240:
033000,000241: 07,2753 01710 10660 2DEC 5.90988980 E-2
033001,000242:
033002,000243: 07,2755 77676 52270 2DEC -4.0085592 E-3
033003,000244:
033004,000245: 07,2757 00004 21652 2DEC 2.781528 E-4
033005,000246:
033006,000247: 07,2761 77777 71323 2DEC -1.25610 E-5
033007,000248:
033008,000249: 07,2763 06115 TC DANZIG
033009,000250:
|
|
Page 730 |
033011,000252: # ROUTINE FOR OBTAINING R AND V, NOW THAT THE PROPER X HAS BEEN FOUND.
033012,000253:
033013,000254:
033014,000255:
033015,000256: 07,2764 40345 GETRANDV DLOAD SETPD
033016,000257: 07,2765 00026 21D # A2 FROM LAST ITERATION
033017,000258: 07,2766 00032 25D
033018,000259: 07,2767 74276 DCOMP VXSC
033019,000260: 07,2770 00001 0 # UNIT OF GIVEN POSITION VECTOR
033020,000261: 07,2771 45325 PDDL DSU
033021,000262: 07,2772 00023 18D # LAST VALUE OF T
033022,000263: 07,2773 00030 23D # LAST VALUE OF A3
033023,000264: 07,2774 74312 SL2 VXSC
033024,000265: 07,2775 02007 VRECT
033025,000266:
033026,000267: 07,2776 76455 VAD VSL1 # ADDITION MUST BE DONE IN THIS ORDER
033027,000268: 07,2777 53206 PUSH VAD
033028,000269: 07,3000 02001 RRECT
033029,000270: 07,3001 45455 VAD STADR
033030,000271: 07,3002 75746 STORE FOUNDR # RESULTING CONIC POSITION
033031,000272: 07,3003 77646 ABVAL
033032,000273: 07,3004 14021 STODL 16D
033033,000274: 07,3005 00013 10D # ALPHA
033034,000275: 07,3006 66405 DMP SL2R
033035,000276: 07,3007 00030 23D # A3
033036,000277:
033037,000278: 07,3010 56225 DSU DDV
033038,000279: 07,3011 02051 XKEP
033039,000280: 07,3012 00021 16D # LENGTH OF FOUND POSITION
033040,000281: 07,3013 72561 VXSC VSL2
033041,000282: 07,3014 00001 0 # UNIT OF RECTIFICATION POSITION
033042,000283: 07,3015 70525 PDDL SR1
033043,000284: 07,3016 00021 16D
033044,000285: 07,3017 56225 DSU DDV
033045,000286: 07,3020 00026 21D
033046,000287: 07,3021 00021 16D
033047,000288: 07,3022 53361 VXSC VAD
033048,000289: 07,3023 02007 VRECT
033049,000290:
033050,000291: 07,3024 77772 VSL1
033051,000292: 07,3025 36037 STCALL FOUNDV # THIS COMPLETES THE CALCULATION
033052,000293: 07,3026 02171 HBRANCH
033053,000294:
|
|
Page 731 |
033055,000296: # THE POSTRUE ROUTINES SET UP THE BETA VECTOR AND OTHER INITIAL CONDITIONS FOR THE NEXT ACCOMP.
033056,000297: 07,3027 77331 POSTRUE SSP VLOAD # TIME STEP CALLS TO KEPLER RETURN HERE
033057,000298: 07,3030 02207 SCALEA
033058,000299: 07,3031 00004 4
033059,000300: 07,3032 02053 ALPHAV
033060,000301: 07,3033 53261 VSR VAD
033061,000302:
033062,000303: 07,3034 20613 10D
033063,000304: 07,3035 02031 RCV # POSITION OUTPUT OF KEPLER
033064,000305: 07,3036 43154 LXA,2 BOF
033065,000306: 07,3037 02205 DIFEQCNT
033066,000307: 07,3040 00340 WMATFLAG
033067,000308: 07,3041 17043 NOSAVE1
033068,000309: 07,3042 12111 STORE VECTAB,2
033069,000310: 07,3043 66331 NOSAVE1 SSP SSP # SETS UP SCALE B AND GMODE
033070,000311: 07,3044 02210 SCALEB
033071,000312: 07,3045 00016 14D
033072,000313: 07,3046 02173 GMODE
033073,000314: 07,3047 00002 2
033074,000315:
033075,000316: 07,3050 02061 STORE BETAV
033076,000317:
|
|
Page 732 |
033078,000319: # AGC ROUTINE TO COMPUTE ACCELERATION COMPONENTS.
033079,000320:
033080,000321:
033081,000322:
033082,000323: 07,3051 53575 ACCOMP VLOAD UNIT
033083,000324: 07,3052 02053 ALPHAV
033084,000325: 07,3053 16053 STODL ALPHAV
033085,000326: 07,3054 00045 36D
033086,000327: 07,3055 02155 STORE ALPHAM
033087,000328: 07,3056 77600 BOV
033088,000329: 07,3057 17060 ACCOMP2 # TURN OFF OVERFLOW INDICATOR
033089,000330: 07,3060 74575 ACCOMP2 VLOAD VSR1
033090,000331: 07,3061 02061 BETAV
033091,000332: 07,3062 40236 VSQ SETPD
033092,000333: 07,3063 00001 0
033093,000334: 07,3064 61501 NORM ROUND
033094,000335: 07,3065 00051 S1
033095,000336: 07,3066 60325 PDDL NORM # NORMED B SQUARED TO PD LIST
033096,000337: 07,3067 02155 ALPHAM # NORMALIZE (LESS ONE) LENGTH OF ALPHA
033097,000338: 07,3070 00047 X1 # SAVING NORM SCALE FACTOR IN X1
033098,000339: 07,3071 63342 SR1 PDVL
033099,000340: 07,3072 02061 BETAV # C(PDL+2) = ALMOST NORMED ALPHA
033100,000341: 07,3073 77656 UNIT
033101,000342: 07,3074 16061 STODL BETAV
033102,000343: 07,3075 00045 36D
033103,000344: 07,3076 02157 STORE BETAM
033104,000345: 07,3077 55301 NORM BDDV # FORM NORMALIZED QUOTIENT ALPHAM/BETAM
033105,000346: 07,3100 00050 X2
033106,000347: 07,3101 41562 SR1R PUSH # C(PDL+2) = ALMOST NORMALIZED RHO.
033107,000348: 07,3102 57144 LXC,2 XAD,2
033108,000349: 07,3103 00047 X2 # C(X2) = -SCALE(RHO) + 1
033109,000350: 07,3104 02206 SCALEA # = -S(B)-N(B)+S(A)+N(A)+1
033110,000351: 07,3105 55074 XAD,2 XSU,2
033111,000352: 07,3106 00046 X1
033112,000353: 07,3107 02207 SCALEB
033113,000354: 07,3110 53714 INCR,2 SR*
033114,000355:
033115,000356: 07,3111 00002 2
033116,000357: 07,3112 57176 0,2
033117,000358: 07,3113 64406 PUSH SR2R # RHO/4 PD+6
033118,000359: 07,3114 50315 PDVL DOT
033119,000360: 07,3115 02053 ALPHAV
033120,000361: 07,3116 02061 BETAV
033121,000362: 07,3117 44372 SL1R BDSU # (RHO/4) - 2(ALPHAV/2.BETAV/2)
033122,000363: 07,3120 57206 PUSH DMPR # TO PDL+6
033123,000364: 07,3121 00005 4
033124,000365:
|
|
Page 733 |
033126,000367: 07,3122 43206 PUSH DAD # Q/4 = RHO(C(PDL+4)) TO PD+8D
033127,000368: 07,3123 17630 DQUARTER # (Q+1)/4 TO PD+10D
033128,000369: 07,3124 75406 PUSH SQRT # 3/2
033129,000370: 07,3125 41475 DMPR PUSH # ((Q+1)/4) TO PD+12D
033130,000371: 07,3126 00013 10D
033131,000372: 07,3127 43352 SL1 DAD
033132,000373: 07,3130 17630 DQUARTER # 3/2
033133,000374: 07,3131 43325 PDDL DAD # (1/4)+2((Q+1)/4) TO PD+14D
033134,000375:
033135,000376: 07,3132 00013 10D
033136,000377: 07,3133 17656 DP1/2
033137,000378: 07,3134 72475 DMPR SL1
033138,000379: 07,3135 00011 8D
033139,000380: 07,3136 56215 DAD DDV
033140,000381: 07,3137 17624 THREE/8
033141,000382: 07,3140 00017 14D
033142,000383: 07,3141 74275 DMPR VXSC
033143,000384: 07,3142 00007 6
033144,000385: 07,3143 02061 BETAV # -
033145,000386: 07,3144 64515 PDVL VSR3 # (G/2)(C(PD+4))B/2 TO PD+16D
033146,000387: 07,3145 02053 ALPHAV
033147,000388: 07,3146 41455 VAD PUSH # A12 + C(PD+16D) TO PD+16D
033148,000389: 07,3147 41345 DLOAD DMP
033149,000390: 07,3150 00001 0
033150,000391: 07,3151 00015 12D # -
033151,000392: 07,3152 61501 NORM ROUND # GAMMA TO PD+22D
033152,000393: 07,3153 00052 S2 # - SCALE(GAMMA)-1 TO X1
033153,000394: 07,3154 70065 BDDV LXC,1
033154,000395: 07,3155 00003 2
033155,000396: 07,3156 00047 X2 # C(X2) = SCALE(RHO)
033156,000397: 07,3157 56070 XAD,1 XAD,1
033157,000398: 07,3160 00051 S2 # C(S2) = N((B.B/4)(....)3/2)
033158,000399:
033159,000400: 07,3161 00050 S1 # C(S1) = N(B.B/4)
033160,000401: 07,3162 56070 XAD,1 XAD,1
033161,000402: 07,3163 02207 SCALEB
033162,000403: 07,3164 02207 SCALEB
033163,000404: 07,3165 74276 DCOMP VXSC
033164,000405: 07,3166 00021 16D # RESULT OF PRECEDING EQUATION
033165,000406: 07,3167 73206 PUSH CGOTO
033166,000407: 07,3170 02173 GMODE
033167,000408: 07,3171 17172 GTABLE
033168,000409: 07,3172 17243 GTABLE CADR GMODE10
033169,000410:
033170,000411: 07,3173 17224 CADR GMODE11
033171,000412: 07,3174 17175 CADR GMODE12
033172,000413:
|
|
Page 734 |
033174,000415: # THE GMODE12 ROUTINE SETS UP THE SECONDARY BODY DISTURBING ACCELERATION FOR ACCOMP.
033175,000416:
033176,000417: 07,3175 77657 GMODE12 VSL* # -SCALE(GAMMA)-1 IS LEFT IN X1.
033177,000418: 07,3176 20240 31D,1 # ADJUST GAMMA TO SCALE OF -32
033178,000419: 07,3177 26103 STOVL FV
033179,000420: 07,3200 02061 BETAV
033180,000421: 07,3201 16053 STODL ALPHAV # BETA VECTOR INTO ALPHA FOR NEXT ACCOMP
033181,000422: 07,3202 02157 BETAM
033182,000423: 07,3203 02155 STORE ALPHAM
033183,000424: 07,3204 45014 BOFF CALL
033184,000425: 07,3205 00342 MIDFLAG
033185,000426: 07,3206 17250 OBLATEST
033186,000427: 07,3207 17663 MOONPOS
033187,000428: 07,3210 02061 STORE BETAV # MOON(EARTH) POSITION WILL BE BETA NEXT
033188,000429:
033189,000430: 07,3211 43150 LXA,1 BOF
033190,000431: 07,3212 02205 DIFEQCNT # SAVE R/QV IN VECTAB FOR W-MATRIX UPDATE
033191,000432: 07,3213 00340 WMATFLAG
033192,000433: 07,3214 17216 NOSAVE2
033193,000434: 07,3215 06117 STORE VECTAB +6,1
033194,000435: 07,3216 65174 NOSAVE2 AXT,2 XCHX,2 # SETUP ALPHAM AND SCALEA
033195,000436: 07,3217 00023 19D # SCALE FOR R/QV
033196,000437: 07,3220 02207 SCALEB # SWAP SCALEB AND X2
033197,000438: 07,3221 52134 SXA,2 GOTO
033198,000439: 07,3222 02206 SCALEA
033199,000440: 07,3223 17060 ACCOMP2 # ENTRY IF UNIT(ALPHAV) AVAILABLE
033200,000441:
|
|
Page 735 |
033202,000443: # THE GMODE11 ROUTINE SETS UP THE SUNS DISTURBING ACCELERATION.
033203,000444:
033204,000445:
033205,000446:
033206,000447: 07,3224 45144 GMODE11 LXC,2 CALL # SET X2 TO TABLE OF PROPER CONSTANTS
033207,000448: 07,3225 02227 PBODY
033208,000449:
033209,000450: 07,3226 17362 ADDTOFV
033210,000451: 07,3227 77624 CALL # BARICENTER-TO-SUN POSITION VECTOR.
033211,000452: 07,3230 17665 SUNPOS # LEAVES VECTOR IN PDL
033212,000453: 07,3231 76750 LXA,1 VLOAD* # COMPUTE R/PS USING CORRECT TABLE FOR
033213,000454: 07,3232 02205 DIFEQCNT # MASS RATIO, ETC.
033214,000455: 07,3233 02117 VECTAB +6,1
033215,000456: 07,3234 53357 VXSC* VAD
033216,000457: 07,3235 77770 6,2 # USE SCALAR AT ENTRY 6 IN THE TABLE
033217,000458: 07,3236 66170 AXT,1 SXA,1
033218,000459: 07,3237 00034 28D
033219,000460: 07,3240 02207 SCALEB # SET SCALEB AND RETURN TO ACCOMP
033220,000461: 07,3241 36061 STCALL BETAV
033221,000462:
033222,000463: 07,3242 17060 ACCOMP2
033223,000464:
033224,000465: # THE GMODE10 ROUTINE ADDS IN THE SUNS PERTURBING ACCELERATION AND COMPUTES THE OBLATENESS CONTRIBUTION
033225,000466: 07,3243 63144 GMODE10 LXC,2 INCR,2
033226,000467: 07,3244 02227 PBODY
033227,000468: 07,3245 77774 DEC -3 B-14
033228,000469: 07,3246 77624 CALL
033229,000470: 07,3247 17362 ADDTOFV
033230,000471: 07,3250 77614 OBLATEST BON
033231,000472: 07,3251 00303 MOONFLAG
033232,000473:
033233,000474: 07,3252 17350 NBRANCH
033234,000475:
|
|
Page 736 |
033236,000477: # THE OBLATE ROUTINE COMPUTES THE ACCELERATION DUE TO THE EARTHS OBLATENESS. 2T USES THE UNIT OF THE VEHICLE
033237,000478: # POSITION VECTOR FOUND IN ALPHAV AND THE DISTANCE TO THE CENTER IN ALPHAM. THIS IS ADDED TO THE SUM OF THE
033238,000479:
033239,000480: # DISTURBING ACCELERATIONS IN FV AND THE PROPER DIFEQ STAGE IS CALLED VIA X1.
033240,000481: 07,3253 77745 OBLATE DLOAD
033241,000482: 07,3254 02057 ALPHAV +4 # Z COMPONENT OF POSITION IS COS PHI
033242,000483: 07,3255 57201 SETPD DMPR
033243,000484: 07,3256 00001 0
033244,000485: 07,3257 17640 3/4
033245,000486: 07,3260 63525 PDDL DSQ # P2:/8 TO REGISTER 0
033246,000487: 07,3261 02057 ALPHAV +4
033247,000488: 07,3262 57252 SL3 DMPR
033248,000489: 07,3263 17636 15/16
033249,000490: 07,3264 41425 DSU PUSH # P3:/4 TO REGISTER 2
033250,000491: 07,3265 17624 3/8
033251,000492: 07,3266 57275 DMPR DMPR
033252,000493: 07,3267 02057 ALPHAV +4
033253,000494: 07,3270 17644 7/12
033254,000495: 07,3271 65352 SL1 PDDL # P4:/16 TO REGISTER 4
033255,000496: 07,3272 00001 0
033256,000497: 07,3273 44275 DMPR BDSU
033257,000498: 07,3274 17662 2/3
033258,000499: 07,3275 57206 PUSH DMPR # BEGIN COMPUTING P5:/128
033259,000500: 07,3276 02057 ALPHAV +4
033260,000501: 07,3277 65275 DMPR PDDL
033261,000502: 07,3300 17646 9/16
033262,000503:
033263,000504: 07,3301 00003 2
033264,000505: 07,3302 44275 DMPR BDSU # FINISH P5:/128 AND TERM USING UNIT
033265,000506: 07,3303 17650 5/128 # POSITION VECTOR AT ALPHA
033266,000507: 07,3304 62405 DMP SL2
033267,000508: 07,3305 17654 J4REQ/J3
033268,000509: 07,3306 43271 DDV DAD
033269,000510: 07,3307 02155 ALPHAM
033270,000511: 07,3310 00005 4
033271,000512: 07,3311 56275 DMPR DDV
033272,000513: 07,3312 17652 2J3RE/J2
033273,000514:
033274,000515: 07,3313 02155 ALPHAM
033275,000516: 07,3314 74215 DAD VXSC
033276,000517: 07,3315 00003 2
033277,000518: 07,3316 02053 ALPHAV
033278,000519: 07,3317 16053 STODL ALPHAV
033279,000520:
|
|
Page 737 |
033281,000522: 07,3320 72405 DMP SL1 # COMPUTE TERM USING IZ
033282,000523: 07,3321 17654 J4REQ/J3
033283,000524: 07,3322 43271 DDV DAD
033284,000525: 07,3323 02155 ALPHAM
033285,000526: 07,3324 50525 PDDL SR3
033286,000527: 07,3325 17652 2J3RE/J2
033287,000528: 07,3326 77675 DMPR
033288,000529: 07,3327 43271 DDV DAD
033289,000530:
033290,000531: 07,3330 02155 ALPHAM
033291,000532: 07,3331 77621 BDSU
033292,000533: 07,3332 02057 ALPHAV +4
033293,000534: 07,3333 16057 STODL ALPHAV +4
033294,000535: 07,3334 02155 ALPHAM
033295,000536: 07,3335 63516 DSQ DSQ
033296,000537: 07,3336 55301 NORM BDDV
033297,000538: 07,3337 00047 X1
033298,000539: 07,3340 17642 J2REQSQ
033299,000540: 07,3341 62161 VXSC INCR,1
033300,000541:
033301,000542: 07,3342 02053 ALPHAV
033302,000543: 07,3343 00004 4
033303,000544: 07,3344 53257 VSL* VAD # SHIFTS LEFT ON +, RIGHT ON -.
033304,000545: 07,3345 20201 0,1
033305,000546: 07,3346 02103 FV
033306,000547: 07,3347 02103 STORE FV
033307,000548: 07,3350 72135 NBRANCH SLOAD LXA,1
033308,000549: 07,3351 02206 DIFEQCNT
033309,000550: 07,3352 00144 MPAC
033310,000551: 07,3353 73205 DMP CGOTO
033311,000552: 07,3354 16725 -1/12
033312,000553: 07,3355 00145 MPAC
033313,000554:
033314,000555: 07,3356 17357 DIFEQTAB
033315,000556: 07,3357 17454 DIFEQTAB CADR DIFEQ+0
033316,000557: 07,3360 17460 CADR DIFEQ+1
033317,000558: 07,3361 17471 CADR DIFEQ+2
033318,000559: 07,3362 77743 ADDTOFV DLOAD* # SETS UP S1 AND S2 PER PRIMARY BODY TABLE
033319,000560: 07,3363 77776 0,2
033320,000561: 07,3364 24051 STOVL S1
033321,000562: 07,3365 00027 22D
033322,000563: 07,3366 73670 XAD,1 VXSC*
033323,000564: 07,3367 00050 S1
033324,000565:
033325,000566: 07,3370 77775 1,2
033326,000567: 07,3371 53257 VSL* VAD
033327,000568: 07,3372 20240 31D,1
033328,000569: 07,3373 02103 FV
033329,000570: 07,3374 02103 STORE FV
033330,000571: 07,3375 77616 ITCQ
033331,000572:
|
|
Page 738 |
033333,000574: # BEGIN INTEGRATION STEP WITH RECTIFICATION TEST.
033334,000575:
033335,000576:
033336,000577:
033337,000578: 07,3376 77775 TIMESTEP VLOAD # MOVE TEMPORARY DELTA AND NU VECTORS
033338,000579: 07,3377 02023 TNUV # TO WORKING STORAGE
033339,000580: 07,3400 26222 STOVL ZV
033340,000581: 07,3401 02015 TDELTAV
033341,000582: 07,3402 02214 STORE YV
033342,000583: 07,3403 45246 ABVAL DSU # RECTIFICATION REQUIRED IF THE LENGTH
033343,000584: 07,3404 17656 DP1/2 # OF DELTA IS GREATER THAN .5 (8KM).
033344,000585: 07,3405 45040 BMN CALL
033345,000586: 07,3406 17410 INTGRATE
033346,000587: 07,3407 17426 RECTIFY
033347,000588:
033348,000589: 07,3410 66331 INTGRATE SSP SSP # INITIALIZE INDICES AND SWITCHES
033349,000590: 07,3411 02171 FBRANCH # EXIT FROM DIFEQCOM
033350,000591: 07,3412 16477 CADR FBR3
033351,000592: 07,3413 02172 HBRANCH # EXIT FROM KEPLER
033352,000593: 07,3414 17027 CADR POSTRUE
033353,000594: 07,3415 77614 CLEAR
033354,000595: 07,3416 00261 JSWITCH # 1 FOR W MATRIX EXTRAPOLATE, 0 OTHERWISE.
033355,000596: 07,3417 77775 DIFEQ0 VLOAD # POSITION DEVIATION INTO ALPHA
033356,000597: 07,3420 02214 YV
033357,000598: 07,3421 16053 STODL ALPHAV
033358,000599: 07,3422 17660 DPZERO
033359,000600: 07,3423 02165 STORE H # START H AT ZERO. GOES O(DELT/2)DELT.
033360,000601: 07,3424 36206 STCALL DIFEQCNT # ZERO DIFEQCNT AND REGISTER FOLLOWING.
033361,000602: 07,3425 02171 HBRANCH # GOES 0(-12D)(-24D).
033362,000603:
|
|
Page 739 |
033364,000605: # THE RECTIFY SUBROUTINE IS CALLED BY THE INTEGRATION PROGRAM AND OCCASIONALLY BY THE MEASUREMENT INCORPORATION
033365,000606: # ROUTINES TO ESTABLISH A NEW CONIC.
033366,000607:
033367,000608:
033368,000609:
033369,000610: 07,3426 40575 RECTIFY VLOAD VSR8 # RECTIFY - FORM TOTAL POSITION AND VEL.
033370,000611: 07,3427 02015 TDELTAV # ADJUST SCALE DIFFERENCE (ASSUMED
033371,000612: 07,3430 53342 VSR2 VAD # CONSTANT HERE.)
033372,000613: 07,3431 02031 RCV
033373,000614: 07,3432 02001 STORE RRECT
033374,000615: 07,3433 26031 STOVL RCV # SET UP CONIC ,ANSWER, FOR TIMESTEP
033375,000616: 07,3434 02023 TNUV
033376,000617: 07,3435 53202 VSR8 VAD # SAME FOR VELOCITY.
033377,000618: 07,3436 02037 VCV
033378,000619: 07,3437 02007 STORE VRECT
033379,000620: 07,3440 66370 AXT,1 SSP
033380,000621: 07,3441 00014 12D # ZERO DELTA, NU, AND TIME SINCE RECT.
033381,000622: 07,3442 00051 S1
033382,000623: 07,3443 00002 2
033383,000624:
033384,000625: 07,3444 16037 STODL VCV
033385,000626: 07,3445 17660 DPZERO
033386,000627: 07,3446 02045 STORE TC
033387,000628: 07,3447 02051 STORE XKEP # ZERO X.
033388,000629: 07,3450 06230 ZEROLOOP STORE YV +12D,1 # INDICES CAUSE LOOP TO ZERO 6 CONSECUTIVE
033389,000630: 07,3451 06031 STORE TDELTAV +12D,1 # DP NUMBERS (DELTA AND NU ARE ADJACENT).
033390,000631: 07,3452 43500 TIX,1 ITCQ # LOOP OR START INTEGRATION STEP IF DONE.
033391,000632: 07,3453 17450 ZEROLOOP
033392,000633:
|
|
Page 740 |
033394,000635: # THE THREE DIFEQ ROUTINES - DIFEQ+0, DIFEQ+12, AND DIFEQ+24 - ARE ENTEREDTO PROCESS THE CONTRIBUTIONS AT THE
033395,000636: # BEGINNING, MIDDLE, AND END OF THE TIMESTEP, RESPECTIVELY. THE UPDATING IS DONE BY THE NYSTROM METHOD.
033396,000637:
033397,000638: 07,3454 64575 DIFEQ+0 VLOAD VSR3
033398,000639: 07,3455 02103 FV
033399,000640: 07,3456 36067 STCALL PHIV
033400,000641: 07,3457 17563 DIFEQCOM
033401,000642: 07,3460 74575 DIFEQ+1 VLOAD VSR1
033402,000643: 07,3461 02103 FV
033403,000644: 07,3462 53206 PUSH VAD
033404,000645: 07,3463 02067 PHIV
033405,000646: 07,3464 26075 STOVL PSIV
033406,000647: 07,3465 53362 VSR1 VAD
033407,000648:
033408,000649: 07,3466 02067 PHIV
033409,000650: 07,3467 36067 STCALL PHIV
033410,000651: 07,3470 17563 DIFEQCOM
033411,000652: 07,3471 57345 DIFEQ+2 DLOAD DMPR
033412,000653: 07,3472 02165 H
033413,000654: 07,3473 17662 DP2/3
033414,000655: 07,3474 74206 PUSH VXSC
033415,000656: 07,3475 02067 PHIV
033416,000657: 07,3476 53372 VSL1 VAD
033417,000658: 07,3477 02222 ZV
033418,000659: 07,3500 53361 VXSC VAD
033419,000660: 07,3501 02165 H
033420,000661:
033421,000662: 07,3502 02214 YV
033422,000663: 07,3503 26214 STOVL YV
033423,000664: 07,3504 02103 FV
033424,000665: 07,3505 53322 VSR3 VAD
033425,000666: 07,3506 02075 PSIV
033426,000667: 07,3507 76561 VXSC VSL1
033427,000668: 07,3510 43055 VAD BOFF # SEE IF THIS IS STATE VECTOR OR W COLUMN.
033428,000669: 07,3511 02222 ZV
033429,000670: 07,3512 00341 JSWITCH
033430,000671: 07,3513 17533 ENDSTATE
033431,000672: 07,3514 77770 AXT,1
033432,000673: 07,3515 00000 0
033433,000674: 07,3516 12341 STORE W +72D,2 # VELOCITY COLUMN VECTOR
033434,000675: 07,3517 77745 DLOAD
033435,000676: 07,3520 02214 YV
033436,000677: 07,3521 12275 STORE W +36D,2 # POSITION COLUMN VECTOR
033437,000678: 07,3522 52104 TIX,2 GOTO # **********
033438,000679: 07,3523 17525 NEXTCOL
033439,000680: 07,3524 02204 STEPEXIT
033440,000681:
|
|
Page 741 |
033442,000683: 07,3525 77773 NEXTCOL VLOAD* # SET UP NEXT COLUMNS OF W MATRIX
033443,000684: 07,3526 75502 W +36D,2
033444,000685: 07,3527 32214 STOVL* YV
033445,000686: 07,3530 75436 W +72D,2
033446,000687: 07,3531 36222 STCALL ZV
033447,000688: 07,3532 17417 DIFEQ0
033448,000689: 07,3533 26023 ENDSTATE STOVL TNUV
033449,000690: 07,3534 02214 YV
033450,000691:
033451,000692: 07,3535 16015 STODL TDELTAV
033452,000693: 07,3536 02165 H
033453,000694: 07,3537 64402 SR4 SR2R
033454,000695: 07,3540 77615 DAD
033455,000696: 07,3541 02045 TC
033456,000697: 07,3542 02045 STORE TC
033457,000698: 07,3543 77014 BOFF AXT,2
033458,000699: 07,3544 00340 WMATFLAG
033459,000700: 07,3545 02204 STEPEXIT
033460,000701: 07,3546 00044 36D
033461,000702:
033462,000703: 07,3547 76131 SSP AXT,1
033463,000704: 07,3550 00052 S2
033464,000705: 07,3551 00006 6
033465,000706: 07,3552 17603 DOW..
033466,000707: 07,3553 66130 SXA,1 SXA,1
033467,000708: 07,3554 02170 FBRANCH
033468,000709: 07,3555 02171 HBRANCH
033469,000710: 07,3556 76014 INVERT AXT,1
033470,000711: 07,3557 00161 JSWITCH
033471,000712: 07,3560 00000 0
033472,000713: 07,3561 77650 GOTO
033473,000714: 07,3562 17525 NEXTCOL
033474,000715:
|
|
Page 742 |
033476,000717: # COMES HERE TO FINISH FIRST TWO DIFEQ COMPUTATIONS.
033477,000718:
033478,000719:
033479,000720:
033480,000721: 07,3563 43345 DIFEQCOM DLOAD DAD # INCREMENT H AND DIFEQCNT.
033481,000722: 07,3564 02163 DT/2
033482,000723: 07,3565 02165 H
033483,000724: 07,3566 66110 INCR,1 SXA,1
033484,000725: 07,3567 77763 DEC -12 B-14
033485,000726: 07,3570 02205 DIFEQCNT # DIFEQCNT SET FOR NEXT ENTRY.
033486,000727: 07,3571 02165 STORE H
033487,000728: 07,3572 74561 VXSC VSR1
033488,000729: 07,3573 02103 FV
033489,000730: 07,3574 74255 VAD VXSC
033490,000731:
033491,000732: 07,3575 02222 ZV
033492,000733: 07,3576 02165 H
033493,000734: 07,3577 77655 VAD
033494,000735: 07,3600 02214 YV
033495,000736: 07,3601 36053 STCALL ALPHAV
033496,000737: 07,3602 02170 FBRANCH
033497,000738:
|
|
Page 743 |
033499,000740: # ORBITAL ROUTINE FOR EXTRAPOLATION OF THE W MATRIX. IT COMPUTES THE SECOND DERIVATIVE OF EACH COLUMN POSITION
033500,000741: # VECTOR OF THE MATRIX AND CALLS THE NYSTROM INTEGRATION ROUTINES TO SOLVETHE DIFFERENTIAL EQUATIONS. THE PROGRAM
033501,000742: # USES A TABLE OF VEHICLE POSITION VECTORS COMPUTED DURING THE INTEGRATION OF THE VEHICLES POSITION AND VELOCITY.
033502,000743: 07,3603 50575 DOW.. VLOAD VSR6
033503,000744: 07,3604 02053 ALPHAV
033504,000745: 07,3605 53513 PDVL* UNIT
033505,000746:
033506,000747: 07,3606 02111 VECTAB,1
033507,000748: 07,3607 46315 PDVL VPROJ
033508,000749: 07,3610 02053 ALPHAV
033509,000750: 07,3611 52361 VXSC VSU
033510,000751: 07,3612 17626 3/16
033511,000752: 07,3613 41325 PDDL DMP
033512,000753: 07,3614 00043 34D
033513,000754: 07,3615 00045 36D
033514,000755: 07,3616 45541 V/SC STADR
033515,000756: 07,3617 41674 STCALL FV
033516,000757: 07,3620 17350 NBRANCH
033517,000758:
|
|
Page 744 |
033519,000760: # CONSTANTS
033520,000761:
033521,000762:
033522,000763:
033523,000764: 07,3621 00000 KEPSILON OCT 00000
033524,000765: 07,3622 00002 OCT 00002
033525,000766: 07,3623 14000 00000 THREE/8 2DEC .375
033526,000767:
033527,000768: 07,3625 06000 00000 3/16 2DEC 3 B -4
033528,000769:
033529,000770: 07,3627 10000 00000 DP1/4 2DEC .25
033530,000771:
033531,000772: 07,3631 12525 12525 DP1/3 2DEC .333333333
033532,000773:
033533,000774: 07,3633 07,3627 DQUARTER EQUALS DP1/4
033534,000775: 07,3633 02000 00000 POS1/16 2DEC .0625
033535,000776:
033536,000777: 07,3635 07,3627 POS1/4 EQUALS DP1/4
033537,000778: 07,3635 07,3623 3/8 EQUALS THREE/8
033538,000779: 07,3635 36000 00000 15/16 2DEC 15. B -4
033539,000780:
033540,000781: 07,3637 30000 00000 3/4 2DEC 3.0 B -2
033541,000782:
033542,000783: 07,3641 12577 24106 J2REQSQ 2DEC .335914874 # SECOND HARMONIC TIMES SQUARE OF RADIUS.
033543,000784:
033544,000785: 07,3643 22525 12525 7/12 2DEC .5833333333
033545,000786:
033546,000787:
033547,000788: 07,3645 22000 00000 9/16 2DEC 9 B -4
033548,000789:
033549,000790: 07,3647 01200 00000 5/128 2DEC 5 B -7
033550,000791:
033551,000792: 07,3651 77711 71033 2J3RE/J2 2DEC -.003309146
033552,000793:
033553,000794: 07,3653 23377 06703 J4REQ/J3 2DEC .60932709
033554,000795:
033555,000796: 07,3655 20000 00000 DP1/2 2DEC .5
033556,000797:
033557,000798: 07,3657 00000 00000 DPZERO 2DEC 0.0
033558,000799:
033559,000800: 07,3661 25252 25253 DP2/3 2DEC .6666666667
033560,000801:
033561,000802: 07,3663 07,3661 2/3 EQUALS DP2/3
033562,000803:
|
|
Page 745 |
033564,000805: # DUMMYMOON POSITION ROUTINE, SUN POSITION ROUTINE, AND PBODY TABLE FOR CHECKOUT OF EARTH-ORBITAL ONLY.
033565,000806: 07,3663 43575 MOONPOS VLOAD ITCQ # LOAD CONSTANT VECTOR INTO A AND EXIT.
033566,000807: 07,3664 17672 MOONVEC
033567,000808:
033568,000809: 07,3665 77201 SUNPOS SETPD VLOAD # RETURNS WITH VECTOR IN VAC AND IN PDL.
033569,000810: 07,3666 00001 0
033570,000811: 07,3667 17675 SUNVEC
033571,000812: 07,3670 77616 ITCQ
033572,000813:
033573,000814:
033574,000815:
033575,000816: 07,3671 00006 EARTHTAB DEC 6 B-14
033576,000817: 07,3672 00000 00000 2DEC 0.0
033577,000818:
033578,000819: 07,3674 00012 DEC 10 B-14
033579,000820: 07,3675 00000 00000 2DEC 0.0
033580,000821:
033581,000822: 07,3677 00000 00000 2DEC 0.0
033582,000823:
033583,000824:
033584,000825: 07,3701 77743 DEC -28 B-14 # 28 3/2
033585,000826: 07,3702 24214 11244 2DEC .6335627 # 400/SQRT(MU)
033586,000827:
033587,000828: 07,3704 07,3671 MOONVEC EQUALS EARTHTAB
033588,000829: 07,3704 07,3674 SUNVEC EQUALS EARTHTAB +3 # ******
End of include-file ORBITAL_INTEGRATION_PROGRAM.agc. Parent file is MAIN.agc