Changeset 4737 for LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Oct 24, 2023, 12:13:53 PM (7 months ago)

Author:

oboucher

Message:

moving the calculation of wfbilo wfrain and wfsnow inside phys_output_write_mod so that it uses the timestep values of rain and snow instead of the values from the previous timestep as it was the case in pbl_surface. Thus there is exact mass balance in the output.

File:

• LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -184 +184 @@
        debut,     lafin,                              &
        rlon,      rlat,      rugoro,   rmu0,          &
-       lwdown_m,  cldt,          &
+       lwdown_m,  cldt,                               &
        rain_f,    snow_f,    bs_f, solsw_m,  solswfdiff_m, sollw_m,       &
        gustiness,                                     &
-       t,         q,        qbs,  u,        v,             &
+       t,         q,        qbs,  u,        v,        &
 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
 !!       t_x,       q_x,       t_w,      q_w,           &
@@ -199 +199 @@
        beta, &
 !>jyg
-       alb_dir_m,    alb_dif_m,  zxsens,   zxevap,  zxsnowerosion,  &
+       alb_dir_m,    alb_dif_m,  zxsens,   zxevap,  zxsnowerosion,      &
        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, & 
+       zxtsol,    zxfluxlat, zt2m,     qsat2m, zn2mout,                 &
+       d_t,       d_q,    d_qbs,    d_u,      d_v, d_t_diss,            &
 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
-       d_t_w,     d_q_w,                              &
-       d_t_x,     d_q_x,                              & 
+       d_t_w,     d_q_w,                             &
+       d_t_x,     d_q_x,                             & 
 !!       d_wake_dlt,d_wake_dlq,                         &
-       zxsens_x,  zxfluxlat_x,zxsens_w,zxfluxlat_w,   &
+       zxsens_x,  zxfluxlat_x,zxsens_w,zxfluxlat_w,  &
 !!!
 !!! nrlmd le 13/06/2011
-       delta_tsurf,wake_dens,cdragh_x,cdragh_w,       &
-       cdragm_x,cdragm_w,kh,kh_x,kh_w,                &
-!!!
-       zcoefh,    zcoefm,    slab_wfbils,             &
-       qsol,    zq2m,      s_pblh,   s_plcl,        &
+       delta_tsurf,wake_dens,cdragh_x,cdragh_w,      &
+       cdragm_x,cdragm_w,kh,kh_x,kh_w,               &
+!!!
+       zcoefh,    zcoefm,    slab_wfbils,            &
+       qsol,    zq2m,      s_pblh,   s_plcl,         &
 !!!
 !!! jyg le 08/02/2012
-       s_pblh_x, s_plcl_x,   s_pblh_w, s_plcl_w,      &
-!!!
-       s_capCL,   s_oliqCL,  s_cteiCL, s_pblT,        &
-       s_therm,   s_trmb1,   s_trmb2,  s_trmb3,       &
-       zustar,zu10m,  zv10m,    fder_print,    &
-       zxqsurf, delta_qsurf,                       &
-       rh2m,      zxfluxu,  zxfluxv,               &
+       s_pblh_x, s_plcl_x,   s_pblh_w, s_plcl_w,     &
+!!!
+       s_capCL,   s_oliqCL,  s_cteiCL, s_pblT,       &
+       s_therm,   s_trmb1,   s_trmb2,  s_trmb3,      &
+       zustar,zu10m,  zv10m,    fder_print,          &
+       zxqsurf, delta_qsurf,                         &
+       rh2m,      zxfluxu,  zxfluxv,                 &
        z0m, z0h,   agesno,  sollw,    solsw,         &
        d_ts,      evap,    fluxlat,   t2m,           &
-       wfbils,    wfbilo, wfevap, wfrain, wfsnow,   & 
-       flux_t,   flux_u, flux_v,                    &
-       dflux_t,   dflux_q,   zxsnow,                  &
+       wfbils,    wfevap,                            & 
+       flux_t,   flux_u, flux_v,                     &
+       dflux_t,   dflux_q,   zxsnow,                 &
 !jyg<
 !!       zxfluxt,   zxfluxq,   q2m,      flux_q, tke,   &
@@ -235 +235 @@
 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012
 !!        tke_x,     tke_w                              &
-       wake_dltke,                                     &
-        treedrg                                   &
+       wake_dltke,                                   &
+        treedrg                                      &
 !FC
 !!!
@@ -503 +503 @@
     REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: t2m        ! temperature at 2 meter height
     REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: wfbils     ! heat balance at surface
-    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: wfbilo     ! water balance at surface weighted by srf
     REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: wfevap     ! water balance (evap) at surface weighted by srf
-    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: wfrain     ! water balance (rain) at surface weighted by srf
-    REAL, DIMENSION(klon, nbsrf), INTENT(OUT)       :: wfsnow     ! water balance (snow) at surface weighted by srf
     REAL, DIMENSION(klon, klev, nbsrf), INTENT(OUT) :: flux_t     ! sensible heat flux (CpT) J/m**2/s (W/m**2)
                                                                   ! positve orientation downwards
@@ -1027 +1024 @@
  snowerosion(:,:)=0. 
  fluxlat(:,:)=0.
- wfbils(:,:)=0. ; wfbilo(:,:)=0.
- wfevap(:,:)=0. ; wfrain(:,:)=0. ; wfsnow(:,:)=0.
+ wfbils(:,:)=0. ; wfevap(:,:)=0. ;
  flux_t(:,:,:)=0. ; flux_q(:,:,:)=0. ; flux_u(:,:,:)=0. ; flux_v(:,:,:)=0.
  flux_qbs(:,:,:)=0.
@@ -3458 +3454 @@
           wfbils(i,nsrf) = ( solsw(i,nsrf) + sollw(i,nsrf) &
                + flux_t(i,1,nsrf) + fluxlat(i,nsrf) ) * pctsrf(i,nsrf)
-          wfbilo(i,nsrf) = (evap(i,nsrf)-(rain_f(i)+snow_f(i)))*pctsrf(i,nsrf)
+
           wfevap(i,nsrf) = evap(i,nsrf)*pctsrf(i,nsrf)
-          wfrain(i,nsrf) = rain_f(i)*pctsrf(i,nsrf)
-          wfsnow(i,nsrf) = snow_f(i)*pctsrf(i,nsrf)
 
           zxtsol(i)    = zxtsol(i)    + ts(i,nsrf)      * pctsrf(i,nsrf)