Learn R Programming

medfate (version 0.8.2)

plot.spwb_day: Displays simulation results for one day

Description

Functions to plot the results of the soil water balance model for one day (see spwb_day) or the transpiration calculations with Sperry (see transp_transpirationSperry).

Usage

# S3 method for spwb_day
plot(x, type="PlantTranspiration", bySpecies = FALSE, xlab = NULL, ylab = NULL, ...)
# S3 method for pwb_day
plot(x, type="PlantTranspiration", bySpecies = FALSE, xlab = NULL, ylab = NULL, ...)

Arguments

x

An object of class spwb_day or pwb_day.

type

The information to be plotted:

  • "LeafPsi":Leaf water potential.

  • "StemPsi":(Upper) stem water potential.

  • "RootPsi":Root crown water potential.

  • "StemPLC":(Average) percentage of loss conductance in the stem conduits.

  • "StemRWC":(Average) relative water content in the stem symplasm.

  • "LeafRWC":Relative water content in the leaf symplasm.

  • "SoilPlantConductance":Overall soil plant conductance (calculated as the derivative of the supply function).

  • "PlantExtraction": Water extracted from each soil layer.

  • "PlantTranspiration": Plant cohort transpiration per ground area.

  • "PlantTranspirationPerLeaf": Plant cohort transpiration per leaf area.

  • "PlantPhotosynthesis": Plant cohort daily net photosynthesis per ground area.

  • "PlantPhotosynthesisPerLeaf": Plant cohort daily net photosynthesis per leaf area.

  • "PlantWaterBalance": Difference between water extraction from the soil and transpired water per ground area.

  • "PlantAbsorbedSWR": Absorbed short wave radiation per ground area (differentiates sunlit and shade leaves).

  • "LeafTranspiration": Instantaneous transpiration per leaf area (differentiates sunlit and shade leaves).

  • "LeafPhotosynthesis": Instantaneous net photosynthesis per leaf area (differentiates sunlit and shade leaves).

  • "LeafAbsorbedSWR": Absorbed short wave radiation per leaf area (differentiates sunlit and shade leaves).

  • "LeafVPD": Leaf vapour pressure deficit (differentiates sunlit and shade leaves).

  • "LeafStomatalConductance": Leaf stomatal conductance (differentiates sunlit and shade leaves).

  • "LeafTemperature": Leaf temperature (differentiates sunlit and shade leaves).

  • "LeafCi": Leaf intercellular CO2 concentration (differentiates sunlit and shade leaves).

  • "LeafIntrinsicWUE": Leaf intrinsic water use efficiency, i.e. the ratio between instantaneous photosynthesis and stomatal conductance (differentiates sunlit and shade leaves).

  • "Temperature": Above-canopy, inside-canopy and soil temperature.

  • "CanopyEnergyBalance": Canopy energy balance components.

  • "SoilEnergyBalance": Soil energy balance components.

bySpecies

Allows aggregating output by species, before drawing plots. Aggregation can involve a sum (as for plant lai or transpiration) or a LAI-weighted mean (as for plant stress or plant water potential).

xlab

x-axis label.

ylab

y-axis label.

...

Additional parameters for function plot.

See Also

spwb_day, plot.spwb

Examples

Run this code
# NOT RUN {
#Load example daily meteorological data
data(examplemeteo)

#Load example plot plant data
data(exampleforest)

#Default species parameterization
data(SpParamsMED)

#Initialize control parameters
control = defaultControl()
control$ndailysteps = 24  

#Initialize soil with default soil params (2 layers)
examplesoil = soil(defaultSoilParams(2), W=c(0.5,0.5))


#Switch to 'Sperry' transpiration mode
control$transpirationMode="Sperry"

#Simulate one day only
x2 = forest2spwbInput(exampleforest,examplesoil, SpParamsMED, control)
d = 100
sd2<-spwb_day(x2, examplesoil, rownames(examplemeteo)[d], 
             examplemeteo$MinTemperature[d], examplemeteo$MaxTemperature[d], 
             examplemeteo$MinRelativeHumidity[d], examplemeteo$MaxRelativeHumidity[d], 
             examplemeteo$Radiation[d], examplemeteo$WindSpeed[d], 
             latitude = 41.82592, elevation = 100, 
             slope= 0, aspect = 0, prec = examplemeteo$Precipitation[d])

#Display transpiration for subdaily steps
plot(sd2, "PlantTranspiration")
# }

Run the code above in your browser using DataLab