1 | ;+ |
---|
2 | ; |
---|
3 | ; @file_comments compute the latitudes of the n80 gaussian grid. See: |
---|
4 | ; http://www.ecmwf.int/products/data/technical/gaussian/n80FIS.html |
---|
5 | ; |
---|
6 | ; @categories |
---|
7 | ; Grid |
---|
8 | ; |
---|
9 | ; @examples |
---|
10 | ; IDL> lat=n80gaussian() |
---|
11 | ; |
---|
12 | ; @returns a 1d array |
---|
13 | ; |
---|
14 | ; @history |
---|
15 | ; Sebastien Masson (smasson\@lodyc.jussieu.fr) |
---|
16 | ; June 2004 |
---|
17 | ; |
---|
18 | ; @version $Id$ |
---|
19 | ; |
---|
20 | ;- |
---|
21 | ; |
---|
22 | FUNCTION n80gaussian |
---|
23 | ; |
---|
24 | ; |
---|
25 | compile_opt idl2, strictarrsubs |
---|
26 | ; |
---|
27 | ; latitude reduced regular latitude |
---|
28 | ; number points points |
---|
29 | ; ------- ------- ------- -------- |
---|
30 | n80 = [ 1, 18, 320, 89.14152, $ |
---|
31 | 2, 25, 320, 88.02943, $ |
---|
32 | 3, 36, 320, 86.91077, $ |
---|
33 | 4, 40, 320, 85.79063, $ |
---|
34 | 5, 45, 320, 84.66992, $ |
---|
35 | 6, 54, 320, 83.54895, $ |
---|
36 | 7, 60, 320, 82.42782, $ |
---|
37 | 8, 64, 320, 81.30659, $ |
---|
38 | 9, 72, 320, 80.18531, $ |
---|
39 | 10, 72, 320, 79.06398, $ |
---|
40 | 11, 80, 320, 77.94262, $ |
---|
41 | 12, 90, 320, 76.82124, $ |
---|
42 | 13, 96, 320, 75.69984, $ |
---|
43 | 14, 100, 320, 74.57843, $ |
---|
44 | 15, 108, 320, 73.45701, $ |
---|
45 | 16, 120, 320, 72.33558, $ |
---|
46 | 17, 120, 320, 71.21414, $ |
---|
47 | 18, 128, 320, 70.09269, $ |
---|
48 | 19, 135, 320, 68.97124, $ |
---|
49 | 20, 144, 320, 67.84978, $ |
---|
50 | 21, 144, 320, 66.72833, $ |
---|
51 | 22, 150, 320, 65.60686, $ |
---|
52 | 23, 160, 320, 64.48540, $ |
---|
53 | 24, 160, 320, 63.36393, $ |
---|
54 | 25, 180, 320, 62.24246, $ |
---|
55 | 26, 180, 320, 61.12099, $ |
---|
56 | 27, 180, 320, 59.99952, $ |
---|
57 | 28, 192, 320, 58.87804, $ |
---|
58 | 29, 192, 320, 57.75657, $ |
---|
59 | 30, 200, 320, 56.63509, $ |
---|
60 | 31, 200, 320, 55.51361, $ |
---|
61 | 32, 216, 320, 54.39214, $ |
---|
62 | 33, 216, 320, 53.27066, $ |
---|
63 | 34, 216, 320, 52.14917, $ |
---|
64 | 35, 225, 320, 51.02769, $ |
---|
65 | 36, 225, 320, 49.90621, $ |
---|
66 | 37, 240, 320, 48.78473, $ |
---|
67 | 38, 240, 320, 47.66325, $ |
---|
68 | 39, 240, 320, 46.54176, $ |
---|
69 | 40, 256, 320, 45.42028, $ |
---|
70 | 41, 256, 320, 44.29879, $ |
---|
71 | 42, 256, 320, 43.17731, $ |
---|
72 | 43, 256, 320, 42.05582, $ |
---|
73 | 44, 288, 320, 40.93434, $ |
---|
74 | 45, 288, 320, 39.81285, $ |
---|
75 | 46, 288, 320, 38.69137, $ |
---|
76 | 47, 288, 320, 37.56988, $ |
---|
77 | 48, 288, 320, 36.44839, $ |
---|
78 | 49, 288, 320, 35.32691, $ |
---|
79 | 50, 288, 320, 34.20542, $ |
---|
80 | 51, 288, 320, 33.08393, $ |
---|
81 | 52, 288, 320, 31.96244, $ |
---|
82 | 53, 300, 320, 30.84096, $ |
---|
83 | 54, 300, 320, 29.71947, $ |
---|
84 | 55, 300, 320, 28.59798, $ |
---|
85 | 56, 300, 320, 27.47649, $ |
---|
86 | 57, 320, 320, 26.35500, $ |
---|
87 | 58, 320, 320, 25.23351, $ |
---|
88 | 59, 320, 320, 24.11203, $ |
---|
89 | 60, 320, 320, 22.99054, $ |
---|
90 | 61, 320, 320, 21.86905, $ |
---|
91 | 62, 320, 320, 20.74756, $ |
---|
92 | 63, 320, 320, 19.62607, $ |
---|
93 | 64, 320, 320, 18.50458, $ |
---|
94 | 65, 320, 320, 17.38309, $ |
---|
95 | 66, 320, 320, 16.26160, $ |
---|
96 | 67, 320, 320, 15.14011, $ |
---|
97 | 68, 320, 320, 14.01862, $ |
---|
98 | 69, 320, 320, 12.89713, $ |
---|
99 | 70, 320, 320, 11.77564, $ |
---|
100 | 71, 320, 320, 10.65415, $ |
---|
101 | 72, 320, 320, 9.53266, $ |
---|
102 | 73, 320, 320, 8.41117, $ |
---|
103 | 74, 320, 320, 7.28968, $ |
---|
104 | 75, 320, 320, 6.16819, $ |
---|
105 | 76, 320, 320, 5.04670, $ |
---|
106 | 77, 320, 320, 3.92521, $ |
---|
107 | 78, 320, 320, 2.80372, $ |
---|
108 | 79, 320, 320, 1.68223, $ |
---|
109 | 80, 320, 320, 0.56074] |
---|
110 | ; |
---|
111 | n80 = reform(n80, 4, 80, /over) |
---|
112 | n80 = reform(n80[3, *], /over) |
---|
113 | n80 = [n80, -reverse(n80)] |
---|
114 | |
---|
115 | return, n80 |
---|
116 | end |
---|