Hi Jia,
thank you for the fix. I run some test inputs and it seems working as expected. I have however another question/suggestion/discussion related to the use of CCM in particular and the use of 'react' in general.
The system I am modelling can be in general described as a system containing mineral surfaces with a constant capacitance. With your fix now I am able to run it with a given capacitance value. However, the system capacitance is a function of solution ionic strength, I, so I need to give the C value calculated by hand for the given ionic strength. A better way, and a more general way, is to be able to let the code calculate the capacitance as a function of ionic strength automatically so the model can be run for varying ionic strength but not only for one ionic strength per input. I hope I explained myself clearly. In other words, is there a way to input capacitance as a function of I, but not as a fixed numerical value as it is now. In general, this is like what the 'field variable' is allowed in Xt codes there one may define a field variable which changes in the course of the run. To my understanding there is no such field variables can be used in 'react' now. Is there a way around to allow a field variable in 'react' code?
Once again thx a lot for the fix, it helped a lot.
best regards,
Lian