1 | ;+ |
---|
2 | ; |
---|
3 | ; @file_comments |
---|
4 | ; compute the latitudes of the n128 gaussian grid. See: |
---|
5 | ; <a href="http://www.ecmwf.int/products/data/technical/gaussian/n128FIS.html">Gaussian 128</a> |
---|
6 | ; |
---|
7 | ; @categories |
---|
8 | ; Grid |
---|
9 | ; |
---|
10 | ; @examples |
---|
11 | ; IDL> lat=n128gaussian() |
---|
12 | ; |
---|
13 | ; @returns |
---|
14 | ; a 1d array |
---|
15 | ; |
---|
16 | ; @history |
---|
17 | ; Sebastien Masson (smasson\@lodyc.jussieu.fr) |
---|
18 | ; June 2004 |
---|
19 | ; |
---|
20 | ; @version |
---|
21 | ; $Id$ |
---|
22 | ; |
---|
23 | ;- |
---|
24 | ; |
---|
25 | FUNCTION n128gaussian |
---|
26 | ; |
---|
27 | compile_opt idl2, strictarrsubs |
---|
28 | ; |
---|
29 | ; latitude reduced regular latitude |
---|
30 | ; number points points |
---|
31 | ; ------- ------- ------- -------- |
---|
32 | n128 = [ 1, 18, 512, 89.46282, $ |
---|
33 | 2, 25, 512, 88.76695, $ |
---|
34 | 3, 36, 512, 88.06697, $ |
---|
35 | 4, 40, 512, 87.36606, $ |
---|
36 | 5, 45, 512, 86.66480, $ |
---|
37 | 6, 50, 512, 85.96337, $ |
---|
38 | 7, 60, 512, 85.26184, $ |
---|
39 | 8, 64, 512, 84.56026, $ |
---|
40 | 9, 72, 512, 83.85863, $ |
---|
41 | 10, 72, 512, 83.15698, $ |
---|
42 | 11, 80, 512, 82.45531, $ |
---|
43 | 12, 90, 512, 81.75363, $ |
---|
44 | 13, 90, 512, 81.05194, $ |
---|
45 | 14, 100, 512, 80.35023, $ |
---|
46 | 15, 108, 512, 79.64852, $ |
---|
47 | 16, 120, 512, 78.94681, $ |
---|
48 | 17, 120, 512, 78.24509, $ |
---|
49 | 18, 125, 512, 77.54336, $ |
---|
50 | 19, 128, 512, 76.84163, $ |
---|
51 | 20, 144, 512, 76.13990, $ |
---|
52 | 21, 144, 512, 75.43817, $ |
---|
53 | 22, 150, 512, 74.73644, $ |
---|
54 | 23, 160, 512, 74.03470, $ |
---|
55 | 24, 160, 512, 73.33296, $ |
---|
56 | 25, 180, 512, 72.63123, $ |
---|
57 | 26, 180, 512, 71.92949, $ |
---|
58 | 27, 180, 512, 71.22774, $ |
---|
59 | 28, 192, 512, 70.52600, $ |
---|
60 | 29, 192, 512, 69.82426, $ |
---|
61 | 30, 200, 512, 69.12252, $ |
---|
62 | 31, 216, 512, 68.42077, $ |
---|
63 | 32, 216, 512, 67.71903, $ |
---|
64 | 33, 216, 512, 67.01728, $ |
---|
65 | 34, 225, 512, 66.31554, $ |
---|
66 | 35, 240, 512, 65.61379, $ |
---|
67 | 36, 240, 512, 64.91204, $ |
---|
68 | 37, 240, 512, 64.21030, $ |
---|
69 | 38, 250, 512, 63.50855, $ |
---|
70 | 39, 250, 512, 62.80680, $ |
---|
71 | 40, 256, 512, 62.10505, $ |
---|
72 | 41, 270, 512, 61.40330, $ |
---|
73 | 42, 270, 512, 60.70156, $ |
---|
74 | 43, 288, 512, 59.99981, $ |
---|
75 | 44, 288, 512, 59.29806, $ |
---|
76 | 45, 288, 512, 58.59631, $ |
---|
77 | 46, 300, 512, 57.89456, $ |
---|
78 | 47, 300, 512, 57.19281, $ |
---|
79 | 48, 320, 512, 56.49106, $ |
---|
80 | 49, 320, 512, 55.78931, $ |
---|
81 | 50, 320, 512, 55.08756, $ |
---|
82 | 51, 320, 512, 54.38581, $ |
---|
83 | 52, 324, 512, 53.68406, $ |
---|
84 | 53, 360, 512, 52.98231, $ |
---|
85 | 54, 360, 512, 52.28056, $ |
---|
86 | 55, 360, 512, 51.57881, $ |
---|
87 | 56, 360, 512, 50.87705, $ |
---|
88 | 57, 360, 512, 50.17530, $ |
---|
89 | 58, 360, 512, 49.47355, $ |
---|
90 | 59, 360, 512, 48.77180, $ |
---|
91 | 60, 375, 512, 48.07005, $ |
---|
92 | 61, 375, 512, 47.36830, $ |
---|
93 | 62, 375, 512, 46.66655, $ |
---|
94 | 63, 375, 512, 45.96479, $ |
---|
95 | 64, 384, 512, 45.26304, $ |
---|
96 | 65, 384, 512, 44.56129, $ |
---|
97 | 66, 400, 512, 43.85954, $ |
---|
98 | 67, 400, 512, 43.15779, $ |
---|
99 | 68, 400, 512, 42.45604, $ |
---|
100 | 69, 400, 512, 41.75428, $ |
---|
101 | 70, 405, 512, 41.05253, $ |
---|
102 | 71, 432, 512, 40.35078, $ |
---|
103 | 72, 432, 512, 39.64903, $ |
---|
104 | 73, 432, 512, 38.94728, $ |
---|
105 | 74, 432, 512, 38.24552, $ |
---|
106 | 75, 432, 512, 37.54377, $ |
---|
107 | 76, 432, 512, 36.84202, $ |
---|
108 | 77, 432, 512, 36.14027, $ |
---|
109 | 78, 450, 512, 35.43851, $ |
---|
110 | 79, 450, 512, 34.73676, $ |
---|
111 | 80, 450, 512, 34.03501] |
---|
112 | n128 = [n128, $ |
---|
113 | 81, 450, 512, 33.33326, $ |
---|
114 | 82, 450, 512, 32.63150, $ |
---|
115 | 83, 480, 512, 31.92975, $ |
---|
116 | 84, 480, 512, 31.22800, $ |
---|
117 | 85, 480, 512, 30.52625, $ |
---|
118 | 86, 480, 512, 29.82449, $ |
---|
119 | 87, 480, 512, 29.12274, $ |
---|
120 | 88, 480, 512, 28.42099, $ |
---|
121 | 89, 480, 512, 27.71924, $ |
---|
122 | 90, 480, 512, 27.01748, $ |
---|
123 | 91, 480, 512, 26.31573, $ |
---|
124 | 92, 480, 512, 25.61398, $ |
---|
125 | 93, 486, 512, 24.91223, $ |
---|
126 | 94, 486, 512, 24.21047, $ |
---|
127 | 95, 486, 512, 23.50872, $ |
---|
128 | 96, 500, 512, 22.80697, $ |
---|
129 | 97, 500, 512, 22.10521, $ |
---|
130 | 98, 500, 512, 21.40346, $ |
---|
131 | 99, 500, 512, 20.70171, $ |
---|
132 | 100, 500, 512, 19.99996, $ |
---|
133 | 101, 500, 512, 19.29820, $ |
---|
134 | 102, 500, 512, 18.59645, $ |
---|
135 | 103, 512, 512, 17.89470, $ |
---|
136 | 104, 512, 512, 17.19294, $ |
---|
137 | 105, 512, 512, 16.49119, $ |
---|
138 | 106, 512, 512, 15.78944, $ |
---|
139 | 107, 512, 512, 15.08768, $ |
---|
140 | 108, 512, 512, 14.38593, $ |
---|
141 | 109, 512, 512, 13.68418, $ |
---|
142 | 110, 512, 512, 12.98243, $ |
---|
143 | 111, 512, 512, 12.28067, $ |
---|
144 | 112, 512, 512, 11.57892, $ |
---|
145 | 113, 512, 512, 10.87717, $ |
---|
146 | 114, 512, 512, 10.17541, $ |
---|
147 | 115, 512, 512, 9.47366, $ |
---|
148 | 116, 512, 512, 8.77191, $ |
---|
149 | 117, 512, 512, 8.07016, $ |
---|
150 | 118, 512, 512, 7.36840, $ |
---|
151 | 119, 512, 512, 6.66665, $ |
---|
152 | 120, 512, 512, 5.96490, $ |
---|
153 | 121, 512, 512, 5.26314, $ |
---|
154 | 122, 512, 512, 4.56139, $ |
---|
155 | 123, 512, 512, 3.85964, $ |
---|
156 | 124, 512, 512, 3.15788, $ |
---|
157 | 125, 512, 512, 2.45613, $ |
---|
158 | 126, 512, 512, 1.75438, $ |
---|
159 | 127, 512, 512, 1.05262, $ |
---|
160 | 128, 512, 512, 0.35087] |
---|
161 | ; |
---|
162 | n128 = reform(n128, 4, 128, /over) |
---|
163 | n128 = reform(n128[3, *], /over) |
---|
164 | n128 = [n128, -reverse(n128)] |
---|
165 | |
---|
166 | return, n128 |
---|
167 | end |
---|