Note
Click here to download the full example code
Profiling componentsΒΆ
Out:
SETTING DEBUG MODE
2020-10-15 13:21:27.736 DEBUG ; Propagator:_check_settings
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:epoch_format = mjd
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:epoch_scale = utc
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:time_format = sec
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:time_scale = None
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:heartbeat = False
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:in_frame = Orekit-ITRF
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:out_frame = Orekit-EME
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:frame_tidal_effects = False
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:integrator = DormandPrince853
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:min_step = 0.001
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:max_step = 120.0
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:position_tolerance = 10.0
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:earth_gravity = HolmesFeatherstone
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:gravity_order = (10, 10)
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:solarsystem_perturbers = ['Moon', 'Sun']
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:drag_force = False
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:atmosphere = DTM2000
2020-10-15 13:21:27.736 DEBUG ; Propagator:settings:radiation_pressure = False
2020-10-15 13:21:27.737 DEBUG ; Propagator:settings:solar_activity = Marshall
2020-10-15 13:21:27.737 DEBUG ; Propagator:settings:constants_source = WGS84
2020-10-15 13:21:27.737 DEBUG ; Propagator:settings:solar_activity_strength = WEAK
2020-10-15 13:21:27.737 DEBUG ; sorts.propagator.Orekit:init
2020-10-15 13:21:27.737 DEBUG ; Orekit:init:orekit-data = /home/danielk/IRF/IRF_GITLAB/orekit_build/orekit-data-master.zip
2020-10-15 13:21:27.899 DEBUG ; Orekit:init:_forces:earth_gravity = <class 'org.orekit.forces.gravity.HolmesFeatherstoneAttractionModel'>
2020-10-15 13:21:27.900 DEBUG ; Orekit:init:_forces:perturbation_Moon = <class 'org.orekit.forces.gravity.ThirdBodyAttraction'>
2020-10-15 13:21:27.900 DEBUG ; Orekit:init:_forces:perturbation_Sun = <class 'org.orekit.forces.gravity.ThirdBodyAttraction'>
2020-10-15 13:21:27.900 DEBUG ; Propagator:convert_time
2020-10-15 13:21:27.900 DEBUG ; Propagator:convert_time:completed
2020-10-15 13:21:27.905 DEBUG ; Orekit:propagate:len(t) = 5000
2020-10-15 13:21:27.906 DEBUG ; Orekit:set_forces:A = 1.0
2020-10-15 13:21:27.906 DEBUG ; Orekit:set_forces:cd = 1.0
2020-10-15 13:21:27.906 DEBUG ; Orekit:set_forces:cr = 1.0
2020-10-15 13:21:27.906 DEBUG ; Orekit:set_forces:re_calc = True
2020-10-15 13:21:28.460 DEBUG ; Orekit:propagate:completed
import numpy as np
import sorts
logger = sorts.profiling.get_logger('orekit')
orekit_data = '/home/danielk/IRF/IRF_GITLAB/orekit_build/orekit-data-master.zip'
def run_prop():
prop = sorts.propagator.Orekit(
orekit_data = orekit_data,
settings=dict(
in_frame='Orekit-ITRF',
out_frame='Orekit-EME',
drag_force = False,
radiation_pressure = False,
),
logger = logger,
)
state0 = np.array([-7100297.113,-3897715.442,18568433.707,86.771,-3407.231,2961.571])
t = np.linspace(0,3600*24.0*2,num=5000)
mjd0 = 53005
states = prop.propagate(t, state0, mjd0, A=1.0, C_R = 1.0, C_D = 1.0)
return states
states = run_prop()
print('\n SETTING DEBUG MODE \n')
sorts.profiling.term_level(logger, 'DEBUG')
states = run_prop()
Total running time of the script: ( 0 minutes 1.565 seconds)