I am using the phreeqc module of HYDRUS-1D. I want to write two formulas of nonequilibrium transport model (two-site model) in rate module. I've written a piece of code that works, but I'm not sure it does, exactly as the formulas say. So I am hoping you could help me correct it! Here are two formula which are shown as jpg and code.
Code: Select all
Rates
Sea
-start
20 beta = parm(1)
30 rate = TOT("Cd")^ beta
40 moles = rate * time
50 PUT(rate, 10)
200 SAVE moles
-end
Seb
############################
-start
20 rho = parm(1)
25 f = parm(2)
26 kd = parm(3)
30 rate1 = GET(10)
40 rate2 = f * rho * kd * rate1
50 moles = - rate2 * time
200 SAVE moles
-end
Sk
############################
-start
20 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 * time
200 SAVE moles
-end