medfate (version 0.2.2)

forest2swbInput: Input variables for soil water balance

Description

Function forest2swbInput takes an object of class forest to calculate input variables for swb. Function swbInput does the same from input vectors.

Usage

forest2swbInput(x, SpParams, d, gdd = NA, petMode = "Input", hydraulicMode = "Simple")
swbInput(SP, LAI, H, CR,  V, SpParams,  petMode = "Input", hydraulicMode = "Simple")

Arguments

x

An object of class forest.

SpParams

A data frame with species parameters (see SpParamsMED and SpParamsMED).

d

A vector containing the depth of each soil layer.

gdd

Growth degree days to account for leaf phenology effects (in Celsius).

petMode

Potential evapotranspiration mode (either 'Input' or 'PenmanMonteith').

hydraulicMode

Hydraulic model (either 'Simple' or 'Sperry').

SP

An integer vector of plant cohort species identity (first species is 0).

LAI

A numeric vector (same length as SP) with plant cohort leaf area index (one-side leaf area relative to plot area).

H

A numeric vector (same length as SP) with plant cohort total heights (in cm).

CR

A numeric vector (same length as SP) with plant cohort crown ratio values (between 0 and 1).

V

A numeric matrix (with three columns and as many rows as the length as SP) containing the proportion of roots of each plant in each soil layer.

Value

A data frame with the following columns (additional columns will difer depending on the value of hydraulicMode and petMode):

• LAI: Leaf area index (in m2/m2).

• SP: Species identity (an integer).

• H: Plant height (in cm).

• V.1, V.2, V.3: The proportion of fine roots in each layer.

• RC_min: Minimum canopy (stomatal) resistance (in s/m) (for petMode = "PenmanMonteith").

• psiExtr: Water potential corresponding to 50% relative conductance (in kPa) (for hydraulicMode = "Simple").

• VC_kxmax: Maximum hydraulic conductance (in L/m2/day) (for hydraulicMode = "Sperry").

• VC_c, VC_d: Parameters of the vulnerability curve (for hydraulicMode = "Sperry").

• k: PAR extinction coefficient.

• g: Canopy water retention capacity per LAI unit (mm/LAI).

• CR: Crown ratio (crown length to total height).

• Sgdd: Growth degree days needed for leaf budburst (for winter decideous species).

• Transpiration: Plant cohort transpiration (filled with zeroes before simulations).

Details

Function forest2swbInput extracts height and species identity from plant cohorts of x. For each plant cohort also calculates leaf area index and the distribution of fine roots across soil. Both forest2swbInput and swbInput find parameter values for each plant cohort according to the parameters of its species as specified in SpParams.

See Also

swb, soil, forest, SpParamsMED, defaultSoilParams

Examples

# NOT RUN {
#Load example plot plant data
data(exampleforest)

#Default species parameterization
data(SpParamsMED)

#Initialize soil with default soil params
s = soil(defaultSoilParams())

#Prepare input
swbInput = forest2swbInput(exampleforest,SpParamsMED, s$dVec)

# }