XYZrot(1) =
XYZrot(2) =
XYZrot(3) = SIN( Y(I,J,bi,bj)*deg2rad )
XYZgeo(iA) = XYZgeo(iA) + Ainv(iA,jA)*XYZrot(jA)
_RL XYZgeo(3), XYZrot(3)