Changeset 4296 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Nov 20, 2019 11:35:12 AM (5 years ago)

Author:

maronga

Message:

improved treatment of interception reservoir in land surface model

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: partitioning of LE from liquid water reservoir fixed. Bare soils are
+! now allowed to store liquid water at the surface.
+! 
+! 4261 2019-10-09 17:58:00Z scharf
 ! bugfix for rev. 4258: deallocate temporary arrays
 ! 
@@ -272 +276 @@
               rho_c_soil         = 2.19E6_wp,  & ! volumetric heat capacity of soil (J/m3/K)
               rho_c_water        = 4.20E6_wp,  & ! volumetric heat capacity of water (J/m3/K)
-              m_max_depth        = 0.0002_wp     ! Maximum capacity of the water reservoir (m)
+              m_max_depth        = 0.0002_wp     ! Maximum capacity of the water reservoir on a flat surface (leaf/bare soil) (m)
 
 
@@ -1823 +1827 @@
 !--    assumed, while paved surfaces might hold up 1 mm of water. The 
 !--    liquid water fraction for paved surfaces is calculated after 
-!--    Noilhan & Planton (1989), while the ECMWF formulation is used for
-!--    vegetated surfaces and bare soils.
+!--    Masson (2000) (TEB model) and originates from Noilhan & Planton (1989), 
+!--    while the ECMWF formulation is used for vegetated surfaces and bare soils.
        IF ( surf%pavement_surface(m) )  THEN
           m_liq_max = m_max_depth * 5.0_wp
@@ -1858 +1862 @@
        ELSE
           f_qsws_veg  = rho_lv * surf%c_veg(m) *                               &
-                            ( 1.0_wp        - surf%c_liq(m)    ) /             &
+                            ( 1.0_wp        - surf%c_liq(m)  ) /               &
                             ( surf%r_a(m) + surf%r_canopy(m) )
-          f_qsws_soil = rho_lv * (1.0_wp    - surf%c_veg(m)    ) /             &
+          f_qsws_soil = rho_lv * (1.0_wp    - surf%c_veg(m)  )                 &
+                               * (1.0_wp    - surf%c_liq(m)  ) /               &
                             ( surf%r_a(m) + surf%r_soil(m)   )
-          f_qsws_liq  = rho_lv * surf%c_veg(m) * surf%c_liq(m)   /             &
-                              surf%r_a(m)
+          f_qsws_liq  = rho_lv * surf%c_liq(m) / surf%r_a(m)
        ENDIF