From 39483913e9fcb05bf4ab2b6d30b79b8742b52f6c Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Fri, 29 Aug 2014 09:08:13 -0300 Subject: [PATCH] COORDINATES: Add domain boundary areas. Signed-off-by: Grzegorz Kowal --- src/coordinates.F90 | 12 ++++++++++++ 1 file changed, 12 insertions(+) diff --git a/src/coordinates.F90 b/src/coordinates.F90 index 9281ec2..2ecca6b 100644 --- a/src/coordinates.F90 +++ b/src/coordinates.F90 @@ -88,6 +88,12 @@ module coordinates real(kind=8), save :: vol = 1.0d+00 real(kind=8), save :: voli = 1.0d+00 +! the domain boundary areas +! + real(kind=8), save :: xarea = 1.0d+00 + real(kind=8), save :: yarea = 1.0d+00 + real(kind=8), save :: zarea = 1.0d+00 + ! the block coordinates for all levels of refinement ! real(kind=8), dimension(:,:), allocatable, save :: ax , ay , az @@ -248,6 +254,12 @@ module coordinates vol = xlen * ylen * zlen voli = 1.0d+00 / vol +! calculate the boundary areas +! + xarea = ylen * zlen + yarea = xlen * zlen + zarea = xlen * ylen + ! allocate space for coordinate variables ! allocate(ax (toplev, im))