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