Minor modif to v_i calculation in seabed_stress_factor_prob to preven…

…t rare instability
CICE-Consortium · Jan 19, 2023 · 6ce76f1 · 6ce76f1
1 parent 79d22b5
commit 6ce76f1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/cicecore/cicedyn/dynamics/ice_dyn_shared.F90 b/cicecore/cicedyn/dynamics/ice_dyn_shared.F90
@@ -1482,7 +1482,7 @@ subroutine seabed_stress_factor_prob (nx_block, ny_block,          &
             do n =1, ncat
                v_i = v_i + vcat(n)**2 / (max(acat(n), puny))
             enddo
-            v_i = v_i - m_i**2
+            v_i = max((v_i - m_i**2), puny)
 
             mu_i    = log(m_i/sqrt(c1 + v_i/m_i**2)) ! parameters for the log-normal
             sigma_i = sqrt(log(c1 + v_i/m_i**2))