My mistake, I meant to write
y_0 = (f(x)-f(x_c))sin(\theta_2)cos(\phi_2)+ y_{c_0}
At x=x_c this gives y_0=y_{c_0} as desired.
The errors should be independent in this case. I’m measuring a vector field and estimating a spherical uncertainty from measurements in multiple orientations.