[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 322)

[roland.haas at physics.gatech.edu]

User: rhaas
Date: 2012/05/13 02:08 PM

Modified:
 /trunk/src/
  GRHydro_Con2Prim.F90

Log:
 Error check for C2P_failed mask greater or equal 1, do not abort code if in the
 atmosphere though.
 
 patch by Philipp Moesta

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_Con2Prim.F90
Delta lines: +4 -1
===================================================================
--- trunk/src/GRHydro_Con2Prim.F90	2012-05-13 19:03:31 UTC (rev 321)
+++ trunk/src/GRHydro_Con2Prim.F90	2012-05-13 19:08:58 UTC (rev 322)
@@ -1980,7 +1980,7 @@
     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 @@
               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.')