diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2011-04-30 00:22:45 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2011-04-30 00:22:45 +0000 |
| commit | 6cc99a9de0769e25654a1413e73157bb750aa278 (patch) | |
| tree | 0a952788cfd9b7112b228d854048d0dc5a4705bb | |
| parent | d01d6ce889030ae0322409021d5c004b3fbcd649 (diff) | |
GRHydro_InitialData: apply boundary conditions to Bcons for diagonal shocks
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@128 ac85fae7-cede-4708-beff-ae01c7fa1c26
| -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 |