C     THIS IS THE EXAMPLE PROGRAM DESCRIBING HOW TO USE KDB ROUTINES.
      PROGRAM HC_HVP_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) 
	DIMENSION TEXP(5),HVPEXP(5)
	DATA (TEXP(J),J=1,5) /298.2, 358.2, 411.2, 431.1,499.7/
	DATA (HVPEXP(J),J=1,5) /56.94D3, 51.22D3, 44.37D3,
	1						41.32D3, 31.61D3/
	OPEN(UNIT=5,FILE="RESULT.TXT")

C     NUMBER OF COMPONENT = 1
      NCMP = 1
C     SPECIFY THE ID OF EACH COMPONENT 
      ICMP(1) = 825
	CNAME(1) = "1-PENTANOL"
	ICLASS(1) = 11
C	SPECIFY THE COMPONENT DATA 
	TC(1) =	588.1D0
	PC(1) = 3897.D0
	TB(1) = 411.13D0
C     PRINT COMPONENT NAMES
	WRITE(5,*) CNAME(1)
	WRITE(5,*) 'COMPARISON OF EXPERIMENTAL AND ESTIMATED VALUE ',
	1		'(% ERROR)'
	WRITE(5,*) '-----------------------------------------------------'
	WRITE(5,*) ' T(K)   HVP(EXP)   VETERE     ERROR'
	WRITE(5,*) '-----------------------------------------------------'
	
	DO I = 1, 5
		CALL HC_HVPVET(1,TEXP(I),HVP1,IST1)
		ERR1 = DABS( (HVP1-HVPEXP(I))/HVPEXP(I) )*100.D0
		WRITE(5,5) TEXP(I),HVPEXP(I),HVP1,ERR1
	ENDDO
5	FORMAT(1x,F6.2,1X, 3(1X,G10.5)) 
	WRITE(5,*) '-----------------------------------------------------'
      STOP 
      END