Improving the numerical stability of the MAGIC model.
MAGIC is one of the most widely used models for forecasting long-term acidification. The model's code, however, has been experiencing numerical instability, though this might not be widely known to the public users. The major instability comes from the analytical solution to two cubic equations for calculating SO4(2-) concentration and the exchangeable fraction of Al on the soils. The mathematical algorithm for calculating the concentration of SO4(2-) from a quadratic equation is also found unstable. This paper is aimed at improving the instability above through proved numerical algorithms.
Deng
Y X
YX
Department of Environment Sciences and Engineering, Tsinghua University, Beijing 100084, China. dyx99g@mails.tsinghua.edu.cn
Chen
J N
JN
Hu
X T
XT
D000327:Adsorption; D000465:Algorithms; D017607:Aluminum Compounds / Q000032:analysis / Q000737:chemistry; D005618:Fresh Water / Q000737:chemistry*; D006813:Humidity; D006863:Hydrogen-Ion Concentration; D007474:Ion Exchange; D008956:Models, Chemical; D012987:Soil / Q000032:analysis; D013431:Sulfates / Q000032:analysis / Q000737:chemistry; D014867:Water / Q000032:analysis
D017607:Aluminum Compounds; D012987:Soil; D013431:Sulfates; D014867:Water