Context Navigation

n160gaussian.pro @ 352

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

Line
1	;+
2	;
3	; @file_comments
4	; compute the latitudes of the n160 gaussian grid. See:
5	; <a href="http://www.ecmwf.int/products/data/technical/gaussian/n160FIS.html">Gaussian 160</a>
6	;
7	; @categories
8	; Grid
9	;
10	; @examples
11	; IDL> lat=n160gaussian()
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 n160gaussian
25	;
26	compile_opt idl2, strictarrsubs
27	;
28	; latitude reduced regular latitude
29	; number points points
30	; ------- ------- ------- --------
31	n160 = [ 1, 18, 640, 89.57009, $
32	2, 25, 640, 89.01318, $
33	3, 36, 640, 88.45297, $
34	4, 40, 640, 87.89203, $
35	5, 45, 640, 87.33080, $
36	6, 50, 640, 86.76944, $
37	7, 60, 640, 86.20800, $
38	8, 64, 640, 85.64651, $
39	9, 72, 640, 85.08499, $
40	10, 72, 640, 84.52345, $
41	11, 80, 640, 83.96190, $
42	12, 90, 640, 83.40033, $
43	13, 90, 640, 82.83876, $
44	14, 96, 640, 82.27718, $
45	15, 108, 640, 81.71559, $
46	16, 120, 640, 81.15400, $
47	17, 120, 640, 80.59240, $
48	18, 125, 640, 80.03080, $
49	19, 128, 640, 79.46920, $
50	20, 135, 640, 78.90760, $
51	21, 144, 640, 78.34600, $
52	22, 150, 640, 77.78439, $
53	23, 160, 640, 77.22278, $
54	24, 160, 640, 76.66117, $
55	25, 180, 640, 76.09956, $
56	26, 180, 640, 75.53795, $
57	27, 180, 640, 74.97634, $
58	28, 192, 640, 74.41473, $
59	29, 192, 640, 73.85311, $
60	30, 200, 640, 73.29150, $
61	31, 216, 640, 72.72988, $
62	32, 216, 640, 72.16827, $
63	33, 225, 640, 71.60665, $
64	34, 225, 640, 71.04504, $
65	35, 240, 640, 70.48342, $
66	36, 240, 640, 69.92181, $
67	37, 243, 640, 69.36019, $
68	38, 250, 640, 68.79857, $
69	39, 256, 640, 68.23695, $
70	40, 270, 640, 67.67534, $
71	41, 270, 640, 67.11372, $
72	42, 288, 640, 66.55210, $
73	43, 288, 640, 65.99048, $
74	44, 288, 640, 65.42886, $
75	45, 300, 640, 64.86725, $
76	46, 300, 640, 64.30563, $
77	47, 320, 640, 63.74401, $
78	48, 320, 640, 63.18239, $
79	49, 320, 640, 62.62077, $
80	50, 320, 640, 62.05915, $
81	51, 324, 640, 61.49753, $
82	52, 360, 640, 60.93591, $
83	53, 360, 640, 60.37429, $
84	54, 360, 640, 59.81267, $
85	55, 360, 640, 59.25105, $
86	56, 360, 640, 58.68943, $
87	57, 360, 640, 58.12781, $
88	58, 375, 640, 57.56619, $
89	59, 375, 640, 57.00457, $
90	60, 375, 640, 56.44295, $
91	61, 384, 640, 55.88133, $
92	62, 384, 640, 55.31971, $
93	63, 400, 640, 54.75809, $
94	64, 400, 640, 54.19647, $
95	65, 400, 640, 53.63485, $
96	66, 405, 640, 53.07323, $
97	67, 432, 640, 52.51161, $
98	68, 432, 640, 51.94999, $
99	69, 432, 640, 51.38837, $
100	70, 432, 640, 50.82675, $
101	71, 432, 640, 50.26513, $
102	72, 450, 640, 49.70351, $
103	73, 450, 640, 49.14189, $
104	74, 450, 640, 48.58026, $
105	75, 450, 640, 48.01864, $
106	76, 480, 640, 47.45702, $
107	77, 480, 640, 46.89540, $
108	78, 480, 640, 46.33378, $
109	79, 480, 640, 45.77216, $
110	80, 480, 640, 45.21054]
111	n160 = [n160, $
112	81, 480, 640, 44.64892, $
113	82, 480, 640, 44.08730, $
114	83, 500, 640, 43.52567, $
115	84, 500, 640, 42.96405, $
116	85, 500, 640, 42.40243, $
117	86, 500, 640, 41.84081, $
118	87, 500, 640, 41.27919, $
119	88, 512, 640, 40.71757, $
120	89, 512, 640, 40.15595, $
121	90, 540, 640, 39.59433, $
122	91, 540, 640, 39.03270, $
123	92, 540, 640, 38.47108, $
124	93, 540, 640, 37.90946, $
125	94, 540, 640, 37.34784, $
126	95, 540, 640, 36.78622, $
127	96, 540, 640, 36.22460, $
128	97, 540, 640, 35.66298, $
129	98, 576, 640, 35.10136, $
130	99, 576, 640, 34.53973, $
131	100, 576, 640, 33.97811, $
132	101, 576, 640, 33.41649, $
133	102, 576, 640, 32.85487, $
134	103, 576, 640, 32.29325, $
135	104, 576, 640, 31.73163, $
136	105, 576, 640, 31.17000, $
137	106, 576, 640, 30.60838, $
138	107, 576, 640, 30.04676, $
139	108, 600, 640, 29.48514, $
140	109, 600, 640, 28.92352, $
141	110, 600, 640, 28.36190, $
142	111, 600, 640, 27.80028, $
143	112, 600, 640, 27.23865, $
144	113, 600, 640, 26.67703, $
145	114, 600, 640, 26.11541, $
146	115, 600, 640, 25.55379, $
147	116, 600, 640, 24.99217, $
148	117, 640, 640, 24.43055, $
149	118, 640, 640, 23.86892, $
150	119, 640, 640, 23.30730, $
151	120, 640, 640, 22.74568, $
152	121, 640, 640, 22.18406, $
153	122, 640, 640, 21.62244, $
154	123, 640, 640, 21.06082, $
155	124, 640, 640, 20.49919, $
156	125, 640, 640, 19.93757, $
157	126, 640, 640, 19.37595, $
158	127, 640, 640, 18.81433, $
159	128, 640, 640, 18.25271, $
160	129, 640, 640, 17.69109, $
161	130, 640, 640, 17.12946, $
162	131, 640, 640, 16.56784, $
163	132, 640, 640, 16.00622, $
164	133, 640, 640, 15.44460, $
165	134, 640, 640, 14.88298, $
166	135, 640, 640, 14.32136, $
167	136, 640, 640, 13.75973, $
168	137, 640, 640, 13.19811, $
169	138, 640, 640, 12.63649, $
170	139, 640, 640, 12.07487, $
171	140, 640, 640, 11.51325, $
172	141, 640, 640, 10.95162, $
173	142, 640, 640, 10.39000, $
174	143, 640, 640, 9.82838, $
175	144, 640, 640, 9.26676, $
176	145, 640, 640, 8.70514, $
177	146, 640, 640, 8.14352, $
178	147, 640, 640, 7.58189, $
179	148, 640, 640, 7.02027, $
180	149, 640, 640, 6.45865, $
181	150, 640, 640, 5.89703, $
182	151, 640, 640, 5.33541, $
183	152, 640, 640, 4.77379, $
184	153, 640, 640, 4.21216, $
185	154, 640, 640, 3.65054, $
186	155, 640, 640, 3.08892, $
187	156, 640, 640, 2.52730, $
188	157, 640, 640, 1.96568, $
189	158, 640, 640, 1.40405, $
190	159, 640, 640, 0.84243, $
191	160, 640, 640, 0.28081]
192	;
193	n160 = reform(n160, 4, 160, /over)
194	n160 = reform(n160[3, *], /over)
195	n160 = [n160, -reverse(n160)]
196
197	return, n160
198	end

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

Download in other formats:

Original Format