diff options
Diffstat (limited to 'src/GRHydro_P2C2PM.F90')
| -rw-r--r-- | src/GRHydro_P2C2PM.F90 | 20 |
1 files changed, 15 insertions, 5 deletions
diff --git a/src/GRHydro_P2C2PM.F90 b/src/GRHydro_P2C2PM.F90 index 62c998b..15cacf0 100644 --- a/src/GRHydro_P2C2PM.F90 +++ b/src/GRHydro_P2C2PM.F90 @@ -40,6 +40,7 @@ subroutine p2c2pm(CCTK_ARGUMENTS) CCTK_REAL :: uxx,uxy,uxz,uyy,uyz,uzz CCTK_REAL :: gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send CCTK_REAL :: dens_send,sx_send,sy_send,sz_send,tau_send + CCTK_REAL :: bconsx_send,bconsy_send,bconsz_send CCTK_REAL :: rho_send(1),velx_send,vely_send,velz_send,eps_send(1) CCTK_REAL :: press_send(1),w_lorentz_send,x_send,y_send,z_send,r_send CCTK_REAL :: bvcx_send,bvcy_send,bvcz_send,b2_send @@ -86,6 +87,9 @@ subroutine p2c2pm(CCTK_ARGUMENTS) bvcx_send = Bx_init bvcy_send = By_init bvcz_send = Bz_init + bconsx_send = Bx_init + bconsy_send = By_init + bconsz_send = Bz_init w_lorentz_send = 1.d0/sqrt(1.0d0-velx_send*velx_send-vely_send*vely_send-velz_send*velz_send) @@ -127,8 +131,9 @@ subroutine p2c2pm(CCTK_ARGUMENTS) call Prim2ConGenM(GRHydro_eos_handle,& gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send,det, & dens_send,sx_send,sy_send,sz_send,tau_send, & - bvcx_send,bvcy_send,bvcz_send,rho_send(1), & - velx_send,vely_send,velz_send,eps_send(1),press_send(1),w_lorentz_send) + bconsx_send,bcibsy_send,bconsz_send, & + rho_send(1),velx_send,vely_send,velz_send,eps_send(1),press_send(1), & + bvcx_send,bvcy_send,bvcz_send,w_lorentz_send) write(*,*) 'P2C2PM test: initial values.' write(*,*) ' conservative variables: ' @@ -137,6 +142,9 @@ subroutine p2c2pm(CCTK_ARGUMENTS) write(*,*) ' sy : ',sy_send write(*,*) ' sz : ',sz_send write(*,*) ' tau : ',tau_send + write(*,*) ' Bconsx : ',bconsx_send + write(*,*) ' Bconsy : ',bconsy_send + write(*,*) ' Bconsz : ',bconsz_send write(*,*) ' eps : ',eps_send write(*,*) ' W : ',w_lorentz_send write(*,*) ' Bvecx : ',bvcx_send @@ -146,11 +154,13 @@ subroutine p2c2pm(CCTK_ARGUMENTS) write(*,*) 'P2C2PM test: getting the associated primitive variables.' call Con2PrimGenM(GRHydro_eos_handle,local_gam(1),dens_send,sx_send,sy_send,sz_send, & - tau_send,rho_send(1),velx_send,vely_send,velz_send, & - eps_send(1),press_send(1),w_lorentz_send, & + tau_send,bconsx_send,bconsy_send,bconsz_send,& + rho_send(1),velx_send,vely_send,velz_send, & + eps_send(1),press_send(1), & + bvcx_send,bvcy_send,bvcz_send,b2_send,& + w_lorentz_send, & gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send,& uxx,uxy,uxz,uyy,uyz,uzz,det,& - bvcx_send,bvcy_send,bvcz_send,b2_send,& epsnegative,C2P_failed) write(*,*) 'P2C2PM test: the primitive variables are' |