$DEBUG
C     THIS IS THE EXAMPLE PROGRAM DESCRIBING HOW TO USE KDB ROUTINES.
      PROGRAM HC_THL_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_KTHL/IEQ_THL(50),THL_A(50),THL_B(50),THL_C(50),
     1	THL_D(50),THL_E(50),THL_F(50),THL_G(50),THL_T1(50),THL_T2(50)

	REAL*8 TEXP(3),THLEXP(3)
C	EXPERIMENTAL DATA FOR ETHANEOL AT LOW PRESSURE
	DATA (TEXP(J),J=1,3) /200.D0,300.D0,400.D0/
	DATA (THLEXP(J),J=1,3) /0.196D0,0.166D0,0.143D0/
	OPEN(UNIT=5,FILE="RESULT.TXT")
C     NUMBER OF COMPONENT = 1
      NCMP = 1

C     SPECIFY THE ID OF EACH COMPONENT 
      ICMP(1) = 818		! EHTANOL
	CNAME(1) = "ETHANOL"
	ICLASS(1) = 11	
C	LOW PRESSURE GAS THERMAL CONDUCTIVITY COEFFICIENT FOR ETHANOL
	IEQ_THL(1) = 1		
      THL_A(1)=0.2629D0
	THL_B(1)=-0.0003847D0
	THL_C(1)=2.211D-7
	THL_D(1)=0.D0
	THL_T1(1)=160.D0
	THL_T2(1)=463.D0

C	PURE PROPERTIES FOR ETHANOL
	WT(1)=4.60688D1
	TC(1)=514.0D0
	TB(1)=351.44D0			 

	WRITE(5,*) CNAME(1)
	WRITE(5,*) 'COMPARISON OF EXPERIMENTAL AND ESTIMATED VALUE ',
	1		'(% ERROR)'
	WRITE(5,*) '-----------------------------------------------------'
	1	   ,'----------------'
	WRITE(5,*) ' T(K)   THL(EXP)     KDB      LATINI ',
	1		'     IST1  IST2'
	WRITE(5,*) '-----------------------------------------------------'
	1	   ,'----------------'

	DO I = 1, 3
		CALL HC_THLEQN(1,TEXP(I),THL1,IST1)
		CALL HC_THLLATINI(1,TEXP(I),THL2,IST2)
		ERR1 = DABS( (THL1-THLEXP(I))/THLEXP(I) )*100.D0
		ERR2 = DABS( (THL2-THLEXP(I))/THLEXP(I) )*100.D0
		WRITE(5,5) TEXP(I),THLEXP(I),ERR1,ERR2,IST1,IST2
	ENDDO
5	FORMAT(1x,F6.2,1X,3(1X,G10.5),2I6) 
	WRITE(5,*) '-----------------------------------------------------'
	1	   ,'----------------'


	CLOSE(5) 
      STOP 
      END