Fix mbi_threshold to have sigma in [s}
Add plot_critical_mass to instabilities (ions)
Add sigma(position) to synchrotron to compute the beam size at any ring location