indAlphaM2SD = np.argwhere(np.diff(np.sign(fminus2SD - z))).flatten()