-
Notifications
You must be signed in to change notification settings - Fork 106
/
pymakegrid2d.1
138 lines (138 loc) · 4.68 KB
/
pymakegrid2d.1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
.\" Automatically generated by Pandoc 1.17.2
.\"
.TH "pymakegrid2d" "1" "2016\-07\-27" "Python" "SHTOOLS 3.3"
.hy
.SH MakeGrid2D
.PP
Create a 2D cylindrical map of arbitrary grid spacing from a set of
spherical harmonic coefficients.
.SH Usage
.PP
\f[C]grid\f[] = pyshtools.MakeGrid2D (\f[C]cilm\f[], \f[C]interval\f[],
[\f[C]lmax\f[], \f[C]norm\f[], \f[C]csphase\f[], \f[C]f\f[], \f[C]a\f[],
\f[C]north\f[], \f[C]south\f[], \f[C]east\f[], \f[C]west\f[],
\f[C]dealloc\f[]])
.SH Returns
.TP
.B \f[C]grid\f[] : float, dimension ((\f[C]north\f[]\-\f[C]south\f[])/\f[C]interval\f[]+1, (\f[C]east\f[]\-\f[C]west\f[])/\f[C]interval\f[]+1)
A 2D equally spaced map of the input spherical harmonic coefficients
\f[C]cilm\f[].
The array is in raster format with upper\-left and lower\-right
coordinates of (90 N, 0 E) and (90 S, 360 E), respectively.
.RS
.RE
.SH Parameters
.TP
.B \f[C]cilm\f[] : float, dimension (2, \f[C]lmaxin\f[]+1, \f[C]lmaxin\f[]+1)
The real spherical harmonic coefficients to be expanded in the space
domain.
The coefficients \f[C]C1lm\f[] and \f[C]C2lm\f[] refer to the cosine
(\f[C]Clm\f[]) and sine (\f[C]Slm\f[]) coefficients, respectively, with
\f[C]Clm=cilm[0,l,m]\f[] and \f[C]Slm=cilm[1,l,m]\f[].
.RS
.RE
.TP
.B \f[C]interval\f[] : float
The latitudinal and longitudinal spacing of \f[C]grid\f[].
.RS
.RE
.TP
.B \f[C]lmax\f[] :optional, integer, default = \f[C]lmaxin\f[]
The maximum spherical harmonic degree of the coefficients \f[C]cilm\f[]
used when calculating the grid.
.RS
.RE
.TP
.B \f[C]norm\f[] : optional, integer, default = 1
1 (default) = 4\-pi (geodesy) normalized harmonics; 2 = Schmidt
semi\-normalized harmonics; 3 = unnormalized harmonics; 4 = orthonormal
harmonics.
.RS
.RE
.TP
.B \f[C]csphase\f[] : 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.
.RS
.RE
.TP
.B \f[C]f\f[] : optional, float, default = 0
The flattening of the reference ellipoid that is subtracted from the
function.
This is given by (\f[C]R_equator\-R_pole)/R_equator\f[].
.RS
.RE
.TP
.B \f[C]a\f[] : optional, float, default = 0
The semi\-major axis of the reference ellispoid that is subtracted from
the function.
.RS
.RE
.TP
.B \f[C]north\f[] : optional, float, default = 90
The maximum latitude of the output raster grid, in degrees.
The default is 90 degrees.
.RS
.RE
.TP
.B \f[C]south\f[] : optional, float, default = \-90
The minimum latitude of the output raster grid, in degrees.
The default is \-90 degrees.
.RS
.RE
.TP
.B \f[C]east\f[] : optional, float, default = 360
The maximum longitude of the output raster grid, in degrees.
The default is 360 degrees.
.RS
.RE
.TP
.B \f[C]west\f[] : optional, float, default = 0
The minimum longitude of the output raster grid, in degrees.
The default is 0 degrees.
.RS
.RE
.TP
.B \f[C]dealloc\f[] : 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.
.RS
.RE
.SH Description
.PP
\f[C]MakeGrid2D\f[] 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 \f[C]north\f[], \f[C]south\f[], \f[C]east\f[]
and \f[C]west\f[] are specified, then the output grid will possess
upper\-left and lower\-right coordinates of (\f[C]north\f[],
\f[C]west\f[]) and (\f[C]south\f[], \f[C]east\f[]), repsectively.
The employed spherical harmonic normalization and Condon\-Shortley phase
convention can be set by the optional arguments \f[C]norm\f[] and
\f[C]csphase\f[]; if not set, the default is to use geodesy 4\-pi
normalized harmonics that exclude the Condon\-Shortley phase of (\-1)^m.
.PP
If the optional arguments \f[C]f\f[] and \f[C]a\f[] are specified, the
output function will be referenced to an ellipsoid with flattening
\f[C]f\f[] and semimajor axis \f[C]a\f[].
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.
.SH References
.PP
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.
.SH See also
.PP
\f[C]makegriddh\f[], \f[C]makegriddhc\f[], \f[C]makegridglq\f[],
\f[C]makegridglqc\f[], \f[C]makegravgriddh\f[], \f[C]makemaggriddh\f[]