10000 REM HGK + Pabalan & Pitzer for KCl(aq) 10050 DEFDBL A-H, M-Z 10100 DIM HGKG(40),II(40),JJ(40),BP(10),BQ(10) 10150 DIM ATZ(4),ADZ(4),AAT(4),AAD(4) 10200 DIM BV(10),A(8),C(18) 10250 DIM QR(11),QT(10),QZR(9),QZT(9) 10300 DIM FFD(2),FFP(5),NNT$(2),NND$(2),NNP$(5) 10350 DIM DU(10),B0KCL(15),B1KCL(15),CKCL(15),CPKCL(10),V0KCL(20) 10400 DIM DBDP(20),B0LKCL(15),B1LKCL(15),CLKCL(15),B0JKCL(15),B1JKCL(15),CJKCL(15) 10450 GOSUB *BLOCKDATA 10500 GOSUB *UNIT 10850 GOSUB *PARAMETERS 10900 LPRINT 10950 INPUT"Pressure? If end, input 0";X 11000 IF X=0 THEN GOTO 14600 11050 INPUT"Temperature";TT 11100 T=TT 11150 GOSUB *TTTT 11200 T=TTT 11300 RT=GASCON*T 11350 GOSUB *BBT 11400 INPUT "Molality";MOL 11410 IF MOL=<0 THEN GOTO 11400 12250 PRES=X 12300 PINPUT=PRES/FP 12350 DGSS=PINPUT/(T*.4#) 12360 IF T>=647.126# THEN GOTO 12700 12400 DLL=0 : DVV=0 12450 DL=DLL : DLIQ=0 12500 GOSUB *PCORRTPDLDV 12550 IF ABS((PINPUT-P)/P)=<5D-005 THEN PPP=PINPUT : GOTO 13550 12600 IF PINPUT>P THEN DGSS=DL:GOTO 12700 12650 IF PINPUT

1.9# THEN DD=1.9# 15050 D=DD 15100 GOSUB *QQTD 15150 Q0=Q 15200 GOSUB *BASEDT 15250 PP=RT*DD*BASEF+Q0 15300 DPD=RT*(Z+Y*DZB)+Q5 : DQ=DPD 15350 IF DPD>0 THEN GOTO 15550 15400 IF D>=.2967# THEN DD=DD*1.02 15450 IF D<.2967# THEN DD=DD*.98# 15500 IF LL=<10 GOTO 14900 15550 DPDX=DPD*1.1# 15600 IF DPDX<.1# THEN DPDX=.1# 15650 DP=ABS(1#-PP/PPP) 15700 IF DP<1D-009 THEN GOTO 16100 15750 IF D>.3# AND DP<1D-008 THEN GOTO 16100 15800 IF D>.7# AND DP<1D-007 THEN GOTO 16100 15850 XP=(PPP-PP)/DPDX 15900 IF ABS(XP)>.1# THEN XP=XP*.1#/ABS(XP) 15950 DD=DD+XP 16000 IF DD=<0 THEN DD=1D-008 16050 IF LL=<30 THEN GOTO 14900 16100 DOUT=DD 16150 RETURN 16200 *CORRTPDLDVDELG 16250 IF T>646.3# THEN GOTO 17150 16300 DLIQ=DLL 16350 IF DLL=<0 THEN DLIQ=1.11#-.0004*T 16400 DLL=DLIQ:D=DLIQ 16450 GOSUB *DFINDDOUTPDTDPD 16500 D=DOUT:DL=DOUT 16550 GOSUB *THERMDT 16600 GL=GD 16650 DVAP=DVV 16700 IF DVV=<0 THEN DVAP=PPP/RT 16750 D=DVAP:DVV=DVAP 16800 GOSUB *DFINDDOUTPDTDPD 16850 IF DOUT<5D-007 THEN DOUT=5D-007 16900 D=DOUT:DV=DOUT 16950 GOSUB *THERMDT 17000 GV=GD 17050 DELG=GL-GV 17100 RETURN 17150 PPP=0 17200 IF T>647.126# THEN RETURN 17250 DELG=0 17300 TAUC=.657128#*(1#-T/647.126#)^.325# 17350 DL=.322#+TAUC 17400 DV=.322#-TAUC 17450 D=DV 17500 GOSUB *BASEDT 17550 GOSUB *QQTD 17650 PPP=RT*DV*BASEF+Q 17700 RETURN 17750 *BBT 17800 BV(1)=1# 17850 FOR I=2 TO 10 17900 BV(I)=BV(I-1)*TZ/T 17950 NEXT I 18000 B1=BP(1)+BP(2)*LOG(1#/BV(2)) 18050 B2=BQ(1) 18100 B1T=BP(2)*BV(2)/TZ 18150 B2T=0 18200 B1TT=0 18250 B2TT=0 18300 FOR I=3 TO 10 18350 B1=B1+BP(I)*BV(I-1) 18400 B2=B2+BQ(I)*BV(I-1) 18450 B1T=B1T-CDBL((I-2))*BP(I)*BV(I-1)/T 18500 B2T=B2T-CDBL((I-2))*BQ(I)*BV(I-1)/T 18550 B1TT=B1TT+BP(I)*CDBL((I-2))*CDBL((I-2))*BV(I-1)/(T*T) 18600 B2TT=B2TT+BQ(I)*CDBL((I-2))*CDBL((I-2))*BV(I-1)/(T*T) 18650 NEXT I 18700 B1TT=B1TT-B1T/T 18750 B2TT=B2TT-B2T/T 18800 RETURN 18850 *BASEDT 18900 Y=.25#*B1*D 18950 XX=1#-Y 19000 Z0=(1#+ALPHA*Y+BETA*Y*Y)/(XX*XX*XX) 19050 Z=Z0+4#*Y*(B2/B1-GAMMA) 19100 DZ0=(ALPHA+2#*BETA*Y)/(XX*XX*XX)+3#*(1#+ALPHA*Y+BETA*Y*Y)/(XX*XX*XX*XX) 19150 DZB=DZ0+4#*(B2/B1-GAMMA) 19200 AB=(-1#)*LOG(XX)-(BETA-1#)/XX+28.16666667#/(XX*XX)+4#*Y*(B2/B1-GAMMA)+15.166666667#+LOG(D*RT/.101325) 19300 BASEF=Z 19350 BB2TT=T*T*B2TT 19400 UB=(-1#)*T*B1T*(Z-1#-D*B2)/B1-D*T*B2T 19500 CVB=2#*UB+(Z0-1#)*((T*B1T/B1)*(T*B1T/B1)-T*T*B1TT/B1)-D*(BB2TT-GAMMA*B1TT*T*T)-(T*B1T/B1)*(T*B1T/B1)*Y*DZ0 19550 DPDTB=BASEF/T+BASEF*D/Z*(DZB*B1T/4#+B2T-B2/B1*B1T) 19600 SB=UB-AB 19650 RETURN 19700 *QQTD 19750 QR(1)=0 19800 Q5=0 19850 Q=0 19900 AR=0 19950 DADT=0 20000 CVR=0 20050 DPDTR=0 20100 E=EXP((-1#)*AA*D) 20150 Q10=D*D*E 20200 Q20=1#-E 20250 QR(2)=Q10 20300 QV=TZ/T 20350 QT(1)=T/TZ 20400 FOR I=2 TO 10 20450 QR(I+1)=QR(I)*Q20 20500 QT(I)=QT(I-1)*QV 20550 NEXT I 20600 FOR I=1 TO INC 20650 K=II(I)+1 20700 L=JJ(I) 20750 QK=CDBL(K) : QL=CDBL(L) 20850 QZR(K-1)=QR(K+1):QZT(L)=QT(L+1):QZR(K)=QR(K+2):QZT(L+1)=QT(L+2) 20900 QP=HGKG(I)*AA*QZR(K-1)*QZT(L) 20950 Q=Q+QP 21000 Q5=Q5+AA*(2#/D-AA*(1#-E*(QK-1#)/Q20))*QP 21050 AR=AR+HGKG(I)*QZR(K)*QZT(L)/(Q10*QK*RT) 21100 DFDT=Q20^QK*(1#-QL)*QZT(L+1)/(TZ*QK) 21150 D2F=QL*DFDT 21200 DPT=DFDT*Q10*AA*QK/Q20 21250 DADT=DADT+HGKG(I)*DFDT 21300 DPDTR=DPDTR+HGKG(I)*DPT 21350 CVR=CVR+HGKG(I)*D2F/GASCON 21400 NEXT I 21450 QP=0 21500 Q2A=0 21550 FOR J=37 TO 40 21600 IF HGKG(J)=0 THEN GOTO 22850 21650 K=II(J) 21700 KM=JJ(J) 21750 QK=CDBL(K) : QKM=CDBL(KM) 21800 DDZ=ADZ(J-36) 21850 DEL=D/DDZ-1# 21900 IF ABS(DEL)<1D-010 THEN DEL=1D-010 22000 EX1=(-1#)*AAD(J-36)*DEL^QK 22050 DEX=EXP(EX1)*DEL^QKM 22100 ATT=AAT(J-36) 22150 TX=ATZ(J-36) 22200 TAU=T/TX-1# 22250 EX2=(-1#)*ATT*TAU*TAU 22300 TEX=EXP(EX2) 22350 Q10=DEX*TEX 22400 QM=QKM/DEL-QK*AAD(J-36)*DEL^(QK-1#) 22450 FCT=QM*D*D*Q10/DDZ 22500 Q5T=FCT*(2#/D+QM/DDZ)-(D/DDZ)*(D/DDZ)*Q10*(QKM/(DEL*DEL)+QK*(QK-1#)*AAD(J-36)*DEL^(QK-2#)) 22550 Q5=Q5+Q5T*HGKG(J) 22600 QP=QP+HGKG(J)*FCT 22650 DADT=DADT-2#*HGKG(J)*ATT*TAU*Q10/TX 22700 DPDTR=DPDTR-2#*HGKG(J)*ATT*TAU*FCT/TX 22750 Q2A=Q2A+T*HGKG(J)*(4#*ATT*EX2+2#*ATT)*Q10/(TX*TX) 22800 AR=AR+Q10*HGKG(J)/RT 22850 NEXT J 22900 SR=(-1#)*DADT/GASCON 22950 UR=AR+SR 23000 CVR=CVR+Q2A/GASCON 23050 Q=Q+QP 23100 RETURN 23150 *THERMDT 23200 GOSUB *IDEALT 23250 GOSUB *BASEDT 23300 GOSUB *QQTD 23350 QPQ=Q:QDPQ=Q5 23400 Z=BASEF+QPQ/(RT*D) 23450 DPDD=RT*(BASEF+Y*DZB)+QDPQ 23500 AD=AB+AR+AI-UREF/T+SREF 23550 GD=AD+Z 23600 UD=UB+UR+UI-UREF/T 23650 DPDT=RT*D*DPDTB+DPDTR 23700 CVDX=CVB+CVR+CVIX 23750 CPD=CVDX+T*DPDT*DPDT/(D*D*DPDD*GASCON) 23800 HD=UD+Z 23850 SD=SB+SR+SI-SREF 23900 RETURN 23950 *SECDERIVP 24100 D2PDDDD1=0 : D2PDD2=0 : D2PDTDD=0 : D2PDTDT1=0 : D2PDT2=0 : ZX=0 24150 D2PRESIDDDDD1=0 : D2PRESIDDDDT=0 : D2PRESIDDT2=0 : D2PRESIDDD2=0 24160 D2DDT2A=0 : D2DDT2=0 : DDDT=0 24250 D2PDDDD1=3#+ALPHA+3#*Y+4#*ALPHA*Y+3#*BETA*Y+ALPHA*Y*Y+3#*BETA*Y*Y 24300 D2PDDDD1=D2PDDDD1*B1/(2#*XX*XX*XX*XX*XX)+2#*B1*(B2/B1-GAMMA) 24350 D2PDD2=D2PDDDD1*RT 24400 D2PDTDD=Z0+8#*Y*(B2/B1-GAMMA)+(Y+B1T*D*T/2#)*((ALPHA+2#*BETA*Y)/(XX*XX*XX)+3#*Z0/XX) 24450 ZX=6#+3#*ALPHA+BETA+3#*ALPHA*Y+4#*BETA*Y+BETA*Y*Y 24500 D2PDTDD=D2PDTDD+B1T*D*T*Y*ZX/(2#*XX*XX*XX*XX*XX) 24550 D2PDTDD=D2PDTDD+2#*B2T*D*T-2#*B1T*D*T*GAMMA 24600 D2PDTDD=D2PDTDD*GASCON 24650 D2PDTDT1=(ALPHA+2#*BETA*Y)/(XX*XX*XX)+3#*Z0/XX 24700 D2PDTDT1=D2PDTDT1*(2#*B1T*D+B1TT*D*T)/4# 24750 D2PDTDT1=D2PDTDT1+B1T*B1T*D*D*T*ZX/(8#*XX*XX*XX*XX*XX)+2#*B2T*D 24800 D2PDTDT1=D2PDTDT1-2#*B1T*D*GAMMA+B2TT*D*T-B1TT*D*T*GAMMA 24850 D2PDT2=D2PDTDT1*GASCON*D 24900 FOR I=1 TO INC 24950 K=II(I)+1 25000 L=JJ(I) 25050 QK=CDBL(K) : QL=CDBL(L) 25100 D2PRESIDDDDD1=2#/(D*D)-4#/D+4#*(QK-1#)*E/(D*Q20)+1#-3#*(QK-1#)*E/Q20+(QK-1#)*(QK-2#)*E*E/(Q20*Q20) 25150 D2PRESIDDD2=D2PRESIDDD2+D2PRESIDDDDD1*HGKG(I)*QT(L+1)*QR(K+1) 25200 D2PRESIDDDDT=D2PRESIDDDDT-(QL-1#)*HGKG(I)*QT(L+1)*QR(K+1)*(2#/D-1#+(QK-1#)*E/Q20)/T 25250 D2PRESIDDT2=D2PRESIDDT2+(QL-1#)*QL*HGKG(I)*QR(K+1)*QT(L+1)/(T*T) 25300 NEXT I 25350 D2PRESIDDD2A=0 : D2PRESIDDDDTC=0 25400 FOR J=37 TO 40 25410 IF HGKG(J)=0 THEN GOTO 26750 25450 K=II(J) 25500 KM=JJ(J) 25550 QK=CDBL(K) : QKM=CDBL(KM) 25600 DDZ=ADZ(J-36) 25650 DEL=D/DDZ-1# 25700 IF ABS(DEL)<1D-010 THEN DEL=1D-010 25750 EX1=(-1#)*AAD(J-36)*DEL^QK 25800 DEX=EXP(EX1)*DEL^QKM 25850 ATT=AAT(J-36) 25900 TX=ATZ(J-36) 25950 TAU=T/TX-1# 26000 EX2=(-1#)*ATT*TAU*TAU 26050 TEX=EXP(EX2) 26100 Q10=DEX*TEX 26150 QM=QKM/DEL-QK*AAD(J-36)*DEL^(QK-1#) 26200 D2PRESIDDD2A=QM*(2#/(D*D)+4#*QKM/(D*DDZ*DEL)+4#*QK*EX1/(D*DDZ*DEL)) 26250 D2PRESIDDD2B=QM*(QKM*(QKM-1#)+2#*QK*QKM*EX1+QK*(QK-1#)*EX1+QK*QK*EX1*EX1) 26300 D2PRESIDDD2A=D2PRESIDDD2A/DDZ+D2PRESIDDD2B/(DDZ*DDZ*DDZ*DEL*DEL) 26350 D2PRESIDDD2C=(QKM-QK*(QK-1#)*EX1)*(4#/D+2*QKM/(DDZ*DEL)+2#*QK*EX1/(DDZ*DEL))/(DDZ*DDZ*DEL*DEL) 26400 D2PRESIDDD2A=D2PRESIDDD2A-D2PRESIDDD2C 26450 D2PRESIDDD2A=D2PRESIDDD2A+(2#*QKM+QK*(QK-1#)*(QK-2#)*EX1)/(DDZ*DDZ*DDZ*DEL*DEL*DEL) 26500 D2PRESIDDD2=D2PRESIDDD2+HGKG(J)*Q10*D2PRESIDDD2A*D*D 26550 D2PRESIDDDDTC=2#*QM+D*QKM*QM/(DDZ*DEL)+D*QK*EX1*QM/(DDZ*DEL) 26600 D2PRESIDDDDTC=D2PRESIDDDDTC+D*(QK*(QK-1#)*EX1/(DEL*DEL)-QKM/(DEL*DEL))/DDZ 26650 D2PRESIDDDDT=D2PRESIDDDDT-2#*D*HGKG(J)*ATT*TAU*Q10*D2PRESIDDDDTC/(TX*DDZ) 26700 D2PRESIDDT2=D2PRESIDDT2-2#*D*D*HGKG(J)*ATT*(1#+2#*EX2)*Q10*QM/(DDZ*TX*TX) 26750 NEXT J 26800 D2PDD2=D2PDD2+D2PRESIDDD2 26850 D2PDTDD=D2PDTDD+D2PRESIDDDDT 26900 D2PDT2=D2PDT2+D2PRESIDDT2 26950 D2DDT2A=DPDD*DPDD*D2PDT2-2#*DPDT*DPDD*D2PDTDD+DPDT*DPDT*D2PDD2 27000 D2DDT2=(-1#)*D2DDT2A/(DPDD*DPDD*DPDD) 27050 DDDT=(-1#)*DPDT/DPDD 27100 RETURN 27150 *PST 27200 IF T>314# THEN GOTO 27400 27250 PL=6.3573118#-8858.843#/T+607.56335*T^(-.6#) 27300 PS=.1#*EXP(PL) 27350 RETURN 27400 TR=T/647.25# 27450 W=ABS(1#-TR) 27500 BPST=0 27550 FOR I=1 TO 8 27600 ZPST=CDBL(I) 27650 BPST=BPST+A(I)*W^((ZPST+1#)/2#) 27700 NEXT I 27750 QPST=BPST/TR 27800 PS=22.093*EXP(QPST) 27850 RETURN 27900 *IDEALT 27950 TIDEAL=T/100 28000 TL=LOG(TIDEAL) 28050 GI=(-1#)*(C(1)/TIDEAL+C(2))*TL 28100 HI=(C(2)+C(1)*(1#-TL)/TIDEAL) 28150 CPI=C(2)-C(1)/TIDEAL 28200 FOR I=3 TO 18 28250 GI=GI-C(I)*TIDEAL^CDBL(I-6) 28300 HI=HI+C(I)*CDBL((I-6))*TIDEAL^CDBL(I-6) 28350 CPI=CPI+C(I)*CDBL((I-6))*CDBL((I-5))*TIDEAL^CDBL(I-6) 28400 NEXT I 28450 AI=GI-1# 28500 UI=HI-1# 28550 CVIX=CPI-1# 28600 SI=UI-AI 28650 RETURN 28700 *PCORRTPDLDV 28750 GOSUB *PST 28800 PPP=PS 28850 GOSUB *CORRTPDLDVDELG 28900 DP=0 28950 DP=DELG*RT/(1#/DV-1#/DL) 29000 PPP=PPP+DP 29050 IF ABS(DELG)<1D-005 THEN GOTO 29150 29100 DLL=DL:DVV=DV:GOTO 28850 29150 P=PPP 29200 RETURN 29250 *UNIT 29300 PRINT"*******************" 29350 PRINT"* Enter units *" 29400 PRINT"*******************" 29450 PRINT A1$ 29500 PRINT"Choose from 1=deg K, 2=deg C" 29550 INPUT IT 29600 IF IT<1 OR IT>2 THEN GOTO 29500 29650 NT$=NNT$(IT) 29700 PRINT A2$ 29750 PRINT"Choose from 1=kg/m3, 2=g/cm3" 29800 INPUT ID 29850 IF ID>2 OR ID<1 THEN GOTO 29750 29900 ND$=NND$(ID) 29950 FD=FFD(ID) 30000 PRINT A3$ 30050 PRINT"Choose from 1=MPa, 2=bar " 30100 INPUT IP 30150 IF IP>2 OR IP<1 THEN GOTO 30050 30200 NP$=NNP$(IP) 30250 FP=FFP(IP) 30600 RETURN 30650 *TTTT 30700 ON IT GOTO 30750, 30900 30750 TTT=T 30850 GOTO 31000 30900 TTT=T+273.15# 31000 RETURN 31050 *BLOCKDATA 31100 FOR I=1 TO 4:READ ATZ(I):NEXT I 31150 DATA 640#,640#,641.6#,270# 31200 FOR I=1 TO 4:READ ADZ(I):NEXT I 31250 DATA 0.319#,0.319#,0.319#,1.55# 31300 FOR I=1 TO 4:READ AAT(I):NEXT I 31350 DATA 2.0D+004,2.0D+004,4.0D+004,25.0# 31400 FOR I=1 TO 4:READ AAD(I):NEXT I 31450 DATA 34.0#,40.0#,30.0#,1.05D+003 31500 GASCON=.461522# : TZ=647.073 : AA=1# : INC=36 31550 UREF=-4328.454977# : SREF=7.618072# 31600 ALPHA=11#:BETA=44.333333333333#:GAMMA=3.5# 31650 FOR I=1 TO 10:READ BP(I):NEXT I 31700 DATA 0.7478629#,-0.3540782#,0.0#,0.0#,0.007159876#,0.0#,-0.003528426#,0.0#,0.0#,0.0# 31750 FOR I=1 TO 10:READ BQ(I):NEXT I 31800 DATA 1.1278334#,0.0#,-0.5944001#,-5.010996#,0.0#,0.63684256#,0.0#,0.0#,0.0#,0.0# 31850 FOR I=1 TO 40:READ HGKG(I):NEXT I 31900 DATA -5.3062968529023D+002,2.2744901424408D+003,7.8779333020687D+002 31950 DATA -6.9830527374994D+001,1.7863832875422D+004,-3.9514731563338D+004 32000 DATA 3.3803884280753D+004,-1.3855050202703D+004,-2.5637436613260D+005 32050 DATA 4.8212575981415D+005,-3.4183016969660D+005,1.2223156417448D+005 32100 DATA 1.1797433655832D+006,-2.1734810110373D+006,1.0829952168620D+006 32150 DATA -2.5441998064049D+005,-3.1377774947767D+006,5.2911910757704D+006 32200 DATA -1.3802577177877D+006,-2.5109914369001D+005,4.6561826115608D+006 32250 DATA -7.2752773275387D+006,4.1774246148294D+005,1.4016358244614D+006 32300 DATA -3.1555231392127D+006,4.7929666384584D+006,4.0912664781209D+005 32350 DATA -1.3626369388386D+006,6.9625220862664D+005,-1.0834900096447D+006 32400 DATA -2.2722827401688D+005,3.8365486000660D+005,6.8833257944332D+003 32450 DATA 2.1757245522644D+004,-2.6627944829770D+003,-7.0730418082074D+004 32500 DATA -0.225#,-1.68#,0.055#,-93.0# 32550 FOR I=1 TO 40:READ II(I):NEXT I 32600 DATA 0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,6,8,8,8,8,2,2,0,4,2,2,2,4 32650 FOR I=1 TO 40:READ JJ(I):NEXT I 32700 DATA 2,3,5,7,2,3,5,7,2,3,5,7,2,3,5,7,2,3,5,7,2,3,5,7,2,3,5,7,2,3,5,7,1,4,4,4,0,2,0,0 32750 FOR I=1 TO 8:READ A(I):NEXT I 32800 DATA -7.8889166#,2.5514255#,-6.716169#,33.239495# 32850 DATA -105.38479#,174.35319#,-148.39348#,48.631602# 32900 FOR I=1 TO 18:READ C(I):NEXT I 32950 DATA 1.9730271018D+001,2.09662681977D+001,-4.83429455355D-001,6.05743189245D+000 33000 DATA 2.256023885D+001,-9.87532442D+000,-4.3135538513D+000,4.58155781D-001 33050 DATA -4.7754901883D-002,4.1238460633D-003,-2.7929052852D-004 33100 DATA 1.4481695261D-005,-5.6473658748D-007,1.6200446D-008,-3.303822796D-010 33150 DATA 4.51916067368D-012,-3.70734122708D-014,1.37546068238D-016 33200 FOR I=1 TO 2:READ FFD(I):NEXT I 33250 DATA 1.0D-003,1.0# 33300 FOR I=1 TO 2:READ FFP(I):NEXT I 33350 DATA 1.0#,10.0# 33500 FOR I=1 TO 2:READ NNT$(I):NEXT I 33550 DATA "K","deg C" 33600 FOR I=1 TO 2:READ NND$(I):NEXT I 33650 DATA "kg/m3","g/cm3" 33700 FOR I=1 TO 2:READ NNP$(I):NEXT I 33750 DATA "MPa","bar" 33900 A1$="TEMPERATURE":A2$="DENSITY":A3$="PRESSURE" 33950 RETURN 34000 *PARAMETERS 34050 REM B0,B1,C 34100 FOR I=1 TO 11 : READ B0KCL(I) : READ B1KCL(I) : READ CKCL(I) : NEXT I 34150 DATA -2.10289D-002, 2.20813D-001, 0.0# 34200 DATA 6.03967D-001, -4.61849#, 7.64891D-004 34250 DATA 3.67768D-003, -4.10116D-002, 0.0# 34300 DATA -7.05537D-006, 1.10445D-004, -1.12131D-008 34350 DATA 1.97968D-009, -4.73196D-008, 1.72256D-011 34400 DATA -2.47588D-003, -2.74120D-002, 0.0# 34450 DATA 1.44160D-001, 3.32883D-001, -5.71188D-003 34500 REM K1 34550 DATA -2931.268116#, 6353.355434#, 28.172180# 34600 REM fL(Tr,Pr) 34650 DATA 6.56838D-004, 9.67854D-004, -4.12364D-005 34700 REM K2 34750 DATA -33.953143#, 193.004059#, -0.125567# 34800 REM fG(Tr,Pr) 34850 DATA 5.0038D-002, 2.18752D-001, -3.94D-004 34900 REM B0L,B1L,CL 34950 FOR I=1 TO 9 : B0LKCL(I)=B0KCL(I) : B1LKCL(I)=B1KCL(I) : CLKCL(I)=CKCL(I) : NEXT I 35000 REM B0J,B1J,CJ 35050 FOR I=1 TO 7 : B0JKCL(I)=B0KCL(I) : B1JKCL(I)=B1KCL(I) : CJKCL(I)=CKCL(I) : NEXT I 35100 REM BV(KCl) 35150 FOR I=1 TO 15 : READ DBDP(I) : NEXT I 35200 DATA 0.0#, 0.0#, 9.45015D-008, -2.90741D-010, 3.26205D-003, 8.39662D-007, 0.0#, -4.41638D-009, 6.71235D-012 35250 DATA -4.42327D-005, -7.97437D-010, 0.0#, 4.12771D-012, -6.24996D-015, 4.16221D-008 35300 FOR I=1 TO 15 : READ V0KCL(I) : NEXT I 35350 DATA 1.56152D+003, -1.69234D+005, -4.29918#, 4.59233D-003, -3.25686D+004, -6.86887# 35400 DATA 7.35220D+002, 2.02245D-002, -2.15779D-005, 1.03212D+002, 5.34941D-003 35450 DATA -5.73121D-001, -1.57862D-005, 1.66987D-008, -7.22012D-002 35500 REM CP(KCl) 35550 FOR I=1 TO 7 : READ CPKCL(I) : NEXT I 35600 DATA 3.71110D+004, 0.0#, -7.90247D+003, 3.30367D+001, -1.76733D-002, -2.91950D+004, -5.92362D+006 35650 KCL=74.551# : TTR=298.15# : PR=179# : ATM=1.01325 35700 RVGAS=83.1441# : RGAS=8.31441# : MW=18.0152 35750 FOR I=1 TO 9: READ DU(I) : NEXT I 35800 DATA 3.4279D+002,-5.0866D-003,9.4690D-007,-2.0525#,3.1159D+003 35850 DATA -1.8289D+002,-8.0325D+003,4.2142D+006,2.1417# 35900 EE=4.803242D-010 : BC=1.380662D-016 35950 RETURN 36000 *DEBYEHUCKEL 36050 PRES=(PRES/FP)*(FFP(2)/FFP(1)) 36150 DPDD=DPDD*(FFP(2)/FFP(1)) 36200 DPDT=DPDT*(FFP(2)/FFP(1)) 36350 EPS=DU(1)*EXP(DU(2)*T+DU(3)*T*T) 36400 E=1#+(PRES-1000)/(DU(7)+DU(8)/T+DU(9)*T+1000) 36450 E=LOG(E) 36500 EPS=EPS+(DU(4)+(DU(5)/(DU(6)+T)))*E 36550 APHI=SQR(2#*3.14159265#*6.022045D+023*D/1000)/3# 36600 APHI=APHI*EE*EE*EE/(BC*SQR(BC)*T*SQR(T)*EPS*SQR(EPS)) 36650 DET=DU(7)+DU(8)/T+DU(9)*T+PRES 36700 DEPSDP=DU(4)+DU(5)/(DU(6)+T) 36750 DEPSDP=DEPSDP/DET 36800 DRHODPDD=1#/(D*DPDD) 36850 AV=2#*RVGAS*T*APHI*(3#*DEPSDP/EPS-DRHODPDD) 36900 ALPH=DPDT/(D*DPDD) 36950 DE=DU(7)+DU(8)/T+DU(9)*T 37000 DEPS=DU(1)*(DU(2)+2#*DU(3)*T)*EXP(DU(2)*T+DU(3)*T*T) 37050 DEPS=DEPS-DU(5)*LOG(1#+(PRES-1000)/(DE+1000))/((DU(6)+T)*(DU(6)+T)) 37100 DEPS=DEPS+(DU(4)+DU(5)/(DU(6)+T))*(1000-PRES)*(DU(9)-DU(8)/(T*T))/((DE+PRES)*(DE+1000)) 37150 DEPS=DEPS/EPS 37200 AH=1#+T*DEPS+T*ALPH/3# 37250 AH=AH*(-6#)*APHI*RGAS*T 37300 D2EPS=DU(1)*(DU(2)+2#*DU(3)*T)*(DU(2)+2#*DU(3)*T)*EXP(DU(2)*T+DU(3)*T*T) 37350 D2EPS=D2EPS+2#*DU(1)*DU(3)*EXP(DU(2)*T+DU(3)*T*T) 37400 D2EPS=D2EPS+(2#*DU(5)/((DU(6)+T)*(DU(6)+T)*(DU(6)+T)))*LOG(1#+(PRES-1000)/(DE+1000)) 37450 D2EPS=D2EPS-(2#*DU(5)/((DU(6)+T)*(DU(6)+T)))*(DU(9)-DU(8)/(T*T))*(1#/(DE+PRES)-1#/(DE+1000)) 37500 D2EPS=D2EPS+(DU(4)+DU(5)/(DU(6)+T))*(2#*DU(8)/(T*T*T))*(1#/(DE+PRES)-1#/(DE+1000)) 37550 D2EPSX=(DU(9)-DU(8)/(T*T))*(DU(9)-DU(8)/(T*T))*(1#/((DE+PRES)*(DE+PRES))-1#/((DE+1000)*(DE+1000))) 37600 D2EPS=D2EPS-(DU(4)+DU(5)/(DU(6)+T))*D2EPSX 37650 D2EPS=D2EPS/EPS 37700 DWDDDT=DDDT/D 37750 DWD2DDT2=D2DDT2/D 37800 AJ=2#*DWD2DDT2-DWDDDT*DWDDDT-2#*DWDDDT/T-6#*D2EPS+15#*DEPS*DEPS+6#*DEPS/T 37850 AJ=AJ-6#*DWDDDT*DEPS+3#/(T*T) 37900 AJ=AJ*APHI*RGAS*T*T 37950 RETURN 38000 *KCLAQ 38050 MI=MOL 38051 S=S*MW : H=H*MW : G=G*MW 38052 S=S+SREF*RGAS : H=H+UREF*RGAS : G=G+UREF*RGAS-T*SREF*RGAS 38060 CPW=CPD*RGAS 38400 TL=T-227# 38450 TH=647#-T 38500 HF=LOG(1#+1.2#*SQR(MI))/(2#*1.2#) 38550 BVKCL=DBDP(1)+DBDP(2)/T+DBDP(3)*T+DBDP(4)*T*T+DBDP(5)/TH 38600 BVKCL=BVKCL+(DBDP(6)+DBDP(7)/T+DBDP(8)*T+DBDP(9)*T*T+DBDP(10)/TH)*PRES 38650 BVKCL=BVKCL+(DBDP(11)+DBDP(12)/T+DBDP(13)*T+DBDP(14)*T*T+DBDP(15)/TH)*PRES*PRES 38700 VKCL=V0KCL(1)+V0KCL(2)/T+V0KCL(3)*T+V0KCL(4)*T*T+V0KCL(5)/TH 38750 VKCL=VKCL+(V0KCL(6)+V0KCL(7)/T+V0KCL(8)*T+V0KCL(9)*T*T+V0KCL(10)/TH)*PRES 38800 VKCL=VKCL+(V0KCL(11)+V0KCL(12)/T+V0KCL(13)*T+V0KCL(14)*T*T+V0KCL(15)/TH)*PRES*PRES 38850 VPHI=VKCL+2#*AV*HF+2#*RVGAS*T*MOL*BVKCL 38900 VSOLN=VPHI*MOL+(1000/D) 38950 DSOLN=(1000+MOL*KCL)/VSOLN 39000 LPRINT 39050 LPRINT USING"V(water)=+##.### VKCl=+###.##";MW/D,VKCL 39100 LPRINT 39150 REM 179bar 39200 HKCLTRPR=(V0KCL(1)+2#*V0KCL(2)/TTR-V0KCL(4)*TTR*TTR)*(PR-ATM) 39250 HKCLTRPR=HKCLTRPR+((647#-2#*TTR)*V0KCL(5)/((647#-TTR)*(647#-TTR)))*(PR-ATM) 39300 HKCLTRPR=HKCLTRPR+(V0KCL(6)+2#*V0KCL(7)/TTR-V0KCL(9)*TTR*TTR)*(PR*PR-ATM*ATM)/2# 39350 HKCLTRPR=HKCLTRPR+((647#-2#*TTR)*V0KCL(10)/((647#-TTR)*(647#-TTR)))*(PR*PR-ATM*ATM)/2# 39400 HKCLTRPR=HKCLTRPR+(V0KCL(11)+2#*V0KCL(12)/TTR-V0KCL(14)*TTR*TTR)*(PR*PR*PR-ATM*ATM*ATM)/3# 39450 HKCLTRPR=HKCLTRPR+((647#-2#*TTR)*V0KCL(15)/((647#-TTR)*(647#-TTR)))*(PR*PR*PR-ATM*ATM*ATM)/3# 39500 HKCLTRPR=HKCLTRPR/10 39550 HKCLPR=HKCLTRPR+CPKCL(1)*(T-TTR)+CPKCL(2)*(LOG(T)-LOG(TTR)) 39600 HKCLPR=HKCLPR+CPKCL(3)*(T*LOG(T)-T-TTR*LOG(TTR)+TTR)+CPKCL(4)*(T*T-TTR*TTR)/2# 39650 HKCLPR=HKCLPR+CPKCL(5)*(T*T*T-TTR*TTR*TTR)/3#+CPKCL(6)*(LOG(TL)-LOG(TTR-227#)) 39700 HKCLPR=HKCLPR+CPKCL(7)*(1#/TH-1#/(647#-TTR)) 39750 DHDPP=(V0KCL(1)+2#*V0KCL(2)/T-V0KCL(4)*T*T+(647#-2#*T)*V0KCL(5)/(TH*TH))*(PRES-PR) 39800 DHDPP=DHDPP+(V0KCL(6)+2#*V0KCL(7)/T-V0KCL(9)*T*T+(647#-2#*T)*V0KCL(10)/(TH*TH))*(PRES*PRES-PR*PR)/2# 39850 DHDPP=DHDPP+(V0KCL(11)+2#*V0KCL(12)/T)*(PRES*PRES*PRES-PR*PR*PR)/3# 39900 DHDPP=DHDPP-(V0KCL(14)*T*T-(647#-2#*T)*V0KCL(15)/(TH*TH))*(PRES*PRES*PRES-PR*PR*PR)/3# 39950 DHDPP=DHDPP/10 40000 HKCL=HKCLPR+DHDPP 40050 REM SKCl(298.15K,0.1MPa)=157.9349 after Archer 40150 DS0KCL=-((-1#)*V0KCL(2)/(TTR*TTR)+V0KCL(3)+2#*V0KCL(4)*TTR+V0KCL(5)/((647#-TTR)*(647#-TTR)))*(ATM-1#) 40200 DS0KCL=DS0KCL-((-1#)*V0KCL(7)/(TTR*TTR)+V0KCL(8)+2#*V0KCL(9)*TTR)*(ATM*ATM-1#*1#)/2# 40250 DS0KCL=DS0KCL-(V0KCL(10)/((647#-TTR)*(647#-TTR)))*(ATM*ATM-1#*1#)/2# 40300 DS0KCL=DS0KCL-((-1#)*V0KCL(12)/(TTR*TTR)+V0KCL(13))*(ATM*ATM*ATM-1#*1#*1#)/3# 40350 DS0KCL=DS0KCL-(2#*V0KCL(14)*TTR)*(ATM*ATM*ATM-1#*1#*1#)/3# 40400 DS0KCL=DS0KCL-(V0KCL(15)/((647#-TTR)*(647#-TTR)))*(ATM*ATM*ATM-1#*1#*1#)/3# 40450 DS0KCL=DS0KCL/10 40500 S0KCL=157.9349#+DS0KCL 40550 DSKCLTRPR=-((-1#)*V0KCL(2)/(TTR*TTR)+V0KCL(3)+2#*V0KCL(4)*TTR+V0KCL(5)/((647#-TTR)*(647#-TTR)))*(PR-ATM) 40600 DSKCLTRPR=DSKCLTRPR-((-1#)*V0KCL(7)/(TTR*TTR)+V0KCL(8))*(PR*PR-ATM*ATM)/2# 40650 DSKCLTRPR=DSKCLTRPR-(2#*V0KCL(9)*TTR+V0KCL(10)/((647#-TTR)*(647#-TTR)))*(PR*PR-ATM*ATM)/2# 40700 DSKCLTRPR=DSKCLTRPR-((-1#)*V0KCL(12)/(TTR*TTR)+V0KCL(13))*(PR*PR*PR-ATM*ATM*ATM)/3# 40750 DSKCLTRPR=DSKCLTRPR-(2#*V0KCL(14)*TTR)*(PR*PR*PR-ATM*ATM*ATM)/3# 40800 DSKCLTRPR=DSKCLTRPR-(V0KCL(15)/((647#-TTR)*(647#-TTR)))*(PR*PR*PR-ATM*ATM*ATM)/3# 40850 DSKCLTRPR=DSKCLTRPR/10 40900 SKCLTRPR=S0KCL+DSKCLTRPR 40950 SKCLPR=SKCLTRPR+CPKCL(1)*LOG(T/TTR)-CPKCL(2)*(1#/T-1#/TTR) 41000 SKCLPR=SKCLPR+CPKCL(3)*(LOG(T)*LOG(T)-LOG(TTR)*LOG(TTR))/2# 41050 SKCLPR=SKCLPR+CPKCL(4)*(T-TTR)+CPKCL(5)*(T*T-TTR*TTR)/2#-(CPKCL(6)/227#)*(LOG(T*(TTR-227#)/(TTR*TL))) 41100 SKCLPRX=LOG(T/TTR)+(1294#-T)/TH-LOG(TH)-(1294#-TTR)/(647#-TTR)+LOG(647#-TTR) 41150 SKCLPR=SKCLPR+(CPKCL(7)/(647#*647#))*SKCLPRX 41200 DSKCL=-((-1#)*V0KCL(2)/(T*T)+V0KCL(3)+2#*V0KCL(4)*T+V0KCL(5)/(TH*TH))*(PRES-PR) 41250 DSKCL=DSKCL-((-1#)*V0KCL(7)/(T*T)+V0KCL(8)+2#*V0KCL(9)*T+V0KCL(10)/(TH*TH))*(PRES*PRES-PR*PR)/2# 41300 DSKCL=DSKCL-((-1#)*V0KCL(12)/(T*T)+V0KCL(13))*(PRES*PRES*PRES-PR*PR*PR)/3# 41350 DSKCL=DSKCL-(2#*V0KCL(14)*T+V0KCL(15)/(TH*TH))*(PRES*PRES*PRES-PR*PR*PR)/3# 41400 DSKCL=DSKCL/10 41450 SKCL=SKCLPR+DSKCL 41500 GKCL=HKCL-T*SKCL 41550 CP0=CPKCL(1)+CPKCL(2)/T+CPKCL(3)*LOG(T)+CPKCL(4)*T+CPKCL(5)*T*T+CPKCL(6)/TL+CPKCL(7)/(TH*TH) 41600 TD2VDT2=(2#*V0KCL(2)/(T*T)+2#*V0KCL(4)*T+2#*V0KCL(5)*T/(TH*TH*TH))*(PRES-PR) 41650 TD2VDT2=TD2VDT2+(2#*V0KCL(7)/(T*T)+2#*V0KCL(9)*T+2#*V0KCL(10)*T/(TH*TH*TH))*(PRES*PRES-PR*PR)/2# 41700 TD2VDT2=TD2VDT2+(2#*V0KCL(12)/(T*T)+2#*V0KCL(14)*T)*(PRES*PRES*PRES-PR*PR*PR)/3# 41750 TD2VDT2=TD2VDT2+(2#*V0KCL(15)*T/(TH*TH*TH))*(PRES*PRES*PRES-PR*PR*PR)/3# 41800 TD2VDT2=TD2VDT2/10 41850 CPK=CP0-TD2VDT2 41950 LPRINT USING"G/RT= +##.#### GKCl/RT= +##.###";G/(RGAS*T),GKCL/(RGAS*T) 42000 LPRINT USING"H/RT= +##.#### HKCl/RT= +##.###";H/(RGAS*T),HKCL/(RGAS*T) 42050 LPRINT USING"S/R= +##.#### SKCl/R= +##.###";S/RGAS,SKCL/RGAS 42100 LPRINT USING"Cp/R= +##.### CpKCl/R= +###.##";CPW/RGAS,CPK/RGAS 42150 LPRINT 42250 BM0=B0KCL(1)*T*T/6#+B0KCL(2)*T/2#+B0KCL(3)*T*T*(LOG(T)/2#-5#/12#)/3# 42300 BM0=BM0+B0KCL(4)*T*T*T/12#+B0KCL(5)*T*T*T*T/20 42350 BM0=BM0+B0KCL(6)*(T/2#+3#*227#*227#/(2#*T)+227#*TL*LOG(TL)/T)-B0KCL(7)*(2#*TH*LOG(TH)/T+LOG(TH)) 42400 BM0=BM0-B0KCL(8)/T-B0KCL(9)*TTR*TTR/T+B0KCL(10)+B0KCL(11) 42450 BM1=B1KCL(1)*T*T/6#+B1KCL(2)*T/2#+B1KCL(3)*T*T*(LOG(T)/2#-5#/12#)/3# 42500 BM1=BM1+B1KCL(4)*T*T*T/12#+B1KCL(5)*T*T*T*T/20 42550 BM1=BM1+B1KCL(6)*(T/2#+3#*227#*227#/(2#*T)+227#*TL*LOG(TL)/T)-B1KCL(7)*(2#*TH*LOG(TH)/T+LOG(TH)) 42600 BM1=BM1-B1KCL(8)/T-B1KCL(9)*TTR*TTR/T+B1KCL(10)+B1KCL(11) 42650 CM=CKCL(1)*T*T/6#+CKCL(2)*T/2#+CKCL(3)*T*T*(LOG(T)/2#-5#/12#)/3#+CKCL(4)*T*T*T/12#+CKCL(5)*T*T*T*T/20 42700 CM=CM+CKCL(6)*(T/2#+3#*227#*227#/(2#*T)+227#*TL*LOG(TL)/T)-CKCL(7)*(2#*TH*LOG(TH)/T+LOG(TH)) 42750 CM=CM-CKCL(8)/T-CKCL(9)*TTR*TTR/T+CKCL(10)+CKCL(11) 42800 B0L=B0LKCL(1)*T/3#+B0LKCL(2)/2#+B0LKCL(3)*T*(LOG(T)-1#/3#)/3#+B0LKCL(4)*T*T/4#+B0LKCL(5)*T*T*T/5# 42850 B0L=B0L+(B0LKCL(6)/(T*T))*(TL*TL/2#+454#*TL+227#*227#*LOG(TL)) 42900 B0L=B0L+(B0LKCL(7)/(T*T))*((-1#)*TH+1294#*LOG(TH)+647#*647#/TH) 42950 B0L=B0L+B0LKCL(8)/(T*T)+B0LKCL(9)*TTR*TTR/(T*T) 43000 B1L=B1LKCL(1)*T/3#+B1LKCL(2)/2#+B1LKCL(3)*T*(LOG(T)-1#/3#)/3#+B1LKCL(4)*T*T/4#+B1LKCL(5)*T*T*T/5# 43050 B1L=B1L+(B1LKCL(6)/(T*T))*(TL*TL/2#+454#*TL+227#*227#*LOG(TL)) 43100 B1L=B1L+(B1LKCL(7)/(T*T))*((-1#)*TH+1294#*LOG(TH)+647#*647#/TH) 43150 B1L=B1L+B1LKCL(8)/(T*T)+B1LKCL(9)*TTR*TTR/(T*T) 43200 CL=CLKCL(1)*T/3#+CLKCL(2)/2#+CLKCL(3)*T*(LOG(T)-1#/3#)/3#+CLKCL(4)*T*T/4#+CLKCL(5)*T*T*T/5# 43250 CL=CL+(CLKCL(6)/(T*T))*(TL*TL/2#+454#*TL+227#*227#*LOG(TL)) 43300 CL=CL+(CLKCL(7)/(T*T))*((-1#)*TH+1294#*LOG(TH)+647#*647#/TH) 43350 CL=CL+CLKCL(8)/(T*T)+CLKCL(9)*TTR*TTR/(T*T) 43400 B0J=B0JKCL(1)+B0JKCL(2)/T+B0JKCL(3)*LOG(T)+B0JKCL(4)*T+B0JKCL(5)*T*T+B0JKCL(6)/TL+B0JKCL(7)/(TH*TH) 43450 B1J=B1JKCL(1)+B1JKCL(2)/T+B1JKCL(3)*LOG(T)+B1JKCL(4)*T+B1JKCL(5)*T*T+B1JKCL(6)/TL+B1JKCL(7)/(TH*TH) 43500 CJ=CJKCL(1)+CJKCL(2)/T+CJKCL(3)*LOG(T)+CJKCL(4)*T+CJKCL(5)*T*T+CJKCL(6)/TL+CJKCL(7)/(TH*TH) 43550 OSC=0 43600 GM=0 43700 EAI=EXP((-2#)*SQR(MI)) 43850 GM1=1#-(1#+2#*SQR(MI)-4#*MI/2#)*EAI 44000 BL=B0L+2#*B1L*(1#-(1#+2#*SQR(MI))*EAI)/(4#*MI) 44150 BJ=B0J+B1J*(2#/(4#*MI))*(1#-(1#+2#*SQR(MI))*EAI) 44250 REM Pressure dependence 44300 BITG=(DBDP(1)+DBDP(2)/T+DBDP(3)*T+DBDP(4)*T*T+DBDP(5)/TH)*(PRES-PR) 44350 BITG=BITG+(DBDP(6)+DBDP(7)/T+DBDP(8)*T+DBDP(9)*T*T+DBDP(10)/TH)*(PRES*PRES-PR*PR)/2# 44400 BITG=BITG+(DBDP(11)+DBDP(12)/T+DBDP(13)*T+DBDP(14)*T*T+DBDP(15)/TH)*(PRES*PRES*PRES-PR*PR*PR)/3# 44600 BITGT=((-1#)*DBDP(2)/(T*T)+DBDP(3)+2#*DBDP(4)*T+DBDP(5)/(TH*TH))*(PRES-PR) 44650 BITGT=BITGT+((-1#)*DBDP(7)/(T*T)+DBDP(8)+2#*DBDP(9)*T+DBDP(10)/(TH*TH))*(PRES*PRES-PR*PR)/2# 44700 BITGT=BITGT+((-1#)*DBDP(12)/(T*T)+DBDP(13)+2#*DBDP(14)*T+DBDP(15)/(TH*TH))*(PRES*PRES*PRES-PR*PR*PR)/3# 44800 DBLDTITG=(2#*DBDP(3)/T+6#*DBDP(4)+1294#*DBDP(5)/(T*TH*TH*TH))*(PRES-PR) 44850 DBLDTITG=DBLDTITG+(2#*DBDP(8)/T+6#*DBDP(9)+1294#*DBDP(10)/(T*TH*TH*TH))*(PRES*PRES-PR*PR)/2# 44900 DBLDTITG=DBLDTITG+(2#*DBDP(13)/T+6#*DBDP(14))*(PRES*PRES*PRES-PR*PR*PR)/3# 44950 DBLDTITG=DBLDTITG+(1294#*DBDP(15)/(T*TH*TH*TH))*(PRES*PRES*PRES-PR*PR*PR)/3# 44960 OSC=1#-APHI*SQR(MI)/(1#+1.2#*SQR(MI))+MOL*(BM0+BITG+BM1*EAI)+2#*MOL*MOL*CM 44970 GM=(-1#)*APHI*(SQR(MI)/(1#+1.2#*SQR(MI))+4#*HF)+2#*MOL*(BM0+BITG+(BM1/(4#*MI))*GM1)+3#*MOL*MOL*CM 44980 GM=EXP(GM) 44990 PHIL=2#*AH*HF-2#*MOL*RGAS*T*T*(BL+BITGT+MOL*CL) 45000 PHICP=CPK+2#*AJ*HF-2#*MOL*RGAS*T*T*(BJ+DBLDTITG+MOL*CJ) 45050 CPX=PHICP-CPK 45100 BM=BM0+2#*BM1*(1#-(1#+2#*SQR(MI))*EAI)/(4#*MI) 45150 BM=BM+BITG 45200 GEX=(-8#)*APHI*HF*MI+2#*MOL*MOL*BM+2#*MOL*MOL*MOL*CM 45250 GEX=GEX*RGAS*T 45300 SX=(PHIL-GEX/MOL)/T 45350 SSPEC=S*(1000/MW)+MOL*(SKCL+SX)+2#*MOL*RGAS*(1#-LOG(MOL)) 45400 SSPEC=SSPEC/(1000+MOL*KCL) 45450 HSPEC=H*(1000/MW)+MOL*(HKCL+PHIL) 45500 HSPEC=HSPEC/(1000+MOL*KCL) 45550 CPSPEC=CPW*(1000/MW)+MOL*PHICP 45600 CPSPEC=CPSPEC/(1000+MOL*KCL) 45650 LPRINT USING"m=#.##### Density(g/cm3)= +#.#####";MOL,DSOLN 45700 LPRINT USING" Osmotic coeff= +#.###";OSC 45750 LPRINT USING" Activity coeff= +#.###";GM 45800 LPRINT USING" phiL/RT= +##.###";PHIL/(RGAS*T) 45850 LPRINT USING" Ex entr/R= +##.###";SX/RGAS 45900 LPRINT USING" phiCp/R= +###.##";PHICP/RGAS 45950 LPRINT USING" Hspecific(J/g)= +#.####^^^^^";HSPEC 46000 LPRINT USING" Sspecific(J/g K)= +#.###";SSPEC 46050 LPRINT USING" Cpspecific(J/g K)= +#.###";CPSPEC 46100 LPRINT 46150 RETURN