fplus1SD = z[0] + np.tan(np.deg2rad(-alphaSD[0])) * s