# NOT RUN {
#Load example daily meteorological data
data(examplemeteo)
#Load example plot plant data
data(exampleforestMED)
#Default species parameterization
data(SpParamsMED)
#Initialize soil with default soil params (2 layers)
examplesoil = soil(defaultSoilParams(2))
#Initialize control parameters
control = defaultControl("Granier")
#Initialize input
x = forest2spwbInput(exampleforestMED,examplesoil, SpParamsMED, control)
#Call simulation function
S1<-spwb(x, examplemeteo, latitude = 41.82592, elevation = 100)
#Monthly summary (averages) of soil status
summary(S1, freq="months",FUN=mean, output="Soil")
# }
Run the code above in your browser using DataLab