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