eCanA: Simulates canopy assimilation (von Caemmerer model)

In serbinsh/biocro: Simulation of Energycane and Sugarcane

Description

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

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)

Examples

## Not run: 
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,')')))

## End(Not run)

serbinsh/biocro documentation built on May 29, 2019, 6:57 p.m.