[4bd0f2d] | 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 | 20110310 000000 |
---|
| 14 | YYYYMMDD HHMISS BEGINNING DATE OF SIMULATION |
---|
| 15 | |
---|
| 16 | 3. ________ ______ 3X, I8, 1X, I6 |
---|
| 17 | 20110310 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 | 25. _ 4X, I1 |
---|
| 105 | 0 |
---|
| 106 | SURF_ONLY IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER |
---|
| 107 | |
---|
| 108 | |
---|
| 109 | 1. Simulation direction, 1 for forward, -1 for backward in time |
---|
| 110 | (consult Seibert and Frank, 2004 for backward runs) |
---|
| 111 | |
---|
| 112 | 2. Beginning date and time of simulation. Must be given in format |
---|
| 113 | YYYYMMDD HHMISS, where YYYY is YEAR, MM is MONTH, DD is DAY, HH is HOUR, |
---|
| 114 | MI is MINUTE and SS is SECOND. Current version utilizes UTC. |
---|
| 115 | |
---|
| 116 | 3. Ending date and time of simulation. Same format as 3. |
---|
| 117 | |
---|
| 118 | 4. Average concentrations are calculated every SSSSS seconds. |
---|
| 119 | |
---|
| 120 | 5. The average concentrations are time averages of SSSSS seconds |
---|
| 121 | duration. If SSSSS is 0, instantaneous concentrations are outputted. |
---|
| 122 | |
---|
| 123 | 6. The concentrations are sampled every SSSSS seconds to calculate the time |
---|
| 124 | average concentration. This period must be shorter than the averaging time. |
---|
| 125 | |
---|
| 126 | 7. Time constant for particle splitting. Particles are split into two |
---|
| 127 | after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on. |
---|
| 128 | |
---|
| 129 | 8. All processes are synchronized with this time interval (lsynctime). |
---|
| 130 | Therefore, all other time constants must be multiples of this value. |
---|
| 131 | Output interval and time average of output must be at least twice lsynctime. |
---|
| 132 | |
---|
| 133 | 9. CTL must be >1 for time steps shorter than the Lagrangian time scale |
---|
| 134 | If CTL<0, a purely random walk simulation is done |
---|
| 135 | |
---|
| 136 | 10.IFINE=Reduction factor for time step used for vertical wind |
---|
| 137 | |
---|
| 138 | 11.IOUT determines how the output shall be made: concentration |
---|
| 139 | (ng/m3, Bq/m3), mixing ratio (pptv), or both, or plume trajectory mode, |
---|
| 140 | or concentration + plume trajectory mode. |
---|
| 141 | In plume trajectory mode, output is in the form of average trajectories. |
---|
| 142 | |
---|
| 143 | 12.IPOUT determines whether particle positions are outputted (in addition |
---|
| 144 | to the gridded concentrations or mixing ratios) or not. |
---|
| 145 | 0=no output, 1 output every output interval, 2 only at end of the |
---|
| 146 | simulation |
---|
| 147 | |
---|
| 148 | 13.Switch on/off subgridscale terrain parameterization (increase of |
---|
| 149 | mixing heights due to subgridscale orographic variations) |
---|
| 150 | |
---|
| 151 | 14.Switch on/off the convection parameterization |
---|
| 152 | |
---|
| 153 | 15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES |
---|
| 154 | must be available |
---|
| 155 | |
---|
| 156 | 16. If IPIN=1, a file "partposit_end" from a previous run must be available in |
---|
| 157 | the output directory. Particle positions are read in and previous simulation |
---|
| 158 | is continued. If IPIN=0, no particles from a previous run are used |
---|
| 159 | |
---|
| 160 | 17. IF IOUTPUTFOREACHRELEASE is set to 1, one output field for each location |
---|
| 161 | in the RLEASE file is created. For backward calculation this should be |
---|
| 162 | set to 1. For forward calculation both possibilities are applicable. |
---|
| 163 | |
---|
| 164 | 18. If IFLUX is set to 1, fluxes of each species through each of the output |
---|
| 165 | boxes are calculated. Six fluxes, corresponding to northward, southward, |
---|
| 166 | eastward, westward, upward and downward are calculated for each grid cell of |
---|
| 167 | the output grid. The control surfaces are placed in the middle of each |
---|
| 168 | output grid cell. If IFLUX is set to 0, no fluxes are determined. |
---|
| 169 | |
---|
| 170 | 19. If MDOMAINFILL is set to 1, the first box specified in file RELEASES is used |
---|
| 171 | as the domain where domain-filling trajectory calculations are to be done. |
---|
| 172 | Particles are initialized uniformly distributed (according to the air mass |
---|
| 173 | distribution) in that domain at the beginning of the simulation, and are |
---|
| 174 | created at the boundaries throughout the simulation period. |
---|
| 175 | |
---|
| 176 | 20. IND_SOURCE switches between different units for concentrations at the source |
---|
| 177 | NOTE that in backward simulations the release of computational particles |
---|
| 178 | takes place at the "receptor" and the sampling of particles at the "source". |
---|
| 179 | 1=mass units (for bwd-runs = concentration) |
---|
| 180 | 2=mass mixing ratio units |
---|
| 181 | 21. IND_RECEPTOR switches between different units for concentrations at the receptor |
---|
| 182 | 1=mass units (concentrations) |
---|
| 183 | 2=mass mixing ratio units |
---|
| 184 | |
---|
| 185 | 22. MQUASILAG indicates whether particles shall be numbered consecutively (1) or |
---|
| 186 | with their release location number (0). The first option allows tracking of |
---|
| 187 | individual particles using the partposit output files |
---|
| 188 | |
---|
| 189 | 23. NESTED_OUTPUT decides whether model output shall be made also for a nested |
---|
| 190 | output field (normally with higher resolution) |
---|
| 191 | |
---|
| 192 | 24. LINIT_COND determines whether, for backward runs only, the sensitivity to initial |
---|
| 193 | conditions shall be calculated and written to output files |
---|
| 194 | 0=no output, 1 or 2 determines in which units the initial conditions are provided. |
---|
| 195 | |
---|
| 196 | 25. SURF_ONLY: When set to 1, concentration/emission sensitivity is written out only |
---|
| 197 | for the surface layer; useful for instance when only footprint emission sensitivity is needed |
---|
| 198 | but initial conditions are needed on a full 3-D grid |
---|