Added variable "plotyrs" to select a subset of 9 years to show within season dynamics
#1) change ode solver (now lsodar) so it ends the season when R is numerically 0
2) Fixed up within and between year plots to work with the ode output
3) make g (germination fraction) a nyrs x nsp matrix that is defined in getSpecies.R
germination fraction, g, is now an nyrs x nsp matrix (rather than calculated on the fly each year) this makes it easier to save and plot later
Added a rootfun to call from lsodar (pkg: deSolve) so that the within year simulation ends when R is numerically 0