Changeset 6180

Timestamp:

Apr 20, 2026, 8:53:25 PM (2 weeks ago)

Author:

evignon

Message:

seuil sur hice pour le premier pas de temps (concerne seulement les nouvelles params de banquise de Guillaume G)

File:

• LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90

@@ -584 +584 @@
                f_cond_s = 0.
                ! 1 / ice-layer heat capacity
-               cal(i) = 2.*RCPD/(hice(ki)*ice_den*ice_cap)
+               cal(i) = 2.*RCPD/(max(hice(ki),h_ice_min)*ice_den*ice_cap)
                ! adjustment time-scale of conductive flux
-               dif_grnd(i) = cal(i)*ice_cond/(hice(ki)*ice_den)/RCPD
+               dif_grnd(i) = cal(i)*ice_cond/(max(hice(ki),h_ice_min)*ice_den)/RCPD
                ! penetrative shortwave flux...
                f_swpen = swnet(i)*pen_frac*exp(-pen_ext*hice(ki))
@@ -645 +645 @@
                f_cond_s = sno_cond*(tice(ki) - tsurf_new(i))/snow(i)
                ! 1 / heat capacity and conductive timescale
-               uscap = 2./ice_cap/(snow(i) + hice(ki)*ice_den)
-               ustau = uscap*ice_cond/(hice(ki)*ice_den)
+               uscap = 2./ice_cap/(snow(i) + max(hice(ki),h_ice_min)*ice_den)
+               ustau = uscap*ice_cond/(max(hice(ki),h_ice_min)*ice_den)
                ! update ice temp
                tice(ki) = (tice(ki) + dtime*(ustau*t_freeze - uscap*f_cond_s)) &
@@ -654 +654 @@
             END IF
             ! ice conductive flux (pos up)
-            f_cond_i = ice_cond*(t_freeze - tice(ki))/(hice(ki)*ice_den)
+            f_cond_i = ice_cond*(t_freeze - tice(ki))/(max(hice(ki),h_ice_min)*ice_den)
             f_bot(i) = fsea - f_cond_i
             fcdi(ki) = f_cond_i - fsea
@@ -746 +746 @@
 !               fsic(i) = pctsrf(ki, is_sic)
 !               hice(ki) = (0.166 + 0.911 * (fsic(i)**2))*(1. + 2.*(fsic(i)-0.351))
+!               hice(ki) = max(h_ice_min , hice(ki))
 !            END DO
           END IF