Second-order upwind reconstruction.
Uses a 3-point stencil (half_stencil = 2, stencil_width = 3). For the positive (left-biased) stencil: stencil(:,1) = F^+{i-1}, stencil(:,2) = F^+_i (centre) F_hat^+ = 1.5 * F^+_i - 0.5 * F^+
For the negative (right-biased, reversed) stencil: stencil(:,1) = F^-{i+2}, stencil(:,2) = F^- (centre) F_hat^- = 1.5 * F^-{i+1} - 0.5 * F^- [same formula]
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | f_stencil(:,:) | |||
| real(kind=wp), | intent(out) | :: | f_hat(:) |