wofost {Rwofost} | R Documentation |
WOFOST crop growth model
Description
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.
Usage
wofost(crop, weather, soil, control)
Arguments
crop |
list. Crop parameters |
weather |
data.frame with weather data |
soil |
list. Soil parameters |
control |
list. Model control options |
Details
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
Value
matrix
References
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
See Also
Examples
# 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"])
[Package Rwofost version 0.8-3 Index]