In UpdateState_TopLayerFluxes, h2osoi_ice has gone significantly negative

[zhangguangzhi]

I created the case for CLM45BGCDV in CLM5.35 (1850_DATM%CRU_CLM45%BGCDV_SICE_SOCN_RTM_SGLC_SWAV). However, an error was reported In UpdateState_TopLayerFluxes, h2osoi_ice has gone significantly negative. I have referred to the previous posts (ERROR: In UpdateState_TopLayerFluxes, h2osoi_ice has gone significantly negative)

I changed "custom_rel_epsilon" from 1. E - 12 _r8 to 1. E - 11 _r8 here in the components/CLM/SRC/biogeophys SnowHydrologyMod. F90. But there will still be this mistake.

The following compressed package is the log file

 

[oleson]

Unfortunately, we no longer support the DV model, so all I can suggest is to try 1.e-10. Otherwise, you can look here for a more robust fix that might work:

Need some changes to avoid negative h2osoi_ice in UpdateState_TopLayerFluxes · Issue #1979 · ESCOMP/CTSM

Brief summary of bug In UpdateState_TopLayerFluxes, there is a somewhat arbitrary epsilon used to determine if h2osoi_ice and h2osoi_liq are close enough to zero that they should be truncated to ze...

github.com

which is available in more recent versions of the model.