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.')