Update inject.f90

When injecting particles at rinj with Keplerian velocity, we do not take into account the velocity and position of the central object
danieljprice · Jun 3, 2024 · 9fa15bf · 9fa15bf
1 parent 11749fe
commit 9fa15bf
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/main/inject_keplerian.f90 b/src/main/inject_keplerian.f90
@@ -171,8 +171,8 @@ subroutine inject_particles(time,dtlast,xyzh,vxyzu,xyzmh_ptmass,vxyz_ptmass,&
 
        vphi = vkep*(1. - (zi/rinj)**2)**(-0.75)  ! see Martire et al. (2024)
 
-       xyzi = (/rinj*cosphi,rinj*sinphi,zi/) + x0
-       vxyz = (/-vphi*sinphi, vphi*cosphi, 0./) + v0
+       xyzi = (/rinj*cosphi,rinj*sinphi,zi/) 
+       vxyz = (/-vphi*sinphi, vphi*cosphi, 0./) 
 
        u = 1.5*cs**2