diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-15 16:44:26 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-15 16:44:26 +0000 |
| commit | c1d720daaf5b837471db77629803bd7aa5c92e8c (patch) | |
| tree | 5085634cc564c730e961574d58d7a07397c26511 /src/GRHydro_UpdateMask.F90 | |
| parent | 76ba623bb876010bae337587f20a440531956fe0 (diff) | |
include Y_e_con in atmosphere handling
original commit by Christian Ott (cott)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@268 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_UpdateMask.F90')
| -rw-r--r-- | src/GRHydro_UpdateMask.F90 | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90 index d408533..d6bf540 100644 --- a/src/GRHydro_UpdateMask.F90 +++ b/src/GRHydro_UpdateMask.F90 @@ -41,7 +41,7 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS) frac = CCTK_DELTA_TIME - if(evolve_temper.ne.1) then + if(evolve_temper.ne.1.and.evolve_Y_e.ne.1) then !$OMP PARALLEL DO PRIVATE(k,j,i) do k = 1, cctk_lsh(3) do j = 1, cctk_lsh(2) @@ -72,6 +72,7 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS) (SpaceMask_CheckStateBitsF90(space_mask, i, j, k, \ type_bits, atmosphere)) .or. & (dens(i,j,k) + frac * densrhs(i,j,k) .le. 0.d0) ) then + y_e_con_rhs(i,j,k) = 0.0d0 densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 |