  We present a new, elementary way to obtain axially symmetric Gaussian
processes on the sphere, in order to accommodate for the directional anisotropy
of global climate data in geostatistical analysis.
