Merge pull request #562 from trackow/production_DE_icefct

production_DE fix in ice_fem_fct for occurring Overflows
FESOM · Jan 26, 2024 · f6b392c · f6b392c
2 parents 2fc601a + 8a2bfd3
commit f6b392c
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