A dynamic model of succession on alp Stabelchod in the Swiss Nationl Park using differential equations and numerial integration. 6 species guilds are considered. Space is conceived as a grid of 30 times 40 cells. Typical simulation time is around 500yr.
SNPsm(trange,tsl,diff,r6,...)
SNPsm2(trange=100,tsl=5.0,diff=0.001,r6=NULL)# S3 method for default
SNPsm(trange, tsl, diff, r6, ...)
# S3 method for SNPsm
plot(x, ...,out.seq=1,col=FALSE)
Time range of simulation in yr
Time range of simulation in yr
Time interval (yr) at which maps of the state are printed
A diffusion coefficient driving random spatial propagation
Growth rates of 6 guilds involved, increase in cover percentage per yr
Parameter out.seq, the plotting interval
An object of class "SNPsm"
A logical variable to suppress color printing
An object of class "SNPsm" with at least the following items:
Number of time steps used for numerical integration
Vertical grid count
Horizontal grid count
The time step length in yr
The names of the vegetation types, i.e., the species
A nspecies x nspecies matrix defining composition of vegetation types
The growth rates given upon input
Simulated scores of all species (guilds) during simulation time
The 30 x 40 grid map of types used as initial condition
The same as tmap
A 30 x 40 grid showing initial forest edges, used for printing
Wildi, O. 2002. Modeling succession from pasture to forest in time and space. Community Ecology 3: 181--189.
Wildi, O. 2017. Data Analysis in Vegetation Ecology. 3rd ed. CABI, Oxfordshire, Boston.
# NOT RUN {
r6=NULL # imposes default growth rates
o.stSNP<- SNPsm(trange=100,tsl=10.0,diff=0.001,r6)
plot(o.stSNP,out.seq=50)
# }
Run the code above in your browser using DataLab