MakeSeasLAI: Construct the seasonal course of leaf area index from...

In pschmidtwalter/brook90r: Run the LWF-Brook90 hydrological model from within R.

Description

A daily sequence of leaf area index is derived from maximum and minimum values, dates and shape parameters using different methods.

Usage

MakeSeasLAI(method = "b90", year, maxlai, winlaifrac = 0,
  budburst.doy = 121, leaffall.doy = 279, emerge.dur = 28,
  leaffall.dur = 58, shape.optdoy = 220, shape.budburst = 0.5,
  shape.leaffall = 10, lai.doy = c(1, 121, 150, 280, 320, 365),
  lai.frac = c(0, 0, 0.5, 1, 0.5, 0))

Arguments

method	name of method for generating the sequence. Must be one of "b90", "linear", "Coupmodel"
year	years to be returned
maxlai	maximum value during summer
winlaifrac	fraction of maxlai during winter (ignored when method = 'linear')
budburst.doy	budburst day of year (ignored when method = 'linear')
leaffall.doy	day of year when leaf fall begins (ignored when method = 'linear')
emerge.dur	number of days from budburst until maximum leaf area index is reached
leaffall.dur	number of days until minimum leaf are index is reached
shape.budburst	shape parameter for the growth phase (required when method = "Coupmodel")
shape.leaffall	shape parameter growth cessation (required when method = "Coupmodel")
lai.doy	integer vector of days of years
lai.frac	vector of values of fractional leaf area index corresponding to lai.doy (required when method = "linear")
opt.doy	day of year when optimum value is reached (required when method = "Coupmodel")

Value

a vector of daily lai values covering the years specified

Examples

# leaf area index for a single year

param.b90 <- MakeParam.B90()

plot(MakeSeasLAI(method = "linear",
year = 2001,
maxlai = param.b90$maxlai,
winlaifrac = param.b90$winlaifrac,
budburst.doy = param.b90$budburstdoy,
leaffall.doy = param.b90$leaffalldoy,
emerge.dur = param.b90$emergedur,
leaffall.dur = param.b90$leaffalldur))

plot(MakeSeasLAI(method = "Coupmodel",
                 year = 2001:2003,
                 maxlai = param.b90$maxlai,
                 winlaifrac = param.b90$winlaifrac,
                 budburst.doy = param.b90$budburstdoy,
                 leaffall.doy = param.b90$leaffalldoy,
                 leaffall.dur = param.b90$leaffalldur,
                 shape.optdoy = 180,
                 shape.budburst = 0.5,
                 shape.leaffall = 5))

# leaf area index for multiple years, using multiple values
plot(MakeSeasLAI(method = "linear",
year = 2001:2003,
maxlai = c(6,4,5),
winlaifrac = param.b90$winlaifrac,
budburst.doy = param.b90$budburstdoy,
leaffall.doy = param.b90$leaffalldoy,
emerge.dur = param.b90$emergedur,
leaffall.dur = param.b90$leaffalldur))

pschmidtwalter/brook90r documentation built on April 6, 2020, 6:35 p.m.