C
C
C     LIST OF SUBROUTINES IN THIS FILE 
C
C     HC_CPGEQN : IDEAL GAS HEAT CAPACITY EQUATIONS 
C
C
C*******************************************************************
C
C KDB THEMOPHYSICAL PROPERTIES CALCULATION FORTRAN ROUTINE LIBRARY 
C
C [NAME   ]  HC_CPGEQN 
C
C [TYPE   ]  FORTRAN SUBROUTINE                              
C
C [PURPOSE]  IDEAL GAS HEAT CAPACITY USING INTERNAL CORRELATION EQUATION 
C
C [USAGE  ]  CALL HC_CPGEQN(ICN,T,CPG,IST)
C
C [ARGUMENTS]                                                                  
C     ICN        : COMPONENT NUMBER (1-50) TO CALCULATE VAPOR PRESSURE (INTEGER, INPUT) 
C     T          : TEMPERATURE IN KELVIN (REAL*8, INPUT) 
C
C     CPG        : IDEAL GAS HEAT CAPACITY IN KJ/KG-MOL.K (REAL*8, OUTPUT) 
C     IST        : STATUS OF CALCULATION (INTEGER, OUTPUT) 
C                = 0	    : NORMAL TERMINATION 
C                = 301   	: IDEAL GAS HEAT CAPACITY COEFFICIENT NOT AVAILABLE 
C
C [COMMENTS]
C     
C [REQUIRED COMMON BLOCKS]
C     COMMON /HC_KCPG/ IDEAL GAS HEAT CAPACITY COEFFICIENTS 
C                                                                   
C [REQUIRED SUBROUTINES OR FUNCTIONS] 
C	NONE 
C
C [REFERENCE]
C     NONE  
C	
C [REVISION INFORMATION]
C     1.PROGRMAMMED BY J.W.KANG, KOREA UNIVERSITY, 1996
C     2.REVISED BY J.W.KANG, KOREA UNIVERSITY, 2000                                                              
C     3.REVISED BY J.W.KANG, KOREA UNIVERSITY, 2001  
C*******************************************************************
      SUBROUTINE HC_CPGEQN(ICN,T,CPG,IST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
	INTEGER IST,ICN  
      REAL*8  T,CPG
      COMMON/HC_KCPG/IEQ_CPG(50),CPG_A(50),CPG_B(50),CPG_C(50),
     1	CPG_D(50),CPG_E(50),CPG_F(50),CPG_G(50),CPG_T1(50),CPG_T2(50)
C --- INITIALIZATION 
      IST = 0
	CPG = 0.0
C --- CHECK IF EQUATION COEFF.S AVAILABLE 
	IF(IEQ_CPG(ICN).LE.0) IST = 301 
C --- CHECK IF THE RANGE OF APPLICATION 
	IF(IST.NE.0) RETURN 
C --- LIQUID HEAT CAPAICITY CALCULATION 
	CPG = CPG_A(ICN) + CPG_B(ICN)*T + CPG_C(ICN)*T*T + 
	1	CPG_D(ICN)*T*T*T + CPG_E(ICN)*T*T*T*T 
      RETURN
      END