Changeset 2670 for LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Oct 16, 2016, 11:17:32 AM (8 years ago)

Author:

oboucher

Message:

adding new diagnostics for evap/rain/snow weighted by surface types
this parallels what was done for wbilo (net flux)
new diagnostics are wevap_??? wrain_??? wsnow_???
where ??? is ter, oce, lic or sic.

File:

• LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -179 +179 @@
        z0m, z0h,   agesno,  sollw,    solsw,         &
        d_ts,      evap,    fluxlat,  t2m,           &
-       wfbils,    wfbilo,    flux_t,   flux_u, flux_v,&
+       wfbils,    wfbilo, wfevap, wfrain, wfsnow,   & 
+       flux_t,   flux_u, flux_v,                    &
        dflux_t,   dflux_q,   zxsnow,                  &
 !jyg<
@@ -431 +432 @@
     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
+    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
@@ -865 +869 @@
  fluxlat(:,:)=0.
  wfbils(:,:)=0. ; wfbilo(:,:)=0.
+ wfevap(:,:)=0. ; wfrain(:,:)=0. ; wfsnow(:,:)=0.
  flux_t(:,:,:)=0. ; flux_q(:,:,:)=0. ; flux_u(:,:,:)=0. ; flux_v(:,:,:)=0.
  dflux_t(:)=0. ; dflux_q(:)=0.
@@ -2823 +2828 @@
           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)
+          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)