Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

n256gaussian.pro @ 98

Last change on this file since 98 was 13, checked in by pinsard, 18 years ago
upgrade of GRILLE/Utilities according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files
Property svn:executable set to ``*
File size: 14.0 KB

Line
1	;+
2	; NAME:n256gaussian
3	;
4	; PURPOSE:compute the latitudes of the n256 gaussian grid. See:
5	; http://www.ecmwf.int/products/data/technical/gaussian/n256FIS.html
6	;
7	; CATEGORY:grid
8	;
9	; CALLING SEQUENCE:lat=n256gaussian()
10	;
11	; INPUTS:None
12	;
13	; KEYWORD PARAMETERS:None
14	;
15	; OUTPUTS:a 1d array
16	;
17	; COMMON BLOCKS:None
18	;
19	; SIDE EFFECTS:None
20	;
21	; RESTRICTIONS:None
22	;
23	; EXAMPLE:
24	;
25	; MODIFICATION HISTORY:
26	; Sebastien Masson (smasson@lodyc.jussieu.fr)
27	; June 2004
28	;-
29	;
30	FUNCTION n256gaussian
31	;
32	; latitude reduced regular latitude
33	; number points points
34	; ------- ------- ------- --------
35	n256 = [ 1, 18, 1024, 89.73115, $
36	2, 25, 1024, 89.38287, $
37	3, 32, 1024, 89.03254, $
38	4, 40, 1024, 88.68175, $
39	5, 45, 1024, 88.33077, $
40	6, 50, 1024, 87.97972, $
41	7, 60, 1024, 87.62861, $
42	8, 64, 1024, 87.27748, $
43	9, 72, 1024, 86.92632, $
44	10, 72, 1024, 86.57515, $
45	11, 75, 1024, 86.22398, $
46	12, 81, 1024, 85.87279, $
47	13, 90, 1024, 85.52160, $
48	14, 96, 1024, 85.17041, $
49	15, 100, 1024, 84.81921, $
50	16, 108, 1024, 84.46801, $
51	17, 120, 1024, 84.11681, $
52	18, 120, 1024, 83.76560, $
53	19, 125, 1024, 83.41440, $
54	20, 135, 1024, 83.06319, $
55	21, 144, 1024, 82.71198, $
56	22, 150, 1024, 82.36077, $
57	23, 160, 1024, 82.00956, $
58	24, 160, 1024, 81.65835, $
59	25, 180, 1024, 81.30714, $
60	26, 180, 1024, 80.95593, $
61	27, 180, 1024, 80.60471, $
62	28, 192, 1024, 80.25350, $
63	29, 192, 1024, 79.90229, $
64	30, 200, 1024, 79.55107, $
65	31, 216, 1024, 79.19986, $
66	32, 216, 1024, 78.84864, $
67	33, 216, 1024, 78.49743, $
68	34, 225, 1024, 78.14621, $
69	35, 240, 1024, 77.79500, $
70	36, 240, 1024, 77.44378, $
71	37, 243, 1024, 77.09256, $
72	38, 250, 1024, 76.74135, $
73	39, 256, 1024, 76.39013, $
74	40, 270, 1024, 76.03891, $
75	41, 270, 1024, 75.68770, $
76	42, 288, 1024, 75.33648, $
77	43, 288, 1024, 74.98526, $
78	44, 288, 1024, 74.63405, $
79	45, 300, 1024, 74.28283, $
80	46, 300, 1024, 73.93161, $
81	47, 320, 1024, 73.58040, $
82	48, 320, 1024, 73.22918, $
83	49, 320, 1024, 72.87796, $
84	50, 324, 1024, 72.52674, $
85	51, 360, 1024, 72.17552, $
86	52, 360, 1024, 71.82431, $
87	53, 360, 1024, 71.47309, $
88	54, 360, 1024, 71.12187, $
89	55, 360, 1024, 70.77065, $
90	56, 360, 1024, 70.41944, $
91	57, 375, 1024, 70.06822, $
92	58, 375, 1024, 69.71700, $
93	59, 384, 1024, 69.36578, $
94	60, 384, 1024, 69.01456, $
95	61, 400, 1024, 68.66334, $
96	62, 400, 1024, 68.31213, $
97	63, 400, 1024, 67.96091, $
98	64, 432, 1024, 67.60969, $
99	65, 432, 1024, 67.25847, $
100	66, 432, 1024, 66.90725, $
101	67, 432, 1024, 66.55603, $
102	68, 432, 1024, 66.20482, $
103	69, 450, 1024, 65.85360, $
104	70, 450, 1024, 65.50238, $
105	71, 450, 1024, 65.15116, $
106	72, 480, 1024, 64.79994, $
107	73, 480, 1024, 64.44872, $
108	74, 480, 1024, 64.09750, $
109	75, 480, 1024, 63.74629, $
110	76, 480, 1024, 63.39507, $
111	77, 486, 1024, 63.04385, $
112	78, 500, 1024, 62.69263, $
113	79, 500, 1024, 62.34141, $
114	80, 500, 1024, 61.99019]
115	n256 = [n256, $
116	81, 512, 1024, 61.63897, $
117	82, 512, 1024, 61.28776, $
118	83, 540, 1024, 60.93654, $
119	84, 540, 1024, 60.58532, $
120	85, 540, 1024, 60.23410, $
121	86, 540, 1024, 59.88288, $
122	87, 540, 1024, 59.53166, $
123	88, 576, 1024, 59.18044, $
124	89, 576, 1024, 58.82922, $
125	90, 576, 1024, 58.47800, $
126	91, 576, 1024, 58.12679, $
127	92, 576, 1024, 57.77557, $
128	93, 576, 1024, 57.42435, $
129	94, 600, 1024, 57.07313, $
130	95, 600, 1024, 56.72191, $
131	96, 600, 1024, 56.37069, $
132	97, 600, 1024, 56.01947, $
133	98, 600, 1024, 55.66825, $
134	99, 640, 1024, 55.31703, $
135	100, 640, 1024, 54.96581, $
136	101, 640, 1024, 54.61460, $
137	102, 640, 1024, 54.26338, $
138	103, 640, 1024, 53.91216, $
139	104, 640, 1024, 53.56094, $
140	105, 640, 1024, 53.20972, $
141	106, 640, 1024, 52.85850, $
142	107, 648, 1024, 52.50728, $
143	108, 675, 1024, 52.15606, $
144	109, 675, 1024, 51.80484, $
145	110, 675, 1024, 51.45362, $
146	111, 675, 1024, 51.10241, $
147	112, 675, 1024, 50.75119, $
148	113, 675, 1024, 50.39997, $
149	114, 720, 1024, 50.04875, $
150	115, 720, 1024, 49.69753, $
151	116, 720, 1024, 49.34631, $
152	117, 720, 1024, 48.99509, $
153	118, 720, 1024, 48.64387, $
154	119, 720, 1024, 48.29265, $
155	120, 720, 1024, 47.94143, $
156	121, 720, 1024, 47.59021, $
157	122, 720, 1024, 47.23899, $
158	123, 729, 1024, 46.88778, $
159	124, 729, 1024, 46.53656, $
160	125, 750, 1024, 46.18534, $
161	126, 750, 1024, 45.83412, $
162	127, 750, 1024, 45.48290, $
163	128, 750, 1024, 45.13168, $
164	129, 750, 1024, 44.78046, $
165	130, 768, 1024, 44.42924, $
166	131, 768, 1024, 44.07802, $
167	132, 768, 1024, 43.72680, $
168	133, 768, 1024, 43.37558, $
169	134, 800, 1024, 43.02436, $
170	135, 800, 1024, 42.67315, $
171	136, 800, 1024, 42.32193, $
172	137, 800, 1024, 41.97071, $
173	138, 800, 1024, 41.61949, $
174	139, 800, 1024, 41.26827, $
175	140, 800, 1024, 40.91705, $
176	141, 800, 1024, 40.56583, $
177	142, 810, 1024, 40.21461, $
178	143, 810, 1024, 39.86339, $
179	144, 864, 1024, 39.51217, $
180	145, 864, 1024, 39.16095, $
181	146, 864, 1024, 38.80973, $
182	147, 864, 1024, 38.45851, $
183	148, 864, 1024, 38.10730, $
184	149, 864, 1024, 37.75608, $
185	150, 864, 1024, 37.40486, $
186	151, 864, 1024, 37.05364, $
187	152, 864, 1024, 36.70242, $
188	153, 864, 1024, 36.35120, $
189	154, 864, 1024, 35.99998, $
190	155, 864, 1024, 35.64876, $
191	156, 864, 1024, 35.29754, $
192	157, 864, 1024, 34.94632, $
193	158, 900, 1024, 34.59510, $
194	159, 900, 1024, 34.24388, $
195	160, 900, 1024, 33.89266]
196	n256 = [n256, $
197	161, 900, 1024, 33.54145, $
198	162, 900, 1024, 33.19023, $
199	163, 900, 1024, 32.83901, $
200	164, 900, 1024, 32.48779, $
201	165, 900, 1024, 32.13657, $
202	166, 900, 1024, 31.78535, $
203	167, 900, 1024, 31.43413, $
204	168, 900, 1024, 31.08291, $
205	169, 960, 1024, 30.73169, $
206	170, 960, 1024, 30.38047, $
207	171, 960, 1024, 30.02925, $
208	172, 960, 1024, 29.67803, $
209	173, 960, 1024, 29.32681, $
210	174, 960, 1024, 28.97559, $
211	175, 960, 1024, 28.62438, $
212	176, 960, 1024, 28.27316, $
213	177, 960, 1024, 27.92194, $
214	178, 960, 1024, 27.57072, $
215	179, 960, 1024, 27.21950, $
216	180, 960, 1024, 26.86828, $
217	181, 960, 1024, 26.51706, $
218	182, 960, 1024, 26.16584, $
219	183, 960, 1024, 25.81462, $
220	184, 960, 1024, 25.46340, $
221	185, 960, 1024, 25.11218, $
222	186, 960, 1024, 24.76096, $
223	187, 960, 1024, 24.40974, $
224	188, 960, 1024, 24.05852, $
225	189, 960, 1024, 23.70731, $
226	190, 960, 1024, 23.35609, $
227	191, 972, 1024, 23.00487, $
228	192, 972, 1024, 22.65365, $
229	193, 972, 1024, 22.30243, $
230	194, 972, 1024, 21.95121, $
231	195, 972, 1024, 21.59999, $
232	196, 1000, 1024, 21.24877, $
233	197, 1000, 1024, 20.89755, $
234	198, 1000, 1024, 20.54633, $
235	199, 1000, 1024, 20.19511, $
236	200, 1000, 1024, 19.84389, $
237	201, 1000, 1024, 19.49267, $
238	202, 1000, 1024, 19.14145, $
239	203, 1000, 1024, 18.79023, $
240	204, 1000, 1024, 18.43902, $
241	205, 1000, 1024, 18.08780, $
242	206, 1000, 1024, 17.73658, $
243	207, 1000, 1024, 17.38536, $
244	208, 1000, 1024, 17.03414, $
245	209, 1000, 1024, 16.68292, $
246	210, 1000, 1024, 16.33170, $
247	211, 1000, 1024, 15.98048, $
248	212, 1024, 1024, 15.62926, $
249	213, 1024, 1024, 15.27804, $
250	214, 1024, 1024, 14.92682, $
251	215, 1024, 1024, 14.57560, $
252	216, 1024, 1024, 14.22438, $
253	217, 1024, 1024, 13.87316, $
254	218, 1024, 1024, 13.52194, $
255	219, 1024, 1024, 13.17073, $
256	220, 1024, 1024, 12.81951, $
257	221, 1024, 1024, 12.46829, $
258	222, 1024, 1024, 12.11707, $
259	223, 1024, 1024, 11.76585, $
260	224, 1024, 1024, 11.41463, $
261	225, 1024, 1024, 11.06341, $
262	226, 1024, 1024, 10.71219, $
263	227, 1024, 1024, 10.36097, $
264	228, 1024, 1024, 10.00975, $
265	229, 1024, 1024, 9.65853, $
266	230, 1024, 1024, 9.30731, $
267	231, 1024, 1024, 8.95609, $
268	232, 1024, 1024, 8.60487, $
269	233, 1024, 1024, 8.25365, $
270	234, 1024, 1024, 7.90244, $
271	235, 1024, 1024, 7.55122, $
272	236, 1024, 1024, 7.20000, $
273	237, 1024, 1024, 6.84878, $
274	238, 1024, 1024, 6.49756, $
275	239, 1024, 1024, 6.14634, $
276	240, 1024, 1024, 5.79512]
277	n256 = [n256, $
278	241, 1024, 1024, 5.44390, $
279	242, 1024, 1024, 5.09268, $
280	243, 1024, 1024, 4.74146, $
281	244, 1024, 1024, 4.39024, $
282	245, 1024, 1024, 4.03902, $
283	246, 1024, 1024, 3.68780, $
284	247, 1024, 1024, 3.33658, $
285	248, 1024, 1024, 2.98536, $
286	249, 1024, 1024, 2.63415, $
287	250, 1024, 1024, 2.28293, $
288	251, 1024, 1024, 1.93171, $
289	252, 1024, 1024, 1.58049, $
290	253, 1024, 1024, 1.22927, $
291	254, 1024, 1024, 0.87805, $
292	255, 1024, 1024, 0.52683, $
293	256, 1024, 1024, 0.17561]
294	;
295	n256 = reform(n256, 4, 256, /over)
296	n256 = reform(n256[3, *], /over)
297	n256 = [n256, -reverse(n256)]
298
299	return, n256
300	end

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

Download in other formats:

Original Format