- [x] replace `localSpp.size()` with `J[step]` - [x] update `death` method to know about `step` - [x] in each step, check for change in `J` - [x] if `J` is smaller, subsample local comm - [x] if `J` is bigger, run `birthImm` and `speciation` until get the right size - [x] calculate the contribution of each step toward generations during every step - [ ] unit tests
localSpp.size()withJ[step]deathmethod to know aboutstepJJis smaller, subsample local commJis bigger, runbirthImmandspeciationuntil get the right size