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