Difference between TOUGHREACT 1.2 and 2

Hi all,

I am observing quite different (quantitative) results between 0D and 1D simulations run with the 1.2 version of TOUGHREACT/ECO2N and the 2.0 version. I am struggling figuring out what could cause the difference.

- I am using the exact same thermodynamic database for the two simulations, hence a v1.2 db. As I take from the manual, the AKCOP parameters in sec-3 and mineral-3 are optional in version 2 and therefore, if no molar volume is specified, I assume the pressure dependency of logKs does not get  computed (and this ensures compatibility with older thermodynamic databases). Am I correct there?

- Did the internal calculation of evolving reactive surface areas change between the two versions?

- What else could originate the difference?

Any hint appreciated, cheers,

Marco