diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-13 19:08:58 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-13 19:08:58 +0000 |
| commit | 5ee96de77ed204c78f07cff1da75bdb125f855f0 (patch) | |
| tree | dac5ed2b91d19d011274ea2920e38a844f70fdac /src/GRHydro_Con2Prim.F90 | |
| parent | ccebc7be6a5189a0dd2e45649c2405cec343b4e5 (diff) | |
Error check for C2P_failed mask greater or equal 1, do not abort code if in the
atmosphere though.
patch by Philipp Moesta
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@322 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Con2Prim.F90')
| -rw-r--r-- | src/GRHydro_Con2Prim.F90 | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/src/GRHydro_Con2Prim.F90 b/src/GRHydro_Con2Prim.F90 index e159a50..6ac8cd2 100644 --- a/src/GRHydro_Con2Prim.F90 +++ b/src/GRHydro_Con2Prim.F90 @@ -1980,7 +1980,7 @@ subroutine check_GRHydro_C2P_failed(CCTK_ARGUMENTS) do j = 1, ny do i = 1, nx - if (GRHydro_C2P_failed(i,j,k) == 1) then + if (GRHydro_C2P_failed(i,j,k) .ge. 1) then if(con2prim_oct_hack.ne.0.and.& (x(i,j,k) .lt. -1.0d-12 .or.& @@ -1989,6 +1989,9 @@ subroutine check_GRHydro_C2P_failed(CCTK_ARGUMENTS) cycle endif + if ( rho(i,j,k) .le. GRHydro_rho_min*(1.d0+GRHydro_atmo_tolerance) ) then + cycle + end if !$OMP CRITICAL call CCTK_WARN(1,'Con2Prim failed; stopping the code.') |