Added debugging stuff

git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/CartGrid3D/trunk@15 c78560ca-4b45-4335-b268-5f3340f3cb52

1 files changed, 28 insertions, 10 deletions

diff --git a/src/CartGrid3D.F b/src/CartGrid3D.F
index b04f5a7..336f878 100644
--- a/src/CartGrid3D.F
+++ b/src/CartGrid3D.F
@@ -7,7 +7,7 @@
    @enddesc 
  @@*/
 
-/*#define CARTGRID3D_DEBUG*/
+#define CARTGRID3D_DEBUG
 
 #include "cctk.h"
 #include "declare_arguments.h"
@@ -108,10 +108,15 @@ c     -----------------------------------------------------------
 
       else if (CCTK_Equals(grid,'byspacing') == 1) then
 
+c        Dx, Dy, Dx on the coarsest grid
+         coarse_dx = grid_dx
+         coarse_dy = grid_dy
+         coarse_dz = grid_dz
+
 c        dx, dy, dz on the grid we are on
-         dx = grid_dx*iconv/levfac
-         dy = grid_dy*iconv/levfac
-         dz = grid_dz*iconv/levfac
+         dx = coarse_dx*iconv/levfac
+         dy = coarse_dy*iconv/levfac
+         dz = coarse_dz*iconv/levfac
 
          if (CCTK_Equals(symmetry,'bitant') == 1) then 
 
@@ -202,12 +207,6 @@ c     -----------------------
                x(i,j,k) = dx*(i+lb(1)-1) + x0
                y(i,j,k) = dy*(j+lb(2)-1) + y0
                z(i,j,k) = dz*(k+lb(3)-1) + z0
-
-#ifdef CARTGRID3D_DEBUG
-               write(*,*) "Index: ",i,j,k
-               write(*,*) "Coord: ",x(i,j,k),y(i,j,k),z(i,j,k)
-#endif
-
             end do
          end do
       end do
@@ -218,6 +217,25 @@ c     -----------------------
       delta_space(2) = dy
       delta_space(3) = dz
 
+#ifdef CARTGRID3D_DEBUG
+      write(*,*)
+      write(*,*) "CartGrid3D"
+      write(*,*) "----------"
+      write(*,*) "Dx, Dy, Dz on coarse grid = ",
+     &     coarse_dx,coarse_dy,coarse_dz
+      write(*,*) "Dx, Dy, Dz on this grid = ",delta_space(1),
+     &     delta_space(2),delta_space(3)
+      write(*,*) "Convergence level = ",convlevel
+      write(*,*) "Grid level = ",levfac
+      write(*,*) "Minimum Global Coords = ",x0,y0,z0
+      write(*,*) "Maximum Global Coords = ",x0+dx*(global_sh(1)-1),
+     &     y0+dy*(global_sh(2)-1),z0+dz**(global_sh(3)-1)
+      write(*,*) "Minimum Local Coords = ",x(1,1,1),
+     &     y(1,1,1),z(1,1,1)
+      write(*,*) "Maximum Local Coords = ",x(sh(1),sh(2),sh(3)),
+     &     y(sh(1),sh(2),sh(3)),z(sh(1),sh(2),sh(3))
+#endif
+
       return
       end subroutine CartGrid3D