Compute interstates velocity and pressure.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(riemann_pattern_compressible_pvl), | intent(inout) | :: | self | Riemann (states) pattern solution. |
elemental subroutine compute_up23(self)
!< Compute interstates velocity and pressure.
class(riemann_pattern_compressible_pvl), intent(inout) :: self !< Riemann (states) pattern solution.
real(R8P) :: ram !< Mean value of `r * a`.
ram = 0.25_R8P * (self%r_1 + self%r_4) * (self%a_1 + self%a_4)
self%u23 = 0.5_R8P * ((self%u_1 + self%u_4) - (self%p_4 - self%p_1) / ram)
self%p23 = 0.5_R8P * ((self%p_1 + self%p_4) - (self%u_4 - self%u_1) * ram)
endsubroutine compute_up23