diff options
| author | allen <allen@c78560ca-4b45-4335-b268-5f3340f3cb52> | 1999-03-01 13:06:13 +0000 |
|---|---|---|
| committer | allen <allen@c78560ca-4b45-4335-b268-5f3340f3cb52> | 1999-03-01 13:06:13 +0000 |
| commit | 75e618a69ac9b56f303b8ac98664eaa6ea686d1f (patch) | |
| tree | 3fbb26c3505503376f8b006aaee9394eff9ec372 | |
| parent | 9599406fb54226fe67f0f19c520de7ebe9f6d0af (diff) | |
Setting flesh delta_space variables
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/CartGrid3D/trunk@7 c78560ca-4b45-4335-b268-5f3340f3cb52
| -rw-r--r-- | src/CartGrid3D.F | 23 |
1 files 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 |