NAME
Grid_subcellGeometry
SYNOPSIS
call Grid_subcellGeometry(integer,VALUE(in) :: nsubi,
integer,VALUE(in) :: nsubj,
integer,VALUE(in) :: nsubk,
real(in) :: dvcell,
real, intent(OUT) :: dvsub(nsubI,nsubJ),
real,OPTIONAL(in) :: xl,
real,OPTIONAL(in) :: xr,
real,OPTIONAL(in) :: yl,
real,OPTIONAL(in) :: yr,
integer,OPTIONAL(in) :: pos(*),
integer,OPTIONAL(in) :: blockID)
DESCRIPTION
Geometrically correct computation of the volumes of subcells.
ARGUMENTS
nsubi : Number of subcell lengths per cell length in X direction.
nsubj : Number of subcell lengths per cell length in Y direction.
nsubk : Number of subcell lengths per cell length in Z direction.
dvcell : Volume of the whole cell.
dvsub : Volumes of subcells. Note that this is a 2-dimensional array.
For the geometries that FLASH supports, subcells volumes
never depend on the third coordinate.
xl : X-coordinate of left cell face.
xr : X-coordinate of right cell face.
yl : Y-coordinate of lower cell face.
yr : Y-coordinate of upper cell face.
pos : currently unused.
blockID : ID of block in current processor, currently unused