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