indAlphaP1SD = np.argwhere(np.diff(np.sign(fplus1SD - z))).flatten()