1 | ;+ |
---|
2 | ; |
---|
3 | ; @file_comments |
---|
4 | ; compute the latitudes of the n256 gaussian grid. See: |
---|
5 | ; <a href="http://www.ecmwf.int/products/data/technical/gaussian/n256FIS.html">Gaussian 256</a> |
---|
6 | ; |
---|
7 | ; @categories |
---|
8 | ; Grid |
---|
9 | ; |
---|
10 | ; @examples |
---|
11 | ; IDL> lat=n256gaussian() |
---|
12 | ; |
---|
13 | ; @returns |
---|
14 | ; a 1d array |
---|
15 | ; |
---|
16 | ; @history |
---|
17 | ; Sebastien Masson (smasson\@lodyc.jussieu.fr) |
---|
18 | ; June 2004 |
---|
19 | ; |
---|
20 | ; @version |
---|
21 | ; $Id$ |
---|
22 | ; |
---|
23 | ;- |
---|
24 | ; |
---|
25 | FUNCTION n256gaussian |
---|
26 | ; |
---|
27 | compile_opt idl2, strictarrsubs |
---|
28 | ; |
---|
29 | ; latitude reduced regular latitude |
---|
30 | ; number points points |
---|
31 | ; ------- ------- ------- -------- |
---|
32 | n256 = [ 1, 18, 1024, 89.73115, $ |
---|
33 | 2, 25, 1024, 89.38287, $ |
---|
34 | 3, 32, 1024, 89.03254, $ |
---|
35 | 4, 40, 1024, 88.68175, $ |
---|
36 | 5, 45, 1024, 88.33077, $ |
---|
37 | 6, 50, 1024, 87.97972, $ |
---|
38 | 7, 60, 1024, 87.62861, $ |
---|
39 | 8, 64, 1024, 87.27748, $ |
---|
40 | 9, 72, 1024, 86.92632, $ |
---|
41 | 10, 72, 1024, 86.57515, $ |
---|
42 | 11, 75, 1024, 86.22398, $ |
---|
43 | 12, 81, 1024, 85.87279, $ |
---|
44 | 13, 90, 1024, 85.52160, $ |
---|
45 | 14, 96, 1024, 85.17041, $ |
---|
46 | 15, 100, 1024, 84.81921, $ |
---|
47 | 16, 108, 1024, 84.46801, $ |
---|
48 | 17, 120, 1024, 84.11681, $ |
---|
49 | 18, 120, 1024, 83.76560, $ |
---|
50 | 19, 125, 1024, 83.41440, $ |
---|
51 | 20, 135, 1024, 83.06319, $ |
---|
52 | 21, 144, 1024, 82.71198, $ |
---|
53 | 22, 150, 1024, 82.36077, $ |
---|
54 | 23, 160, 1024, 82.00956, $ |
---|
55 | 24, 160, 1024, 81.65835, $ |
---|
56 | 25, 180, 1024, 81.30714, $ |
---|
57 | 26, 180, 1024, 80.95593, $ |
---|
58 | 27, 180, 1024, 80.60471, $ |
---|
59 | 28, 192, 1024, 80.25350, $ |
---|
60 | 29, 192, 1024, 79.90229, $ |
---|
61 | 30, 200, 1024, 79.55107, $ |
---|
62 | 31, 216, 1024, 79.19986, $ |
---|
63 | 32, 216, 1024, 78.84864, $ |
---|
64 | 33, 216, 1024, 78.49743, $ |
---|
65 | 34, 225, 1024, 78.14621, $ |
---|
66 | 35, 240, 1024, 77.79500, $ |
---|
67 | 36, 240, 1024, 77.44378, $ |
---|
68 | 37, 243, 1024, 77.09256, $ |
---|
69 | 38, 250, 1024, 76.74135, $ |
---|
70 | 39, 256, 1024, 76.39013, $ |
---|
71 | 40, 270, 1024, 76.03891, $ |
---|
72 | 41, 270, 1024, 75.68770, $ |
---|
73 | 42, 288, 1024, 75.33648, $ |
---|
74 | 43, 288, 1024, 74.98526, $ |
---|
75 | 44, 288, 1024, 74.63405, $ |
---|
76 | 45, 300, 1024, 74.28283, $ |
---|
77 | 46, 300, 1024, 73.93161, $ |
---|
78 | 47, 320, 1024, 73.58040, $ |
---|
79 | 48, 320, 1024, 73.22918, $ |
---|
80 | 49, 320, 1024, 72.87796, $ |
---|
81 | 50, 324, 1024, 72.52674, $ |
---|
82 | 51, 360, 1024, 72.17552, $ |
---|
83 | 52, 360, 1024, 71.82431, $ |
---|
84 | 53, 360, 1024, 71.47309, $ |
---|
85 | 54, 360, 1024, 71.12187, $ |
---|
86 | 55, 360, 1024, 70.77065, $ |
---|
87 | 56, 360, 1024, 70.41944, $ |
---|
88 | 57, 375, 1024, 70.06822, $ |
---|
89 | 58, 375, 1024, 69.71700, $ |
---|
90 | 59, 384, 1024, 69.36578, $ |
---|
91 | 60, 384, 1024, 69.01456, $ |
---|
92 | 61, 400, 1024, 68.66334, $ |
---|
93 | 62, 400, 1024, 68.31213, $ |
---|
94 | 63, 400, 1024, 67.96091, $ |
---|
95 | 64, 432, 1024, 67.60969, $ |
---|
96 | 65, 432, 1024, 67.25847, $ |
---|
97 | 66, 432, 1024, 66.90725, $ |
---|
98 | 67, 432, 1024, 66.55603, $ |
---|
99 | 68, 432, 1024, 66.20482, $ |
---|
100 | 69, 450, 1024, 65.85360, $ |
---|
101 | 70, 450, 1024, 65.50238, $ |
---|
102 | 71, 450, 1024, 65.15116, $ |
---|
103 | 72, 480, 1024, 64.79994, $ |
---|
104 | 73, 480, 1024, 64.44872, $ |
---|
105 | 74, 480, 1024, 64.09750, $ |
---|
106 | 75, 480, 1024, 63.74629, $ |
---|
107 | 76, 480, 1024, 63.39507, $ |
---|
108 | 77, 486, 1024, 63.04385, $ |
---|
109 | 78, 500, 1024, 62.69263, $ |
---|
110 | 79, 500, 1024, 62.34141, $ |
---|
111 | 80, 500, 1024, 61.99019] |
---|
112 | n256 = [n256, $ |
---|
113 | 81, 512, 1024, 61.63897, $ |
---|
114 | 82, 512, 1024, 61.28776, $ |
---|
115 | 83, 540, 1024, 60.93654, $ |
---|
116 | 84, 540, 1024, 60.58532, $ |
---|
117 | 85, 540, 1024, 60.23410, $ |
---|
118 | 86, 540, 1024, 59.88288, $ |
---|
119 | 87, 540, 1024, 59.53166, $ |
---|
120 | 88, 576, 1024, 59.18044, $ |
---|
121 | 89, 576, 1024, 58.82922, $ |
---|
122 | 90, 576, 1024, 58.47800, $ |
---|
123 | 91, 576, 1024, 58.12679, $ |
---|
124 | 92, 576, 1024, 57.77557, $ |
---|
125 | 93, 576, 1024, 57.42435, $ |
---|
126 | 94, 600, 1024, 57.07313, $ |
---|
127 | 95, 600, 1024, 56.72191, $ |
---|
128 | 96, 600, 1024, 56.37069, $ |
---|
129 | 97, 600, 1024, 56.01947, $ |
---|
130 | 98, 600, 1024, 55.66825, $ |
---|
131 | 99, 640, 1024, 55.31703, $ |
---|
132 | 100, 640, 1024, 54.96581, $ |
---|
133 | 101, 640, 1024, 54.61460, $ |
---|
134 | 102, 640, 1024, 54.26338, $ |
---|
135 | 103, 640, 1024, 53.91216, $ |
---|
136 | 104, 640, 1024, 53.56094, $ |
---|
137 | 105, 640, 1024, 53.20972, $ |
---|
138 | 106, 640, 1024, 52.85850, $ |
---|
139 | 107, 648, 1024, 52.50728, $ |
---|
140 | 108, 675, 1024, 52.15606, $ |
---|
141 | 109, 675, 1024, 51.80484, $ |
---|
142 | 110, 675, 1024, 51.45362, $ |
---|
143 | 111, 675, 1024, 51.10241, $ |
---|
144 | 112, 675, 1024, 50.75119, $ |
---|
145 | 113, 675, 1024, 50.39997, $ |
---|
146 | 114, 720, 1024, 50.04875, $ |
---|
147 | 115, 720, 1024, 49.69753, $ |
---|
148 | 116, 720, 1024, 49.34631, $ |
---|
149 | 117, 720, 1024, 48.99509, $ |
---|
150 | 118, 720, 1024, 48.64387, $ |
---|
151 | 119, 720, 1024, 48.29265, $ |
---|
152 | 120, 720, 1024, 47.94143, $ |
---|
153 | 121, 720, 1024, 47.59021, $ |
---|
154 | 122, 720, 1024, 47.23899, $ |
---|
155 | 123, 729, 1024, 46.88778, $ |
---|
156 | 124, 729, 1024, 46.53656, $ |
---|
157 | 125, 750, 1024, 46.18534, $ |
---|
158 | 126, 750, 1024, 45.83412, $ |
---|
159 | 127, 750, 1024, 45.48290, $ |
---|
160 | 128, 750, 1024, 45.13168, $ |
---|
161 | 129, 750, 1024, 44.78046, $ |
---|
162 | 130, 768, 1024, 44.42924, $ |
---|
163 | 131, 768, 1024, 44.07802, $ |
---|
164 | 132, 768, 1024, 43.72680, $ |
---|
165 | 133, 768, 1024, 43.37558, $ |
---|
166 | 134, 800, 1024, 43.02436, $ |
---|
167 | 135, 800, 1024, 42.67315, $ |
---|
168 | 136, 800, 1024, 42.32193, $ |
---|
169 | 137, 800, 1024, 41.97071, $ |
---|
170 | 138, 800, 1024, 41.61949, $ |
---|
171 | 139, 800, 1024, 41.26827, $ |
---|
172 | 140, 800, 1024, 40.91705, $ |
---|
173 | 141, 800, 1024, 40.56583, $ |
---|
174 | 142, 810, 1024, 40.21461, $ |
---|
175 | 143, 810, 1024, 39.86339, $ |
---|
176 | 144, 864, 1024, 39.51217, $ |
---|
177 | 145, 864, 1024, 39.16095, $ |
---|
178 | 146, 864, 1024, 38.80973, $ |
---|
179 | 147, 864, 1024, 38.45851, $ |
---|
180 | 148, 864, 1024, 38.10730, $ |
---|
181 | 149, 864, 1024, 37.75608, $ |
---|
182 | 150, 864, 1024, 37.40486, $ |
---|
183 | 151, 864, 1024, 37.05364, $ |
---|
184 | 152, 864, 1024, 36.70242, $ |
---|
185 | 153, 864, 1024, 36.35120, $ |
---|
186 | 154, 864, 1024, 35.99998, $ |
---|
187 | 155, 864, 1024, 35.64876, $ |
---|
188 | 156, 864, 1024, 35.29754, $ |
---|
189 | 157, 864, 1024, 34.94632, $ |
---|
190 | 158, 900, 1024, 34.59510, $ |
---|
191 | 159, 900, 1024, 34.24388, $ |
---|
192 | 160, 900, 1024, 33.89266] |
---|
193 | n256 = [n256, $ |
---|
194 | 161, 900, 1024, 33.54145, $ |
---|
195 | 162, 900, 1024, 33.19023, $ |
---|
196 | 163, 900, 1024, 32.83901, $ |
---|
197 | 164, 900, 1024, 32.48779, $ |
---|
198 | 165, 900, 1024, 32.13657, $ |
---|
199 | 166, 900, 1024, 31.78535, $ |
---|
200 | 167, 900, 1024, 31.43413, $ |
---|
201 | 168, 900, 1024, 31.08291, $ |
---|
202 | 169, 960, 1024, 30.73169, $ |
---|
203 | 170, 960, 1024, 30.38047, $ |
---|
204 | 171, 960, 1024, 30.02925, $ |
---|
205 | 172, 960, 1024, 29.67803, $ |
---|
206 | 173, 960, 1024, 29.32681, $ |
---|
207 | 174, 960, 1024, 28.97559, $ |
---|
208 | 175, 960, 1024, 28.62438, $ |
---|
209 | 176, 960, 1024, 28.27316, $ |
---|
210 | 177, 960, 1024, 27.92194, $ |
---|
211 | 178, 960, 1024, 27.57072, $ |
---|
212 | 179, 960, 1024, 27.21950, $ |
---|
213 | 180, 960, 1024, 26.86828, $ |
---|
214 | 181, 960, 1024, 26.51706, $ |
---|
215 | 182, 960, 1024, 26.16584, $ |
---|
216 | 183, 960, 1024, 25.81462, $ |
---|
217 | 184, 960, 1024, 25.46340, $ |
---|
218 | 185, 960, 1024, 25.11218, $ |
---|
219 | 186, 960, 1024, 24.76096, $ |
---|
220 | 187, 960, 1024, 24.40974, $ |
---|
221 | 188, 960, 1024, 24.05852, $ |
---|
222 | 189, 960, 1024, 23.70731, $ |
---|
223 | 190, 960, 1024, 23.35609, $ |
---|
224 | 191, 972, 1024, 23.00487, $ |
---|
225 | 192, 972, 1024, 22.65365, $ |
---|
226 | 193, 972, 1024, 22.30243, $ |
---|
227 | 194, 972, 1024, 21.95121, $ |
---|
228 | 195, 972, 1024, 21.59999, $ |
---|
229 | 196, 1000, 1024, 21.24877, $ |
---|
230 | 197, 1000, 1024, 20.89755, $ |
---|
231 | 198, 1000, 1024, 20.54633, $ |
---|
232 | 199, 1000, 1024, 20.19511, $ |
---|
233 | 200, 1000, 1024, 19.84389, $ |
---|
234 | 201, 1000, 1024, 19.49267, $ |
---|
235 | 202, 1000, 1024, 19.14145, $ |
---|
236 | 203, 1000, 1024, 18.79023, $ |
---|
237 | 204, 1000, 1024, 18.43902, $ |
---|
238 | 205, 1000, 1024, 18.08780, $ |
---|
239 | 206, 1000, 1024, 17.73658, $ |
---|
240 | 207, 1000, 1024, 17.38536, $ |
---|
241 | 208, 1000, 1024, 17.03414, $ |
---|
242 | 209, 1000, 1024, 16.68292, $ |
---|
243 | 210, 1000, 1024, 16.33170, $ |
---|
244 | 211, 1000, 1024, 15.98048, $ |
---|
245 | 212, 1024, 1024, 15.62926, $ |
---|
246 | 213, 1024, 1024, 15.27804, $ |
---|
247 | 214, 1024, 1024, 14.92682, $ |
---|
248 | 215, 1024, 1024, 14.57560, $ |
---|
249 | 216, 1024, 1024, 14.22438, $ |
---|
250 | 217, 1024, 1024, 13.87316, $ |
---|
251 | 218, 1024, 1024, 13.52194, $ |
---|
252 | 219, 1024, 1024, 13.17073, $ |
---|
253 | 220, 1024, 1024, 12.81951, $ |
---|
254 | 221, 1024, 1024, 12.46829, $ |
---|
255 | 222, 1024, 1024, 12.11707, $ |
---|
256 | 223, 1024, 1024, 11.76585, $ |
---|
257 | 224, 1024, 1024, 11.41463, $ |
---|
258 | 225, 1024, 1024, 11.06341, $ |
---|
259 | 226, 1024, 1024, 10.71219, $ |
---|
260 | 227, 1024, 1024, 10.36097, $ |
---|
261 | 228, 1024, 1024, 10.00975, $ |
---|
262 | 229, 1024, 1024, 9.65853, $ |
---|
263 | 230, 1024, 1024, 9.30731, $ |
---|
264 | 231, 1024, 1024, 8.95609, $ |
---|
265 | 232, 1024, 1024, 8.60487, $ |
---|
266 | 233, 1024, 1024, 8.25365, $ |
---|
267 | 234, 1024, 1024, 7.90244, $ |
---|
268 | 235, 1024, 1024, 7.55122, $ |
---|
269 | 236, 1024, 1024, 7.20000, $ |
---|
270 | 237, 1024, 1024, 6.84878, $ |
---|
271 | 238, 1024, 1024, 6.49756, $ |
---|
272 | 239, 1024, 1024, 6.14634, $ |
---|
273 | 240, 1024, 1024, 5.79512] |
---|
274 | n256 = [n256, $ |
---|
275 | 241, 1024, 1024, 5.44390, $ |
---|
276 | 242, 1024, 1024, 5.09268, $ |
---|
277 | 243, 1024, 1024, 4.74146, $ |
---|
278 | 244, 1024, 1024, 4.39024, $ |
---|
279 | 245, 1024, 1024, 4.03902, $ |
---|
280 | 246, 1024, 1024, 3.68780, $ |
---|
281 | 247, 1024, 1024, 3.33658, $ |
---|
282 | 248, 1024, 1024, 2.98536, $ |
---|
283 | 249, 1024, 1024, 2.63415, $ |
---|
284 | 250, 1024, 1024, 2.28293, $ |
---|
285 | 251, 1024, 1024, 1.93171, $ |
---|
286 | 252, 1024, 1024, 1.58049, $ |
---|
287 | 253, 1024, 1024, 1.22927, $ |
---|
288 | 254, 1024, 1024, 0.87805, $ |
---|
289 | 255, 1024, 1024, 0.52683, $ |
---|
290 | 256, 1024, 1024, 0.17561] |
---|
291 | ; |
---|
292 | n256 = reform(n256, 4, 256, /over) |
---|
293 | n256 = reform(n256[3, *], /over) |
---|
294 | n256 = [n256, -reverse(n256)] |
---|
295 | |
---|
296 | return, n256 |
---|
297 | end |
---|