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