forrie Posted February 4, 2013 Share Posted February 4, 2013 Hi I am trying to model calcite dissolution under both neutral and acid conditions as the model start with the former and end with the latter. However, GWB totally ignores dissolution when I assign a script to the model (calcite reacted = 0). The script is given below. Is my script wrong that GWB ignores it? I also assigned a rate const and a surface area in the interface for the neutral conditions. I assumed that neutral/acid conditions change at 6.3. Thanks in advance. Regards Johan IF pH > 6.3 THEN 20 ELSE 40 20: rate = rate_const * surface * (1 - Q/K) GOTO 60 40: rate = 3.235936569E+05 * rate_const * surface * (1 - Q/K) * activity("H+")^powerD(1) 60: RETURN rate Quote Link to comment Share on other sites More sharing options...
Join the conversation
You can post now and register later. If you have an account, sign in now to post with your account.