Context Navigation

n160gaussian.pro @ 118

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

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

Download in other formats:

Original Format