diff options
| author | tbode <tbode@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-04-28 17:55:57 +0000 |
|---|---|---|
| committer | tbode <tbode@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-04-28 17:55:57 +0000 |
| commit | e0baa8110196356ee96faa0788d950126fddd8ee (patch) | |
| tree | b1d6b50cc3728ab696edcf304f767d92e88ee9a9 /src/GRHydro_Boundaries.F90 | |
| parent | 4adfcb611a12689096a604e2a510b7ff12c145f1 (diff) | |
MERGE divergence cleaning feature into trunk.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@244 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Boundaries.F90')
| -rw-r--r-- | src/GRHydro_Boundaries.F90 | 25 |
1 files changed, 20 insertions, 5 deletions
diff --git a/src/GRHydro_Boundaries.F90 b/src/GRHydro_Boundaries.F90 index 9e95adc..765755b 100644 --- a/src/GRHydro_Boundaries.F90 +++ b/src/GRHydro_Boundaries.F90 @@ -24,6 +24,9 @@ #define Bvecx(i,j,k) Bvec(i,j,k,1) #define Bvecy(i,j,k) Bvec(i,j,k,2) #define Bvecz(i,j,k) Bvec(i,j,k,3) +#define Bconsx(i,j,k) Bcons(i,j,k,1) +#define Bconsy(i,j,k) Bcons(i,j,k,2) +#define Bconsz(i,j,k) Bcons(i,j,k,3) /*@@ @routine GRHydro_InitSymBound @@ -93,8 +96,10 @@ subroutine GRHydro_InitSymBound(CCTK_ARGUMENTS) call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[0]") - if(evolve_mhd.ne.0)call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[0]") - + if(evolve_mhd.ne.0) then + call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[0]") + call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::Bcons[0]") + endif sym(1) = 1 sym(2) = -1 @@ -102,15 +107,21 @@ subroutine GRHydro_InitSymBound(CCTK_ARGUMENTS) call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[1]") - if(evolve_mhd.ne.0)call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]") - + if(evolve_mhd.ne.0) then + call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]") + call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::Bcons[1]") + endif + sym(1) = 1 sym(2) = 1 sym(3) = -1 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]") call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[2]") - if(evolve_mhd.ne.0)call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[2]") + if(evolve_mhd.ne.0) then + call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[2]") + call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::Bcons[2]") + endif end subroutine GRHydro_InitSymBound @@ -169,6 +180,8 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS) if(evolve_mhd.ne.0) then ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::Bvec", "Flat") + ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + "GRHydro::Bcons", "Flat") if(clean_divergence.ne.0) then ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::psidc", "Flat") @@ -219,6 +232,8 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS) if(evolve_mhd.ne.0) then ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::Bvec", "None") + ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + "GRHydro::Bcons", "None") if(clean_divergence.ne.0) then ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::psidc", "None") |