Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=r64), | intent(in) | :: | FluidTemp | |||
real(kind=r64), | intent(out) | :: | HLiqWater |
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed arrows point from an interface to procedures which implement that interface. This could include the module procedures in a generic interface or the implementation in a submodule of an interface in a parent module. Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.
SUBROUTINE FigureLiquidWaterEnthalpy(FluidTemp, HLiqWater)
! SUBROUTINE INFORMATION:
! AUTHOR B griffith
! DATE WRITTEN December 2005
! MODIFIED na
! RE-ENGINEERED na
! PURPOSE OF THIS SUBROUTINE:
! calculate Enthalpy from Shomate equations for liquid water
! No enthalpy of formation in this one
! METHODOLOGY EMPLOYED:
! REFERENCES:
! NIST Webbook on gas phase thermochemistry
! USE STATEMENTS:
! na
IMPLICIT NONE ! Enforce explicit typing of all variables in this routine
! SUBROUTINE ARGUMENT DEFINITIONS:
REAL(r64), INTENT(IN) :: FluidTemp ! degree C
REAL(r64), INTENT(OUT) :: HLiqWater ! kJ/mol
! SUBROUTINE PARAMETER DEFINITIONS:
! na
! INTERFACE BLOCK SPECIFICATIONS:
! na
! DERIVED TYPE DEFINITIONS:
! na
! SUBROUTINE LOCAL VARIABLE DECLARATIONS:
REAL(r64) :: Tsho ! temp for Shomate eq in (Kelvin/1000)
REAL(r64) :: A ! shomate coeff
REAL(r64) :: B ! shomate coeff
REAL(r64) :: C ! shomate coeff
REAL(r64) :: D ! shomate coeff
REAL(r64) :: E ! shomate coeff
REAL(r64) :: F ! shomate coeff
REAL(r64) :: H ! shomate coeff
Tsho = (FluidTemp +KelvinConv) / 1000.0d0
A = -203.606d0
B = 1523.29d0
C = -3196.413d0
D = 2474.455d0
E = 3.85533d0
F = -256.5478d0
H = -285.8304d0
HLiqWater = A*Tsho + B*(Tsho**2)/2.0d0 + C*(Tsho**3)/3.0d0 + D*(Tsho**4)/4.0d0 - E/Tsho + F !- H
RETURN
END SUBROUTINE FigureLiquidWaterEnthalpy