CHANGES:

1. The posterior NGM is now not saved, only the sum over all destinations to give $R_{it}$
2. Replace Rt as the dominant eigenvalue of the NGM with weighted average by population size.
3. ngm.nc now contains both $R_{it}$ and $R_t$ posterior predictive values.