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