An integral equation theory for water around globular solutes is described that gives the pair distribution function with spatial and orientational dependence, g(r, Ω). After a decomposition of the position and orientation dependence, an HNC-OZ theory is used for the position-dependent part and a probability function approximation for the orientation-dependent part. For a model system of water around frozen water clusters, the theory yields results in good agreement with those from Monte Carlo computer simulations for the dipolar hard-sphere plus sticky potential model for water by Bratko, Blum and Luzar. Extensions to other water models are also described.