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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/11/08 07:53 PM

Modified:
 /trunk/src/
  GRHydro_Con2Prim.F90, GRHydro_Prim2ConM.F90

Log:
 GRHydro: add extensive information output to check_C2P_failed
 
 From: Roland Haas <roland.haas at physics.gatech.edu>

File Changes:

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

File [modified]: GRHydro_Con2Prim.F90
Delta lines: +55 -17
===================================================================
--- trunk/src/GRHydro_Con2Prim.F90	2012-11-09 01:53:00 UTC (rev 431)
+++ trunk/src/GRHydro_Con2Prim.F90	2012-11-09 01:53:05 UTC (rev 432)
@@ -2231,18 +2231,8 @@
   DECLARE_CCTK_ARGUMENTS
 
   integer :: i, j, k, nx, ny, nz
-  character(len=100) warnline
+  character(len=300) warnline
 
-  ! save memory when MP is not used
-  CCTK_INT :: GRHydro_UseGeneralCoordinates
-  CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup
-
-  if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then
-    vup => lvel
-  else
-    vup => vel
-  end if
-
 !  call CCTK_INFO("Checking the C2P failure mask.")
 
   nx = cctk_lsh(1)
@@ -2273,16 +2263,64 @@
           !$OMP CRITICAL
           call CCTK_WARN(1,'Con2Prim failed; stopping the code.')
           call CCTK_WARN(1,'Even with mesh refinement, this point is not restricted from a finer level, so this is really an error')
-          write(warnline,'(a28,i2)') 'on carpet reflevel: ',GRHydro_reflevel
+          write(warnline,'(a32,i2)') 'on carpet reflevel: ',GRHydro_reflevel
           call CCTK_WARN(1,warnline)
-          write(warnline,'(a20,3g16.7)') 'xyz location: ',&
+          write(warnline,'(a32,3g16.7)') 'xyz location: ',&
                x(i,j,k),y(i,j,k),z(i,j,k)
           call CCTK_WARN(1,warnline)
-          write(warnline,'(a20,g16.7)') 'radius: ',r(i,j,k)
+          write(warnline,'(a32,g16.7)') 'radius: ',r(i,j,k)
           call CCTK_WARN(1,warnline)
-          write(warnline,'(a20,3g16.7)') 'velocities: ',&
-               vup(i,j,k,1),vup(i,j,k,2),vup(i,j,k,3)
-          call CCTK_WARN(0,warnline)
+          write(warnline,'(a32,i3,1g16.7)') 'hydro_excision_mask, C2Pfailed: ', hydro_excision_mask(i,j,k), GRHydro_C2P_failed(i,j,k)
+          call CCTK_WARN(1,warnline)
+          write(warnline,'(a32,4g16.7)') 'rho eps press w_lorentz: ',&
+               rho(i,j,k),eps(i,j,k),press(i,j,k),w_lorentz(i,j,k)
+          call CCTK_WARN(1,warnline)
+          write(warnline,'(a32,3g16.7)') 'velocities: ',&
+               vel(i,j,k,1),vel(i,j,k,2),vel(i,j,k,3)
+          call CCTK_WARN(1,warnline)
+          write(warnline,'(a32,2g16.7)') 'dens tau: ',&
+               dens(i,j,k),tau(i,j,k)
+          call CCTK_WARN(1,warnline)
+          write(warnline,'(a32,3g16.7)') 'scon: ',&
+               scon(i,j,k,1),scon(i,j,k,2),scon(i,j,k,3)
+          call CCTK_WARN(1,warnline)
+          if (evolve_mhd.ne.0) then
+            write(warnline,'(a32,3g16.7)') 'magnetic field: ',&
+                 Bvec(i,j,k,1),Bvec(i,j,k,2),Bvec(i,j,k,3)
+            call CCTK_WARN(1,warnline)
+            write(warnline,'(a32,3g16.7)') 'cons magnetic field: ',&
+                 Bcons(i,j,k,1),Bcons(i,j,k,2),Bcons(i,j,k,3)
+            call CCTK_WARN(1,warnline)
+            if (calculate_bcom.ne.0) then
+              write(warnline,'(a32,6g16.7)') 'com magnetic field, b0, b^2: ', bcom(i,j,k,1), bcom(i,j,k,2), bcom(i,j,k,3), bcom0(i,j,k), bcom_sq(i,j,k)
+              call CCTK_WARN(1,warnline)
+            end if
+            if (clean_divergence.ne.0) then
+              write(warnline,'(a32,g16.7)') 'psidc: ',psidc(i,j,k)
+              call CCTK_WARN(1,warnline)
+            end if
+            if (track_divB.ne.0) then
+              write(warnline,'(a32,g16.7)') 'divB: ',divB(i,j,k)
+              call CCTK_WARN(1,warnline)
+            end if
+          end if
+          if (evolve_temper.ne.0) then
+            write(warnline,'(a32,2g16.7)') 'temperature, entropy: ',&
+                  temperature(i,j,k),entropy(i,j,k)
+            call CCTK_WARN(1,warnline)
+          end if
+          if (evolve_Y_e.ne.0) then
+            write(warnline,'(a32,2g16.7)') 'Y_e, Y_e_con: ',&
+                  Y_e(i,j,k),Y_e_con(i,j,k)
+            call CCTK_WARN(1,warnline)
+          end if
+          write(warnline,'(a32,6g16.7)') '3-metric: ',&
+               gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)
+          call CCTK_WARN(1,warnline)
+          write(warnline,'(a32,4g16.7)') 'lapse, shift: ',&
+               alp(i,j,k),betax(i,j,k),betay(i,j,k),betaz(i,j,k)
+          call CCTK_WARN(1,warnline)
+          call CCTK_WARN(0,"Aborting.")
           !$OMP END CRITICAL
 
         end if

File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +2 -0
===================================================================
--- trunk/src/GRHydro_Prim2ConM.F90	2012-11-09 01:53:00 UTC (rev 431)
+++ trunk/src/GRHydro_Prim2ConM.F90	2012-11-09 01:53:05 UTC (rev 432)
@@ -386,6 +386,8 @@
 !!$ BEGIN: Check for NaN value
   if (w .ne. w) then
     !$OMP CRITICAL
+    write(NaN_WarnLine,'(a100,6g15.6)') 'NaN produced in sqrt(): (gxx,gxy,gxz,gyy,gyz,gzz)', gxx, gxy, gxz, gyy, gyz, gzz
+    call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine)
     write(NaN_WarnLine,'(a100,3g15.6)') 'NaN produced in sqrt(): (dvelx,dvely,dvelz)', dvelx, dvely, dvelz
     call CCTK_WARN(GRHydro_NaN_verbose, NaN_WarnLine)
     !$OMP END CRITICAL