Ticket #12: COMMAND

File COMMAND, 7.3 KB (added by victoria.sinclair@…, 9 years ago)
Line 
1********************************************************************************
2*                                                                              *
3*      Input file for the Lagrangian particle dispersion model FLEXPART        *
4*                           Please select your options                         *
5*                                                                              *
6********************************************************************************
7
81. __                3X, I2
9    1       
10   LDIRECT           1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION
11
122. ________ ______   3X, I8, 1X, I6
13   20040720 000000
14   YYYYMMDD HHMISS   BEGINNING DATE OF SIMULATION
15
163. ________ ______   3X, I8, 1X, I6
17   20040721 120000
18   YYYYMMDD HHMISS   ENDING DATE OF SIMULATION
19
204. _____             3X, I5
21   10800
22   SSSSS             OUTPUT EVERY SSSSS SECONDS
23
245. _____             3X, I5
25   10800
26   SSSSS             TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
27
286. _____             3X, I5
29     900
30   SSSSS             SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
31
327. _________         3X, I9
33   999999999
34   SSSSSSSSS         TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)
35
368. _____             3X, I5
37     900
38   SSSSS             SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
39
409.  ---.--           4X, F6.4
41     -5.0
42    CTL              FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL
43
4410. ---              4X, I3
45      4
46    IFINE            DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE
47
4811. -                4X, I1
49    3 
50    IOUT             1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4
51
5212. -                4X, I1
53    0 
54    IPOUT            PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END
55
5613. _                4X, I1
57    1
58    LSUBGRID         SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO
59
6014. _                4X, I1
61    1
62    LCONVECTION      CONVECTION: 1 YES, 0 NO
63
6415. _                4X, I1
65    0
66    LAGESPECTRA      AGE SPECTRA: 1 YES, 0 NO
67
6816. _                4X, I1
69    0
70    IPIN             CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO
71
7217. _               
73    0                4X,I1
74    IOFR             IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
75
7618. _                4X, I1
77    0
78    IFLUX            CALCULATE FLUXES: 1 YES, 0 NO
79
8019. _                4X, I1
81    0
82    MDOMAINFILL      DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT. O3 TRACER
83
8420. _                4X, I1
85    1
86    IND_SOURCE       1=MASS UNIT , 2=MASS MIXING RATIO UNIT
87
8821. _                4X, I1
89    1
90    IND_RECEPTOR     1=MASS UNIT , 2=MASS MIXING RATIO UNIT
91
9222. _                4X, I1
93    0
94    MQUASILAG        QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO
95
9623. _                4X, I1
97    0
98    NESTED_OUTPUT    SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO
99
10024. _                4X, I1
101    2
102    LINIT_COND       INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
103
104
1051. Simulation direction, 1 for forward, -1 for backward in time
106        (consult Seibert and Frank, 2004 for backward runs)
107
1082. Beginning date and time of simulation. Must be given in format
109   YYYYMMDD HHMISS, where YYYY is YEAR, MM is MONTH, DD is DAY, HH is HOUR,
110   MI is MINUTE and SS is SECOND. Current version utilizes UTC.
111
1123. Ending date and time of simulation. Same format as 3.
113
1144. Average concentrations are calculated every SSSSS seconds.
115
1165. The average concentrations are time averages of SSSSS seconds
117   duration. If SSSSS is 0, instantaneous concentrations are outputted.
118
1196. The concentrations are sampled every SSSSS seconds to calculate the time
120   average concentration. This period must be shorter than the averaging time.
121
1227. Time constant for particle splitting. Particles are split into two
123   after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on.
124
1258. All processes are synchronized with this time interval (lsynctime).
126   Therefore, all other time constants must be multiples of this value.
127   Output interval and time average of output must be at least twice lsynctime.
128
1299. CTL must be >1 for time steps shorter than the  Lagrangian time scale
130   If CTL<0, a purely random walk simulation is done
131
13210.IFINE=Reduction factor for time step used for vertical wind
133
13411.IOUT determines how the output shall be made: concentration
135   (ng/m3, Bq/m3), mixing ratio (pptv), or both, or plume trajectory mode,
136   or concentration + plume trajectory mode.
137   In plume trajectory mode, output is in the form of average trajectories.
138
13912.IPOUT determines whether particle positions are outputted (in addition
140   to the gridded concentrations or mixing ratios) or not.
141   0=no output, 1 output every output interval, 2 only at end of the
142   simulation
143
14413.Switch on/off subgridscale terrain parameterization (increase of
145   mixing heights due to subgridscale orographic variations)
146
14714.Switch on/off the convection parameterization
148
14915.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES
150   must be available
151
15216. If IPIN=1, a file "partposit_end" from a previous run must be available in
153    the output directory. Particle positions are read in and previous simulation
154    is continued. If IPIN=0, no particles from a previous run are used
155
15617. IF IOUTPUTFOREACHRELEASE is set to 1, one output field for each location
157    in the RLEASE file is created. For backward calculation this should be
158    set to 1. For forward calculation both possibilities are applicable.
159
16018. If IFLUX is set to 1, fluxes of each species through each of the output
161    boxes are calculated. Six fluxes, corresponding to northward, southward,
162    eastward, westward, upward and downward are calculated for each grid cell of
163    the output grid. The control surfaces are placed in the middle of each
164    output grid cell. If IFLUX is set to 0, no fluxes are determined.
165
16619. If MDOMAINFILL is set to 1, the first box specified in file RELEASES is used
167    as the domain where domain-filling trajectory calculations are to be done.
168    Particles are initialized uniformly distributed (according to the air mass
169    distribution) in that domain at the beginning of the simulation, and are
170    created at the boundaries throughout the simulation period.
171
17220. IND_SOURCE switches between different units for concentrations at the source
173    NOTE that in backward simulations the release of computational particles
174    takes place at the "receptor" and the sampling of particles at the "source".
175          1=mass units (for bwd-runs = concentration)
176          2=mass mixing ratio units
17721. IND_RECEPTOR switches between different units for concentrations at the receptor
178          1=mass units (concentrations)
179          2=mass mixing ratio units
180
18122. MQUASILAG indicates whether particles shall be numbered consecutively (1) or
182    with their release location number (0). The first option allows tracking of
183    individual particles using the partposit output files
184
18523. NESTED_OUTPUT decides whether model output shall be made also for a nested
186    output field (normally with higher resolution)
187
18824. LINIT_COND determines whether, for backward runs only, the sensitivity to initial
189    conditions shall be calculated and written to output files
190    0=no output, 1 or 2 determines in which units the initial conditions are provided.
hosted by ZAMG