Context Navigation

n128gaussian.pro @ 325

Last change on this file since 325 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: 7.4 KB

Line
1	;+
2	;
3	; @file_comments
4	; compute the latitudes of the n128 gaussian grid. See:
5	; <a href="http://www.ecmwf.int/products/data/technical/gaussian/n128FIS.html">Gaussian 128</a>
6	;
7	; @categories
8	; Grid
9	;
10	; @examples
11	; IDL> lat=n128gaussian()
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 n128gaussian
25	;
26	compile_opt idl2, strictarrsubs
27	;
28	; latitude reduced regular latitude
29	; number points points
30	; ------- ------- ------- --------
31	n128 = [ 1, 18, 512, 89.46282, $
32	2, 25, 512, 88.76695, $
33	3, 36, 512, 88.06697, $
34	4, 40, 512, 87.36606, $
35	5, 45, 512, 86.66480, $
36	6, 50, 512, 85.96337, $
37	7, 60, 512, 85.26184, $
38	8, 64, 512, 84.56026, $
39	9, 72, 512, 83.85863, $
40	10, 72, 512, 83.15698, $
41	11, 80, 512, 82.45531, $
42	12, 90, 512, 81.75363, $
43	13, 90, 512, 81.05194, $
44	14, 100, 512, 80.35023, $
45	15, 108, 512, 79.64852, $
46	16, 120, 512, 78.94681, $
47	17, 120, 512, 78.24509, $
48	18, 125, 512, 77.54336, $
49	19, 128, 512, 76.84163, $
50	20, 144, 512, 76.13990, $
51	21, 144, 512, 75.43817, $
52	22, 150, 512, 74.73644, $
53	23, 160, 512, 74.03470, $
54	24, 160, 512, 73.33296, $
55	25, 180, 512, 72.63123, $
56	26, 180, 512, 71.92949, $
57	27, 180, 512, 71.22774, $
58	28, 192, 512, 70.52600, $
59	29, 192, 512, 69.82426, $
60	30, 200, 512, 69.12252, $
61	31, 216, 512, 68.42077, $
62	32, 216, 512, 67.71903, $
63	33, 216, 512, 67.01728, $
64	34, 225, 512, 66.31554, $
65	35, 240, 512, 65.61379, $
66	36, 240, 512, 64.91204, $
67	37, 240, 512, 64.21030, $
68	38, 250, 512, 63.50855, $
69	39, 250, 512, 62.80680, $
70	40, 256, 512, 62.10505, $
71	41, 270, 512, 61.40330, $
72	42, 270, 512, 60.70156, $
73	43, 288, 512, 59.99981, $
74	44, 288, 512, 59.29806, $
75	45, 288, 512, 58.59631, $
76	46, 300, 512, 57.89456, $
77	47, 300, 512, 57.19281, $
78	48, 320, 512, 56.49106, $
79	49, 320, 512, 55.78931, $
80	50, 320, 512, 55.08756, $
81	51, 320, 512, 54.38581, $
82	52, 324, 512, 53.68406, $
83	53, 360, 512, 52.98231, $
84	54, 360, 512, 52.28056, $
85	55, 360, 512, 51.57881, $
86	56, 360, 512, 50.87705, $
87	57, 360, 512, 50.17530, $
88	58, 360, 512, 49.47355, $
89	59, 360, 512, 48.77180, $
90	60, 375, 512, 48.07005, $
91	61, 375, 512, 47.36830, $
92	62, 375, 512, 46.66655, $
93	63, 375, 512, 45.96479, $
94	64, 384, 512, 45.26304, $
95	65, 384, 512, 44.56129, $
96	66, 400, 512, 43.85954, $
97	67, 400, 512, 43.15779, $
98	68, 400, 512, 42.45604, $
99	69, 400, 512, 41.75428, $
100	70, 405, 512, 41.05253, $
101	71, 432, 512, 40.35078, $
102	72, 432, 512, 39.64903, $
103	73, 432, 512, 38.94728, $
104	74, 432, 512, 38.24552, $
105	75, 432, 512, 37.54377, $
106	76, 432, 512, 36.84202, $
107	77, 432, 512, 36.14027, $
108	78, 450, 512, 35.43851, $
109	79, 450, 512, 34.73676, $
110	80, 450, 512, 34.03501]
111	n128 = [n128, $
112	81, 450, 512, 33.33326, $
113	82, 450, 512, 32.63150, $
114	83, 480, 512, 31.92975, $
115	84, 480, 512, 31.22800, $
116	85, 480, 512, 30.52625, $
117	86, 480, 512, 29.82449, $
118	87, 480, 512, 29.12274, $
119	88, 480, 512, 28.42099, $
120	89, 480, 512, 27.71924, $
121	90, 480, 512, 27.01748, $
122	91, 480, 512, 26.31573, $
123	92, 480, 512, 25.61398, $
124	93, 486, 512, 24.91223, $
125	94, 486, 512, 24.21047, $
126	95, 486, 512, 23.50872, $
127	96, 500, 512, 22.80697, $
128	97, 500, 512, 22.10521, $
129	98, 500, 512, 21.40346, $
130	99, 500, 512, 20.70171, $
131	100, 500, 512, 19.99996, $
132	101, 500, 512, 19.29820, $
133	102, 500, 512, 18.59645, $
134	103, 512, 512, 17.89470, $
135	104, 512, 512, 17.19294, $
136	105, 512, 512, 16.49119, $
137	106, 512, 512, 15.78944, $
138	107, 512, 512, 15.08768, $
139	108, 512, 512, 14.38593, $
140	109, 512, 512, 13.68418, $
141	110, 512, 512, 12.98243, $
142	111, 512, 512, 12.28067, $
143	112, 512, 512, 11.57892, $
144	113, 512, 512, 10.87717, $
145	114, 512, 512, 10.17541, $
146	115, 512, 512, 9.47366, $
147	116, 512, 512, 8.77191, $
148	117, 512, 512, 8.07016, $
149	118, 512, 512, 7.36840, $
150	119, 512, 512, 6.66665, $
151	120, 512, 512, 5.96490, $
152	121, 512, 512, 5.26314, $
153	122, 512, 512, 4.56139, $
154	123, 512, 512, 3.85964, $
155	124, 512, 512, 3.15788, $
156	125, 512, 512, 2.45613, $
157	126, 512, 512, 1.75438, $
158	127, 512, 512, 1.05262, $
159	128, 512, 512, 0.35087]
160	;
161	n128 = reform(n128, 4, 128, /over)
162	n128 = reform(n128[3, *], /over)
163	n128 = [n128, -reverse(n128)]
164
165	return, n128
166	end

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

Download in other formats:

Original Format