In a DEM case, these code are written to prevent divergence. But it is shown that the bounded temperatures are abnormally low.
I don’t know wyh the temperatures can be nearly zero. I am gratful for any suggestions or discussions. usr_rates_des.f is shared.
usr_rates_des.f (8.9 KB)