32 |
|
|
33 |
COMMON /CARBON_CHEM/ |
COMMON /CARBON_CHEM/ |
34 |
& ak0,ak1,ak2,akw,akb,aks,akf, |
& ak0,ak1,ak2,akw,akb,aks,akf, |
35 |
& ak1p,ak2p,ak3p,aksi, |
& ak1p,ak2p,ak3p,aksi, fugf, |
36 |
& ff,ft,st,bt, Ksp_TP_Calc |
& ff,ft,st,bt, Ksp_TP_Calc |
37 |
_RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
38 |
_RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
45 |
_RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
46 |
_RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
47 |
_RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
48 |
|
C Fugacity Factor added following Val Bennington |
49 |
|
_RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
50 |
_RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
51 |
_RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
52 |
_RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |