CpCurveFits.f90

Go to the documentation of this file.

module cpcurvefits
!
!       This module contains the curve fit data for Cp, or better
!       Cp/R, as a function of the temperature. The temperature is
!       assumed to be in Kelvin.
!
    use constants, only: inttype, realtype
    implicit none
    save
!
!       The definition of the derived data type CpTempFitType, which
!       stores the curve fit data for a certain temperature range.
!
    type cptempfittype
 
        ! nterm:     Number of terms in the polynomial expansion.
        ! exponents: The powers of the polynomial fit.
        ! constants: The constants in front of each term.
 
        integer(kind=intType) :: nterm
        integer(kind=intType), dimension(:), pointer :: exponents
        real(kind=realtype), dimension(:), pointer :: constants
 
        ! Additional constant to compute the internal energy.
 
        real(kind=realtype) :: eint0
 
        ! Values of the integrand of Cp/(R*T) at the lower (_1) and
        ! upper (_2) curve fit boundary. Needed to compute the total
        ! pressure.
 
        real(kind=realtype) :: intcpovrt_1, intcpovrt_2
 
    end type cptempfittype
 
    ! CpNParts:             Number of temperature ranges.
    ! CpTRange(0:CpNParts): The temperature ranges.
    ! CpEint(0:CpNParts):   Internal energies at the curve fit
    !                       boundaries.
    ! CpHint(0:CpNParts):   Internal enthalpies at the curve fit
    !                       boundaries.
    ! CpTempFit(CpNParts):  The actual curve fit data.
 
    integer(kind=intType) :: cpnparts
 
    real(kind=realtype), dimension(:), allocatable :: cptrange
    real(kind=realtype), dimension(:), allocatable :: cpeint
    real(kind=realtype), dimension(:), allocatable :: cphint
    type(cptempfittype), dimension(:), allocatable :: cptempfit
 
    ! cv0:   The cv value at the the temperature CpTRange(0). If the
    !        temperature is lower than the lowest curve fit boundary
    !        this value is needed to extrapolate the energy (assuming
    !        constant cv).
    ! cvn:   Idem, but than at CpTRange(CpNParts). Used when the
    !        temperature is higher than the highest curve fit boundary.
 
    real(kind=realtype) :: cv0, cvn
 
end module cpcurvefits