Create a 2D cylindrical map of arbitrary grid spacing from a set of spherical harmonic coefficients.
call MakeGrid2D (
grid: output, real*8, dimension (180/
- 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*8, dimension (2,
- The real spherical harmonic coefficients to be expanded in the space domain. The coefficients
C2lmrefer to the cosine (
Clm) and sine (
Slm) coefficients, respectively, with
lmax: input, integer
- The maximum spherical harmonic degree of the coefficients
cilmused when calculating the grid.
interval: input, real*8
- The latitudinal and longitudinal spacing of
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*8
- 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
R_equator) must be specified for this calculation.
a: input, optional, real*8
- The semi-major axis of the reference ellispoid that is subtracted from the function. The flattening
fmust 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*8, default = -90
- The minimum latitude of the output raster grid, in degrees. The default is -90 degrees.
east: input, optional, real*8, default = 360
- The maximum longitude of the output raster grid, in degrees. The default is 360 degrees.
west: input, optional, real*8, 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.
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
west are specified, then the output grid will possess upper-left and lower-right coordinates of (
west) and (
east), repsectively. The employed spherical harmonic normalization and Condon-Shortley phase convention can be set by the optional arguments
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
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.
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.