diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:10:18 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-07-06 18:10:18 +0000 |
| commit | 5e766fded84f4e5e50f26864cfe628d9fb4c27be (patch) | |
| tree | 62c5da53ecf207a85821ce34c448e6d076fa3377 | |
| parent | a25126b7eea2576c4cae5075befae6d67e09b7c6 (diff) | |
GRHydro: Fixed a problem with excision mask and atmopshere treatment in hot Con2Prim.
From: Christian Reisswig <reisswig@scriwalker.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@543 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_Con2PrimHot.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/GRHydro_Con2PrimHot.F90 b/src/GRHydro_Con2PrimHot.F90 index a9e37f4..0d7015d 100644 --- a/src/GRHydro_Con2PrimHot.F90 +++ b/src/GRHydro_Con2PrimHot.F90 @@ -83,7 +83,7 @@ subroutine Conservative2PrimitiveHot(CCTK_ARGUMENTS) !$OMP PARALLEL DO PRIVATE(i,j,k,itracer,& !$OMP uxx, uxy, uxz, uyy, uyz, uzz, det, epsnegative, anyerr, keyerr, keytemp,& - !$OMP warnline, dummy1, dummy2) + !$OMP warnline, dummy1, dummy2,reset_to_atmo) do k = 1, nz do j = 1, ny do i = 1, nx @@ -125,7 +125,7 @@ subroutine Conservative2PrimitiveHot(CCTK_ARGUMENTS) reset_to_atmo = 1 endif - if (reset_to_atmo .gt. 0 .or. hydro_excision_mask(i,j,k) .gt. 0) then + if (reset_to_atmo .gt. 0 .or. (GRHydro_enable_internal_excision /= 0 .and. hydro_excision_mask(i,j,k) .gt. 0)) then SET_ATMO_MIN(dens(i,j,k), sqrt(det)*GRHydro_rho_min, r(i,j,k)) SET_ATMO_MIN(rho(i,j,k), GRHydro_rho_min, r(i,j,k)) scon(i,j,k,:) = 0.d0 |