Gauss2lats

Estimates the latitudes of a Gaussian grid.

Spectral atmospheric models provide data on a Gaussian grid that has regular longitude spacing but irregular latitude spacing. Gauss2lats provides the latitudes for such data. This information is needed, for example, when making a contour map.

The function determines the Gaussian latitudes by finding the roots of the ordinary Legendre polynomial of degree NLAT using Newton's iteration method. The only input required is the number of latitude lines in the Gaussian grid. The outputs are Gaussian latitudes, latitude spacing, and cosine of latitude.