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

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