diff options
-rw-r--r-- | src/CartGrid3D.F | 38 |
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 |