DESCRIPTION OF CHANGES:
Added a collection of changes to the energy balance at the surface for thin snow on lakes and to enforce an upper bound snow temperature at 280K(!) this is mostly a numerical stability issue, but the excess heat may be required to melt the snow, so didn't want to enforce 273.16K
TESTS CONDUCTED:
Ran tests with CMIP5 model simulations that crashed, runs more stably now.
NOTES:
Ideally the numerical stability should be addressed with better numerics; however, this is a very rare case (1x in ~1000 years of simulation).
TYPE: stability enhancements
KEYWORDS: lake model, numerical stability
SOURCE: Ethan Gutmann, NCAR
DESCRIPTION OF CHANGES: Added a collection of changes to the energy balance at the surface for thin snow on lakes and to enforce an upper bound snow temperature at 280K(!) this is mostly a numerical stability issue, but the excess heat may be required to melt the snow, so didn't want to enforce 273.16K
TESTS CONDUCTED: Ran tests with CMIP5 model simulations that crashed, runs more stably now.
NOTES: Ideally the numerical stability should be addressed with better numerics; however, this is a very rare case (1x in ~1000 years of simulation).
Checklist