diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:30 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:30 +0000 |
commit | d6e45e79d92f3d2adf3eb4ebcc7641db07119f30 (patch) | |
tree | bcfd7e5e200638f34bd270a438555950de5d34ed /src/GRHydro_ParamCheck.F90 | |
parent | c141ed811d4a5b84ecc077999d1fec0a3c598b95 (diff) |
GRHydro: optionally use CarpetEvolutionMask::evolution_mask to skip validity tests in regions that are not needed for evolution
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@454 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_ParamCheck.F90')
-rw-r--r-- | src/GRHydro_ParamCheck.F90 | 9 |
1 files changed, 8 insertions, 1 deletions
diff --git a/src/GRHydro_ParamCheck.F90 b/src/GRHydro_ParamCheck.F90 index b9b3997..085fac1 100644 --- a/src/GRHydro_ParamCheck.F90 +++ b/src/GRHydro_ParamCheck.F90 @@ -35,7 +35,7 @@ subroutine GRHydro_ParamCheck(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS - integer :: coordinates_is_active + integer :: coordinates_is_active, evolution_mask_idx if (GRHydro_stencil > minval(cctk_nghostzones)) then call CCTK_PARAMWARN("The stencil is larger than the number of ghost zones. Answer will be dependent on processor number...") @@ -138,6 +138,13 @@ subroutine GRHydro_ParamCheck(CCTK_ARGUMENTS) end if end if + if(use_evolution_mask.ne.0) then + call CCTK_VarIndex(evolution_mask_idx, "CarpetEvolutionMask::evolution_mask") + if(evolution_mask_idx .lt. 0) then + call CCTK_PARAMWARN("You activated use_evolution_mask but I cannot find 'CarpetEvolutionMask::evolution_mask'. If you use Carpet, then you should activate thorn 'CarpetEvolutionMask', if using PUGH then evolution_mask makes no sense and you should disable this option.") + end if + end if + if(clean_divergence.ne.0 .and. evolve_MHD.eq.0) then call CCTK_WARN(1, "You activated divergence cleaning but do not evolve a magnetic field. At best, nothing should happen, be preapered for a segfault otherwise.") end if |