Description Usage Arguments Details Value References See Also Examples
Run the WOFOST crop growth model. Through this interface, you provide weather data, and crop, soil and control parameters to run the model once. For multiple runs it might be preferable to use wofost_model
instead.
1 | wofost(crop, weather, soil, control)
|
crop |
list. Crop parameters |
weather |
data.frame with weather data |
soil |
list. Soil parameters |
control |
list. Model control options |
The weather data must be passed as a data.frame with the following variables and units.
variable | description | class / unit | |
date | "Date" class variable | - | |
srad | Solar radiation | kJ m-2 day-1 | |
tmin | Minimum temperature | degrees C | |
tmax | Maximum temperature | degrees C | |
vapr | Vapor pressure | kPa | |
wind | Wind speed | m s-1 | |
prec | Precipitation | mm day-1 | |
Note that there should not be any time gaps between the days in the data.frame
matrix
Van Diepen, C.A., J. Wolf, and H van Keulen, 1989. WOFOST: a simulation model of crop production. Soil Use and Management, 5: 16-24
Van Keulen, H. and J. Wolf, 1986. Modelling of agricultural production : weather, soils and crops. http://edepot.wur.nl/168025
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 | # weather data
f <- system.file("extdata/Netherlands_Swifterbant.csv", package="meteor")
w <- read.csv(f)
w$date <- as.Date(w$date)
head(w)
# crop and soil parameters
crop <- wofost_crop("barley")
soil <- wofost_soil("ec1")
# "control" parameters
contr <- wofost_control()
contr$modelstart <- as.Date("1980-02-06")
contr$latitude=52.57
contr$elevation=50
# run model
d <- wofost(crop, w, soil, contr)
# output
head(d)
tail(d)
plot(d[,"step"], d[, "LAI"])
## Another example
crop <- wofost_crop("rapeseed_1001")
soil <- wofost_soil("soil_5")
contr$modelstart <- as.Date("1977-01-01")
rp <- wofost(crop, w, soil, contr)
plot(rp[,"step"], rp[, "LAI"])
# yield
plot(rp[, 1], rp[,"WSO"])
## water limited
contr$water_limited <- TRUE
contr$modelstart <- as.Date("1985-01-01")
crop <- wofost_crop("maize_1")
f <- system.file("extdata/Philippines_IRRI.csv", package="meteor")
wth <- read.csv(f)
wth$date <- as.Date(wth$date)
contr$elevation <- 21
contr$latitude <- 14.18
ma <- wofost(crop, wth, soil, contr)
plot(ma[,"step"], ma[, "LAI"])
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.