LINTUL1: LINTUL1 crop growth model

In Rlintul: R implementation of LINTUL crop growth simulation models

Description

LINTUL1 is a relatively simple crop growth simulation model. It computes potential production.

Usage

lintul1(crop, control, weather)
lintul1_crop()
## S4 replacement method for signature 'Rcpp_Lintul1Model'
crop(x)<-value
## S4 replacement method for signature 'Rcpp_Lintul1Model'
control(x)<-value
## S4 replacement method for signature 'Rcpp_Lintul1Model'
weather(x)<-value
## S4 method for signature 'Rcpp_Lintul1Model'
run(x, ...)

Arguments

x	Rcpp_Lintul1Model object
crop	list with named crop parameters. See Details. An example is returned by lintul1_crop()
control	list with named control parameters
weather	data.frame with daily weather data
value	a list with crop or control parameters, or a data.frame with weather data. As above
...	additional arguments. None implemented

Value

matrix

References

Spitters, C.J.T. and A.H.C.M. Schapendonk, 1990. Evaluation of breeding strategies for drought tolerance in potato by means of crop growth simulation. Plant and Soil 123: 193-203.

Examples

library(Rlintul)
w <- example_weather()
crop <- lintul1_crop()
contr <- list(emergence=as.Date('2012-04-01'), maxdur=200)

m <- lintul1(crop, contr, w)
# or
m <- lintul1()
crop(m) <- crop
control(m) <- contr
weather(m) <- w

x <- run(m)
head(x)

# or 
m$run()
m


# change parameters
m$crop$RUE <- 3.5
m$control$emergence=as.Date('2012-04-01')
m$run()

out <- m$out
out


date <- as.Date(m$control$emergence, origin="1970-01-01") + out$step
par(mfrow=c(1,2))
plot(date, out$LAI)

prod <- cbind(out$WRT, out$WLV, out$WST, out$WSO)

wtot <- rowSums(prod)
plot(date, wtot, type='l', lwd=2)
lines(date, prod[,2], col='green', lwd=2)
lines(date, prod[,4], col='red', lwd=2)


# test LINTUL1

library(meteor)
library(Rlintul)
f1 <- system.file("lintul/test/1/res.dat", package="Rlintul")
s <- readLIN1output(f1)
f2 <- system.file("extdata/Netherlands_Wageningen.csv", package="meteor")
wth <- read.csv(f2)
wth$date <- as.Date(wth$date)
wth$srad <- wth$srad / 1000

crop <- lintul1_crop()
sdate <- dateFromDoy(58, 1971)
contr <- list(emergence=sdate+2, start = sdate, maxdur=365)
x <- run(lintul1(crop, contr, wth))
x$DOY <- doyFromDate(x$date)

plot(s[,'TIME'], s[,'LAI'], type='l')
points(x[, 'DOY'], x[, 'LAI'])

Rlintul documentation built on May 2, 2019, 5:34 p.m.