Context Navigation

n256gaussian.pro @ 371

Last change on this file since 371 was 371, checked in by pinsard, 16 years ago
improvements of headers (alignments of IDL prompt in examples)
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	;
12	; IDL> lat=n256gaussian()
13	;
14	; @returns
15	; a 1d array
16	;
17	; @history
18	; Sebastien Masson (smasson\@lodyc.jussieu.fr)
19	; June 2004
20	;
21	; @version
22	; $Id$
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

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

Download in other formats:

Original Format