From 5e766fded84f4e5e50f26864cfe628d9fb4c27be Mon Sep 17 00:00:00 2001 From: rhaas Date: Sat, 6 Jul 2013 18:10:18 +0000 Subject: GRHydro: Fixed a problem with excision mask and atmopshere treatment in hot Con2Prim. From: Christian Reisswig git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@543 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Con2PrimHot.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src') 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 -- cgit v1.2.3