Context Navigation

n256gaussian.pro @ 368

Last change on this file since 368 was 325, checked in by pinsard, 17 years ago
modification of some headers (+some corrections) to prepare usage of the new idldoc
Property svn:keywords set to `Id`
File size: 13.9 KB

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

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

Download in other formats:

Original Format