Use adr(:) to get dr in shape_binaries().

This commit is contained in:
Grzegorz Kowal 2011-05-12 08:55:56 -03:00
parent 0055e86ac2
commit 355d7c900f

View File

@ -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
!