calc_h.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Function name: CALC_H                                               C
!  Purpose: Calculate specific enthalpy of species N in phase M        C
!                                                                      C
!  Author: M. Syamlal                                Date: 27-DEC-2007 C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      DOUBLE PRECISION FUNCTION calc_h(refT, M, NN)

! Modules
!---------------------------------------------------------------------//
      USE physprop, only: mw_g, mw_s, hfrefor
      USE constant, only: rgas => gas_const_cal
      USE read_thermochemical, only: calc_icpor
      IMPLICIT NONE

! Dummy arguments
!---------------------------------------------------------------------//
! cell, phase and species indices
      DOUBLE PRECISION, INTENT(IN) :: refT   ! Temperature
      INTEGER, INTENT(IN) :: M ! Phase index
      INTEGER, INTENT(IN) :: NN ! Species index

! Local variables
!---------------------------------------------------------------------//
      DOUBLE PRECISION :: ICpoR
      DOUBLE PRECISION :: lMW
      INTEGER :: IER
!---------------------------------------------------------------------//

      ier = 0

      if(m == 0)then
         lmw = mw_g(nn)
      else
         lmw = mw_s(m,nn)
      endif

! Integrate the specific heat from zero to refT
      icpor = calc_icpor(reft, m, nn, ier)

! Evaluate the enthalpy of species N at refT
      calc_h = (hfrefor(m,nn)  + icpor) * (rgas / lmw)

      RETURN
      END FUNCTION calc_h