PROGRAM MAP_3 c c Givn two uniform grids defined by the
locations of the lefthand c zone interfaces XL and XNUL. The
grids are offset by an amount c DX, that is XNUL(I) = XL(I) +
DX. c This program creates piecewise parabolic interpolations c
from zone averaged data, defined on the XL grid. c Monotonicity
constraints are applied. For those zones c considered to be in a
discontinuity, an interpolation more c appropriate for
discontinuous functions is used, producing c a better
representation of unresolved shocks. c Regions of over lap
between the two grid are determined and c advective fluxes found
by integrating the PPM interpolations c over these regions. The
zone averages on the grid XNUL are found c by differeing these
fluxes. c c Note, the interpolations require 3 boundary zones at
each end, and c the integrations will require an additional
zone. c PARAMETER ( NBDY = 4 ) ALLOCATABLE :: XL(:), A_AVG(:),
AL(:), AR(:), DA(:), A6(:), ^ SMALLDA(:), UNSMOOTH(:),
DISCONT(:), DX(:), DA_AVL(:), ^ DVOLL(:), XF(:),