source: trunk/SRC/Grid/n128gaussian.pro @ 74

Last change on this file since 74 was 13, checked in by pinsard, 18 years ago

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