diff options
Diffstat (limited to 'src/GRHydro_Con2Prim.F90')
-rw-r--r-- | src/GRHydro_Con2Prim.F90 | 26 |
1 files changed, 14 insertions, 12 deletions
diff --git a/src/GRHydro_Con2Prim.F90 b/src/GRHydro_Con2Prim.F90 index 550e2bd..e159a50 100644 --- a/src/GRHydro_Con2Prim.F90 +++ b/src/GRHydro_Con2Prim.F90 @@ -1026,6 +1026,8 @@ subroutine Conservative2PrimitiveBounds(CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS + CCTK_REAL, parameter :: half = 0.5d0 + integer :: i, j, k, itracer, nx, ny, nz CCTK_REAL :: uxxl, uxyl, uxzl, uyyl, uyzl, uzzl,& uxxr, uxyr, uxzr, uyyr, uyzr, uzzr, pmin, epsmin @@ -1135,9 +1137,9 @@ subroutine Conservative2PrimitiveBounds(CCTK_ARGUMENTS) velyminus(i,j,k),velzminus(i,j,k),epsminus(i,j,k),& pressminus(i,j,k),w_lorentzminus(i,j,k),& uxxl,uxyl,uxzl,uyyl,uyzl,uzzl,avg_detl,& - x(i,j,k)-0.5d0*CCTK_DELTA_SPACE(1),& - y(i,j,k)-0.5d0*CCTK_DELTA_SPACE(2), & - z(i,j,k)-0.5d0*CCTK_DELTA_SPACE(3),r(i,j,k),& + x(i,j,k)-half*CCTK_DELTA_SPACE(1),& + y(i,j,k)-half*CCTK_DELTA_SPACE(2), & + z(i,j,k)-half*CCTK_DELTA_SPACE(3),r(i,j,k),& epsnegative,GRHydro_rho_min,pmin, epsmin, GRHydro_reflevel, GRHydro_C2P_failed(i,j,k)) if (epsnegative) then !$OMP CRITICAL @@ -1149,9 +1151,9 @@ subroutine Conservative2PrimitiveBounds(CCTK_ARGUMENTS) velyminus(i,j,k),velzminus(i,j,k),epsminus(i,j,k),& pressminus(i,j,k),w_lorentzminus(i,j,k),& uxxl,uxyl,uxzl,uyyl,uyzl,uzzl,avg_detl,& - x(i,j,k)-0.5d0*CCTK_DELTA_SPACE(1),& - y(i,j,k)-0.5d0*CCTK_DELTA_SPACE(2), & - z(i,j,k)-0.5d0*CCTK_DELTA_SPACE(3),r(i,j,k),GRHydro_rho_min, GRHydro_reflevel, GRHydro_C2P_failed(i,j,k)) + x(i,j,k)-half*CCTK_DELTA_SPACE(1),& + y(i,j,k)-half*CCTK_DELTA_SPACE(2), & + z(i,j,k)-half*CCTK_DELTA_SPACE(3),r(i,j,k),GRHydro_rho_min, GRHydro_reflevel, GRHydro_C2P_failed(i,j,k)) end if if (epsminus(i,j,k) .lt. 0.0d0) then @@ -1182,9 +1184,9 @@ subroutine Conservative2PrimitiveBounds(CCTK_ARGUMENTS) velyplus(i,j,k),velzplus(i,j,k),epsplus(i,j,k),& pressplus(i,j,k),w_lorentzplus(i,j,k),& uxxr,uxyr,uxzr,uyyr,uyzr,uzzr,avg_detr,& - x(i,j,k)+0.5d0*CCTK_DELTA_SPACE(1),& - y(i,j,k)+0.5d0*CCTK_DELTA_SPACE(2), & - z(i,j,k)+0.5d0*CCTK_DELTA_SPACE(3),r(i,j,k),& + x(i,j,k)+half*CCTK_DELTA_SPACE(1),& + y(i,j,k)+half*CCTK_DELTA_SPACE(2), & + z(i,j,k)+half*CCTK_DELTA_SPACE(3),r(i,j,k),& epsnegative,GRHydro_rho_min,pmin, epsmin, GRHydro_reflevel,GRHydro_C2P_failed(i,j,k)) if (epsnegative) then !$OMP CRITICAL @@ -1196,9 +1198,9 @@ subroutine Conservative2PrimitiveBounds(CCTK_ARGUMENTS) velyplus(i,j,k),velzplus(i,j,k),epsplus(i,j,k),& pressplus(i,j,k),w_lorentzplus(i,j,k),& uxxr,uxyr,uxzr,uyyr,uyzr,uzzr,avg_detr,& - x(i,j,k)+0.5d0*CCTK_DELTA_SPACE(1),& - y(i,j,k)+0.5d0*CCTK_DELTA_SPACE(2), & - z(i,j,k)+0.5d0*CCTK_DELTA_SPACE(3),r(i,j,k),GRHydro_rho_min, GRHydro_reflevel, GRHydro_C2P_failed(i,j,k)) + x(i,j,k)+half*CCTK_DELTA_SPACE(1),& + y(i,j,k)+half*CCTK_DELTA_SPACE(2), & + z(i,j,k)+half*CCTK_DELTA_SPACE(3),r(i,j,k),GRHydro_rho_min, GRHydro_reflevel, GRHydro_C2P_failed(i,j,k)) end if if (epsplus(i,j,k) .lt. 0.0d0) then |