Changeset 5489 for LMDZ6/branches/contrails/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Jan 17, 2025, 6:16:25 PM (13 days ago)

Author:

aborella

Message:

Merge with trunk

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/pbl_surface_mod.F90 (modified) (7 diffs)

LMDZ6/branches/contrails/libf/phylmd/pbl_surface_mod.F90

@@ -277 +277 @@
 !>jyg
        alb_dir_m,    alb_dif_m,  zxsens,   zxevap,  zxsnowerosion,      &
-       alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
+       icesub_lic, alb3_lic,  runoff,    snowhgt,   qsnow,     to_ice,    sissnow,  &
        zxtsol,    zxfluxlat, zt2m,     qsat2m, zn2mout,                 &
        d_t,       d_q,    d_qbs,    d_u,      d_v, d_t_diss,            &
@@ -522 +522 @@
     REAL, DIMENSION(klon),        INTENT(OUT)       :: zxevap     ! water vapour flux at surface, positiv upwards
     REAL, DIMENSION(klon),        INTENT(OUT)       :: zxsnowerosion     ! blowing snow flux at surface
+    REAL, DIMENSION(klon),        INTENT(OUT)       :: icesub_lic ! ice (no snow!) sublimation over ice sheet 
     REAL, DIMENSION(klon),        INTENT(OUT)       :: zxtsol     ! temperature at surface, mean for each grid point
 !!! jyg le ???
@@ -745 +746 @@
     REAL, DIMENSION(klon)              :: AcoefQBS, BcoefQBS
     REAL, DIMENSION(klon)              :: ypsref
-    REAL, DIMENSION(klon)              :: yevap, yevap_pot, ytsurf_new, yalb3_new
+    REAL, DIMENSION(klon)              :: yevap, yevap_pot, ytsurf_new, yalb3_new, yicesub_lic
 !albedo SB >>>
     REAL, DIMENSION(klon,nsw)          :: yalb_dir_new, yalb_dif_new
@@ -1246 +1247 @@
  zxfluxt(:,:)=0. ; zxfluxq(:,:)=0.; zxfluxqbs(:,:)=0.
  qsnow(:)=0. ; snowhgt(:)=0. ; to_ice(:)=0. ; sissnow(:)=0.
- runoff(:)=0.
+ runoff(:)=0. ; icesub_lic(:)=0.
 #ifdef ISO
 zxxtevap(:,:)=0.
@@ -2498 +2499 @@
                   ypsref, yu1, yv1, ygustiness, yrugoro, pctsrf, &
                   ysnow, yqsurf, yqsol,yqbs1, yagesno, &
-                  ytsoil, yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap,yfluxsens,yfluxlat, &
+                  ytsoil, yz0m, yz0h, SFRWL, yalb_dir_new, yalb_dif_new, yevap, yicesub_lic, yfluxsens,yfluxlat, &
                   yfluxbs, ytsurf_new, y_dflux_t, y_dflux_q, &
                   yzmea, yzsig, ycldt, &
@@ -2521 +2522 @@
                 sissnow(i)   = ysissnow(j)
                 runoff(i)    = yrunoff(j)
+                icesub_lic(i) = yicesub_lic(j)*ypct(j)
              ENDDO
              ! Martin
@@ -3225 +3227 @@
 
        ENDIF  ! (iflag_split .eq.0)
-!!!
+
 
        ! tendencies of blowing snow