tsirSpat function

Function to simulate the spatially-extended seasonally-forced TSIR