From 75e618a69ac9b56f303b8ac98664eaa6ea686d1f Mon Sep 17 00:00:00 2001 From: allen Date: Mon, 1 Mar 1999 13:06:13 +0000 Subject: Setting flesh delta_space variables git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/CartGrid3D/trunk@7 c78560ca-4b45-4335-b268-5f3340f3cb52 --- src/CartGrid3D.F | 23 ++++++++++++++--------- 1 file changed, 14 insertions(+), 9 deletions(-) diff --git a/src/CartGrid3D.F b/src/CartGrid3D.F index b209ec5..ebfdc5f 100644 --- a/src/CartGrid3D.F +++ b/src/CartGrid3D.F @@ -20,17 +20,12 @@ integer :: CCTK_Equals integer :: iflag, iconv, i, j, k - REAL :: x0,y0,z0, dx, dy, dz + REAL :: x0,y0,z0,dx,dy,dz iflag = 0 iconv = 2**convlevel - print *,"GH = ",GH - print *,"Grid = ",grid - print *,"grid_dx = ",grid_dx - print *,"cctk_itfirst = ",cctk_itlast - c -------------------------------------------------------------- c c 1: grids that ignore dx,dy,dz @@ -88,9 +83,14 @@ c Coordinates are all -0.5 to 0.5 z0 = -0.5 c dx,dy,dz on the coarsest grid of each GH - dx = 1.d0/max(global_sh(1)-1,1) - dy = 1.d0/max(global_sh(1)-1,1) - dz = 1.d0/max(global_sh(1)-1,1) + coarse_dx = 1.d0/max(global_sh(1)-1,1) + coarse_dy = 1.d0/max(global_sh(1)-1,1) + coarse_dz = 1.d0/max(global_sh(1)-1,1) + +c dx,dy,dz on the grid we are on + dx = coarse_dx/levfac + dy = coarse_dy/levfac + dz = coarse_dz/levfac c Special cases if (global_sh(1) == 1) x0 = 0.0D0 @@ -186,6 +186,11 @@ c ----------------------- r = sqrt(x**2 + y**2 + z**2) + delta_space(1) = dx + delta_space(2) = dy + delta_space(3) = dz + + print *,"delta_space",delta_space(1) return end subroutine CartGrid3D -- cgit v1.2.3