experimental condition range T=25-70C, pH=1-3

#### -start

1 rem unit should be mol, kgw-1 and second-1 2 rem parm(1) is surface area in the unit of m2/kgw 3 rem calculation of surface area can be found in the note 4 rem M is current moles of minerals 5 rem M0 is the initial moles of minerals 6 rem parm(2) is a correction factor 40 SR\_mineral=SR("Galena") 41 if (M<0) then goto 200 42 if (M=0 and SR\_mineral<1) then goto 200 43 if (M0<=0) then SA=PARM(1) else SA=PARM(1)\*(M/M0)ˆ0.67 50 if (SA<=0) then SA=1 60 R=8.31451 70 J=10ˆ-5.7\*exp(-23000/R/TK)\*ACT("H+")ˆ0.43 90 Rate=J\*(1-Sr\_mineral)\*SA\*parm(2) 100 moles=Rate\*Time 200 save moles -end ########### #As2S3(a) ########### Orpiment
