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