Project

General

Profile

afterburner.txt

Uwe Schulzweida, 2010-02-16 14:18

 
1
NAME
2

    
3
     after - transforms and selects data from T-model GRIB or netCDF files
4

    
5
SYNOPSIS
6

    
7
     after  [ -options ]  infiles  resultfile  <  namelist
8

    
9
DESCRIPTION
10

    
11
     The program "after" reads, unpacks, transforms and selects data
12
     from T-model GRIB/netCDF files, and writes them to
13
     SERVICE/GRIB/netCDF format for further processing (i.e. plotting).
14
     Additionally, pressure or height interpolation is optional
15
     (refer to parameter description of "UNITSEL").
16

    
17
     Possible combinations of TYPE, CODE and MEAN
18

    
19
     TYPE =  0 : CODE = 130  temperature           MEAN = 0
20
                        131  u-velocity
21
                        132  v-velocity
22
                        133  specific humidity
23
                        138  vorticity
24
                        148  streamfunction
25
                        149  velocity potential
26
                        152  LnPs
27
                        155  divergence
28

    
29
     TYPE = 10: CODE = (same as TYPE 0)            MEAN = 0
30

    
31
     TYPE = 11: CODE = (same as TYPE 0)            MEAN = 0
32

    
33
     TYPE > 11: CODE = (all codes)                 MEAN = 0/1
34

    
35
     Vorticity, divergence, streamfunction and velocity potential
36
     need special treatment in the vertical transformation.
37
     They are not available as types 30, 40 and 41.
38
     If you select one of these combinations, type is automatically
39
     switched to the equivalent types 70, 60 and 61. The type of all
40
     other variables will be switched too, because the type is a
41
     global parameter.
42

    
43
OPTIONS
44

    
45
     -P <nthreads> : Set number of OpenMP threads
46
     -c            : Print available codes and names
47
     -d            : Debug mode
48

    
49
NAMELIST
50

    
51
  &SELECT
52
     TYPE = , CODE = , LEVEL = , MULTI = , MEAN = , INTERVAL = , 
53
     TIMESEL = , UNITSEL = , FORMAT = , PRECISION = , SZIP =
54
  &END
55

    
56
  default ->  TYPE  =  0 : Hybrid             level spectral coefficients
57
              TYPE  = 10 : Hybrid             level fourier  coefficients
58
              TYPE  = 11 : Hybrid             level zonal mean sections
59
              TYPE  = 20 : Hybrid             level gauss grids
60
              TYPE  = 30 : Pressure or Height level gauss grids
61
              TYPE  = 40 : Pressure or Height level fourier  coefficients
62
              TYPE  = 41 : Pressure or Height level zonal mean sections
63
              TYPE  = 50 : Pressure or Height level spectral coefficients
64
              TYPE  = 60 : Pressure or Height level fourier  coefficients
65
              TYPE  = 61 : Pressure or Height level zonal mean sections
66
              TYPE  = 70 : Pressure or Height level gauss grids
67

    
68
  default ->  CODE  = -1 : process all detected codes
69
              CODE  =  1 - 255 ECMWF field codes (130=Temperature, etc.)
70
                     256 - 258 : internal
71
                           259 : sqrt(u*u + v*v)
72
                           260 : largescale + convective rain
73
                           261 : total top    radiation
74
                           262 : total bottom radiation
75
                           263 : net heat  flux
76
                           264 : net water flux
77

    
78
  default ->  LEVEL = -1 : process all detected levels
79
              LEVEL =: Up to NN + 1 hybrid, pressure or height levels
80
                       (depending on TYPE and UNITSEL !!!)
81
                       NN depends on model type (e.g. 19, 31 or 39)
82

    
83
  default ->  MULTI =  0 : Process only one input file
84
              MULTI =  n : Process "n" input files - each containing
85
                           one month. Put only the name of the first
86
                           input file on the command line. All subsequent
87
                           files are expected to be in the same directory.
88
                           The filenames must follow the syntax
89
                           experiment_YYMM, where YY denotes the year
90
                           and MM the month.
91

    
92
  default ->  DAYIN = 30 : Day step for incrementing input file names
93
              DAYIN = nn : Day step is nn days
94
                           Used only in conjunction with MULTI != 0.
95
                           Since ECHAM4 output files follow the
96
                           naming convention
97
                              experiment_YYMM.DD
