Context Navigation

← Previous Revision
Next Revision →
Blame
Revision Log

n160gaussian.pro

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

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

Download in other formats:

Original Format