There appears to be a typo in solve_cool.F and has caused problems when using precision-32.
parameter (tolerance = 1.e-52_RKIND)
parameter (tolerance = 1.e-5_RKIND)
This looks correct to me. I'm not 100% convinced that it will require a regeneration of the gold standard, as it is a tolerance that also only applies when RadiativeCooling == 1 and MultiSpecies == 0 and precision == 32. It's definitely the right fix, though. Can you check if this requires a regeneration of the gold standard? Thanks, Yuan!
Thanks indeed! I just checked the code prior to my "upgrades" to the precision handling, and I did add in a typo when updating this value. I'm going to accept the PR.
The gold standard was generated for precision == 64, and this change has no effect on that, so it does not require a regeneration of the gold standard. A related issue is: does precision == 32 pass the gold-standard? I don't know the answer to that -- I plan to find out, but I don't think it should hold up the acceptance of this pull request.