if indAlphaM2SD:
        plt.plot(s[indAlphaM2SD], z[indAlphaM2SD], 'x', markersize=8,