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