diff --git a/src/problem.F90 b/src/problem.F90 index e787696..b7a4a85 100644 --- a/src/problem.F90 +++ b/src/problem.F90 @@ -1370,7 +1370,7 @@ module problem use config , only : gamma, gammam1i #endif /* ADI */ use constants, only : dpi - use coords , only : ax, ay, az, adx, ady, adz + use coords , only : ax, ay, az, adr use variables, only : idn, imx, imy, imz #ifdef ADI use variables, only : ien @@ -1390,7 +1390,7 @@ module problem ! local variables ! integer :: i, j, k - real :: dx, dy, dz, dr + real :: dr real :: dnamb, pramb real :: dnstar, vrstar, prstar, rc real :: dnsat , vrsat , prsat , rs, xs, ys @@ -1446,17 +1446,9 @@ module problem xvl = - asat * om * sn yvl = bsat * om * cs -! obtain cell sizes +! obtain the diagonal size of a cell at the current level ! - dx = adx(pdata%meta%level) - dy = ady(pdata%meta%level) -#if NDIMS == 3 - dz = adz(pdata%meta%level) - dr = sqrt(dx * dx + dy * dy + dz * dz) -#else /* NDIMS == 3 */ - dz = 1.0d0 - dr = sqrt(dx * dx + dy * dy) -#endif /* NDIMS == 3 */ + dr = adr(pdata%meta%level) ! obtain block coordinates !