Return pressure value.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(conservative_compressible), | intent(in) | :: | self | Conservative. |
||
| class(eos_object), | intent(in) | :: | eos | Equation of state. |
Pressure value.
elemental function pressure(self, eos) result(pressure_)
!< Return pressure value.
class(conservative_compressible), intent(in) :: self !< Conservative.
class(eos_object), intent(in) :: eos !< Equation of state.
real(R8P) :: pressure_ !< Pressure value.
type(vector) :: velocity_ !< Velocity vector.
velocity_ = self%velocity()
pressure_ = (eos%g() - 1._R8P) * (self%energy - 0.5_R8P * self%density * velocity_%sq_norm())
endfunction pressure