Context Navigation

n128gaussian.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: 7.4 KB

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

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

Download in other formats:

Original Format