Context Navigation

n160gaussian.pro @ 157

Last change on this file since 157 was 157, checked in by navarro, 18 years ago
header improvements + xxx doc
Property svn:keywords set to `Id`
File size: 9.0 KB

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

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

Download in other formats:

Original Format