work around overflow in ice_fem_fct as discussed with Sergey

FESOM · Jan 25, 2024 · 8a2bfd3 · 8a2bfd3
1 parent 2fc601a
commit 8a2bfd3
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/ice_fct.F90 b/src/ice_fct.F90
@@ -817,14 +817,14 @@ subroutine ice_fem_fct(tr_array_id, ice, partit, mesh)
 
         flux=icepplus(n)
         if (abs(flux)>0) then
-            icepplus(n)=min(1.0_WP,tmax(n)/flux)
+            icepplus(n)=min(1.0_WP,tmax(n)/max(flux,1.e-12))
         else
             icepplus(n)=0._WP
         end if
 
         flux=icepminus(n)
         if (abs(flux)>0) then
-            icepminus(n)=min(1.0_WP,tmin(n)/flux)
+            icepminus(n)=min(1.0_WP,tmin(n)/min(flux,-1.e-12))
         else
             icepminus(n)=0._WP
         end if