merge branch hot_and_MHD_temp_dev into branch at revision r185

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@186 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 35 insertions, 0 deletions

diff --git a/src/Utils.F90 b/src/Utils.F90
index 0ef9af2..0de6b74 100644
--- a/src/Utils.F90
+++ b/src/Utils.F90
@@ -27,6 +27,41 @@
    @endhistory 
 @@*/
 
+subroutine GRHydro_Debug(CCTK_ARGUMENTS)
+
+  implicit none
+  DECLARE_CCTK_ARGUMENTS
+  integer i,j,k
+  integer nx, ny, nz
+
+  GRHydro_reflevel = aint(log10(dble(cctk_levfac(1)))/log10(2.0d0))
+  
+  nx = cctk_lsh(1)
+  ny = cctk_lsh(2)
+  nz = cctk_lsh(3)
+
+  if (GRHydro_reflevel .lt. 3) return
+ 
+  do i=4,nx
+     do j=4,ny
+        do k=4,4
+           if( r(i,j,k)-4.0d0 .lt. 96.0d0 .and. &
+                r(i,j,k)+4.0d0 .gt. 96.0d0) then
+              write(6,"(4i4,1P10E15.6)") GRHydro_reflevel, i,j,k,&
+                   r(i,j,k),rho(i,j,k),eps(i,j,k),y_e(i,j,k),temperature(i,j,k)
+           endif
+        enddo
+     enddo
+  enddo
+
+!  call flush(6)
+!  if(GRHydro_reflevel.eq.4.and.temperature(1,1,1).lt.0.1d0) then
+!     call CCTK_WARN(0,"stop")
+!  endif
+
+end subroutine GRHydro_Debug
+
+
 subroutine GRHydro_RefinementLevel(CCTK_ARGUMENTS)
 
   implicit none