Context Navigation

n256gaussian.pro @ 163

Last change on this file since 163 was 157, checked in by navarro, 18 years ago
header improvements + xxx doc
Property svn:keywords set to `Id`
File size: 13.9 KB

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format