98
                           DAYIN is the time step in days between
99
                           subsequent input files for MULTI != 0.
100

    
101
  default ->  MEAN  =  0 : Write out all terms
102
              MEAN  =  1 : Compute and write monthly or daily mean fields
103
                           !!! Only for TYPE >= 20 !!!
104
              MEAN  =  2 : Compute and write monthly or daily deviations.
105
                           !!! Only for TYPE = 20 or TYPE = 30 !!!
106
              MEAN  =  3 : Combination of MEAN=1 and MEAN=2.
107
                           Each mean field is followed by a deviation
108
                           field with an identical header record.
109
                           !!! Only for TYPE = 20 or TYPE = 30 !!!
110

    
111
  default ->  INTERVAL = 0 : Process data on monthly intervals
112
              INTERVAL = 1 : Process data on daily intervals
113
                             e.g. To compute daily averages
114
                                  set MEAN = 1 and INTERVAL = 1
115
              INTERVAL = 2 : Process data over all inputfiles
116
                             e.g. To compute an average over
117
                                  all inputfiles
118
                                  set MEAN = 1 and INTERVAL = 2
119

    
120
  default ->TIMESEL = -1 : The resultfile contains the same timesteps
121
                           per day as the input file.
122
                           Multiple values may be specified, i.e.
123
                              TIMESEL = 0, 12, 18
124
                           Note that only hours existing in the input
125
                           file can be specified.
126
                           (No time interpolation or average is done)
127

    
128
  default ->UNITSEL =  0 : Pressure interpolation (TYPE >= 30)
129
                           or hybrid level output (TYPE <= 20).
130
                           Unlike before, pressure levels less than 100 Pa
131
                           can be processed (except 0 and negative values).
132
            UNITSEL =  1 : Height   interpolation 
133
                           Values specified for LEVEL are interpreted
134
                           as height-meters       (TYPE >= 30).
135
                           (Maximum LEVEL: 65535 m)
136
            UNITSEL =  2 : Height   interpolation 
137
                           Values specified for LEVEL are interpreted
138
                           as height-kilometers   (TYPE >= 30).
139
                           (Maximum LEVEL: 65   km)
140

    
141
  default -> FORMAT =  0 : SERVICE
142
             FORMAT =  1 : GRIB
143
             FORMAT =  2 : netCDF
144
             FORMAT =  3 : EXTRA
145
	     FORMAT =  4 : netCDF version 2
146
	     FORMAT =  6 : netCDF version 4
147

    
148
  default -> PRECISION =  0 : use input data precision for output
149
             PRECISION =  8 :  8 bit packed output data precision
150
             PRECISION = 16 : 16 bit packed output data precision
151
             PRECISION = 24 : 24 bit packed output data precision
152
             PRECISION = 32 : 32 bit floating point output data precision
153
             PRECISION = 64 : 64 bit floating point output data precision
154

    
155
  default -> SZIP = 0 : No SZIP compression of GRIB records
156
             SZIP = 1 : SZIP compression of GRIB records
157

    
158
             EXTRAPOLATE = 0 : No extrapolation of missing values during the
159
                               interpolation from model to pressure or height level
160
                               (only available with MEAN = 0 and TYPE = 30)
161
  default -> EXTRAPOLATE = 1 : Extrapolate missing values
162

    
163

    
164
EXAMPLES
165

    
166
     To create monthly mean over all codes and hybrid levels in GRIB
167
     format, try
168

    
169
       after  Model-file  GRIB-file.avg  << EOF
170
          &SELECT TYPE=20, MEAN=1, FORMAT=1 &END
171
       EOF
172

    
173
     To select code 129 and 130 and interpolate to pressure levels
174
     100, 250, 500, 850, 925 and 1000 hPa, try
175

    
176
       after  Model-file  GRIB-file.sel  << EOF
177
          &SELECT
178
            TYPE=70, FORMAT=1,
179
            CODE=129,130,
180
            LEVEL=10000, 25000, 50000, 85000, 92500, 100000
181
          &END
182
       EOF
183

    
184
     To calculate u and v velocity from vorticity and divergence
185
     on model level, try
186

    
187
       after  Model-file  GRIB-file.sel  << EOF
188
          &SELECT
189
            TYPE=20, FORMAT=1,
190
            CODE=131,132
191
          &END
192
       EOF
193