User: knarf Date: 2010/07/11 09:54 PM
Modified: /trunk/src/ GRHydro_Prim2Con.F90
Log: improve debug output for NaNs
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_Prim2Con.F90 Delta lines: +5 -1 =================================================================== --- trunk/src/GRHydro_Prim2Con.F90 2010-07-02 22:14:28 UTC (rev 135) +++ trunk/src/GRHydro_Prim2Con.F90 2010-07-12 02:54:26 UTC (rev 136) @@ -261,6 +261,7 @@ integer :: i, j, k CCTK_REAL :: gxxl,gxyl,gxzl,gyyl,gyzl,gzzl,avg_detl,& gxxr,gxyr,gxzr,gyyr,gyzr,gzzr,avg_detr,psi4r,psi4l + character(len=200) warnline
!$OMP PARALLEL DO PRIVATE(i, j, gxxl,gxyl,gxzl,gyyl,gyzl,gzzl,avg_detl,& !$OMP gxxr,gxyr,gxzr,gyyr,gyzr,gzzr,avg_detr,psi4r,psi4l) @@ -317,7 +318,10 @@ w_lorentzplus(i,j,k)) if (densminus(i,j,k) .ne. densminus(i,j,k)) then !$OMP CRITICAL - call CCTK_WARN(1, "NaN in densminus(i,j,k) (Prim2Con)") + call CCTK_WARN(1, "NaN in densminus(i,j,k) (Prim2Con)") + write(warnline,' (a27,3g16.7)') 'avg_detl, velx, rhominus: ', & + avg_detl, velx(i,j,k), rhominus(i,j,k) + call CCTK_WARN(1, warnline) !$OMP END CRITICAL endif end do