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

This commit is contained in:
Grzegorz Kowal 2011-05-12 08:58:13 -03:00
parent 355d7c900f
commit 62c1aabba4

View File

@ -701,7 +701,7 @@ module problem
use config , only : dnfac, dnrat use config , only : dnfac, dnrat
use config , only : rstar, vstar, rsat, dsat, vsat, tsat use config , only : rstar, vstar, rsat, dsat, vsat, tsat
use constants, only : dpi use constants, only : dpi
use coords , only : ax, ay, az, adx, ady, adz use coords , only : ax, ay, az, adr
use scheme , only : prim2cons use scheme , only : prim2cons
use variables, only : nqt use variables, only : nqt
use variables, only : idn, ivx, ivy, ivz use variables, only : idn, ivx, ivy, ivz
@ -722,7 +722,7 @@ module problem
! local variables ! local variables
! !
integer :: i, j, k integer :: i, j, k
real :: dx, dy, dz, dr real :: dr
real :: dnamb, pramb real :: dnamb, pramb
real :: dnstar, prstar, vrstar, rc real :: dnstar, prstar, vrstar, rc
real :: dnsat , prsat , vrsat , rs, xs, ys real :: dnsat , prsat , vrsat , rs, xs, ys
@ -766,17 +766,9 @@ module problem
xvl = - asat * om * sn xvl = - asat * om * sn
yvl = bsat * om * cs yvl = bsat * om * cs
! obtain cell sizes ! obtain the diagonal size of a cell at the current level
! !
dx = adx(pdata%meta%level) dr = adr(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 */
! obtain block coordinates ! obtain block coordinates
! !