Skip to content

Instantly share code, notes, and snippets.

@daniestevez
Created May 3, 2020 15:33
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save daniestevez/7f6c77293b2adb7761ddeb4abdfa5710 to your computer and use it in GitHub Desktop.
Save daniestevez/7f6c77293b2adb7761ddeb4abdfa5710 to your computer and use it in GitHub Desktop.
GMAT's delta-range simulation debug output
Moderator is updating data files...
Moderator is creating core engine...
libmx.so: cannot open shared object file: No such file or directory
*** Library "../plugins/libMatlabInterface" did not open.
../plugins/proprietary/libMsise86.so: cannot open shared object file: No such file or directory
*** Library "../plugins/proprietary/libMsise86" did not open.
../plugins/proprietary/libNRLMsise00.so: cannot open shared object file: No such file or directory
*** Library "../plugins/proprietary/libNRLMsise00" did not open.
../plugins/proprietary/libSNOptimizer.so: cannot open shared object file: No such file or directory
*** Library "../plugins/proprietary/libSNOptimizer" did not open.
../plugins/proprietary/libVF13Optimizer.so: cannot open shared object file: No such file or directory
*** Library "../plugins/proprietary/libVF13Optimizer" did not open.
Successfully set Planetary Source to use: DE405
Successfully set Planetary Source to use: DE405
Successfully set Planetary Source to use: DE405
Setting nutation file to ../data/planetary_coeff/NUTATION.DAT
Setting leap seconds file to ../data/time/tai-utc.dat
2020-05-03 17:32:48 GMAT Moderator successfully created core engine
Successfully set Planetary Source to use: DE405
Successfully set Planetary Source to use: DE405
(GMAT-R2019aBeta1:3593): Gtk-CRITICAL **: 17:32:48.294: IA__gtk_widget_set_size_request: assertion 'width >= -1' failed
2020-05-03 17:32:48 GMAT GUI successfully launched.
*** Error loading "libMatlabInterface": GMAT could not load the MATLAB interface. Please check that MATLAB is installed and configured correctly. See the User Guide for instructions, or disable the appropriate PLUGIN line in bin/gmat_startup_file.txt to remove this error.
*** Unable to load the dynamic library "../plugins/proprietary/libMsise86"
*** Unable to load the dynamic library "../plugins/proprietary/libNRLMsise00"
*** Unable to load the dynamic library "../plugins/proprietary/libSNOptimizer"
*** Unable to load the dynamic library "../plugins/proprietary/libVF13Optimizer"
Interpreting scripts from the file.
***** file: /home/daniel/gmat-dslwp/scripts/test-deltarange.script
Successfully set Planetary Source to use: DE405
Successfully set Planetary Source to use: DE405
Successfully interpreted the script
Running mission...
Successfully set Planetary Source to use: DE405
Successfully set Planetary Source to use: DE405
Kernel ../data/planetary_ephem/spk/DE405AllPlanets.bsp has been loaded.
Kernel ../data/planetary_coeff/SPICEPlanetaryConstantsKernel.tpc has been loaded.
Kernel ../data/time/SPICELeapSecondKernel.tls has been loaded.
Kernel ../data/planetary_coeff/SPICEEarthPredictedKernel.bpc has been loaded.
Kernel ../data/planetary_coeff/SPICEEarthCurrentKernel.bpc has been loaded.
Kernel ../data/planetary_coeff/earth_latest_high_prec.bpc has been loaded.
Kernel ../data/planetary_coeff/SPICELunaCurrentKernel.bpc has been loaded.
Kernel ../data/planetary_coeff/SPICELunaFrameKernel.tf has been loaded.
Using default potential file "../data/gravity/luna/LP165P.cof" for GravityField object "GravityField.Luna"
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Hardware delay calculation for leg from DSLWP_B to PI9CAM :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
. DSLWP_B's hardware delay : 0.000000000000e+00 sec
. PI9CAM's hardware delay : 0.000000000000e+00 sec
. Total hardware delay for this signal leg: 0.000000000000e+00 sec
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Range, relativity correction, and ET-TAI correction calculation for leg from DSLWP_B to PI9CAM at fixed time 28662.7087619720104166697604 A1Mjd:
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1. Compute Range Vector before light time correction for the Leg from <TNode = DSLWP_B> to <RNode = PI9CAM>:
. DSLWP_B state in LunaMJ2000 cs : (8094.881620068918, -7125.732028219994, -9809.286749897787, 0.124798852556, 0.305404411570, -0.022405811802)km at time tT= 28662.7087619720104166697604
. PI9CAM state in EarthMJ2000 cs : (3868.032782534150, -191.721482350998, 5050.752208421422, 0.013973142331, 0.281375212237, -0.000020373168)km at time tR = 28662.7087619720104166697604
. State of Luna's origin SSBMJ2000 cs : (15823062.953505944461, -137497715.996215969324, -59644480.069209396839, 28.474022222022, 3.416219609543, 1.521426940172)km at time tR = 28662.7087619720104166697604
. State of Earth's origin SSBMJ2000 cs : (15515968.399858726189, -137730327.214818507433, -59708640.580756992102, 29.134741651851, 2.737218180597, 1.186422830232)km at time tR = 28662.7087619720104166697604
. DSLWP_B position in SSBMJ2000 : (15831157.835126012564, -137504841.728244185448, -59654289.355959296227)km at time tT = 28662.7087619720104166697604
. PI9CAM position in SSBMJ2000 : (15519836.432641260326, -137730518.936300843954, -59703589.828548572958)km at time tR = 28662.7087619720104166697604
. Range vector in SSBMJ2000: (-311321.402484753053, -225677.208056669100, -49300.472589276054)km
. Range vector in obs CS : (-208573.189838595514, -240837.852848393901, -220853.135804727150)km
. R_Obs_j2k matrix : (0.794549659492 -0.039504839067 -0.605912705174)
(0.049601818299 0.998769069445 -0.000074438963)
(0.605169809411 -0.029995226454 0.795531135888)
DeltaT for light travel over distance 387661.907 km = -1.293101e+00
Starting: dEpoch = 0.000000000000e+00 second, dR = 387661.907 km, dT = -1.293100933975e+00 second
Initial x Positions: DSLWP_B 8094.882 --> PI9CAM 3868.033
Loop iteration 0
Positions in SSBMJ2000 cs: DSLWP_B (15831120.853950809687 -137504846.540667116642 -59654291.294330336154)km --> PI9CAM (15519836.432641260326 -137730518.936300843954 -59703589.828548572958)km
Light Time range = 387629.161154257366 km
Relativity correction = 0.007573161548 km
===> dEpoch = -1.293100933975e+00 second, dR = 387629.16872742 km, dT = -1.292991729390e+00 second, trigger = -1.092046e-04 second
Loop iteration 1
Positions in SSBMJ2000 cs: DSLWP_B (15831120.857073932886 -137504846.540260702372 -59654291.294166646898)km --> PI9CAM (15519836.432641260326 -137730518.936300843954 -59703589.828548572958)km
Light Time range = 387629.163919698738 km
Relativity correction = 0.007573161601 km
===> dEpoch = -1.292991729390e+00 second, dR = 387629.17149286 km, dT = -1.292991738614e+00 second, trigger = 9.224520e-09 second
Loop iteration 2
Positions in SSBMJ2000 cs: DSLWP_B (15831120.857073670253 -137504846.540260732174 -59654291.294166654348)km --> PI9CAM (15519836.432641260326 -137730518.936300843954 -59703589.828548572958)km
Light Time range = 387629.163919469516 km
Relativity correction = 0.007573161601 km
===> dEpoch = -1.292991738614e+00 second, dR = 387629.17149263 km, dT = -1.292991738613e+00 second, trigger = -7.644996e-13 second
2. Compute Range Vector after light time correction for the Leg from <TNode = DSLWP_B> to <RNode = PI9CAM>:
. DSLWP_B state in LunaMJ2000 cs : (8094.720245441031, -7126.126904128575, -9809.257766336139, 0.124815469025, 0.305389759390, -0.022425969778)km at time tT= 28662.7087470068282567825742
. PI9CAM state in EarthMJ2000 cs : (3868.032782534150, -191.721482350998, 5050.752208421422, 0.013973142331, 0.281375212237, -0.000020373168)km at time tR = 28662.7087619720104166697604
. State of Luna's origin SSBMJ2000 cs : (15823026.136828228831, -137497720.413356602192, -59644482.036400318146, 28.474025670240, 3.416214836396, 1.521424551933)km at time tR = 28662.7087470068282567825742
. State of Earth's origin SSBMJ2000 cs : (15515968.399858726189, -137730327.214818507433, -59708640.580756992102, 29.134741651851, 2.737218180597, 1.186422830232)km at time tR = 28662.7087619720104166697604
. DSLWP_B position in SSBMJ2000 : (15831120.857073670253, -137504846.540260732174, -59654291.294166654348)km at time tT = 28662.7087470068282567825742
. PI9CAM position in SSBMJ2000 : (15519836.432641260326, -137730518.936300843954, -59703589.828548572958)km at time tR = 28662.7087619720104166697604
. Range vector in SSBMJ2000: (-311284.424432409520, -225672.396040127642, -49298.534381916397)km
. Range vector in obs CS : (-208545.173422100837, -240831.212720755342, -220829.360237061454)km
. R_Obs_j2k matrix : (0.794549659492 -0.039504839067 -0.605912705174)
(0.049601818299 0.998769069445 -0.000074438963)
(0.605169809411 -0.029995226454 0.795531135888)
3. Summary for signal leg from DSLWP_B to PI9CAM:
. Geometric range = 387661.907438495546 km
. Light time solution range = 387629.163919478131 km
. Relativity correction = 0.007573161601 km
. ET-TAI correction = 32.978179363250 km
. Feasibility = true
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Signal Frequency calculation for leg from DSLWP_B to PI9CAM :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
. Arrival frequency : -1.000000000000e+00 Mhz
. Transmit frequency : 4.354000000000e+02 Mhz
. Doppler shift frequency: 4.353999891583e+02 Mhz
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Media corrections calculation for leg from DSLWP_B to PI9CAM :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
.Frequency : 4.353999891583e+02 Mhz
.Troposphere range correction : 0.004284564790 m
.Troposphere elevation correction : 0.000471377966 rad
===================================================================
==== GNRangeAdapter (DSNsimData_{DSLWP_B,PI9CAM}_Range): Range Calculation for Measurement Data 0th
===================================================================
. Path : DSLWP_B, PI9CAM,
. Measurement epoch : 28662.708761972011
. Measurement type : <DeltaRange>
. C-value w/o noise and bias : 387662.153956567810 km
. Noise adding option : true
. Bias adding option : true
. C-value with noise and bias : 387662.153956567810 km
. Measurement epoch A1Mjd : 28662.708761972011
. Transmit frequency at receive epoch : 4.354000000000e+08 Hz
. Transmit frequency at transmit epoch : 4.354000000000e+08 Hz
. Measurement is feasible
. Feasibility reason : N
. Elevation angle : 34.728773536118 degree
. Covariance matrix : <0x55bf39780460>
. Covariance matrix size = 1
[ 1.000000e+00 ]
===================================================================
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Hardware delay calculation for leg from DSLWP_B to Shahe :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
. DSLWP_B's hardware delay : 0.000000000000e+00 sec
. Shahe's hardware delay : 0.000000000000e+00 sec
. Total hardware delay for this signal leg: 0.000000000000e+00 sec
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Range, relativity correction, and ET-TAI correction calculation for leg from DSLWP_B to Shahe at fixed time 28662.7087470068282567825742 A1Mjd:
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1. Compute Range Vector before light time correction for the Leg from <TNode = DSLWP_B> to <RNode = Shahe>:
. DSLWP_B state in LunaMJ2000 cs : (8094.720245441029, -7126.126904128576, -9809.257766336139, 0.124815469025, 0.305389759390, -0.022425969778)km at time tT= 28662.7087470068282567825742
. Shahe state in EarthMJ2000 cs : (-1419.229842423156, 4671.222940501376, 4090.741087709160, -0.340636346202, -0.104047064570, 0.000632102328)km at time tR = 28662.7087470068282567825742
. State of Luna's origin SSBMJ2000 cs : (15823026.136828228831, -137497720.413356602192, -59644482.036400318146, 28.474025670240, 3.416214836396, 1.521424551933)km at time tR = 28662.7087470068282567825742
. State of Earth's origin SSBMJ2000 cs : (15515930.728877997026, -137730330.754014611244, -59708642.114790007472, 29.134742392363, 2.737211386347, 1.186419888852)km at time tR = 28662.7087470068282567825742
. DSLWP_B position in SSBMJ2000 : (15831120.857073670253, -137504846.540260732174, -59654291.294166654348)km at time tT = 28662.7087470068282567825742
. Shahe position in SSBMJ2000 : (15514511.499035574496, -137725659.531074106693, -59704551.373702295125)km at time tR = 28662.7087470068282567825742
. Range vector in SSBMJ2000: (-316609.358038095990, -220812.990813379089, -50260.079535644029)km
. Range vector in obs CS : (-37615.405169619153, 367214.186377098318, -123550.640311999639)km
. R_Obs_j2k matrix : (-0.189656098441 0.616264375063 -0.764361684252)
(-0.956378457367 -0.292125137387 0.001774934109)
(-0.222195433366 0.731355675533 0.644785286167)
DeltaT for light travel over distance 389263.327 km = 1.298443e+00
Starting: dEpoch = 0.000000000000e+00 second, dR = 389263.327 km, dT = 1.298442696176e+00 second
Initial x Positions: DSLWP_B 8094.720 --> Shahe -1419.230
Loop iteration 0
Positions in SSBMJ2000 cs: DSLWP_B (15831120.857073670253 -137504846.540260732174 -59654291.294166654348)km --> Shahe (15514548.886538181454 -137725656.112077653408 -59704549.832381397486)km
Light Time range = 389230.780126877013 km
Relativity correction = 0.007605963202 km
===> dEpoch = 1.298442696176e+00 second, dR = 389230.78773284 km, dT = 1.298334155334e+00 second, trigger = 1.085408e-04 second
Loop iteration 1
Positions in SSBMJ2000 cs: DSLWP_B (15831120.857073670253 -137504846.540260732174 -59654291.294166654348)km --> Shahe (15514548.883412841707 -137725656.112363427877 -59704549.832510247827)km
Light Time range = 389230.782847557333 km
Relativity correction = 0.007605963255 km
===> dEpoch = 1.298334155334e+00 second, dR = 389230.79045352 km, dT = 1.298334164409e+00 second, trigger = -9.075213e-09 second
Loop iteration 2
Positions in SSBMJ2000 cs: DSLWP_B (15831120.857073670253 -137504846.540260732174 -59654291.294166654348)km --> Shahe (15514548.883413106203 -137725656.112363398075 -59704549.832510240376)km
Light Time range = 389230.782847324386 km
Relativity correction = 0.007605963255 km
===> dEpoch = 1.298334164409e+00 second, dR = 389230.79045329 km, dT = 1.298334164408e+00 second, trigger = 7.769341e-13 second
2. Compute Range Vector after light time correction for the Leg from <TNode = DSLWP_B> to <RNode = Shahe>:
. DSLWP_B state in LunaMJ2000 cs : (8094.720245441029, -7126.126904128576, -9809.257766336139, 0.124815469025, 0.305389759390, -0.022425969778)km at time tT= 28662.7087470068282567825742
. Shahe state in EarthMJ2000 cs : (-1419.672095818080, 4671.087831707016, 4090.741908382495, -0.340626493932, -0.104079314253, 0.000632083341)km at time tR = 28662.7087620338440485522649
. State of Luna's origin SSBMJ2000 cs : (15823026.136828228831, -137497720.413356602192, -59644482.036400318146, 28.474025670240, 3.416214836396, 1.521424551933)km at time tR = 28662.7087470068282567825742
. State of Earth's origin SSBMJ2000 cs : (15515968.555508924648, -137730327.200195103884, -59708640.574418619275, 29.134741648791, 2.737218208669, 1.186422842386)km at time tR = 28662.7087620338440485522649
. DSLWP_B position in SSBMJ2000 : (15831120.857073670253, -137504846.540260732174, -59654291.294166654348)km at time tT = 28662.7087470068282567825742
. Shahe position in SSBMJ2000 : (15514548.883413106203, -137725656.112363398075, -59704549.832510240376)km at time tR = 28662.7087620338440485522649
. Range vector in SSBMJ2000: (-316571.973660563293, -220809.572102666105, -50258.538343582499)km
. Range vector in obs CS : (-37599.166523781809, 367188.675828624808, -123528.868431141673)km
. R_Obs_j2k matrix : (-0.189714441499 0.616246551181 -0.764361575986)
(-0.956350795890 -0.292215682327 0.001774880798)
(-0.222264675301 0.731334522062 0.644785414656)
3. Summary for signal leg from DSLWP_B to Shahe:
. Geometric range = 389263.327458642016 km
. Light time solution range = 389230.782847323397 km
. Relativity correction = 0.007605963255 km
. ET-TAI correction = 33.451455045338 km
. Feasibility = true
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Signal Frequency calculation for leg from DSLWP_B to Shahe :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
. Arrival frequency : -1.000000000000e+00 Mhz
. Transmit frequency : 4.354000000000e+02 Mhz
. Doppler shift frequency: 4.353992753690e+02 Mhz
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
++++ Media corrections calculation for leg from DSLWP_B to Shahe :
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
.Frequency : 4.353992753690e+02 Mhz
.Troposphere range correction : 0.007635640939 m
.Troposphere elevation correction : 0.000968731245 rad
===================================================================
==== GNRangeAdapter (DSNsimData_{DSLWP_B,Shahe}_Range): Range Calculation for Measurement Data 0th
===================================================================
. Path : DSLWP_B, Shahe,
. Measurement epoch : 28662.708747006829
. Measurement type : <DeltaRange>
. C-value w/o noise and bias : 389264.249543972895 km
. Noise adding option : true
. Bias adding option : true
. C-value with noise and bias : 389264.249543972895 km
. Measurement epoch A1Mjd : 28662.708747006829
. Transmit frequency at receive epoch : 4.354000000000e+08 Hz
. Transmit frequency at transmit epoch : 4.354000000000e+08 Hz
. Measurement is feasible
. Feasibility reason : N
. Elevation angle : 18.503745381468 degree
. Covariance matrix : <0x55bf39781f30>
. Covariance matrix size = 1
[ 1.000000e+00 ]
===================================================================
===================================================================
==== DeltaRangeAdapter (DSNsimData_{PI9CAM,DSLWP_B,Shahe}_DeltaRange): DeltaRange Calculation
===================================================================
. Path : [ DSLWP_B -> PI9CAM ] - [ DSLWP_B -> Shahe ]
. Measurement epoch : 28662.708761972011
. Measurement type : <DeltaRange>
. Reference leg C-value w/o noise and bias : 387662.153956567810 km
. Other leg C-value w/o noise and bias : 389264.249543972895 km
. C-value w/o noise and bias : -1602.095587405085 km
. Reference leg corrections : 32.990037089641 km
. Other leg corrections : 33.466696649532 km
. Corrections : -0.476659559891 km
. Noise adding option : true
. Bias adding option : true
. Range noise sigma : 1.100000000000 km
. Range bias : -3.000000000000 km
. Multiplier : 0.000000000000
. C-value with noise and bias : -1605.728116406530 km
. Measurement epoch A1Mjd : 28662.708761972011
. Transmit frequency at receive epoch : 0.000000000000e+00 Hz
. Transmit frequency at transmit epoch : 0.000000000000e+00 Hz
. Measurement is feasible
. Feasibility reason : N
. Elevation angle : 18.503745381468 degree
. Covariance matrix : <0x55bf3976b0d0>
. Covariance matrix size = 1
[ 1.210000e+00 ]
===================================================================
Mission run completed.
===> Total Run Time: 0.785 seconds
========================================
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment