diff options
| -rw-r--r-- | src/GRHydro_ShockTubeM.F90 | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/src/GRHydro_ShockTubeM.F90 b/src/GRHydro_ShockTubeM.F90 index a87c528..a864054 100644 --- a/src/GRHydro_ShockTubeM.F90 +++ b/src/GRHydro_ShockTubeM.F90 @@ -708,9 +708,9 @@ subroutine GRHydro_Diagshock_BoundaryM(CCTK_ARGUMENTS) sy(i,j,k) = sy(inew,jnew,knew) sz(i,j,k) = sz(inew,jnew,knew) tau(i,j,k) = tau(inew,jnew,knew) - Bvecx(i,j,k)=Bvecx(inew,jnew,knew) - Bvecy(i,j,k)=Bvecy(inew,jnew,knew) - Bvecz(i,j,k)=Bvecz(inew,jnew,knew) + Bconsx(i,j,k)=Bconsx(inew,jnew,knew) + Bconsy(i,j,k)=Bconsy(inew,jnew,knew) + Bconsz(i,j,k)=Bconsz(inew,jnew,knew) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(inew,jnew,knew) endif @@ -789,9 +789,9 @@ subroutine GRHydro_Diagshock2D_BoundaryM(CCTK_ARGUMENTS) sy(i,j,k) = sy(inew,jnew,k) sz(i,j,k) = sz(inew,jnew,k) tau(i,j,k) = tau(inew,jnew,k) - Bvecx(i,j,k)=Bvecx(inew,jnew,k) - Bvecy(i,j,k)=Bvecy(inew,jnew,k) - Bvecz(i,j,k)=Bvecz(inew,jnew,k) + Bconsx(i,j,k)=Bconsx(inew,jnew,k) + Bconsy(i,j,k)=Bconsy(inew,jnew,k) + Bconsz(i,j,k)=Bconsz(inew,jnew,k) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(inew,jnew,k) endif @@ -805,9 +805,9 @@ subroutine GRHydro_Diagshock2D_BoundaryM(CCTK_ARGUMENTS) sy(i,j,k) = sy(i,j,kc) sz(i,j,k) = sz(i,j,kc) tau(i,j,k) = tau(i,j,kc) - Bvecx(i,j,k)=Bvecx(i,j,kc) - Bvecy(i,j,k)=Bvecy(i,j,kc) - Bvecz(i,j,k)=Bvecz(i,j,kc) + Bconsx(i,j,k)=Bconsx(i,j,kc) + Bconsy(i,j,k)=Bconsy(i,j,kc) + Bconsz(i,j,k)=Bconsz(i,j,kc) if(clean_divergence.ne.0) then psidc(i,j,k)=psidc(i,j,kc) endif |