indAlphaM1SD = np.argwhere(np.diff(np.sign(fminus1SD - z))).flatten()