Changeset 6053 for LMDZ6/trunk/libf/phylmd/pbl_surface_subsrf_mod.F90

Timestamp:

Jan 23, 2026, 1:11:53 PM (29 hours ago)

Author:

evignon

Message:

corrections dans la routine simplehydrol et ajout de diagnostiques
pour quantifier les flux d'eau liés à la fonte de la neige sur les surfaces
englacees

File:

• LMDZ6/trunk/libf/phylmd/pbl_surface_subsrf_mod.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/pbl_surface_subsrf_mod.F90

@@ -54 +54 @@
        cdragh,    cdragm,                             &
        beta, &
-       icesub_lic, alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
+       icesub_ice, icemelt_ice, alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
        qsat2m,                 &
        d_t,       d_q,    d_qbs,    d_u,      d_v, d_t_diss,            &
@@ -277 +277 @@
     REAL, DIMENSION(klon),        INTENT(INOUT)       :: cdragm     ! drag coefficient for wind
     REAL, DIMENSION(klon),        INTENT(INOUT)       :: alb3_lic
-    REAL, DIMENSION(klon),        INTENT(INOUT)       :: icesub_lic ! ice (no snow!) sublimation over ice sheet 
+    REAL, DIMENSION(klon),        INTENT(INOUT)       :: icesub_ice ! ice (no snow!) sublimation flux over ice sheet and sea ice 
+    REAL, DIMENSION(klon),        INTENT(INOUT)       :: icemelt_ice ! ice (no snow!) melting flux over ice sheet and sea ice 
     REAL, DIMENSION(klon,klev),   INTENT(INOUT)       :: d_t_w      !   !
     REAL, DIMENSION(klon,klev),   INTENT(INOUT)       :: d_q_w      !      !  Tendances dans les poches
@@ -529 +530 @@
     REAL, DIMENSION(knon)              :: AcoefQBS, BcoefQBS
     REAL, DIMENSION(knon)              :: ypsref
-    REAL, DIMENSION(knon)              :: yevap, yevap_pot, ytsurf_new, yalb3_new, yicesub_lic
+    REAL, DIMENSION(knon)              :: yevap, yevap_pot, ytsurf_new, yalb3_new, yicesub, yicemelt
     REAL, DIMENSION(knon,nsw)          :: yalb_dir_new, yalb_dif_new
     REAL, DIMENSION(knon,klev)         :: y_d_t, y_d_q, y_d_t_diss, y_d_qbs
@@ -1637 +1638 @@
                   ypsref, yu1, yv1, ygustiness, yrugoro, pctsrf, &
                   ysnow, yqsurf, yqsol,yqbs1, yagesno, &
-                  ytsoil, yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yicesub_lic, yfluxsens,yfluxlat, &
+                  ytsoil, yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yicesub, yicemelt, yfluxsens,yfluxlat, &
                   yfluxbs, ytsurf_new, y_dflux_t, y_dflux_q, &
                   yzmea, yzsig, ycldt, &
@@ -1657 +1658 @@
                 sissnow(i)   = ysissnow(j)
                 runoff(i)    = yrunoff(j)
-                icesub_lic(i) = yicesub_lic(j)*ypct(j)
+                icesub_ice(i) = icesub_ice(i) + yicesub(j)*ypct(j)
+                icemelt_ice(i) = icemelt_ice(i) + yicemelt(j)*ypct(j)
              ENDDO
              ! Martin
@@ -1783 +1785 @@
                ypsref, yu1, yv1, ygustiness, pctsrf, &
                ysnow, yqsurf, yqsol, yagesno, ytsoil, &
-               yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yfluxsens,yfluxlat,&
+               yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yicesub, yicemelt, yfluxsens,yfluxlat,&
                ytsurf_new, y_dflux_t, y_dflux_q, &
                y_flux_u1, y_flux_v1, &
@@ -1794 +1796 @@
 #endif               
          &      )
-          
+
+             DO j = 1, knon
+                i = ni(j)
+                icesub_ice(i) = icesub_ice(i) + yicesub(j)*ypct(j)
+                icemelt_ice(i) = icemelt_ice(i) + yicemelt(j)*ypct(j)
+             ENDDO
 ! Special DICE MPL 05082013 puis BOMEX MPL 20150410
        IF (ok_prescr_ust) THEN