Create a 2D cylindrical map of arbitrary grid spacing from a set of spherical harmonic coefficients.

## Usage

call MakeGrid2D (grid, cilm, lmax, interval, nlat, nlong, norm, csphase, f, a, north, south, east, west, dealloc, exitstatus)

## Parameters

grid : output, real(dp), dimension (180/interval+1, 360/interval+1)
A 2D equally spaced map of the input spherical harmonic coefficients cilm. The array is in raster format with upper-left and lower-right coordinates of (90 N, 0 E) and (90 S, 360 E), respectively.
cilm : input, real(dp), dimension (2, lmax+1, lmax+1)
The real spherical harmonic coefficients to be expanded in the space domain. The coefficients C1lm and C2lm refer to the cosine (Clm) and sine (Slm) coefficients, respectively, with Clm=cilm(1,l+1,m+1) and Slm=cilm(2,l+1,m+1).
lmax : input, integer
The maximum spherical harmonic degree of the coefficients cilm used when calculating the grid.
interval : input, real(dp)
The latitudinal and longitudinal spacing of grid.
nlat : output, integer
The number of latitudinal samples. Both 90 N and 90 S are included.
nlong : output, integer
The number of longitudinal samples. Both 0 and 360 E are included.
norm : input, optional, integer, default = 1
1 (default) = 4-pi (geodesy) normalized harmonics; 2 = Schmidt semi-normalized harmonics; 3 = unnormalized harmonics; 4 = orthonormal harmonics.
csphase : input, optional, integer, default = 1
1 (default) = do not apply the Condon-Shortley phase factor to the associated Legendre functions; -1 = append the Condon-Shortley phase factor of (-1)^m to the associated Legendre functions.
f : input, optional, real(dp)
The flattening of the reference ellipoid that is subtracted from the function. This is given by (R_equator-R_pole)/R_equator. The semi-major axis a (i.e., R_equator) must be specified for this calculation.
a : input, optional, real(dp)
The semi-major axis of the reference ellispoid that is subtracted from the function. The flattening f must be specified for this calculation.
north : input, real*8, optional, default = 90
The maximum latitude of the output raster grid, in degrees. The default is 90 degrees.
south : input, optional, real(dp), default = -90
The minimum latitude of the output raster grid, in degrees. The default is -90 degrees.
east : input, optional, real(dp), default = 360
The maximum longitude of the output raster grid, in degrees. The default is 360 degrees.
west : input, optional, real(dp), default = 0
The minimum longitude of the output raster grid, in degrees. The default is 0 degrees.
dealloc : input, optional, integer, default = 0
0 (default) = Save variables used in the external Legendre function calls. (1) Deallocate this memory at the end of the funcion call.
exitstatus : output, optional, integer
If present, instead of executing a STOP when an error is encountered, the variable exitstatus will be returned describing the error. 0 = No errors; 1 = Improper dimensions of input array; 2 = Improper bounds for input variable; 3 = Error allocating memory; 4 = File IO error.

## Description

MakeGrid2D will create a 2-dimensional cylindrical map, equally spaced in (geocentric) latitude and longitude, from a set of input spherical harmonic coefficients. The output grid is in raster format possessing upper-left and lower-right coordinates of (90 N, 0 E) and (90 S, 360 E), repsectively. If the optional parameters north, south, east and west are specified, then the output grid will possess upper-left and lower-right coordinates of (north, west) and (south, east), repsectively. The employed spherical harmonic normalization and Condon-Shortley phase convention can be set by the optional arguments norm and csphase; if not set, the default is to use geodesy 4-pi normalized harmonics that exclude the Condon-Shortley phase of (-1)^m.

If the optional arguments f and a are specified, the output function will be referenced to an ellipsoid with flattening f and semimajor axis a. The normalized legendre functions are calculated in this routine using the scaling algorithm of Holmes and Featherstone (2002), which are accurate to about degree 2800. The unnormalized functions are accurate only to about degree 15.

## References

Holmes, S. A., and W. E. Featherstone, A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions, J. Geodesy, 76, 279-299, 2002.