diff options
author | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-03-28 00:59:39 +0000 |
---|---|---|
committer | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2010-03-28 00:59:39 +0000 |
commit | f909bf63ef761417261f5d729464625243ed61a2 (patch) | |
tree | 124d265811ea0751192e3565b9ac1fd834a24ffa | |
parent | c53825d6a37bec894be09536dff341a6a62c40fc (diff) |
* comment out fallback to polytrope in non-polytropic old-EOS Con2Prim and
rather set C2P fail mask to one where this happens. This change is necessary
to (a) catch problems that should not be present in real simulation data and
(b) prevent error messages from being generated in buffer zones that contain crap.
- Christian
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@90 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r-- | src/Whisky_Con2Prim.F90 | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/src/Whisky_Con2Prim.F90 b/src/Whisky_Con2Prim.F90 index 0fd930c..d001433 100644 --- a/src/Whisky_Con2Prim.F90 +++ b/src/Whisky_Con2Prim.F90 @@ -241,6 +241,12 @@ subroutine Conservative2Primitive(CCTK_ARGUMENTS) if (epsnegative) then + whisky_C2P_failed(i,j,k) = 1 + +#if 0 ! cott 2010/03/27: + ! Honestly, this should never happen. We need to flag the point where + ! this happened as having led to failing con2prim. + !$OMP CRITICAL call CCTK_WARN(1, 'Specific internal energy just went below 0, trying polytype.') !$OMP END CRITICAL @@ -250,7 +256,9 @@ subroutine Conservative2Primitive(CCTK_ARGUMENTS) vel(i,j,k,3), eps(i,j,k), press(i,j,k), w_lorentz(i,j,k), & uxx, uxy, uxz, uyy, uyz, uzz, det, x(i,j,k), y(i,j,k), & z(i,j,k), r(i,j,k),whisky_rho_min, whisky_reflevel, whisky_C2P_failed(i,j,k)) - +#endif + + end if if (eps(i,j,k) .lt. 0.0d0) then |