  We study orientational order, subject to thermal fluctuations, on a fixed
curved surface. We derive, in particular, the average density of zeros of
Gaussian distributed vector fields on a closed Riemannian manifold. Results are
compared with the density of disclination charges obtained from a Coulomb gas
model. Our model describes the disordered state of two dimensional objects with
orientational degrees of freedom, such as vector ordering in Langmuir
monolayers and lipid bilayers above the hexatic to fluid transition.
