lintul3 {Rlintul} | R Documentation |
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.
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, ...)
x |
Rcpp_Lintul3Model object |
crop |
list with named crop parameters. See Details. An example is returned by |
soil |
list with named soil parameters. See Details. An example is returned by |
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 |
matrix
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)