Click here to donate to keep PhreeqcUsers open

#formula 1: Se(Sea,Seb)=-rho*f*kd*(d(conc^beta)/dt)#formula 2: Sk=-rho*omega*((1-f)*kd*(conc^beta)-Sk)RatesSea-start20 beta = parm(1)30 rate = TOT("Cd")^ beta40 moles = rate * time50 PUT(rate, 10)200 SAVE moles-endSeb############################-start20 rho = parm(1)25 f = parm(2)26 kd = parm(3)30 rate1 = GET(10)40 rate2 = f * rho * kd * rate150 moles = - rate2 * time51 if (moles > m) then moles = m200 SAVE moles-endSk############################-start20 rho = parm(1)25 f = parm(2)26 kd = parm(3)27 beta = parm(4)28 w = parm(5)30 rate = rho * w * (1 - f) * kd * TOT("Cd") ^ beta - rho * w * kin("Sk")50 moles = - rate * time51 if (moles > m) then moles = m200 SAVE moles-end