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