Rmiscanmod(data, RUE = 2.4, LER = 0.01, Tb = 10, k = 0.67, LAIdrd = 0.8, LAIStop = 1.8, RUEdrd = 1.3, RUEStop = 2.5, SMDdrd = -30, SMDStop = -120, FieldC = 45, iWatCont = 45, a = 6682.2, b = -0.33, soildepth = 0.6)
Simple model to calculate crop growth and yield based on MISCANMOD (see references).
The data.frame or matrix should contain
column 1: year column 2: month column 3: day column 4: JD column 5: max T (Celsius) column 6: min T (Celsius) column 7: PPFD or solar radiation divided by 2 (MJ/m2) column 8: Potential evaporation column 9: precip (mm)
Clifton-Brown, J. C.; Neilson, B.; Lewandowski, I. and Jones, M. B. The modelled productivity of Miscanthus x giganteus (GREEF et DEU) in Ireland. Industrial Crops and Products, 2000, 12, 97-109.
Clifton-brown, J. C.; Stampfl, P. F. and Jones, M. B. Miscanthus biomass production for energy in Europe and its potential contribution to decreasing fossil fuel carbon emissions. Global Change Biology, 2004, 10, 509-518.
## Need to get an example data set and then run it ## <strong>Not run</strong>: # data(WD1979) # # res <- Rmiscanmod(WD1979) # # ## convert to Mg/ha # # Yld <- res$Yield / 100 # # xyplot(Yld ~ 1:365 , # xlab='doy', # ylab='Dry biomass (Mg/ha)') # # ## although the default value for Field Capacity is 45 # ## a more reasonable value is closer to 27 # # # ## <strong>End(Not run)</strong>