Changeset 4296 for palm/trunk


Ignore:
Timestamp:
Nov 20, 2019 11:35:12 AM (5 years ago)
Author:
maronga
Message:

improved treatment of interception reservoir in land surface model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r4261 r4296  
    2525! -----------------
    2626! $Id$
     27! Bugfix: partitioning of LE from liquid water reservoir fixed. Bare soils are
     28! now allowed to store liquid water at the surface.
     29!
     30! 4261 2019-10-09 17:58:00Z scharf
    2731! bugfix for rev. 4258: deallocate temporary arrays
    2832!
     
    272276              rho_c_soil         = 2.19E6_wp,  & ! volumetric heat capacity of soil (J/m3/K)
    273277              rho_c_water        = 4.20E6_wp,  & ! volumetric heat capacity of water (J/m3/K)
    274               m_max_depth        = 0.0002_wp     ! Maximum capacity of the water reservoir (m)
     278              m_max_depth        = 0.0002_wp     ! Maximum capacity of the water reservoir on a flat surface (leaf/bare soil) (m)
    275279
    276280
     
    18231827!--    assumed, while paved surfaces might hold up 1 mm of water. The
    18241828!--    liquid water fraction for paved surfaces is calculated after
    1825 !--    Noilhan & Planton (1989), while the ECMWF formulation is used for
    1826 !--    vegetated surfaces and bare soils.
     1829!--    Masson (2000) (TEB model) and originates from Noilhan & Planton (1989),
     1830!--    while the ECMWF formulation is used for vegetated surfaces and bare soils.
    18271831       IF ( surf%pavement_surface(m) )  THEN
    18281832          m_liq_max = m_max_depth * 5.0_wp
     
    18581862       ELSE
    18591863          f_qsws_veg  = rho_lv * surf%c_veg(m) *                               &
    1860                             ( 1.0_wp        - surf%c_liq(m)    ) /             &
     1864                            ( 1.0_wp        - surf%c_liq(m)  ) /               &
    18611865                            ( surf%r_a(m) + surf%r_canopy(m) )
    1862           f_qsws_soil = rho_lv * (1.0_wp    - surf%c_veg(m)    ) /             &
     1866          f_qsws_soil = rho_lv * (1.0_wp    - surf%c_veg(m)  )                 &
     1867                               * (1.0_wp    - surf%c_liq(m)  ) /               &
    18631868                            ( surf%r_a(m) + surf%r_soil(m)   )
    1864           f_qsws_liq  = rho_lv * surf%c_veg(m) * surf%c_liq(m)   /             &
    1865                               surf%r_a(m)
     1869          f_qsws_liq  = rho_lv * surf%c_liq(m) / surf%r_a(m)
    18661870       ENDIF
    18671871
Note: See TracChangeset for help on using the changeset viewer.