mHM
The mesoscale Hydrological Model
|
Canopy interception. More...
Functions/Subroutines | |
elemental pure subroutine, public | canopy_interc (pet, interc_max, precip, interc, throughfall, evap_canopy) |
Canopy interception. | |
Canopy interception.
This module deals with processes related to canopy interception, evaporation and throughfall.
COPYING
and COPYING.LESSER
provided with this software. The complete GNU license text can also be found at http://www.gnu.org/licenses/. elemental pure subroutine, public mo_canopy_interc::canopy_interc | ( | real(dp), intent(in) | pet, |
real(dp), intent(in) | interc_max, | ||
real(dp), intent(in) | precip, | ||
real(dp), intent(inout) | interc, | ||
real(dp), intent(out) | throughfall, | ||
real(dp), intent(out) | evap_canopy | ||
) |
Canopy interception.
Calculates throughfall. Updates interception and evaporation intensity from canopy. Throughfall ( \(F\)) is estimated as a function of the incoming precipitation ( \(P\)), the current status of the canopy water content ( \(C\)), and the max. water
\[ F = Max( (P + C - C_{max}), 0) \]
Evaporation ( \(E\)) from canopy is estimated as a fraction of the potential evapotranspiration( \(E_{p}\)) depending on the current status of the canopy water content ( \(C\)) and the max. water content( \(C_{max}\)) that can be intecepted by the vegetation.
\[ E = E_{p}(C/C_{max})^{2/3} \]
ADDITIONAL INFORMATION content( \(C_{max}\)) that can be intecepted by the vegetation. canopy_interc(pet, interc_month_max, interc_max, precip, throughfall, evap_canopy, interc)
[in] | REAL(dp) :: pet | Potential evapotranspiration [mm TS-1] |
[in] | REAL(dp) :: interc_max | Maximum interception [mm] |
[in] | REAL(dp) :: precip | Daily mean precipitation [mm] |
[in,out] | REAL(dp) :: interc | Interception [mm] |
[out] | REAL(dp) :: throughfall | Throughfall [mm TS-1] |
[out] | REAL(dp) :: evap_canopy | Real evaporation intensity from canopy[mm TS-1] |
Definition at line 75 of file mo_canopy_interc.f90.
References mo_common_constants::eps_dp.
Referenced by mo_mhm::mhm().