[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 432)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Thu Nov 8 19:53:06 CST 2012
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
More information about the Commits
mailing list