Context Navigation

n160gaussian.pro @ 103

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

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

Download in other formats:

Original Format