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

[knarf at cct.lsu.edu]

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