Bug fix for divide by zero in hail size

joeolson42 · Nov 4, 2023 · 571cb59 · 571cb59
1 parent 6ea2bce
commit 571cb59
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/physics/module_mp_thompson.F90 b/physics/module_mp_thompson.F90
@@ -6527,6 +6527,8 @@ function hail_mass_99th_percentile(kts, kte, qg, temperature, pressure, qv) resu
       rho(k) = 0.622*pressure(k)/(R*temperature(k)*(max(1.e-10, qv(k))+0.622))
       if (qg(k) .gt. R1) then
          rg(k) = qg(k)*rho(k)
+      else
+         rg(k) = R1
       endif 
    enddo