kjodle/Apollo-11
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Proof CONIC_SUBROUTINES

Browse Source
This commit is contained in:
Thomas Danner 2018-12-28 12:33:22 -06:00
parent 2e30258105
commit f7f8dbc625
1 changed files with 45 additions and 45 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

90
Comanche055/CONIC_SUBROUTINES.agc
View File

@ -41,7 +41,7 @@
# RATHER THAN OPTIMIZING EACH FOR A PARTICULAR USE. THEREFORE, MULTIPLE USAGE CAN BE MADE OF 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 # 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 # 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 # THE INHERENT SCALE CHANGE REQUIRED IN POSITION, VELOCITY, MU, AND TIME, AS OUTLINED IN MISSION PROGRAMMING
# DEFINITION MEMO NO. 10. THIS CAN BE ACCOMPLISHED BY SIMPLY ADDING TO THE MUTABLE AND INITIALIZING THE SUBROUTINES # DEFINITION MEMO NO. 10. THIS CAN BE ACCOMPLISHED BY SIMPLY ADDING TO THE MUTABLE AND INITIALIZING THE SUBROUTINES
# APPROPRIATELY. # APPROPRIATELY.
# #
@ -77,8 +77,8 @@
# IF A NEGATIVE TIME-OF-FLIGHT IS INPUT, THE PROGRAM WILL SOLVE FOR THE STATE WHICH WOULD BE PRODUCED BY # 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. # EXTRAPOLATING THE POSITION BACKWARD IN TIME.
# #
# IF THE ABSOLUTE VALUE DESIRED TRANSFER TIME EXCEEDS THE ORBITAL PERIOD, THE SUBROUTINE, THROUGH A # IF THE ABSOLUTE VALUE OF THE DESIRED TRANSFER TIME EXCEEDS THE ORBITAL PERIOD, THE SUBROUTINE, THROUGH A
# MODULAR TECHNIQUE, WILL COMPUTE THE STATE CORRESPONDING TO THE DESIRED TIME AS USUAL. # MODULAR TECHNIQUE, WILL COMPUTE THE STATE CORRESPONDING TO THE DESIRED TIME (WHETHER POSITIVE OR NEGATIVE).
# #
# THE RESTRICTIONS ARE -- # THE RESTRICTIONS ARE --
# 1. (PREVIOUS RESTRICTION ON THE NEGATIVE DESIRED TRANSFER TIME IS NOW DELETED.) # 1. (PREVIOUS RESTRICTION ON THE NEGATIVE DESIRED TRANSFER TIME IS NOW DELETED.)
@ -86,7 +86,7 @@
# ANY OF THESE LIMITS ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS. # 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 # 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 # GUESS, XKEPNEW. THE AGC COMPUTATION TIME IS APPROXIMATELY .061 SECONDS FOR INITIALIZATION, .065 SECONDS FOR THE
# FINAL COMPUTATIONS, PLUS .083 SECONDS FOR EACH ITERATION. # FINAL COMPUTATIONS, PLUS .083 SECONDS FOR EACH ITERATION.
# #
# REFERENCES -- # REFERENCES --
@ -126,7 +126,7 @@
# THE INTERRUPTED JOB. THEREFORE THE USER MUST CALL CSMCONIC OR LEMCONIC WHICH GUARANTEES NO INTERRUPTS AND WHICH # 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. # ALSO CALLS KEPPREP TO COMPUTE A GUESS OF XKEPNEW.
# #
# ABORT EXIT MODE -- # ABORT EXIT MODES --
# NONE # NONE
# #
# OUTPUT -- # OUTPUT --
@ -151,13 +151,13 @@
# URRECT +1 DP UNIT VECTOR OF INITIAL POSITION # URRECT +1 DP UNIT VECTOR OF INITIAL POSITION
# R1 +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION IN METERS # R1 +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION IN METERS
# +27 FOR MOON # +27 FOR MOON
# ALPHA -22 FOR EARTH DP INVERSE OF SEMI-MAJOR AXIS IN 1/METERS # ALPHA -22 FOR EARTH DP INVERSE OF SEMIMAJOR AXIS IN 1/METERS
# -20 FOR MOON # -20 FOR MOON
# TMODULO +28 DP INTEGRAL NUMBER OF PERIODS IN CENTISECS, WHICH WAS SUBTRACTED FROM TAU. TO PRODUCE A # TMODULO +28 DP INTEGRAL NUMBER OF PERIODS IN CENTISECS, WHICH WAS SUBTRACTED FROM TAU. TO PRODUCE A
# TAU. OF LESS THAN ONE PERIOD. # TAU. OF LESS THAN ONE PERIOD.
# #
# PARAMETERS OF NO USE -- # PARAMETERS OF NO USE --
# DP PARAMETERS -- FPSILENT, DELX, DELT, RCNORM, XMODULO, PLUS PUSHLIST REGISTERS 0 THROUGH 39D. # DP PARAMETERS -- EPSILONT, DELX, DELT, RCNORM, XMODULO, PLUS PUSHLIST REGISTERS 0 THROUGH 39D.
# Page 1266 # Page 1266
# PROGRAM DESCRIPTION -- LAMBERT SUBROUTINE DATE -- 1 SEPTEMBER 1967 # PROGRAM DESCRIPTION -- LAMBERT SUBROUTINE DATE -- 1 SEPTEMBER 1967
@ -181,7 +181,7 @@
# 5. THE PARAMETERS IN THE PROBLEM MUST NOT EXCEED THEIR SCALING LIMITS SPECIFIED IN THE GSOP. IF THE # 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. # 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 # THE NUMBER OF ITERATIONS AND, THEREFORE, THE COMPUTATIONS SPEED IS DEPENDENT ON THE ACCURACY OF THE FIRST
# GUESS OF THE INDEPENDENT VARIABLE, COGA. THE AGC COMPUTATION TIME IS APPROXIMATELY # 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. # .105 SECONDS FOR INITIALIZATION, .069 SECONDS FOR FINAL COMPUTATIONS, PLUS .205 SECONDS FOR EACH ITERATION.
# #
@ -204,13 +204,13 @@
# GUESSW NONE AN INTERPRETER SWITCH TO BE SET IF NO GUESS OF COGA IS AVAILABLE, CLEAR IF A GUESS OF # GUESSW NONE AN INTERPRETER SWITCH TO BE SET IF NO GUESS OF COGA IS AVAILABLE, CLEAR IF A GUESS OF
# Page 1267 # Page 1267
# COGA IS TO BE USED BY LAMBERT # COGA IS TO BE USED BY LAMBERT
# COGA +5 DP GUESS OF COTANGENT OF FLIGHT PATH ANGLE (MEASURED FROM VERTICAL). THIS WILL BE # COGA +5 DP GUESS OF COTANGNT OF FLIGHT PATH ANGLE (MEASURED FROM VERTICAL). THIS WILL BE
# IGNORED IF GUESSW IS SET. # 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 # 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). # 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 # 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. # 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 # VTARGTAG NONE A S.P. TAG TO BE SET TO ZERO IF LAMBERT IS TO COMPUTE THE VELOCITY AT R2VEC AS WELL AS
# AT R1VEC. # AT R1VEC.
# #
# SUBROUTINES CALLED -- # SUBROUTINES CALLED --
@ -224,16 +224,16 @@
# L+3 SOLNSW # L+3 SOLNSW
# L+4 LAMABORT # L+4 LAMABORT
# #
# IF A LAMBER RESULT IS TO BE A FIRST GUESS FOR THE NEXT LAMBERT CALCULATION, COGA MUST BE PRESERVED AND # IF A LAMBERT 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. # GUESSW MUST BE CLEAR FOR EACH SUCCEEDING LAMBERT CALL.
# #
# ABORT EXIT MODES -- # 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 # 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 # WAS TOO SMALL TO PRODUCE A REALISTIC TRANSFER BETWEEN R1VEC AND R2VEC. IN EITHER CASE THE FIX MUST BE MADE
# ACCORDING TO THE NEEDS OF THE PARTICULAR USER. THE ABORT EXIT MODE MAY BE CODED AS ... # 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 # LAMBERT DLOAD ABS # A MEASURE OF THE PROXIMITY TO 0 OR
# 1-CHTH # 360 DEGREES. # 1-CSTH # 360 DEGREES.
# DSU BWM # DSU BMN
# ONEBIT # ONEBIT
# CHANGER2 # CHANGE R2VEC DIRECTION SLIGHTLY. # CHANGER2 # CHANGE R2VEC DIRECTION SLIGHTLY.
# DLOAD DAD # DLOAD DAD
@ -266,7 +266,7 @@
# 1-CSTH +2 DP 1-CSTH # 1-CSTH +2 DP 1-CSTH
# COGA +5 DP COTAN OF INITIAL REQUIRED FLIGHT PATH ANGLE MEASURED FROM VERTICAL # COGA +5 DP COTAN OF INITIAL REQUIRED FLIGHT PATH ANGLE MEASURED FROM VERTICAL
# P +4 DP RATIO OF SEMILATUS RECTUM TO INITIAL RADIUS # P +4 DP RATIO OF SEMILATUS RECTUM TO INITIAL RADIUS
# R1A +6 DP RATIO OF INITIAL RADIUS TO SEMI-MAJOR AXIS # R1A +6 DP RATIO OF INITIAL RADIUS TO SEMIMAJOR AXIS
# R1 (32D) +29 FOR EARTH DP INITIAL RADIUS IN METERS # R1 (32D) +29 FOR EARTH DP INITIAL RADIUS IN METERS
# +27 FOR MOON # +27 FOR MOON
# UR1 +1 DP UNIT VECTOR OF R1VEC # UR1 +1 DP UNIT VECTOR OF R1VEC
@ -292,7 +292,7 @@
# THE RESTRICTIONS ARE -- # THE RESTRICTIONS ARE --
# 1. THE ANGLE BETWEEN ANY POSITION VECTOR AND ITS VELOCITY VECTOR MUST BE GREATER THAN 1 DEGREE 47.5 MINUTES # 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. # 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 # 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. # ARE EXCEEDED, THE RESULTING SOLUTION WILL BE MEANINGLESS.
# #
# THE AGC COMPUTATION TIME IS APPROXIMATELY .292 SECONDS. # THE AGC COMPUTATION TIME IS APPROXIMATELY .292 SECONDS.
@ -308,7 +308,7 @@
# +27 FOR MOON # +27 FOR MOON
# VVEC +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND # VVEC +7 FOR EARTH DP INITIAL VELOCITY VECTOR IN METERS/CENTISECOND
# +5 FOR MOON # +5 FOR MOON
# SNTH +1 ` DP SINE OF THE TRUE-ANOMALY-DIFFERENCE THROUGH WHICH THE STATE IS TO BE UPDATED # SNTH +1 DP SINE OF TRUE-ANOMALY-DIFFERENCE THROUGH WHICH THE STATE IS TO BE UPDATED
# CSTH +1 DP COSINE OF THE ANGLE # 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 # 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. # IS TO BE COMPUTED ALSO.
@ -333,7 +333,7 @@
# # THE INITIAL VELOCITY VECTOR IN MPAC. # # THE INITIAL VELOCITY VECTOR IN MPAC.
# L+3 STOVL NEWVVEC # L+3 STOVL NEWVVEC
# L+4 STADR # L+4 STADR
# L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USER'S LOCATIONS. # L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USERS LOCATIONS.
# L+6 ... # CONTINUE. # L+6 ... # CONTINUE.
# #
# ABORT EXIT MODES -- # ABORT EXIT MODES --
@ -344,7 +344,7 @@
# SCALE FACTOR # SCALE FACTOR
# VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS # VARIABLE IN POWERS OF 2 DESCRIPTION AND REMARKS
# -------- -------------- ----------------------- # -------- -------------- -----------------------
# T(30D) +28 DP TRANSFER TIME IN CENTISECONDS # T (30D) +28 DP TRANSFER TIME IN CENTISECONDS
# INFINFLG NONE AN INTERPRETIVE SWITCH WHICH IS SET IF THE TRANSFER ANGLE REQUIRES CLOSURE THROUGH # 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. # 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), # COGAFLAG NONE AN INTERPRETIVE SWITCH WHICH IS SET IF RESTRICTION 1 HAS BEEN VIOLATED (NO SOLUTION),
@ -366,7 +366,7 @@
# -------- -------------- ----------------------- # -------- -------------- -----------------------
# R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS # R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS
# +27 FOR MOON # +27 FOR MOON
# R1A +6 DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES) # R1A +6 DP RATIO OF R1 TO SEMIMAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES)
# P +4 DP RATIO OF SEMILATUS RECTUM TO R1 # P +4 DP RATIO OF SEMILATUS RECTUM TO R1
# COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC # COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC
# UR1 +1 DP UNIT VECTOR OF RVEC # UR1 +1 DP UNIT VECTOR OF RVEC
@ -443,7 +443,7 @@
# # THE INITIAL VELOCITY VECTOR IN MPAC. # # THE INITIAL VELOCITY VECTOR IN MPAC.
# L+3 STOVL NEWVVEC # L+3 STOVL NEWVVEC
# L+4 STADR # L+4 STADR
# L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USER'S LOCATIONS. # L+5 STORE NEWRVEC # NEWVVEC AND NEWRVEC ARE SYMBOLIC REPRESENTATIONS OF THE USERS LOCATIONS.
# L+6 ... # CONTINUE # L+6 ... # CONTINUE
# #
# ABORT EXIT MODES -- # ABORT EXIT MODES --
@ -483,7 +483,7 @@
# -------- -------------- ----------------------- # -------- -------------- -----------------------
# R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS # R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS
# +27 FOR MOON # +27 FOR MOON
# R1A +6 DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES) # R1A +6 DP RATIO OF R1 TO SEMIMAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES)
# P +4 DP RATIO OF SEMILATUS RECTUM TO R1 # P +4 DP RATIO OF SEMILATUS RECTUM TO R1
# COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC # COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC
# UR1 +1 DP UNIT VECTOR OF RVEC # UR1 +1 DP UNIT VECTOR OF RVEC
@ -493,7 +493,7 @@
# SNTH +1 DP SINE OF TRUE ANOMALY DIFFERENCE. # SNTH +1 DP SINE OF TRUE ANOMALY DIFFERENCE.
# #
# PARAMETERS OF NO USE -- # PARAMETERS OF NO USE --
# SP PARAMETERS -- RTNTT, GEOMSGN, RTNPRM, MAGVEC2*R2 (DP), PLUS PUSHLIST LOCATIONS 0-11D, 14D-21D, 24D-39D, 41D # SP PARAMETERS -- RTNTT, GEOMSGN, RTNPRM, MAGVEC2=R2 (DP), PLUS PUSHLIST LOCATIONS 0-11D, 14D-21D, 24D-39D, 41D
# ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW, 360SW # ADDITIONAL INTERPRETIVE SWITCHES USED -- NORMSW, 360SW
# #
@ -541,7 +541,7 @@
# L+1 APSIDES # RETURNS WITH PL AT 0, RADIUS OF APOCENTER IN MPAC AND RADIUS OF PERICENTER IN 0D # L+1 APSIDES # RETURNS WITH PL AT 0, RADIUS OF APOCENTER IN MPAC AND RADIUS OF PERICENTER IN 0D
# L+2 STODL APOAPSE # L+2 STODL APOAPSE
# L+3 0D # L+3 0D
# L+4 STORE PERIAPSE # APOAPSE AND PERIAPSE ARE SYMBOLIC REPRESENTATIONS OF THE USER'S LOCATIONS # L+4 STORE PERIAPSE # APOAPSE AND PERIAPSE ARE SYMBOLIC REPRESENTATIONS OF THE USERS LOCATIONS
# L+5 ... # CONTINUE # L+5 ... # CONTINUE
# #
# OUTPUT -- # OUTPUT --
@ -563,7 +563,7 @@
# -------- -------------- ----------------------- # -------- -------------- -----------------------
# R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS # R1 (32D) +29 FOR EARTH DP MAGNITUDE OF INITIAL POSITION VECTOR, RVEC, IN METERS
# +27 FOR MOON # +27 FOR MOON
# R1A +6 DP RATIO OF R1 TO SEMI-MAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES) # R1A +6 DP RATIO OF R1 TO SEMIMAJOR AXIS (NEG. FOR HYPERBOLIC TRAJECTORIES)
# P +4 DP RATIO OF SEMILATUS RECTUM TO R1 # P +4 DP RATIO OF SEMILATUS RECTUM TO R1
# COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC # COGA +5 DP COTAN OF ANGLE BETWEEN RVEC AND VVEC
# UR1 +1 DP UNIT VECTOR OF RVEC # UR1 +1 DP UNIT VECTOR OF RVEC
@ -810,7 +810,7 @@ KEPCONVG DLOAD SR4R
T T
SL1 VXSC SL1 VXSC
VRECT VRECT
VSL1 VAD VSL1 VAD # PL AT 0
VSL4 VSL4
STORE RCV # RCV (+29 OR +27) STORE RCV # RCV (+29 OR +27)
@ -975,7 +975,7 @@ MODNGDEL DLOAD DSU # TRIAL DELINDEP WOULD EXCEED MIN BOUND
FIRSTIME DLOAD DMP FIRSTIME DLOAD DMP
MIN MIN
TWEEKIT # DLOAD TWEEKIT(40D) SENSITIVE TO CHANGE. TWEEKIT # DLOAD TWEEKIT(40D) SENSITIVE TO CHANGE.
PDDL DMP # S2(41D) SHOULDN'T CONTAIN HI ORDER ONES PDDL DMP # S2(41D) SHOULDNT CONTAIN HI ORDER ONES
# Page 1286 # Page 1286
MAX MAX
TWEEKIT TWEEKIT
@ -1060,8 +1060,8 @@ LAMENTER NORM
R2 R2
BDSU BDSU
D1/256 D1/256
VXSC VAD VXSC VAD # PL AT 6
VVEC # PL AT 6 VVEC
VSL8 RVQ VSL8 RVQ
# Page 1288 # Page 1288
@ -1118,8 +1118,8 @@ TIMETHET STQ SETPD # PL AT 0
VVEC VVEC
CALL CALL
PARAM PARAM
BOV CALL BOV CALL # PL AT 0
COGAOVFL # PL AT 0 COGAOVFL
GETX GETX
COMMNOUT DLOAD BON COMMNOUT DLOAD BON
XI XI
@ -1349,15 +1349,15 @@ INVRSEQN DLOAD SQRT
R1A R1A
DMP SR4 DMP SR4
34D 34D
TAD TAD # PL AT 4
BMN SQRT # PL AT 4 BMN SQRT
INFINITY INFINITY
DAD # PL AT 2 DAD # PL AT 2
TIX,2 NORM TIX,2 NORM
1/WLOOP 1/WLOOP
X1 X1
BDDV BDDV # PL AT 0
SLR* GOTO # PL AT 0 SLR* GOTO
0 -7,1 0 -7,1
POLYCOEF POLYCOEF
@ -1549,8 +1549,8 @@ LOENERGY SETPD DLOAD # LOW ENERGY TRAJECTORY RESULTED
SUFFCHEK DLOAD ABS SUFFCHEK DLOAD ABS
TERRLAMB TERRLAMB
PDDL DMP PDDL DMP # PL AT 2D
TDESIRED # PL AT 2D TDESIRED
D1/4 D1/4
DAD DSU # PL AT 0D DAD DSU # PL AT 0D
ONEBIT ONEBIT
@ -1598,8 +1598,8 @@ INITV DLOAD NORM
PDDL # XCH WITH 0D PL AT 0,6 PDDL # XCH WITH 0D PL AT 0,6
VXSC VSL1 VXSC VSL1
UN UN
VXV VAD VXV VAD # PL AT 0
UR1 # PL AT 0 UR1
VSL1 CLEAR VSL1 CLEAR
SOLNSW SOLNSW
STORE VVEC STORE VVEC
@ -1626,8 +1626,8 @@ TIMERAD STQ SETPD # PL AT 0
VVEC VVEC
CALL CALL
PARAM PARAM
BOV DLOAD BOV DLOAD # PL AT 0
COGAOVFL # PL AT 0 COGAOVFL
D1/32 D1/32
DSU DMP DSU DMP
R1A R1A
@ -1645,7 +1645,7 @@ TIMERAD STQ SETPD # PL AT 0
BOV BOV
CIRCULAR CIRCULAR
PDDL NORM # 0D=UNIT(ECC) (+3) PL AT 6 PDDL NORM # 0D=UNIT(ECC) (+3) PL AT 6
RDESIRED # 35D=ECC (+3) RDESIRED # 36D=ECC (+3)
X1 X1
PDDL DMP # PL AT 8 PDDL DMP # PL AT 8
R1 R1
@ -1826,7 +1826,7 @@ KEPC2 EQUALS 36D
# TDESIRED ERASE +1 # TDESIRED ERASE +1
# GEOMSGN ERASE +0 # GEOMSGN ERASE +0
# GUESSW # 0 IF COGA GUESS AVAILABLE, 1 IF NOT # GUESSW # 0 IF COGA GUESS AVAILABLE, 1 IF NOT
# COGA ERASE +1 # INPUT ONLY IF GUESS IS ZERO. # COGA ERASE +1 # INPUT ONLY IF GUESSW IS ZERO.
# NORMSW # 0 IF UN TO BE COMPUTED, 1 IF UN INPUT # NORMSW # 0 IF UN TO BE COMPUTED, 1 IF UN INPUT
# UN ERASE +5 # ONLY USED IF NORMSW IS 1 # UN ERASE +5 # ONLY USED IF NORMSW IS 1
# VTARGTAG ERASE +0 # VTARGTAG ERASE +0
@ -1865,7 +1865,7 @@ DCOGA EQUALS 12D
# T EQUALS 30D # T EQUALS 30D
# KEPC1 EQUALS 34D # KEPC1 EQUALS 34D
# KEPC2 EQUALS 36D # KEPC2 EQUALS 36D
# SLOPSW # SLOPESW
# SOLNSW # SOLNSW
# OTHERS -- # OTHERS --
@ -1885,7 +1885,7 @@ DCOGA EQUALS 12D
COSF EQUALS 24D COSF EQUALS 24D
# RTNPRM ERASE +0 # RTNPRM ERASE +0
# SCNRDOT ERASE +0 # SGNRDOT ERASE +0
# RDESIRED ERASE +1 # RDESIRED ERASE +1