Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
fixed typos in init_cond.py and added XVEL case to headers
Browse files Browse the repository at this point in the history
  • Loading branch information
cwishard committed Aug 26, 2021
1 parent 60eb897 commit 1e5fdff
Showing 1 changed file with 28 additions and 2 deletions.
30 changes: 28 additions & 2 deletions python/swiftest/swiftest/init_cond.py
Original file line number Diff line number Diff line change
Expand Up @@ -129,7 +129,7 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
ephemerides_step = '1d'

clab, plab, tlab = swiftest.io.make_swiftest_labels(param)
if param['OUT_FORM'] == 'XV':
if param['OUT_FORM'] == 'EL':
plab.append('a')
plab.append('e')
plab.append('inc')
Expand All @@ -142,7 +142,7 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
tlab.append('capom')
tlab.append('omega')
tlab.append('capm')
elif param['OUT_FORM'] == 'EL':
elif param['OUT_FORM'] == 'XV':
plab.append('xhx')
plab.append('xhy')
plab.append('xhz')
Expand All @@ -156,6 +156,32 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
tlab.append('vhy')
tlab.append('vhz')

elif param['OUT_FORM'] == 'XVEL':
plab.append('xhx')
plab.append('xhy')
plab.append('xhz')
plab.append('vhx')
plab.append('vhy')
plab.append('vhz')
tlab.append('xhx')
tlab.append('xhy')
tlab.append('xhz')
tlab.append('vhx')
tlab.append('vhy')
tlab.append('vhz')
plab.append('a')
plab.append('e')
plab.append('inc')
plab.append('capom')
plab.append('omega')
plab.append('capm')
tlab.append('a')
tlab.append('e')
tlab.append('inc')
tlab.append('capom')
tlab.append('omega')
tlab.append('capm')

dims = ['time', 'id', 'vec']
t = np.array([0.0])

Expand Down

0 comments on commit 1e5fdff

Please sign in to comment.