source: trunk/SRC/Grid/n80gaussian.pro @ 114

Last change on this file since 114 was 114, checked in by smasson, 18 years ago

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