AquaticEcoDynamics / libaed-water

Code for the AED water quality model
GNU General Public License v3.0
6 stars 6 forks source link

Deducing threshold so pH-TA converges #58

Closed rqthomas closed 2 years ago

rqthomas commented 2 years ago

alk_mode = 0 and co2_model = 1 does not converge on the pH-TA calculation. This seems to be because the error threshold of 1e-10 is close to machine precision. Increasing the threshold of 1e-5 resulted in convergence. 1e-10 is changed to 1e-5 in two places in aed_carbon.F90