C     THIS IS THE EXAMPLE PROGRAM DESCRIBING HOW TO USE KDB ROUTINES.
      PROGRAM HC_CPL_E
      IMPLICIT DOUBLE PRECISION (A-H,O-Z) 
      INTEGER NCMP
      CHARACTER*50 CNAME
      CHARACTER*20 FORMULA 
      COMMON/HC_NAME/ICMP(50),CNAME(50),FORMULA(50),ICLASS(50) 
      COMMON/HC_PROP/WT(50),TB(50),TF(50),TC(50),PC(50),VC(50),ZC(50)
	1    ,ACCF(50),WSRK(50),VEST(50),ZRA(50),SOLP(50),VOLP(50)
     2    ,QI(50),RI(50),DM(50) 
      COMMON/HC_KCPL/IEQ_CPL(50),CPL_A(50),CPL_B(50),CPL_C(50),
     1	CPL_D(50),CPL_E(50),CPL_F(50),CPL_G(50),CPL_T1(50),CPL_T2(50)
      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)
	DIMENSION TEXP(5),CPLEXP(5)
	DATA (TEXP(J),J=1,5) /100.D0,200.D0,230.D0,300.D0,360.D0/
	DATA (CPLEXP(J),J=1,5) /85.2D0,93.5D0,98.88D0,122.D0,260.1D0/
	OPEN(UNIT=5,FILE="RESULT.TXT")
C     NUMBER OF COMPONENT = 1
      NCMP = 1
C     SPECIFY THE ID OF COMPONENT 
      ICMP(1) = 3
	CNAME(1) = "PROPANE"
C     SPECIFY THE COEFFICIENT DATA FOR HC_CPGEQN
	IEQ_CPG(1) = 1		
      CPG_A(1)=2.95952D1
	CPG_B(1)=8.379912D-2
	CPG_C(1)=3.255759D-4
	CPG_D(1)=-3.957572D-7
	CPG_E(1)=1.312889D-10
C	SPECIFY THE COEFFICIENT DATA FOR HC_CPLEQN
	IEQ_CPL(1) = 1		
      CPL_A(1)=61.3642D0
	CPL_B(1)=0.4302505D0
	CPL_C(1)=-0.0024649D0
	CPL_D(1)=5.564475D-6
	CPL_T1(1)=85.45D0
	CPL_T2(1)=353.15D0
C	PURE PROPERTIES FOR 1-OCTANOL
	TC(1)=369.83D0			 
	ACCF(1)=0.153D0

	WRITE(5,*) CNAME(1)
	WRITE(5,*) 'COMPARISON OF EXPERIMENTAL AND ESTIMATED VALUE ',
	1		'(% ERROR)'
	WRITE(5,*) '-----------------------------------------------------'
	WRITE(5,*) ' T(K)   CPL(EXP)   KDBEQN     CPLCSP      IST1  IST2'
	WRITE(5,*) '-----------------------------------------------------'
	
	DO I = 1, 5
		CALL HC_CPLEQN(1,TEXP(I),CPL1,IST1)
		CALL HC_CPLCSP(1,TEXP(I),CPL2,IST2)
		ERR1 = DABS( (CPL1-CPLEXP(I))/CPLEXP(I) )*100.D0
		ERR2 = DABS( (CPL2-CPLEXP(I))/CPLEXP(I) )*100.D0
		WRITE(5,5) TEXP(I),CPLEXP(I),ERR1,ERR2,IST1,IST2
	ENDDO
5	FORMAT(1x,F6.2,1X, 3(1X,G10.5),2I6) 
	WRITE(5,*) '-----------------------------------------------------'
	CLOSE(5) 
      STOP 
      END