LINTUL3: LINTUL3 crop growth model

In Rlintul: R implementation of LINTUL crop growth simulation models

Description

LINTUL3 simulates the growth of a crop as function of intercepted radiation, temperature and light use efficiency. Soil water (free drainage) and simple nitrogen, phosphorus and potassium balances are simulated and also the effects of water and N, P and K supplies on crop growth.

LINTUL3 was developed by Joost Wolf at Wageningen University.

Usage

lintul3(crop, soil, control, weather)
lintul3_crop(name)
lintul3_soil(name='p1')
lintul3_control()
## S4 replacement method for signature 'Rcpp_Lintul3Model'
crop(x)<-value
## S4 replacement method for signature 'Rcpp_Lintul3Model'
soil(x)<-value
## S4 replacement method for signature 'Rcpp_Lintul3Model'
control(x)<-value
## S4 replacement method for signature 'Rcpp_Lintul3Model'
weather(x)<-value
## S4 method for signature 'Rcpp_Lintul3Model'
run(x, ...)

Arguments

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

Value

matrix

Examples

library(Rlintul)
wth <- example_weather()
crop <- lintul3_crop('POTATO')
soil <- lintul3_soil(name='p1')

contr <- lintul3_control()
contr$emergence <- as.Date('2012-04-01')
contr$start <- contr$emergence - 10
contr$PL <- FALSE
contr$DAYPL <- contr$emergence
contr$IOPT <- 2
# IOPT 1=potential, 2=water limited, 3=water and N limited and 4=water and N, P and K limited

x <- lintul3(crop, soil, contr, wth)


## Not run: 

z <- run(x)

#or 
x$run()
x$out


# test LINTUL3

wf <- system.file("extdata/Netherlands_Wageningen.csv", package="meteor")
wth <- read.csv(wf)
wth$date <- as.Date(wth$date)
wth$srad <- wth$srad / 1000
crop <- lintul3_crop('WWHEAT')
soil <- lintul3_soil('p1')
sdate <- dateFromDoy(1, 1969)
contr <- lintul3_control()
contr$emergence=sdate
contr$start = sdate
contr$long_output=TRUE
contr$PL = FALSE


contr$IOPT = 1
contr.IRRI = 1
f1 <- system.file("lintul/test/3/WAGI1PTS.OUT", package="Rlintul")
s1 <- readLIN3output(f1)
x <- lintul3(crop, wth, soil, contr)
plot(s1[,'DAY'], s1[,'LAI'], type='l')
points(x[, 'step']+1, x[, 'LAI'])

contr$IOPT = 2
contr.IRRI = 1
f2 <- system.file("lintul/test/3/WAGW2PTS.OUT", package="Rlintul")
s2 <- readLIN3output(f2)
x2 <- lintul3(crop, wth, soil, contr)
plot(s2[,'DAY'], s2[,'LAI'], type='l')
points(x2[, 'step']+1, x2[, 'LAI'])

contr$IOPT = 3
f3 <- system.file("lintul/test/3/WAGW3PTS.OUT", package="Rlintul")
s3 <- readLIN3output(f3)
x3 <- lintul3(crop, wth, soil, contr)
plot(s3[,'DAY'], s3[,'LAI'], type='l')
points(x3[, 'step']+1, x3[, 'LAI'])

contr$IOPT = 4
f4 <- system.file("lintul/test/3/WAGW4PTS.OUT", package="Rlintul")
s4 <- readLIN3output(f4)
x4 <- lintul3(crop, wth, soil, contr)
plot(s4[,'DAY'], s4[,'LAI'], type='l')
points(x4[, 'step']+1, x4[, 'LAI'])

## End(Not run)

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