Simulation models

Statistical vs mechanistic models

Models

LINTUL

WOFOST

Here is an example

Input data

library(Rwofost)
## Loading required package: meteor
##
## Attaching package: 'meteor'
## The following objects are masked from 'package:luna':
##
##     dateFromDoy, monthFromDate, yearFromDate
w <- example_weather()
## Error in example_weather(): could not find function "example_weather"
# radiation should be in kJ
w$srad <- w$srad * 1000
## Error in eval(expr, envir, enclos): object 'w' not found
head(w)
## Error in head(w): object 'w' not found

Crop parameters

crop <- wofost_crop('barley')

Soil parameters

soil <- wofost_soil('ec1')
contr <- wofost_control()
contr$modelstart <- as.Date("2012-02-06")

Run model

d <- wofost(crop, w, soil, contr)
## Error in wofost(crop, w, soil, contr): object 'w' not found

Inspect output

head(d)
## [1] "statistical"
tail(d)
## [1] "statistical"
plot(d[,'step'], d[, 'LAI'])
## Error in d[, "step"]: incorrect number of dimensions

APSIM

DSSAT