# Copyright: Public domain. # Filename: CONIC_SUBROUTINES.agc # Purpose: Part of the source code for Colossus 2A, AKA Comanche 055. # It is part of the source code for the Command Module's (CM) # Apollo Guidance Computer (AGC), for Apollo 11. # Assembler: yaYUL # Contact: Ron Burkey . # Website: www.ibiblio.org/apollo. # Pages: 1262-1308 # Mod history: 2009-05-08 RSB Adapted from the Colossus249/ file of the # same name, using Comanche055 page images. # 2009-05-20 RSB Corrected: Fixed four interpreter # instructions. # # This source code has been transcribed or otherwise adapted from digitized # images of a hardcopy from the MIT Museum. The digitization was performed # by Paul Fjeld, and arranged for by Deborah Douglas of the Museum. Many # thanks to both. The images (with suitable reduction in storage size and # consequent reduction in image quality as well) are available online at # www.ibiblio.org/apollo. If for some reason you find that the images are # illegible, contact me at info@sandroid.org about getting access to the # (much) higher-quality images which Paul actually created. # # Notations on the hardcopy document read, in part: # # Assemble revision 055 of AGC program Comanche by NASA # 2021113-051. 10:28 APR. 1, 1969 # # This AGC program shall also be referred to as # Colossus 2A # Page 1262 # PROGRAM DESCRIPTION -- ENTIRE CONIC SUBROUTINE LOG SECTION DATE -- 1 SEPTEMBER 1967 # MOD NO. -- 0 LOG SECTION -- CONIC SUBROUTINES # MOD BY KRAUSE ASSEMBLY -- COLOSSUS REVISION 88 # # FUNCTIONAL DESCRIPTION -- # THE FOLLOWING SET OF SUBROUTINES SOLVE VARIOUS PROBLEMS INVOLVING THE TRAJECTORY PRODUCED BY A CENTRAL # INVERSE-SQUARE FORCE ACTING ON A POINT MASS, AS OUTLINED IN THE CMC AND LGC LUNAR LANDING MISSION GSOP, SECTION # 5.5.1.2. A GENERAL USAGE POINT-OF-VIEW WAS TAKEN IN FORMULATING, MECHANIZING, AND SCALING THE SUBROUTINES, # RATHER THAN OPTIMIZING EACH FOR A PARTICULAR USE. THEREFORE, MULTIPLE USAGE CAN BE MADE OF THE SUBROUTINES # INVOLVING ANY REALISTIC SET OF CONSTRAINTS. IT SHOULD BE NOTED THAT ONLY ONE SET OF CODING IS USED, WHETHER THE # EARTH, MOON, OR ANY OTHER CELESTIAL BODY IS SPECIFIED AS THE CENTRAL BODY OF THE PROBLEM, PROVIDED ONE OBSERVES # THE INHERENT SCALE CHANGE REQUIRED IN POSITION, VELOCITY, MU, AND TIME, AS OUTLINES IN MISSION PROGRAMMING # DEFINITION MEMO NO. 10. THIS CAN BE ACCOMPLISHED BY SIMPLY ADDING TO THE MUTABLE AND INITIALIZING THE SUBROUTINES # APPROPRIATELY. # # DUE TO THE UNIFORMITY OF THE EQUATIONS INVOLVED, CODING WAS MINIMIZED BY TREATING INDIVIDUAL EQUATIONS AND # BLOCKS OF EQUATIONS AS SUBROUTINES OF LOWER RANK WHENEVER POSSIBLE. AS A RESULT, THREE BY-PRODUCTS SUBROUTINES, # DIRECTLY USABLE AS INDEPENDENT SUBROUTINES, WERE GENERATED. # # RESTRICTIONS -- # THE ONLY LIMITATION IN THE SCOPE OF THE PROBLEM WHICH CAN BE SOLVED BY A PARTICULAR SUBROUTINE IS THE SCALING # LIMIT OF EACH PARAMETER AS SPECIFIED IN THE GSOP. THESE SCALING LIMITS WERE CHOSEN SO THAT ALL FEASIBLE TRAJECTORIES # COULD BE HANDLED. # # SINCE THE SUBROUTINES (EXCEPT KEPLER) USE COMMON SUBROUTINES OF LOWER RANK WHICH USE ERASABLE OTHER THAN # THE PUSHLIST (DUE TO ITS LIMITED SIZE) AND COMMON INTERPRETIVE SWITCHES, THE CONIC SUBROUTINES CANNOT BE ALLOWED # TO INTERRUPT EACH OTHER. IT IS UP TO THE USER TO GUARANTEE THIS CONDITION. # Page 1263 # PROGRAM DESCRIPTION -- KEPLER SUBROUTINE DATE -- 11 OCTOBER 1967 # MOD NO. -- 1 LOG SECTION -- CONIC SUBROUTINES # MOD BY KRAUSE ASSEMBLY -- COLOSSUS 103 AND SUNDANCE 222 # MOD NO. -- 2 (AUGUST 1968) BY ROBERTSON: TO PERMIT BACKDATING BY MORE THAN ONE ORBITAL PERIOD. # MOD NO. -- 3 (DEC 1968) BY ROBERTSON: SUPPRESSION OF X-MODULO-ING # MOD NO. -- 4 (JAN 1969) BY ROBERTSON: CLEAR OVFIND AT KEPLER ENTRY # # FUNCTIONAL DESCRIPTION -- # THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR AND THE DESIRED TRANSFER TIME THROUGH WHICH THE STATE IS TO # BE UPDATED ALONG A CONIC TRAJECTORY, COMPUTES THE NEW, UPDATED STATE VECTOR. THE TRAJECTORY MAY BE ANY CONIC # SECTION -- CIRCULAR, ELLIPTIC, PARABOLIC, HYPERBOLIC, OR RECTILINEAR WITH RESPECT TO THE EARTH OR THE MOON. THE # USE OF THE SUBROUTINE CAN BE EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT # INTRODUCING ANY CODING CHANGES, ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY. AN ITERATION # TECHNIQUE IS UTILIZED IN THE COMPUTATION. # # IF A NEGATIVE TIME-OF-FLIGHT IS INPUT, THE PROGRAM WILL SOLVE FOR THE STATE WHICH WOULD BE PRODUCED BY # EXTRAPOLATING THE POSITION BACKWARD IN TIME. # # IF THE ABSOLUTE VALUE DESIRED TRANSFER TIME EXCEEDS THE ORBITAL PERIOD, THE SUBROUTINE, THROUGH A # MODULAR TECHNIQUE, WILL COMPUTE THE STATE CORRESPONDING TO THE DESIRED TIME AS USUAL. # # THE RESTRICTIONS ARE -- # 1. (PREVIOUS RESTRICTION ON THE NEGATIVE DESIRED TRANSFER TIME IS NOW DELETED.) # 2. THE PARAMETERS IN THE PROBLEM CANNOT EXCEED THEIR SCALING LIMITS AS SPECIFIED IN THE GSOP. IF # ANY OF THESE LIMITS ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS. # # THE NUMBER OF ITERATIONS AND, THEREFORE, THE COMPUTATION SPEED IS DEPENDENT ON THE ACCURACY OF THE # GUESS, XKFPNEW. THE AGC COMPUTATION TIME IS APPROXIMATELY .061 SECONDS FOR INITIALIZATION, .065 SECONDS FOR THE # FINAL COMPUTATIONS, PLUS .083 SECONDS FOR EACH ITERATION. # # REFERENCES -- # R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP, SECTION 5.5, SGA # MEMO 67-4. # # INPUT -- ERASABLE INITIALIZATION REQUIRED # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # RRECT +29 FOR EARTH DP INITIAL POSITION VECTOR IN METERS # +27 FOR MOON # Page 1264 # VRECT +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND # +5 FOR MOON # X1 (38D) NONE INDEX REGISTER SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON, # RESPECTIVELY, IS THE CENTRAL BODY # TAU +28 DESIRED TRANSFER TIME IN CENTISECONDS (DP) # MAY BE POS OR NEG AND ABSOLUTE VALUE MAY BE GREATER OR LESS THAN ONE ORBITAL PERIOD. # XKEPNEW +17 FOR EARTH DP GUESS OF ROOT X OF KEPLERS EQN IN SQRT(METERS).SIGN SHOULD AGREE WITH THAT OF TAU. # +16 FOR MOON AND ABS VALUE SHOULD BE LESS THAN THAT CORRESPONDING TO A PERIOD, VIZ, 2PI SQRT(SEMI- # MAJOR AXIS), FOR SPEED OF CONVERGENCE, BUT IF EITHER CONDITION FAILS, XKEPNEW IS RESET # BY KEPLER TO A POOR BUT VALID GUESS. # TC +28 DP PREV. VALUE OF TIME IN CENTISECS. MUST BE LESS THAN ONE ORBITAL PERIOD. # XPREV +17 FOR EARTH PREVIOUS VALUE OF X IN SQRT(METERS). MUST BE LESS THAN AN X CORRESPONDING TO ONE # +16 FOR MOON ORBITAL PERIOD, VIZ, 2PI SQRT(SEMI-MAJOR AXIS) # # SUBROUTINES CALLED -- # DELTIME # # CALLING SEQUENCE AND NORMAL EXIT MODES -- # KEPRTN-2 GOTO # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # KEPRTN-1 KEPLER # RETURNS WITH XPREV IN MPAC. PL IS AT 0. # KEPRTN ... # CONTINUE # # KEPLER MUST NOT BE CALLED DIRECTLY SINCE AN INTERRUPTION OF IT WOULD DESTROY THE ERASABLES IT NEEDS TO COMPLETE # THE INTERRUPTED JOB. THEREFORE THE USER MUST CALL CSMCONIC OR LEMCONIC WHICH GUARANTEES NO INTERRUPTS AND WHICH # ALSO CALLS KEPPREP TO COMPUTE A GUESS OF XKEPNEW. # # ABORT EXIT MODE -- # NONE # # OUTPUT -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # RCV +29 FOR EARTH DP TERMINAL POSITION VECTOR IN METERS # +27 FOR MOON # VCV +7 FOR EARTH DP TERMINAL VELOCITY VECTOR IN METERS/CENTISEC # +5 FOR MOON # TC +28 DP TRANSFER TIME IN CENTISECS TO WHICH KEPLER CONVERGED. ALWAYS LESS THAN ONE PERIOD. # XPREV +17 FOR EARTH DP VALUE OF X IN SQRT(METERS) TO WHICH KEPLER CONVERGED. ALWAYS LESS THAN THE X # +16 FOR MOON CORRESPONDING TO ONE PERIOD. # Page 1265 # FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS. # # DEBRIS -- # PARAMETERS WHICH MAY BE OF USE -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # URRECT +1 DP UNIT VECTOR OF INITIAL POSITION # R1 +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION IN METERS # +27 FOR MOON # ALPHA -22 FOR EARTH DP INVERSE OF SEMI-MAJOR AXIS IN 1/METERS # -20 FOR MOON # TMODULO +28 DP INTEGRAL NUMBER OF PERIODS IN CENTISECS, WHICH WAS SUBTRACTED FROM TAU. TO PRODUCE A # TAU. OF LESS THAN ONE PERIOD. # # PARAMETERS OF NO USE -- # DP PARAMETERS -- FPSILENT, DELX, DELT, RCNORM, XMODULO, PLUS PUSHLIST REGISTERS 0 THROUGH 39D. # Page 1266 # PROGRAM DESCRIPTION -- LAMBERT SUBROUTINE DATE -- 1 SEPTEMBER 1967 # MOD NO. -- 0 LOG SECTION -- CONIC SUBROUTINES # MOD BY KRAUSE ASSEMBLY -- COLOSSUS REVISION 88 # # FUNCTIONAL DESCRIPTION -- # THIS SUBROUTINE CALCULATES THE INITIAL VELOCITY REQUIRED TO TRANSFER A POINT-MASS ALONG A CONIC TRAJECTORY # FROM AN INITIAL POSITION TO A TERMINAL POSITION IN A PRESCRIBED TIME INTERVAL. THE RESULTING TRAJECTORY MAY BE # A SECTION OF A CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON. THE USE OF THE # SUBROUTINE CAN BE EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY # CODING CHANGES, ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY. AN ITERATION TECHNIQUE IS # UTILIZED IN THE COMPUTATION. # # THE RESTRICTIONS ARE: -- # 1. RECTILINEAR TRAJECTORIES CANNOT BE COMPUTED. # 2. AN ACCURACY DEGRADATION OCCURS AS THE COSINE OF THE TRUE ANOMALY DIFFERENCE APPROACHES +1.0. # 3. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES # AND LESS THAN 178 DEGREES 12.5 MINUTES. # 4. NEGATIVE TRANSFER TIME IS AMBIGUOUS AND WILL RESULT IN NO SOLUTION. # 5. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP. IF THE # LIMITS ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS. # # THE NUMBER OF ITERATIONS AND, THEREFORE, THE COMPUTATION'S SPEED IS DEPENDENT ON THE ACCURACY OF THE FIRST # GUESS OF THE INDEPENDENT VARIABLE, COGA. THE AGC COMPUTATION TIME IS APPROXIMATELY # .105 SECONDS FOR INITIALIZATION, .069 SECONDS FOR FINAL COMPUTATIONS, PLUS .205 SECONDS FOR EACH ITERATION. # # REFERENCES -- # R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP -- SECTION 5.5, SGA MEMO 67-8, # SGA MEMO 67-4. # # INPUT -- ERASABLE INITIALIZATION REQUIRED # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # R1VEC +29 FOR EARTH DP INITIAL POSITION VECTOR IN METERS # +27 FOR MOON # R2VEC +29 FOR EARTH DP TARGET OR TERMINAL POSITION VECTOR IN METERS # +27 FOR MOON # TDESIRED +28 DP DESIRED TRANSFER TIME IN CENTISECONDS # X1 (38D) NONE INDEX REGISTER SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON, # RESPECTIVELY, IS THE CENTRAL BODY # GEOMSGN NONE SP +.5 IF DESIRED TRANSFER ANGLE IS LESS THAN 180 DEGREES, -.5 IF GREATER THAN 180 DEG. # GUESSW NONE AN INTERPRETER SWITCH TO BE SET IF NO GUESS OF COGA IS AVAILABLE, CLEAR IF A GUESS OF # Page 1267 # COGA IS TO BE USED BY LAMBERT # COGA +5 DP GUESS OF COTANGENT OF FLIGHT PATH ANGLE (MEASURED FROM VERTICAL). THIS WILL BE # IGNORED IF GUESSW IS SET. # NORMSW NONE AN INTERPRETER SWITCH TO BE SET IF UN IS TO BE AN INPUT TO THE SUBROUTINE, CLEAR IF # LAMBERT IS TO COMPUTE ITS OWN NORMAL (UN). # UN +1 DP UNIT NORMAL TO THE DESIRED ORBIT PLANE IN THE DIRECTION OF THE RESULTING ANGULAR # MOMENTUM VECTOR. THIS WILL BE IGNORED IF NORMSW IS CLEAR. # VTARGTAG NONE A S.P. TAG TO BE SET TO ZERO IF LAMBERT IS TO COMPUTE THE VELOCITY OF R2VEC AS WELL AS # AT R1VEC. # # SUBROUTINES CALLED -- # GEOM, GETX, DELTIME, ITERATOR, LAMENTER (PART OF NEWSTATE) # # CALLING SEQUENCE AND NORMAL EXIT MODES -- # L CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # L+1 LAMBERT # RETURNS WITH PL AT 0 AND WITH VVEC IN MPAC IF VTARGTAG WAS WAS NON-ZERO OR VTARGET # # IN MPAC IF VTARGTAG WAS ZERO # L+2 BON # CONTINUE IF SOLNSW CLEAR SINCE SOLUTION IS ACCEPTABLE # L+3 SOLNSW # L+4 LAMABORT # # IF A LAMBER RESULT IS TO BE A FIRST GUESS FOR THE NEXT LAMBERT CALCULATION, COGA MUST BE PRESERVED AND # GUESSW MUST BE CLEAR FOR EACH SUCCEEDING LAMBERT CALL. # # ABORT EXIT MODES -- # IF SOLNSW WAS SET UPON EXITING, EITHER LAMBERT WAS ASKED TO COMPUTE A TRANSFER TOO NEAR 0 OR 360 DEG, OR T # WAS TOO SMALL TO PRODUCE A REALISTIC TRANSFER BETWEEN R1VEC AND R2FEC. IN EITHER CASE THE FIX MUST BE MADE # ACCORDING TO THE NEEDS OF THE PARTICULAR USER. THE ABORT EXIT MODE MAY BE CODED AS ... # LAMBERT DLOAD ABS # A MEASURE OF THE PROXIMITY TO 0 OR # 1-CHTH # 360 DEGREES. # DSU BWM # ONEBIT # CHANGER2 # CHANGE R2VEC DIRECTION SLIGHTLY. # DLOAD DAD # TDESIRED # SOMETIME # STCALL TDESIRED # INCRESE TDESIRED # LAMBERT # # OUTPUT -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # Page 1268 # -------- -------------- ----------------------- # VVEC +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND REQUIRED TO SATISFY THE BOUNDARY VALUE # +5 FOR MOON PROBLEM. # VTARGET +7 FOR EARTH DP RESULTANT VELOCITY VECTOR AT R2VEC IN METERS/CENTISECOND. # +5 FOR MOON # SOLNSW NONE INTERPRETER SWITCH WHICH IS SET IF THE SUBROUTINE CANNOT SOLVE THE PROBLEM, CLEAR IF THE # SOLUTION EXISTS. # # FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS. # # DEBRIS -- # PARAMETERS WHICH MAY BE OF USE -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # SNTH +1 DP SIN OF ANGLE BETWEEN R1VEC AND R2VEC # CSTH +1 DP COSINE OF ANGLE # 1-CSTH +2 DP 1-CSTH # COGA +5 DP COTAN OF INITIAL REQUIRED FLIGHT PATH ANGLE MEASURED FROM VERTICAL # P +4 DP RATIO OF SEMILATUS RECTUM TO INITIAL RADIUS # R1A +6 DP RATIO OF INITIAL RADIUS TO SEMI-MAJOR AXIS # R1 (32D) +29 FOR EARTH DP INITIAL RADIUS IN METERS # +27 FOR MOON # UR1 +1 DP UNIT VECTOR OF R1VEC # U2 +1 DP UNIT VECTOR OF R2VEC # # PARAMETERS OF NO USE -- # DP PARAMETERS -- EPSILONL, CSTH-RHO, TPREV, TERRLAMB, R2, RTNLAMB (SP), PLUS PUSHLIST REGISTER 0 THROUGH 41D # ADDITIONAL INTERPRETIVE SWITCHES USED -- INFINFLG, 360SW, SLOPESW, ORDERSW # Page 1269 # PROGRAM DESCRIPTION -- TIME-THETA SUBROUTINE DATE -- 1 SEPTEMBER 1967 # MOD NO. -- 0 LOG SECTION -- CONIC SUBROUTINES # MOD BY KRAUSE ASSEMBLY -- COLOSSUS REVISION 88 # # FUNCTIONAL DESCRIPTION -- # THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR AND A DESIRED TRUE-ANOMALY-DIFFERENCE THROUGH WHICH THE # STATE IS TO BE UPDATED ALONG A CONIC TRAJECTORY, CALCULATES THE CORRESPONDING TIME-OF-FLIGHT AND, IN ADDITION, # PROVIDES THE OPTION OF COMPUTING THE NEW UPDATED STATE VECTOR. THE RESULTING TRAJECTORY MAY BE A SECTION OF A # CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON. THE USE OF THE SUBROUTINE CAN BE # EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY CODING CHANGES, # ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY. # # THE RESTRICTIONS ARE -- # 1. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES # AND LESS THAN 178 DEGREES 12.5 MINUTES. # 2. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSCP. IF THE LIMITS # ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS. # # THE AGC COMPUTATION TIME IS APPROXIMATELY .292 SECONDS. # # REFERENCES -- # R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP-SECTION 5.5, SGA MEMO 67-8. # # INPUT -- ERASABLE INITIALIZATION REQUIRED # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # RVEC +29 FOR EARTH DP INITIAL POSITION VECTOR IN METERS # +27 FOR MOON # VVEC +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND # +5 FOR MOON # SNTH +1 ` DP SINE OF THE TRUE-ANOMALY-DIFFERENCE THROUGH WHICH THE STATE IS TO BE UPDATED # CSTH +1 DP COSINE OF THE ANGLE # RVSW NONE AN INTERPRETIVE SWITCH TO BE SET IF ONLY TIME IS TO BE AN OUTPUT, CLEAR IF THE NEW STATE # IS TO BE COMPUTED ALSO. # X1 (38D) NONE INDEX REGISTER TO BE SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON, # RESPECTIVELY, IS THE CENTRAL BODY. # # SUBROUTINES CALLED -- # Page 1270 # PARAM, GEOM, GETX, DELTIME, NEWSTATE # # CALLING SEQUENCE AND NORMAL EXIT MODES -- # IF ONLY TIME IS DESIRED AS OUTPUT -- # L SET CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # L+1 RVSW # L+2 TIMETHET # RETURN WITH PL AT 0 AND T IN MPAC # L+3 ... # CONTINUE # # IF THE UPDATE STATE VECTOR IS DESIRED AS WELL -- # L CLEAR CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # L+1 RVSW # L+2 TIMETHET # RETURNS WITH PL AT 6. THE INITIAL POSITION VECTOR IS IN 0D OF THE PUSHLIST AND # # THE INITIAL VELOCITY VECTOR IN MPAC. # L+3 STOVL NEWVVEC # L+4 STADR # L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USER'S LOCATIONS. # L+6 ... # CONTINUE. # # ABORT EXIT MODES -- # IF COGAFLAG AND/OR INFINFLG IS SET AT THE EXIT TO TIME-THETA, TIME-THETA WILL TRANSFER TO P00DOO WITH # AN ALARM CODE (ORIGINALLY 00607), AND NOT RETURN TO THE CALLING PROGRAM. (PCR 692 AND 721). # # OUTPUT -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # T(30D) +28 DP TRANSFER TIME IN CENTISECONDS # INFINFLG NONE AN INTERPRETIVE SWITCH WHICH IS SET IF THE TRANSFER ANGLE REQUIRES CLOSURE THROUGH # INFINITY (NO SOLUTION), CLEAR IF A PHYSICAL SOLUTION IS POSSIBLE. # COGAFLAG NONE AN INTERPRETIVE SWITCH WHICH IS SET IF RESTRICTION 1 HAS BEEN VIOLATED (NO SOLUTION), # CLEAR IF NOT. # # IN ADDITION, IF RVSW IS CLEAR, THE FOLLOWING ARE OUTPUT -- # MPAC - +7 FOR EARTH DP TERMINAL VELOCITY VECTOR IN METERS/CENTISEC. # MPAC+5 +5 FOR MOON # 0D - 5D +29 FOR EARTH DP TERMINAL POSITION VECTOR IN METERS (PL AT 6D) # +27 FOR MOON # # FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS. # Page 1271 # # DEBRIS -- # PARAMETERS WHICH MAY BE OF USE -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS # +27 FOR MOON # R1A +6 DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES) # P +4 DP RATIO OF SEMILATUS RECTUM TO R1 # COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC # UR1 +1 DP UNIT VECTOR OF RVEC # U2 +1 DP UNIT VECTOR OF VVEC # UN +1 DP UNIT VECTOR OF UR1*U2 # # PARAMETERS OF NO USE -- # SP PARAMETERS -- RTNTT, GEOMSGN, RTNPRM, MAGVEC2=R2 (DP), PLUS PUSHLIST LOCATIONS 0-11D, 14D-21D, 24D-39D, 41D # ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW, 360SW # Page 1272 # PROGRAM DESCRIPTION -- TIME-RADIUS SUBROUTINE DATE -- 11 OCTOBER 1967 # MOD NO. -1 LOG SECTION -- CONIC SUBROUTINES # MOD BY KRAUSE ASSEMBLY -- COLOSSUS REVISION 88 # # FUNCTIONAL DESCRIPTION -- # THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR AND A DESIRED RADIUS TO WHICH THE # STATE IS TO BE UPDATED ALONG A CONIC TRAJECTORY, CALCULATES THE CORRESPONDING TIME-OF-FLIGHT AND, IN ADDITION, # PROVIDES THE OPTION OF COMPUTING THE NEW UPDATED STATE VECTOR. THE RESULTING TRAJECTORY MAY BE A SECTION OF A # CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON. THE USE OF THE SUBROUTINE CAN BE # EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY CODING CHANGES, # ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY. # # IF THE DESIRED RADIUS IS BEYOND THE RADIUS OF APOCENTER OF THE CONIC OR BELOW THE RADIUS OF PERICENTER, # APSESW WILL BE SET AND THE SUBROUTINE WILL RETURN THE APOCENTER OR PERICENTER SOLUTION, RESPECTIVELY. # # THE RESTRICTIONS ARE -- # 1. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES # AND LESS THAN 178 DEGREES 12.5 MINUTES. # 2. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP. IF THE LIMITS # EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS. # 3. AN ACCURACY DEGRADATION OCCURS AS THE SENSITIVITIES OF TIME AND UPDATED STATE VECTOR TO CHANGES IN # RDESIRED INCREASE. THIS WILL OCCUR NEAR EITHER APSIS OF THE CONIC AND WHEN THE CONIC IS NEARLY CIRCULAR. IN # PARTICULAR, IF THE CONIC IS AN EXACT CIRCLE, THE PROBLEM IS UNDEFINED AND THE SUBROUTINE WILL ABORT. # # THE AGC COMPUTATION TIME IS APPROXIMATELY .363 SECONDS. # # REFERENCES -- # R-479, MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP-SECTION 5.5, SGA MEMO 67-8. # # INPUT -- ERASABLE INITIALIZATION REQUIRED. # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # RVEC +29 FOR EARTH DP INITIAL POSITION VECTOR IN METERS # +27 FOR MOON # VVEC +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND # +5 FOR MOON # RDESIRED +29 FOR EARTH DP TERMINAL RADIAL DISTANCE ON CONIC TRAJECTORY FOR WHICH TRANSFER TIME IS TO BE # +27 FOR MOON COMPUTED # SGNRDOT NONE SP TAG SET TO +.5 OR -.5 ACCORDING TO WHETHER THE RADIAL VELOCITY AT RDESIRED IS TO BE # POSITIVE OR NEGATIVE, RESPECTIVELY. THIS TAG REDUCES THE DOUBLE-VALUED PROBLEM TO A # Page 1273 # SINGLE-VALUED PROBLEM. # X1 (38D) NONE INDEX REGISTER TO BE SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON, # RESPECTIVELY, IS THE CENTRAL BODY. # RVSW NONE AN INTERPRETIVE SWITCH TO BE SET IF ONLY TIME IS TO BE AN OUTPUT, CLEAR IF THE NEW STATE # IS TO BE COMPUTED ALSO. # # SUBROUTINES CALLED -- # PARAM, GEOM, GETX, DELTIME, NEWSTATE # # CALLING SEQUENCE AND NORMAL EXIT MODES -- # IF ONLY TIME IS DESIRED AS OUTPUT -- # L SET CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # L+1 RVSW # L+2 TIMERAD # RETURN WITH PL AT 0 AND T IN MPAC # L+3 ... # CONTINUE # # IF THE UPDATE STATE VECTOR IS DESIRED AS WELL -- # L CLEAR CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # L+1 RVSW # L+2 TIMERAD # RETURNS WITH PL AT 6. THE INITIAL POSITION VECTOR IS IN 0D OF THE PUSHLIST AND # # THE INITIAL VELOCITY VECTOR IN MPAC. # L+3 STOVL NEWVVEC # L+4 STADR # L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USER'S LOCATIONS. # L+6 ... # CONTINUE # # ABORT EXIT MODES -- # IF SOLNSW AND/OR COGAFLAG AND/OR INFINFLG IS SET AT THE EXIT TO TIME-RADIUS, TIME-RADIUS WILL TRANSFER # TO P00DOO WITH AN ALARM CODE (ORIGINALLY 00607), AND NOT RETURN TO THE CALLING PROGRAM. (PCR 692 & 721). # # OUTPUT -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # T (30D) +28 DP TRANSFER TIME IN CENTISECONDS. # INFINFLG NONE AN INTERPRETIVE SWITCH WHICH IS SET IF RDESIRED AND SGNRDOT REQUIRE CLOSURE THROUGH # INFINITY (NO SOLUTION), CLEAR IF A PHYSICAL SOLUTION IS POSSIBLE. # COGAFLAG NONE AN INTERPRETIVE SWITCH WHICH IS SET IF RESTRICTION 1 HAS BEEN VIOLATED (NO SOLUTION), # CLEAR IF NOT. # APSESW NONE AN INTERPRETIVE SWITCH WHICH IS SET IF RDESIRED WAS GREATER THAN RADIUS OF APOCENTER OR # Page 1274 # LESS THAN RADIUS OF PERICENTER. THE APOCENTER OR PERICENTER SOLUTION, RESPECTIVELY, # WILL THEN BE RETURNED. THE SWITCH IS CLEAR IF RDESIRED WAS BETWEEN PERICENTER AND # APOCENTER. # SOLNSW NONE AN INTERPRETIVE SWITCH WHICH IS SET IF THE CONIC IS SO CLOSE TO A CIRCLE THAT THE TERMIN # POINT IS AMBIGUOUS, VIOLATING RESTRICTION 3. IF ECCENTRICITY IS GREATER THAN 2-TO-THE- # MINUS-18, THE SWITCH IS CLEAR. # # IN ADDITION, IF RVSW IS CLEAR, THE FOLLOWING ARE OUTPUT -- # MPAC - +7 FOR EARTH DP TERMINAL VELOCITY VECTOR IN METERS/CENTISEC. # MPAC+5 +5 FOR MOON # 0D - 5D +29 FOR EARTH DP TERMINAL POSITION VECTOR IN METERS (PL AT 6D) # +27 FOR MOON # # FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS. # # DEBRIS -- # PARAMETERS WHICH MAY BE OF USE -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS # +27 FOR MOON # R1A +6 DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES) # P +4 DP RATIO OF SEMILATUS RECTUM TO R1 # COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC # UR1 +1 DP UNIT VECTOR OF RVEC # U2 +1 DP UNIT VECTOR OF VVEC # UN +1 DP UNIT VECTOR OF UR1*U2 # CSTH +1 DP COSINE OF TRUE ANOMALY DIFFERENCE BETWEEN RVEC AND RDESIRED. # SNTH +1 DP SINE OF TRUE ANOMALY DIFFERENCE. # # PARAMETERS OF NO USE -- # SP PARAMETERS -- RTNTT, GEOMSGN, RTNPRM, MAGVEC2*R2 (DP), PLUS PUSHLIST LOCATIONS 0-11D, 14D-21D, 24D-39D, 41D # ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW, 360SW # # Page 1275 # PROGRAM DESCRIPTION -- APSIDES SUBROUTINE DATE -- 1 SEPTEMBER 1967 # MOD NO. -- 0 LOG SECTION -- CONIC SUBROUTINES # MOD BY KRAUSE ASSEMBLY -- COLOSSUS REVISION 88 # # FUNCTIONAL DESCRIPTION -- # THIS SUBROUTINE, GIVEN AN INITIAL STATE VECTOR, CALCULATES THE RADIUS OF PERICENTER AND OF APOCENTER AND THE # ECCENTRICITY OF THE RESULTING CONIC TRAJECTORY, WHICH MAY BE A STRAIGHT LINE, # CIRCLE, ELLIPSE, PARABOLA, OR HYPERBOLA WITH RESPECT TO THE EARTH OR THE MOON. THE USE OF THE SUBROUTINE CAN # BE EXTENDED USING OTHER PRIMARY BODIES BY SIMPLE ADDITIONS TO THE MUTABLE WITHOUT INTRODUCING ANY CODING CHANGES, # ACCEPTING THE INHERENT SCALE FACTOR CHANGES IN POSITION AND VELOCITY. # # THE RESTRICTIONS ARE -- # 1. IF APOCENTER IS BEYOND THE SCALING OF POSITION, THE SCALE FACTOR LIMIT (536,870,910 METERS WITH RESPECT # TO THE EARTH OR 134,217,727.5 METERS WITH RESPECT TO THE MOON) WILL BE RETURNED. # 2. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP. IF THE LIMITS # ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS. # # THE AGC COMPUTATION TIME IS APPROXIMATELY .103 SECONDS. # # REFERENCES -- # MISSION PROGRAMMING DEFINITION MEMO NO. 10, LUNAR LANDING MISSION GSOP-SECTION 5.5. # # INPUT -- ERASABLE INITIALIZATION REQUIRED # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # RVEC +29 FOR EARTH DP INITIAL POSITION VECTOR IN METERS # +27 FOR MOON # VVEC +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND # +5 FOR MOON # X1 (38D) NONE INDEX REGISTER TO BE SET TO -2D OR -10D ACCORDING TO WHETHER THE EARTH OR MOON, # RESPECTIVELY, IS THE CENTRAL BODY. # # SUBROUTINES CALLED -- # PARAM, GEOM # # CALLING SEQUENCE AND NORMAL EXIT MODES -- # Page 1276 # IF ONLY TIME IS DESIRED AS OUTPUT -- # L CALL # MUST BE IN INTERPRETIVE MODE BUT OVFIND ARBITRARY. # L+1 APSIDES # RETURNS WITH PL AT 0, RADIUS OF APOCENTER IN MPAC AND RADIUS OF PERICENTER IN 0D # L+2 STODL APOAPSE # L+3 0D # L+4 STORE PERIAPSE # APOAPSE AND PERIAPSE ARE SYMBOLIC REPRESENTATIONS OF THE USER'S LOCATIONS # L+5 ... # CONTINUE # # OUTPUT -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # MPAC +29 FOR EARTH DP RADIUS OF APOCENTER IN METERS # +27 FOR MOON # 0D-1D +29 FOR EARTH DP RADIUS OF PERICENTER IN METERS # +27 FOR MOON # ECC +3 DP ECCENTRICITY OF CONIC TRAJECTORY # # FOR OTHER OUTPUT WHICH MAY BE OF USE, SEE DEBRIS. # # DEBRIS -- # PARAMETERS WHICH MAY BE OF USE -- # SCALE FACTOR # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # -------- -------------- ----------------------- # R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS # +27 FOR MOON # R1A +6 DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES) # P +4 DP RATIO OF SEMILATUS RECTUM TO R1 # COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC # UR1 +1 DP UNIT VECTOR OF RVEC # U2 +1 DP UNIT VECTOR OF VVEC # UN +1 DP UNIT VECTOR OF UR1*U2 # MAGVEC2 +7 FOR EARTH DP MAGNITUDE OF VVEC # +5 FOR MOON # # PARAMETERS OF NO USE -- # SP PARAMETERS -- RTNAPSE, GOMSGN, RTNPRM, PLUS PUSHLIST LOCATIONS 0-5, 10D-11D, 14D-21D, 31D-38D. # ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW SETLOC CONICS # Page 1277 BANK COUNT 12/CONIC EBANK= UR1 KEPLERN SETPD BOV 0 +1 VLOAD* MUTABLE,1 STOVL 14D RRECT UNIT SSP ITERCTR 20D STODL URRECT 36D STOVL R1 RRECT DOT SL1R VRECT DMP SL1R 1/ROOTMU # 1/ROOTMU (-17 OR -14) STOVL KEPC1 # C1=R.V/ROOTMU (+17 OR +16) VRECT VSQ DMPR 1/MU # 1/MU (-34 OR -28) DMP SL3 R1 DSU ROUND D1/64 STORE KEPC2 # C2=RV.V/MU -1 (+6) BDSU SR1R D1/64 DDV R1 STORE ALPHA # ALPHA=(1-C2)/R1 (-22 OR -20) BPL DLOAD # MAXIMUM X DEPENDS ON TYPE OF CONIC 1REV -50SC # -50SC (+12) DDV BOV ALPHA STOREMAX SQRT GOTO STOREMAX # Page 1278 1REV SQRT BDDV 2PISC # 2PISC (+6) BOV STOREMAX STOREMAX STORE XMAX DMP PDDL 1/ROOTMU ALPHA NORM PDDL X1 SL* DDV 0 -6,1 BOV BMN MODDONE MODDONE # MPAC=PERIOD PERIODCH PDDL ABS # 0D=PERIOD TAU. DSU BMN 0D MODDONE SIGN TAU. STODL TAU. GOTO PERIODCH MODDONE SETPD DLOAD 0 XKEPNEW STORE X SIGN BZE TAU. BADX BMN ABS BADX DSU BPL XMAX BADX STORBNDS DLOAD BPL TAU. STOREMIN DLOAD DCOMP XMAX STODL XMIN KEPZERO STCALL XMAX DXCOMP STOREMIN DLOAD KEPZERO STORE XMIN DXCOMP DLOAD DMPR # Page 1279 TAU. BEE22 ABS STODL EPSILONT XPREV XDIFF BDSU X STORE DELX KEPLOOP DLOAD DSQ X # X=XKEP NORM PUSH # 0D=XSQ (+34 OR +32 -N1) PL AT 2 X1 DMP SRR* ALPHA 0 -6,1 STCALL XI # XI=ALPHA XSQ (+6) DELTIME BOV BDSU TIMEOVFL # UNLIKELY TAU. STORE DELT # DELT=DELINDEP ABS BDSU EPSILONT BPL DLOAD KEPCONVG T DSU NORM TC X1 PDDL NORM DELX X2 XSU,1 DMP X2 DELT SLR* DDV 1,1 SR1 PUSH # 0D=TRIAL DELX PL AT 2 BPL DLOAD POSDELX X STORE XMAX # MOVE MAX BOUND IN BDSU DSU # PL AT 0 XMIN BOV BPL NDXCHNGE NDXCHNGE DLOAD GOTO # Page 1280 0D NEWDELX NDXCHNGE DLOAD DSU XMIN X DMPR GOTO # TO FORCE MPAC +2 TO ZERO DP9/10 NEWDELX POSDELX DLOAD X STORE XMIN # MOVE MIN BOUND IN BDSU DSU # PL AT 0 XMAX BOV BMN PDXCHNGE PDXCHNGE DLOAD 0D NEWDELX STORE DELX BZE DAD KEPCONVG X STODL X T STORE TC BRNCHCTR RTB BHIZ CHECKCTR KEPCONVG GOTO KEPLOOP # ITERATE PDXCHNGE DLOAD DSU XMAX X DMPR GOTO # TO FORCE MPAC +2 TO ZERO DP9/10 NEWDELX BADX DLOAD SR1 XMAX SIGN TAU. STCALL X STORBNDS # Page 1281 TIMEOVFL DLOAD BMN # X WAS TOO BIG X NEGTOVFL STORE XMAX CMNTOVFL DLOAD SR1 DELX STORE DELX BZE BDSU KEPRTN X STODL X TC STCALL T BRNCHCTR NEGTOVFL STCALL XMIN CMNTOVFL KEPCONVG DLOAD SR4R R1 DSU VXSC XSQC(XI) URRECT VSL1 PDDL # 0D=(R1-XSQC(XI))URRECT (+33 OR +31) X DSQ NORM X1 DMPR DMPR 1/ROOTMU X DMP SRR* S(XI) 0 -7,1 BDSU T SL1 VXSC VRECT VSL1 VAD VSL4 STORE RCV # RCV (+29 OR +27) ABVAL NORM X2 STODL RCNORM XI DMPR DSU S(XI) D1/128 DMP SL1R ROOTMU DMP SLR* # Page 1282 X 0 -3,2 DDV VXSC RCNORM URRECT VSL1 PDDL # 0D=URRECT(XI S(XI)-1)X ROOTMU/RCV (+15 XSQC(XI) # OR +13) PL AT 6 SLR* DDV 0 -4,2 RCNORM BDSU VXSC D1/256 VRECT VAD VSL8 STADR # PL AT 0 STODL VCV # VCV (+7 OR +5) T STODL TC X STCALL XPREV KEPRTN # Page 1283 DELTIME EXIT # MPAC=XI (+6), 0D=XSQ (+34 OR +32 -N1) TC POLY DEC 8 2DEC .083333334 2DEC -.266666684 2DEC .406349155 2DEC -.361198675 2DEC .210153242 2DEC -.086221951 2DEC .026268812 2DEC -.006163316 2DEC .001177342 2DEC -.000199055 TC INTPRET STODL S(XI) XI EXIT TC POLY DEC 8 2DEC .031250001 2DEC -.166666719 2DEC .355555413 2DEC -.406347410 2DEC .288962094 2DEC -.140117894 2DEC .049247387 2DEC -.013081923 2DEC .002806389 2DEC -.000529414 TC INTPRET # Page 1284 DMP SRR* # PL AT 0 0D 0 -5,1 STORE XSQC(XI) # XSQC(XI) (+33 OR +31) DMP SL1 KEPC1 RTB PDDL # XCH WITH PL. 0D=C1 XSQ C(XI) (+49 OR +46) TPMODE # PL AT 0,3 DMP SRR* S(XI) 0 -5,1 DMP SL1 KEPC2 RTB PDDL # 3D=C2 XSQ S(XI) (+35 OR +33) PL AT 6 TPMODE R1 SR TAD # PL AT 3 6 NORM DMP # TO PRESERVE SIGNIF. X1 X SR* TAD # X(C2 XSQ S(XI) +R1) (+49 OR +46) PL AT 0 0 -3,1 SL4R DMPR 1/ROOTMU STORE T RVQ # Page 1285 ITERATOR BONCLR DLOAD SLOPESW FIRSTIME DEP DSU NORM DEPREV X1 PDDL NORM DELINDEP X2 XSU,1 DMP X2 DELDEP SLR* DDV # PL UP 2 1,1 SR1 BOFF ORDERSW SGNCHECK ABS SIGN # IN CASE 2ND DERIV. CHANGED SIGN, MUST DELDEP # DISREGARD IT TO FIND MIN. SGNCHECK PUSH BPL # TRIAL DELINDEP PL DOWN 2 POSDEL DLOAD BON INDEP ORDERSW MINCHECK STORE MAX # IF NOT 2ND ORDER, CAN MOVE MAX BOUND IN. MINCHECK BDSU DSU MIN BOV BPL MODNGDEL MODNGDEL GOTO DELOK MODNGDEL DLOAD DSU # TRIAL DELINDEP WOULD EXCEED MIN BOUND MIN INDEP DMP GOTO DP9/10 NEWDEL FIRSTIME DLOAD DMP MIN TWEEKIT # DLOAD TWEEKIT(40D) SENSITIVE TO CHANGE. PDDL DMP # S2(41D) SHOULDN'T CONTAIN HI ORDER ONES # Page 1286 MAX TWEEKIT DSU SIGN GOTO DELDEP SGNCHECK POSDEL DLOAD BON INDEP ORDERSW MAXCHECK STORE MIN # IF NOT 2ND ORDER, CAN MOVE MIN BOUND IN. MAXCHECK BDSU DSU MAX BOV BMN MODPSDEL MODPSDEL DELOK DLOAD 0D NEWDEL STORE DELINDEP RVQ MODPSDEL DLOAD DSU MAX INDEP DMP GOTO DP9/10 NEWDEL CHECKCTR CS ONE INDEX FIXLOC AD ITERCTR INDEX FIXLOC TS ITERCTR TS MPAC TC DANZIG # Page 1287 NEWSTATE DLOAD SR4R R1 DSU VXSC XSQC(XI) UR1 VSL1 PDDL # 0D=(R1-XSQC(XI))UR1 (+33 OR 31) PL AT 6 X DSQ NORM X1 DMPR DMPR 1/ROOTMU X DMP SRR* S(XI) 0 -7,1 BDSU T SL1 VXSC VVEC VSL1 VAD # PL AT 0 VSL4 PUSH ABVAL LAMENTER NORM X1 STODL R2 XI DMP DSU S(XI) D1/128 DMP SL1R ROOTMU DMP SLR* X 0 -3,1 DDV VXSC R2 UR1 VSL1 PDDL # 6D=V2VEC PART (+15 OR 13) PL AT 12 XSQC(XI) SLR* DDV 0 -4,1 R2 BDSU D1/256 VXSC VAD VVEC # PL AT 6 VSL8 RVQ # Page 1288 SETLOC CONICS1 BANK COUNT 04/CONIC # DO NOT DISTURB THE ORDER OF THESE CDS, OVERLAYS HAVE BEEN MADE. BEE17 DEC 0 # KEEP WITH D1/8 2DEC 1.0B-17 (0000004000) D1/8 2DEC 1.0 B-3 D1/128 2DEC 1.0 B-7 D1/64 2DEC 1.0 B-6 D1/4 2DEC 1.0 B-2 D1/16 2DEC 1.0 B-4 D1/32 2DEC 1.0 B-5 D1/1024 2DEC 1.0 B-10 D1/256 2DEC 1.0 B-8 DP9/10 2DEC .9 KEPZERO EQUALS LO6ZEROS -50SC 2DEC -50.0 B-12 2PISC 2DEC 6.28318530 B-6 BEE19 EQUALS D1/32 -1 # 2DEC 1.0 B-19 (00000 01000) BEE22 EQUALS D1/256 -1 # 2DEC 1.0 B-22 (00000 00100) ONEBIT 2DEC 1.0 B-28 COGUPLIM 2DEC .999511597 COGLOLIM 2DEC -.999511597 # Page 1289 SETLOC CONICS BANK COUNT 12/CONIC TIMETHET STQ SETPD # PL AT 0 RTNTT 0 BOV +1 VLOAD PDVL # SETUP FOR PARAM CALL PL AT 6 RVEC VVEC CALL PARAM BOV CALL COGAOVFL # PL AT 0 GETX COMMNOUT DLOAD BON XI INFINFLG ABTCONIC CLEAR CALL COGAFLAG DELTIME BON CALL RVSW RTNTT NEWSTATE GOTO RTNTT COGAOVFL SETGO COGAFLAG ABTCONIC BANK 4 SETLOC CONICS1 BANK COUNT* $$/CONIC PARAM STQ CLEAR # MPAC=V1VEC, 0D=R1VEC PL AT 6 RTNPRM NORMSW CLEAR COGAFLAG SSP CALL GEOMSGN 37777 # GAMMA ALWAYS LESS THAN 180DEG GEOM # MPAC=SNGA (+1), 0D=CSGA (+1) PL AT 2 STODL 36D # 36D=SIN GAMMA (+1) PL AT 0 # Page 1290 SR DDV 5 36D STOVL* COGA MUTABLE,1 STODL 1/MU MAGVEC2 DSQ NORM X1 DMPR DMP 1/MU R1 SRR* 0 -3,1 PUSH BDSU # 0D=R1 V1SQ/MU (+6) PL AT 2 D1/32 STODL R1A # R1A (+6) PL AT 0 DMP NORM 36D X1 DMP SR* 36D 0 -4,1 STCALL P # P (+4) RTNPRM # Page 1291 GEOM UNIT # MPAC=V2VEC, 0D=R1VEC PL AT 6 STODL U2 # U2 (+1) 36D STOVL MAGVEC2 # PL AT 0 UNIT STORE UR1 # UR1 (+1) DOT SL1 U2 PDDL # 0D=CSTH (+1) PL AT 2 36D STOVL R1 # R1 (+29 OR +27) UR1 VXV VSL1 U2 BON SIGN NORMSW HAVENORM GEOMSGN UNIT BOV COLINEAR UNITNORM STODL UN # UN (+1) 36D SIGN RVQ # MPAC=SNTH (+1), 34D=SNTH.SNTH (+2) GEOMSGN COLINEAR VSR1 GOTO UNITNORM HAVENORM ABVAL SIGN GEOMSGN RVQ # MPAC=SNTH (+1), 34D=SNTH.SNTH (+2) # Page 1292 BANK 12 SETLOC CONICS BANK COUNT 12/CONIC GETX AXT,2 SSP # ASSUMES P (+4) IN MPAC 3 S2 1 CLEAR 360SW SQRT PDDL # 0D=SQRT(P) PL AT 2 CSTH SR1 BDSU D1/4 PDDL SRR # PL AT 4D SNTH 6 DDV # PL AT 2 BOV 360CHECK DSU DMP COGA # PL AT 0 SL2R BOV 360CHECK WLOOP PUSH DSQ # 0D=W (+5) PL AT 2 TLOAD PDDL # 2D=WSQ (+10) PL AT 5 MPAC R1A SR4 TAD # PL AT 2 BMN SQRT INFINITY ROUND DAD # PL AT 0D BOV TIX,2 RESETX2 WLOOP BDDV BOV D1/128 INFINITY POLYCOEF BMN PUSH # 0D=1/W (+2) OR 16/W (+6) PL AT 2 INFINITY DSQ NORM DMP X1 R1A SRR* EXIT 0 -10D,1 # Page 1293 TC POLY DEC 5 2DEC .5 2DEC -.166666770 2DEC .100000392 2DEC -.071401086 2DEC .055503292 2DEC -.047264098 2DEC .040694204 TC INTPRET DMP SL1R # PL AT 0D PUSH BON 360SW TRUE360X XCOMMON DSQ NORM X1 DMP SRR* R1A 0 -12D,1 STODL XI # XI (+6) R1 SR1 SQRT ROUND DMP SL4R # PL AT 0 STORE X # X (+17 OR +16) DSQ NORM X1 PDDL DMP # 0D=XSQ (+34 OR +32 -N1) PL AT 2 P R1 SL3 SQRT DMP SL3R COGA STODL KEPC1 R1A BDSU CLEAR D1/64 INFINFLG STORE KEPC2 # Page 1294 RVQ RESETX2 AXT,2 3 360CHECK SETPD BPL 0D INVRSEQN SET 360SW INVRSEQN DLOAD SQRT P PDDL DMP # 0D=SQRT(P) (+2) PL AT 2 SNTH COGA SL1 PDDL # 2D=SNTH COGA (+5) PL AT 4 CSTH SR4 DAD D1/32 DSU DMP # PL AT 2,0 NORM BDDV X1 SNTH SLR* ABS # NOTE: NEAR 360 CASE TREATED DIFFERENTLY 0 -5,1 PUSH DSQ # 0D=1/W (-1) PL AT 2 STODL 34D D1/16 1/WLOOP PUSH DSQ # 2D=G (+4) PL AT 4 RTB PDDL # PL AT 7 TPMODE R1A DMP SR4 34D TAD BMN SQRT # PL AT 4 INFINITY DAD # PL AT 2 TIX,2 NORM 1/WLOOP X1 BDDV SLR* GOTO # PL AT 0 0 -7,1 POLYCOEF # Page 1295 TRUE360X DLOAD BMN R1A INFINITY SQRT NORM X1 BDDV SL* 2PISC 0 -3,1 DSU PUSH # 0D=2PI/SQRT(R1A) -X PL AT 0,2 GOTO XCOMMON INFINITY SETPD BOV # NO SOLUTION EXISTS SINCE CLOSURE THROUGH 0 # INFINITY IS REQUIRED OVFLCLR OVFLCLR SET RVQ INFINFLG # Page 1296 LAMBERT STQ SETPD RTNLAMB 0D BOV +1 SSP VLOAD* ITERCTR 20D MUTABLE,1 STODL 1/MU TDESIRED DMPR BEE19 STORE EPSILONL SET VLOAD SLOPESW R1VEC PDVL CALL # 0D=R1VEC (+29 OR +27) PL AT 6 R2VEC # MPAC=R2VEC (+29 OR +27) GEOM STODL SNTH # 0D=CSTH (+1) PL AT 2 MAGVEC2 NORM PDDL # PL AT 4 X1 R1 SR1 DDV # PL AT 2 SL* PDDL # DXCH WITH 0D, 0D=R1/R2 (+7) PL AT 0,2 0 -6,1 STADR STORE CSTH # CSTH (+1) SR1 BDSU D1/4 STORE 1-CSTH # 1-CSTH (+2) ROUND BZE 360LAMB NORM PDDL # PL AT 4 X1 0D SR1 DDV # PL AT 2 SL* SQRT 0 -3,1 PDDL SR # 2D=SQRT(2R1/R2(1-CSTH)) (+5) PL AT 4 SNTH 6 DDV DAD # PL AT 2 1-CSTH STADR STORE COGAMAX BOV BMN # IF OVFL, COGAMAX=COGUPLIM # Page 1297 UPLIM # IF NEG, USE EVEN IF LT COGLOLIM, SINCE MAXCOGA # THIS WOULD BE RESET IN LAMBLOOP DSU BMN # IF COGAMAX GT COGUPLIM, COGAMAX=COGUPLIM COGUPLIM MAXCOGA # OTHERWISE OK, SO GO TO MAXCOGA UPLIM DLOAD COGUPLIM # COGUPLIM=.999511597 = MAX VALUE OF COGA STORE COGAMAX # NOT CAUSING OVFL IN R1A CALCULATION MAXCOGA DLOAD CSTH SR DSU # PL AT 0 6 STADR STODL CSTH-RHO GEOMSGN BMN DLOAD LOLIM CSTH-RHO SL1 DDV SNTH BOV LOLIM MINCOGA STORE COGAMIN # COGAMIN (+5) BON SSP GUESSW NOGUESS TWEEKIT 00001 DLOAD COGA LAMBLOOP DMP SNTH SR1 DSU CSTH-RHO NORM PDDL # 0D=SNTH COGA-(CSTH-RHO) (+7+C(XI)) PL=2 X1 1-CSTH SL* DDV # 1-CSTH (+2) PL AT 0 0 -9D,1 BMN BZE NEGP NEGP STODL P # P=(1-CSTH)/(SNTH COGA-(CSTH-RHO)) (+4) COGA DSQ DAD D1/1024 NORM DMP X1 # Page 1298 P SR* BDSU 0 -8D,1 D1/32 STODL R1A # R1A=2-P(1+COGA COGA) (+6) P BOV CALL HIENERGY GETX DLOAD T STODL TPREV XI BON CALL INFINFLG NEGP # HAVE EXCEEDED THEORETICAL BOUNDS DELTIME BOV BDSU BIGTIME TDESIRED STORE TERRLAMB ABS BDSU EPSILONL BPL RTB INITV CHECKCTR BHIZ CALL SUFFCHEK ITERATOR DLOAD BZE MPAC SUFFCHEK DAD COGA STCALL COGA LAMBLOOP NEGP DLOAD BPL # IMPOSSIBLE TRAJECTORY DUE TO INACCURATE DCOGA # BOUND CALCULATION. TRY NEW COGA. LOENERGY HIENERGY SETPD DLOAD # HIGH ENERGY TRAJECTORY RESULTED. 0 COGA # IN OVFL OF P OR R1A, OR XI EXCEEDING 50. STORE COGAMIN # THIS IS THE NEW BOUND. COMMONLM DLOAD SR1 DCOGA # Page 1299 STORE DCOGA # USE DCOGA/2 AS DECREMENT BZE BDSU SUFFCHEK COGA STCALL COGA LAMBLOOP BIGTIME DLOAD TPREV STORE T LOENERGY SETPD DLOAD # LOW ENERGY TRAJECTORY RESULTED 0 COGA # IN OVERFLOW OF TIME. STCALL COGAMAX # THIS IS THE NEW BOUND. COMMONLM SUFFCHEK DLOAD ABS TERRLAMB PDDL DMP TDESIRED # PL AT 2D D1/4 DAD DSU # PL AT 0D ONEBIT BPL SETGO INITV SOLNSW RTNLAMB 360LAMB SETPD SETGO # LAMBERT CANNOT HANDLE CSTH=1 0 SOLNSW RTNLAMB NOGUESS SSP DLOAD TWEEKIT 20000 COGAMIN SR1 PDDL # PL AT 2 COGAMAX SR1 DAD STADR # PL AT 0 STORE COGA STCALL DCOGA LAMBLOOP # Page 1300 LOLIM DLOAD GOTO COGLOLIM # COGLOLIM=-.999511597 MINCOGA INITV DLOAD NORM R1 X1 PDDL SR1 # PL AT 2 P DDV # PL AT 0 SL* SQRT 0 -4,1 DMP SL1 ROOTMU PUSH DMP # 0D=VTAN (+7) PL AT 2 COGA SL VXSC 5 UR1 PDDL # XCH WITH 0D PL AT 0,6 VXSC VSL1 UN VXV VAD UR1 # PL AT 0 VSL1 CLEAR SOLNSW STORE VVEC SLOAD BZE VTARGTAG TARGETV GOTO RTNLAMB TARGETV DLOAD CALL MAGVEC2 LAMENTER STCALL VTARGET RTNLAMB # Page 1301 TIMERAD STQ SETPD # PL AT 0 RTNTR 0 BOV +1 VLOAD PDVL # PL AT 6 RVEC VVEC CALL PARAM BOV DLOAD COGAOVFL # PL AT 0 D1/32 DSU DMP R1A P SQRT DMP COGA SL4 VXSC U2 PDDL DSU # PL AT 6 D1/64 R1A VXSC VSU # PL AT 0 UR1 VSL4 UNIT BOV CIRCULAR PDDL NORM # 0D=UNIT(ECC) (+3) PL AT 6 RDESIRED # 35D=ECC (+3) X1 PDDL DMP # PL AT 8 R1 P SL* DDV # PL AT 6 0,1 DSU DDV D1/16 36D # 36D=ECC (+3) STORE COSF BOV DSQ BADR2 BDSU BMN D1/4 BADR2 SQRT SIGN SGNRDOT CLEAR APSESW # Page 1302 TERMNVEC VXSC VSL1 UN VXV PDVL # VXCH WITH 0D PL AT 0,6 0D VXSC VAD # PL AT 0 COSF VSL1 PUSH # 0D=U2 PL AT 6 DOT DDV # LIMITS RESULT TO POSMAX OR NEGMAX UR1 DP1/4 SR1 BOV # SCALE BACK DOWN TO NORMAL +1 # CLEAR OVFIND IF SET STOVL CSTH # CSTH (+1) UR1 VXV VSL1 DOT SL1 UN STODL SNTH # SNTH (+1) P CALL GETX CLRGO SOLNSW COMMNOUT CIRCULAR SETPD SETGO 0 SOLNSW ABTCONIC BADR2 DLOAD SIGN LODPHALF COSF STODL COSF KEPZERO SETGO APSESW TERMNVEC # Page 1303 APSIDES STQ SETPD # PL AT 0 RTNAPSE 0D BOV +1 VLOAD PDVL # PL AT 6 RVEC VVEC CALL PARAM BOV # PL AT 0 GETECC GETECC DMP SL4 R1A BDSU SQRT D1/64 STORE ECC DAD PDDL # PL AT 2 D1/8 R1 DMP SL1 P DDV # PL AT 0 PDDL NORM # 0D=RP (+29 OR +27) PL AT 2 R1A X1 PDDL SL* # PL AT 4 R1 0 -5,1 DDV DSU # PL AT 2,0 BOV BMN INFINAPO INFINAPO GOTO RTNAPSE INFINAPO DLOAD GOTO # RETURNS WITH APOAPSIS IN MPAC, PERIAPSIS LDPOSMAX RTNAPSE # THAT PL IS AT 0. # Page 1304 ABTCONIC EXIT TC P00DOO OCT 00607 # Page 1305 SETLOC CONICS1 BANK COUNT 04/CONIC MUTABLE 2DEC* 3.986032 E10 B-36* # MUE 2DEC* .25087606 E-10 B+34* # 1/MUE 2DEC* 1.99650495 E5 B-18* # SQRT(MUE) 2DEC* .50087529 E-5 B+17* # 1/SQRT(MUE) 2DEC 4.902778 E8 B-30 # MUM 2DEC .203966 E-8 B+28 # 1/MUM 2DEC* 2.21422176 E4 B-15* # SQRT(MUM) 2DEC* .45162595 E-4 B+14* # 1/SQRT(MUM) LDPOSMAX EQUALS LODPMAX # DPPOSMAX IN LOW MEMORY. # ERASABLE ASSIGNMENTS # KEPLER SUBROUTINE # INPUT -- # RRECT ERASE +5 # VRECT ERASE +5 # TAU. ERASE +1 # XKEP ERASE +1 # TC ERASE +1 # XPREV ERASE +1 1/MU EQUALS 14D ROOTMU EQUALS 16D 1/ROOTMU EQUALS 18D # OUTPUT -- # RCV ERASE +5 # VCV ERASE +5 # RC ERASE +1 # XPREV ERASE +1 # DEBRIS -- ALPHA EQUALS 8D XMAX EQUALS 10D # Page 1306 XMIN EQUALS 12D X EQUALS 20D XI EQUALS 24D S(XI) EQUALS 26D XSQC(XI) EQUALS 28D T EQUALS 30D R1 EQUALS 32D KEPC1 EQUALS 34D KEPC2 EQUALS 36D # DELX ERASE +1 # DELT ERASE +1 # URRECT ERASE +5 # RCNORM ERASE +1 # XPREV EQUALS XKEP # LAMBERT SUBROUTINE # # INPUT -- # R1VEC ERASE +5 # R2VEC ERASE +5 # TDESIRED ERASE +1 # GEOMSGN ERASE +0 # GUESSW # 0 IF COGA GUESS AVAILABLE, 1 IF NOT # COGA ERASE +1 # INPUT ONLY IF GUESS IS ZERO. # NORMSW # 0 IF UN TO BE COMPUTED, 1 IF UN INPUT # UN ERASE +5 # ONLY USED IF NORMSW IS 1 # VTARGTAG ERASE +0 # TWEEKIT EQUALS 40D # ONLY USED IF GUESSW IS 0 # OUTPUT -- # VTARGET ERASE +5 # AVAILABLE ONLY IF VTARGTAG IS ZERO. # V1VEC EQUALS MPAC # DEBRIS -- # RTNLAMB ERASE +0 # U2 ERASE +5 # MAGVEC2 ERASE +1 # UR1 ERASE +5 # R1 EQUALS 31D # UN ERASE +5 # SNTH ERASE +1 # CSTH ERASE +1 # 1-CSTH ERASE +1 # CSTH-RHO ERASE +1 COGAMAX EQUALS 14D # CLOBBERS 1/MU COGAMIN EQUALS 8D DCOGA EQUALS 12D # TWEEKIT EQUALS 40D # P ERASE +1 # Page 1307 # COGA ERASE +1 # R1A ERASE +1 # X EQUALS 20D # XSQ EQUALS 22D # XI EQUALS 24D # S(XI) EQUALS 26D # XSQC(XI) EQUALS 28D # T EQUALS 30D # KEPC1 EQUALS 34D # KEPC2 EQUALS 36D # SLOPSW # SOLNSW # OTHERS -- # RVEC EQUALS R1VEC # VVEC ERASE +5 # COGAFLAG # RVSW # INFINFLG # APSESW # 360SW # RTNTT EQUALS RTNLAMB # ECC ERASE +1 # RTNTR EQUALS RTNLAMB # RTNAPSE EQUALS RTNLAMB # R2 EQUALS MAGVEC2 COSF EQUALS 24D # RTNPRM ERASE +0 # SCNRDOT ERASE +0 # RDESIRED ERASE +1 # ITERATOR SUBROUTINE # ORDERSW MAX EQUALS 14D # CLOBBERS 1/MU MIN EQUALS 8D # INDEP ERASE +1 DELINDEP EQUALS 12D ITERCTR EQUALS 22D DEP EQUALS 30D # DELDEP ERASE +1 # DEPREV ERASE +1 TWEEKIT EQUALS 40D # MORE KEPLER # EPSILONT ERASE +1 # Page 1308 # MORE LAMBERT # TERRLAMB EQUALS DELDEP # TPREV EQUALS DEPREV # EPSILONL EQUALS EPSILONT +2 # DOUBLE PRECISION WORD