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

Last change on this file since 378 was 371, checked in by pinsard, 16 years ago

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