Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
def galpy_nfw_orbit():
# Setting up the potential
nfw = NFWPotential(conc=C, mvir=M_200, H=70.0, wrtcrit=True, overdens=200)
nfw.turn_physical_on()
vxvv = [
8.0 * units.kpc,
0.0 * units.km / units.s,
240.0 * units.km / units.s,
0.0 * units.pc,
5.0 * units.km / units.s,
]
# Calculating the orbit
ts = np.linspace(0.0, 0.58, 1000) * units.Gyr
o = Orbit(vxvv=vxvv)
o.integrate(ts, nfw, method="odeint")
return o