lintul1 {Rlintul}R Documentation

LINTUL1 crop growth model

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'])


[Package Rlintul version 0.1-4 Index]