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

Commit

Permalink
Fixed bug in example problem where escape velocity was miscalculated …
Browse files Browse the repository at this point in the history
…when setting test particle initial conditions
  • Loading branch information
daminton committed May 20, 2023
1 parent a2c5e54 commit b6e61d3
Showing 1 changed file with 1 addition and 1 deletion.
Original file line number Diff line number Diff line change
Expand Up @@ -17,7 +17,7 @@
for i,n in enumerate(pl.name):
pli = pl.sel(name=n)
rstart = 2 * pli['radius'].data[0] # Start the test particles at a multiple of the planet radius away
vstart = 1.5 * np.sqrt(2 * pli['Gmass'].data[0]) / rstart # Start the test particle velocities at a multiple of the escape speed
vstart = 1.5 * np.sqrt(2 * pli['Gmass'].data[0] / rstart) # Start the test particle velocities at a multiple of the escape speed
rstart_vec = np.array([rstart / np.sqrt(2.0), rstart / np.sqrt(2.0), 0.0])
vstart_vec = np.array([vstart, 0.0, 0.0])
rp = pli['rh'].data[0]
Expand Down

0 comments on commit b6e61d3

Please sign in to comment.