filrecur

Signature:

recursive subroutine filrecur(level, nvar, valbig, aux, naux, time, mitot, mjtot, nrowst, ncolst, ilo, ihi, jlo, jhi)

Arguments:

integer, intent(in):

• level: Level of the input region.
• nvar: Number of solution variables.
• naux: Number of aux variables.
• mitot: Horizontal dimension of the input region.
• mjtot: Vertical dimension of the input region.
• nrowst: Starting i-index of the patch, in indices relative to valbig.
• ncolst: Starting j-index of the patch, in indices relative to valbig.
• ilo: Lower horizontal index of the patch to fill, in global indices.
• ihi: Upper horizontal index of the patch to fill, in global indices.
• jlo: Lower vertical index of the patch to fill, in global indices.
• jhi: Upper vertical index of the patch to fill, in global indices.

double precision, intent(in):

• time: Time at which values are needed.

double precision, intent(inout):

• valbig(mitot,mjtot,nvar): Solution values on the input region.
• aux(mitot,mjtot,naux): Aux variable values on the input region.

Description:

Fills in a rectangular patch of valbig and aux, which contain solution and aux values on an input region. Values are needed at time time, and at level level.

In indices relative to valbig and aux, the lower-left corner of the patch is at (nrowst,ncolst). In global indices, its lower-left corner is at (ilo,jlo), and its upper-right corner is at (ihi,jhi).

First, the patch is filled with values obtainable from the level level grids. If any values are left unfilled, the remaining rectangle of unfilled values is enlarged by 1 (for later linear interpolation), and remaining values are recursively obtained from coarser levels.