Context Navigation

n256gaussian.pro @ 114

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

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

Download in other formats:

Original Format