this assumes that the backward difference for hte first cell has been set the the same as the forward difference for the final cell, i.e. x%dx3all(lbnd)==x%dx3all(ubnd+1). In general when doing periodic grids it is probably best to hard code all of the differences outside the domain to be equal (or to use uniform meshes)
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(:,:,:) | :: | f | ||
type(curvmesh), | intent(in) | :: | x | |||
integer, | intent(in) | :: | lbnd | |||
integer, | intent(in) | :: | ubnd |