Context Navigation

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

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

Download in other formats:

Original Format