Simulates canopy assimilation (von Caemmerer model)

Usage

eCanA(LAI, doy, hour, solarR, AirTemp, RH, WindS, Vcmax, Vpmax, Vpr, Jmax, Ca = 380, Oa = 210, StomataWS = 1)

Arguments

LAI
leaf area index.
doy
day of the year, (1--365).
hour
hour of the day, (0--23).
solarR
solar radiation (micro mol /m2/s).
AirTemp
temperature (Celsius).
RH
relative humidity (0--1).
WindS
wind speed (m/s).
Vcmax
Maximum rubisco activity (micro mol /m2/s).
Vpmax
Maximum PEP carboxylase activity (micro mol /m2/s).
Vpr
PEP regeneration rate (micro mol /m2/s).
Jmax
Maximal electron transport rate (micromol electrons m-2 s-1).
Ca
atmospheric carbon dioxide concentration (ppm or microbar) (e.g. 380).
Oa
atmospheric oxygen concentration (mbar) (e.g. 210).
StomataWS
Effect of water stress on assimilation.

Value

numeric

returns a single value which is hourly canopy assimilation (mol m-2 ground s-1)

Description

It represents an integration of the photosynthesis function eC4photo, canopy evapo/transpiration and the multilayer canopy model sunML.

Examples

## <strong>Not run</strong>: # data(doy124) # tmp1 <- numeric(24) # for(i in 1:24){ # lai <- doy124[i,1] # doy <- doy124[i,3] # hr <- doy124[i,4] # solar <- doy124[i,5] # temp <- doy124[i,6] # rh <- doy124[i,7]/100 # ws <- doy124[i,8] # # tmp1[i] <- CanA(lai,doy,hr,solar,temp,rh,ws) # } # # plot(c(0:23),tmp1, # type='l',lwd=2, # xlab='Hour', # ylab=expression(paste('Canopy assimilation (mol ', # m^-2,' ',s^-1,')'))) # ## <strong>End(Not run)</strong>