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

Commit

Permalink
changed the names of velocity labels to match horizons database
Browse files Browse the repository at this point in the history
  • Loading branch information
cwishard committed Aug 26, 2021
1 parent 6b7d3ba commit 4e1e44b
Showing 1 changed file with 7 additions and 6 deletions.
13 changes: 7 additions & 6 deletions python/swiftest/swiftest/init_cond.py
Original file line number Diff line number Diff line change
Expand Up @@ -210,9 +210,9 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
p1.append(pldata[key].vectors()['x'][0] * DCONV)
p2.append(pldata[key].vectors()['y'][0] * DCONV)
p3.append(pldata[key].vectors()['z'][0] * DCONV)
p4.append(pldata[key].vectors()['vhx'][0] * VCONV)
p5.append(pldata[key].vectors()['vhy'][0] * VCONV)
p6.append(pldata[key].vectors()['vhz'][0] * VCONV)
p4.append(pldata[key].vectors()['vx'][0] * VCONV)
p5.append(pldata[key].vectors()['vy'][0] * VCONV)
p6.append(pldata[key].vectors()['vz'][0] * VCONV)
p7.append(pldata[key].elements()['a'][0] * DCONV)
p8.append(pldata[key].elements()['e'][0])
p9.append(pldata[key].elements()['incl'][0])
Expand All @@ -229,10 +229,11 @@ def solar_system_horizons(plname, idval, param, ephemerides_start_date, ds):
p7.append(pldata[key].vectors()['x'][0] * DCONV)
p8.append(pldata[key].vectors()['y'][0] * DCONV)
p9.append(pldata[key].vectors()['z'][0] * DCONV)
p10.append(pldata[key].vectors()['vhx'][0] * VCONV)
p11.append(pldata[key].vectors()['vhy'][0] * VCONV)
p12.append(pldata[key].vectors()['vhz'][0] * VCONV)
p10.append(pldata[key].vectors()['vx'][0] * VCONV)
p11.append(pldata[key].vectors()['vy'][0] * VCONV)
p12.append(pldata[key].vectors()['vz'][0] * VCONV)
pvec = np.vstack([p1, p2, p3, p4, p5, p6])

if ispl:
Rpl.append(planetradius[key] * DCONV)
GMpl.append(GMcb[0] / MSun_over_Mpl[key])
Expand Down

0 comments on commit 4e1e44b

Please sign in to comment.