Hi All, In one of our offline simulations we encountered an improbable, but possible situation that causes CLM to crash because F90 returns NaN for the cube of a negative number. F90 then returns NaN and stops. This should not affect climates. Dai and Xubin, can you please check my logic here and confirm my solution. I'd like to suggest to Keith to update their code and WWW page with this change.
Relavant code follows.
We could add the computation of dthv in the loop. Keith, please update the code and WWW page with this change. We will look at these two suggested fixes and implement one of them. The code is currently undergoing some reordering that is non-answer changing. The fix will be implemented once that is complete.