the pointing constants:

  Equation 4 fits tilt and sin and cos az and 2az terms
         in a different order:

         daz = v(1)*cosel + sinel*(v(2)*sin(az) + v(3)*cos(az))
             + v(4) + v(5)*sinel + cosel*(v(6)*sin(az) + v(7)*cos(az))
                             + cosel*(v(8)*sin(2*az) + v(9)*cos(2*az))

         del = v(1) + v(2)*sin(az)) + v(3)*cos(az) + v(4)*sin(el)
             + v(5)*cos(el) + v(6)*sin(2*az) + v(7)*cos(2az)
             + v(8)*tan(1.5708-el)



    if bd != 0: bd.setPointingModelCoefs( \
            [-6.00, -0.52, -0.24, 31.47, -1.09, 0.06, 0.05, -0.08,
0.10], \
            [-83.29, 0.10, -0.46, -0.44, -0.29, 0.04, 0.01, 0.60, 0])
    bd=Bima(3).drive()
    if bd != 0: bd.setPointingModelCoefs( \
            [-7.47, 0.46, 0.16, 5.05, 0.28, 0.01, -0.04, 0.04, -0.01 ],
\
            [-0.24, -0.02, 0.53, 0.13, 0.56, 0.09, 0.0, 0.70, 0.0])

    bd=Bima(9).drive()
    if bd != 0: bd.setPointingModelCoefs( \
            [-10.19, 0.24, -0.20, -1.22, 0.01, 0.07, 0.00, -0.04,
0.06 ], \
            [-7.23, -0.06, 0.27, -0.90, -0.08, 0.03, 0.01, 0.70, 0.0])

First line of a pair is the az coefficients; bottom is el.