1 | ******************************************************************************** |
---|
2 | * * |
---|
3 | * Input file for the Lagrangian particle dispersion model FLEXPART * |
---|
4 | * Please select your options * |
---|
5 | * * |
---|
6 | ******************************************************************************** |
---|
7 | |
---|
8 | 1. __ 3X, I2 |
---|
9 | 1 |
---|
10 | LDIRECT 1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION |
---|
11 | |
---|
12 | 2. ________ ______ 3X, I8, 1X, I6 |
---|
13 | 20040720 000000 |
---|
14 | YYYYMMDD HHMISS BEGINNING DATE OF SIMULATION |
---|
15 | |
---|
16 | 3. ________ ______ 3X, I8, 1X, I6 |
---|
17 | 20040721 120000 |
---|
18 | YYYYMMDD HHMISS ENDING DATE OF SIMULATION |
---|
19 | |
---|
20 | 4. _____ 3X, I5 |
---|
21 | 10800 |
---|
22 | SSSSS OUTPUT EVERY SSSSS SECONDS |
---|
23 | |
---|
24 | 5. _____ 3X, I5 |
---|
25 | 10800 |
---|
26 | SSSSS TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS) |
---|
27 | |
---|
28 | 6. _____ 3X, I5 |
---|
29 | 900 |
---|
30 | SSSSS SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS) |
---|
31 | |
---|
32 | 7. _________ 3X, I9 |
---|
33 | 999999999 |
---|
34 | SSSSSSSSS TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS) |
---|
35 | |
---|
36 | 8. _____ 3X, I5 |
---|
37 | 900 |
---|
38 | SSSSS SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS) |
---|
39 | |
---|
40 | 9. ---.-- 4X, F6.4 |
---|
41 | -5.0 |
---|
42 | CTL FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL |
---|
43 | |
---|
44 | 10. --- 4X, I3 |
---|
45 | 4 |
---|
46 | IFINE DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE |
---|
47 | |
---|
48 | 11. - 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 | |
---|
52 | 12. - 4X, I1 |
---|
53 | 0 |
---|
54 | IPOUT PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END |
---|
55 | |
---|
56 | 13. _ 4X, I1 |
---|
57 | 1 |
---|
58 | LSUBGRID SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO |
---|
59 | |
---|
60 | 14. _ 4X, I1 |
---|
61 | 1 |
---|
62 | LCONVECTION CONVECTION: 1 YES, 0 NO |
---|
63 | |
---|
64 | 15. _ 4X, I1 |
---|
65 | 0 |
---|
66 | LAGESPECTRA AGE SPECTRA: 1 YES, 0 NO |
---|
67 | |
---|
68 | 16. _ 4X, I1 |
---|
69 | 0 |
---|
70 | IPIN CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO |
---|
71 | |
---|
72 | 17. _ |
---|
73 | 0 4X,I1 |
---|
74 | IOFR IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO |
---|
75 | |
---|
76 | 18. _ 4X, I1 |
---|
77 | 0 |
---|
78 | IFLUX CALCULATE FLUXES: 1 YES, 0 NO |
---|
79 | |
---|
80 | 19. _ 4X, I1 |
---|
81 | 0 |
---|
82 | MDOMAINFILL DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT. O3 TRACER |
---|
83 | |
---|
84 | 20. _ 4X, I1 |
---|
85 | 1 |
---|
86 | IND_SOURCE 1=MASS UNIT , 2=MASS MIXING RATIO UNIT |
---|
87 | |
---|
88 | 21. _ 4X, I1 |
---|
89 | 1 |
---|
90 | IND_RECEPTOR 1=MASS UNIT , 2=MASS MIXING RATIO UNIT |
---|
91 | |
---|
92 | 22. _ 4X, I1 |
---|
93 | 0 |
---|
94 | MQUASILAG QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO |
---|
95 | |
---|
96 | 23. _ 4X, I1 |
---|
97 | 0 |
---|
98 | NESTED_OUTPUT SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO |
---|
99 | |
---|
100 | 24. _ 4X, I1 |
---|
101 | 2 |
---|
102 | LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT |
---|
103 | |
---|
104 | |
---|
105 | 1. Simulation direction, 1 for forward, -1 for backward in time |
---|
106 | (consult Seibert and Frank, 2004 for backward runs) |
---|
107 | |
---|
108 | 2. 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 | |
---|
112 | 3. Ending date and time of simulation. Same format as 3. |
---|
113 | |
---|
114 | 4. Average concentrations are calculated every SSSSS seconds. |
---|
115 | |
---|
116 | 5. The average concentrations are time averages of SSSSS seconds |
---|
117 | duration. If SSSSS is 0, instantaneous concentrations are outputted. |
---|
118 | |
---|
119 | 6. The concentrations are sampled every SSSSS seconds to calculate the time |
---|
120 | average concentration. This period must be shorter than the averaging time. |
---|
121 | |
---|
122 | 7. Time constant for particle splitting. Particles are split into two |
---|
123 | after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on. |
---|
124 | |
---|
125 | 8. 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 | |
---|
129 | 9. 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 | |
---|
132 | 10.IFINE=Reduction factor for time step used for vertical wind |
---|
133 | |
---|
134 | 11.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 | |
---|
139 | 12.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 | |
---|
144 | 13.Switch on/off subgridscale terrain parameterization (increase of |
---|
145 | mixing heights due to subgridscale orographic variations) |
---|
146 | |
---|
147 | 14.Switch on/off the convection parameterization |
---|
148 | |
---|
149 | 15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES |
---|
150 | must be available |
---|
151 | |
---|
152 | 16. 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 | |
---|
156 | 17. 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 | |
---|
160 | 18. 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 | |
---|
166 | 19. 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 | |
---|
172 | 20. 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 |
---|
177 | 21. IND_RECEPTOR switches between different units for concentrations at the receptor |
---|
178 | 1=mass units (concentrations) |
---|
179 | 2=mass mixing ratio units |
---|
180 | |
---|
181 | 22. 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 | |
---|
185 | 23. NESTED_OUTPUT decides whether model output shall be made also for a nested |
---|
186 | output field (normally with higher resolution) |
---|
187 | |
---|
188 | 24. 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. |
---|