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

Last change on this file since 157 was 157, checked in by navarro, 18 years ago

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