compute_resid_2d_curvilinear Subroutine

private subroutine compute_resid_2d_curvilinear(state)

Curvilinear residual: R = -(1/V)[ (Fhat_{i+1}-Fhat_i) + (Ghat_{j+1}-Ghat_j) ] with contravariant face fluxes Fhat = |S_xi| * f_n(qL,qR; n_xi), Ghat = |S_eta| * f_n(qL,qR; n_eta). Reconstruction is identical to the uniform path (uniform computational stencil); geometry enters via s_xi/s_eta/vol.

Arguments

Type IntentOptional Attributes Name
type(solver_state_2d_t), intent(inout) :: state