weno5_substencils Subroutine

public pure subroutine weno5_substencils(f_stencil, vp0, vp1, vp2)

Compute the three 3rd-order substencil reconstructions for WENO5.

Operates component-wise on a 5-point stencil. The reconstructed values are the same for WENO5-JS, WENO5-Z, and TENO5.

Substencil approximations (Jiang & Shu 1996, Eq. 16a-16c): vp0 = (1/3)v1 - (7/6)v2 + (11/6)v3 vp1 = -(1/6)v2 + (5/6)v3 + (1/3)v4 vp2 = (1/3)v3 + (5/6)v4 - (1/6)*v5

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: f_stencil(:,:)
real(kind=wp), intent(out) :: vp0(:)
real(kind=wp), intent(out) :: vp1(:)
real(kind=wp), intent(out) :: vp2(:)