I am using FORTRAN to implement Gillespie's stochastic simulation algorithm. I would be running many simulations in parallel (both parallel instances with different seed and parallel functions); if I have to store every time interval then a lot of memory would be consumed. So I thought I can have some fixed number of time points where the data would be saved (Allocate space beforehand either dynamically or statically).How do I decide the size of the time intervals? Should the intervals have a fixed size or it should be based on local dynamics(some regions would have more dynamics)? Is there any way to predict that based on the deterministic function without actually running a deterministic simulation? I have this vague idea that the tau-leap formula can be used to predict time intervals where n number of a major reaction would occur. But this still doesn't help me predict the number of steps beforehand. 
How to fix time intervals to store data in a stochastic simulation (continous time markov chain)
algorithms;fortran;stochastic;memory management