| 1 | ******************************************************************************** |
|---|
| 2 | * * |
|---|
| 3 | * Input file for the Lagrangian particle dispersion model FLEXPART * |
|---|
| 4 | * Please select your options * |
|---|
| 5 | * * |
|---|
| 6 | ******************************************************************************** |
|---|
| 7 | |
|---|
| 8 | 1 LDIRECT 1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION |
|---|
| 9 | 20040720 000000 YYYYMMDD HHMISS BEGINNING DATE OF SIMULATION |
|---|
| 10 | 20040721 120000 YYYYMMDD HHMISS ENDING DATE OF SIMULATION |
|---|
| 11 | 10800 SSSSS OUTPUT EVERY SSSSS SECONDS |
|---|
| 12 | 10800 SSSSS TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS) |
|---|
| 13 | 900 SSSSS SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS) |
|---|
| 14 | 9999999 SSSSSSS TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS) |
|---|
| 15 | 900 SSSSS SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS) |
|---|
| 16 | -5.0 CTL FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL |
|---|
| 17 | 4 IFINE DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE |
|---|
| 18 | 3 IOUT 1 CONC. (RESID. TIME FOR BACKWARD RUNS) OUTPUT,2 MIX. RATIO OUTPUT,3 BOTH,4 PLUME TRAJECT.,5=1+4 |
|---|
| 19 | 0 IPOUT PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END |
|---|
| 20 | 1 LSUBGRID SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO |
|---|
| 21 | 1 LCONVECTION CONVECTION: 1 YES, 0 NO |
|---|
| 22 | 0 LAGESPECTRA AGE SPECTRA: 1 YES, 0 NO |
|---|
| 23 | 0 IPIN CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO |
|---|
| 24 | 0 IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO |
|---|
| 25 | 0 IFLUX CALCULATE FLUXES: 1 YES, 0 NO |
|---|
| 26 | 0 MDOMAINFILL DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO |
|---|
| 27 | 1 IND_SOURCE 1=MASS UNIT , 2=MASS MIXING RATIO UNIT |
|---|
| 28 | 1 IND_RECEPTOR 1=MASS UNIT , 2=MASS MIXING RATIO UNIT |
|---|
| 29 | 0 MQUASILAG QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO |
|---|
| 30 | 0 NESTED_OUTPUT SHALL NESTED OUTPUT BE USED? YES, 0 NO |
|---|
| 31 | 2 LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT |
|---|
| 32 | 0 SURF_ONLY IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER |
|---|
| 33 | .f. LDEP_INCR IF .TRUE., INCREMENTAL DEPOSITION, ELSE ACCUMULATED DEPOSITION (=DEFAULT) |
|---|
| 34 | |
|---|
| 35 | |
|---|
| 36 | 1. Simulation direction, 1 for forward, -1 for backward in time |
|---|
| 37 | |
|---|
| 38 | 2. Beginning date and time of simulation. Must be given in format |
|---|
| 39 | YYYYMMDD HHMISS, where YYYY is YEAR, MM is MONTH, DD is DAY, HH is HOUR, |
|---|
| 40 | MI is MINUTE and SS is SECOND. Current version utilizes UTC. |
|---|
| 41 | |
|---|
| 42 | 3. Ending date and time of simulation. Same format as 3. |
|---|
| 43 | |
|---|
| 44 | 4. Average concentrations are calculated every SSSSS seconds. |
|---|
| 45 | |
|---|
| 46 | 5. The average concentrations are time averages of SSSSS seconds |
|---|
| 47 | duration. If SSSSS is 0, instantaneous concentrations are outputted. |
|---|
| 48 | |
|---|
| 49 | 6. The concentrations are sampled every SSSSS seconds to calculate the time |
|---|
| 50 | average concentration. This period must be shorter than the averaging time. |
|---|
| 51 | |
|---|
| 52 | 7. Time constant for particle splitting. Particles are split into two |
|---|
| 53 | after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on. |
|---|
| 54 | |
|---|
| 55 | 8. All processes are synchronized with this time interval (lsynctime). |
|---|
| 56 | Therefore, all other time constants must be multiples of this value. |
|---|
| 57 | Output interval and time average of output must be at least twice lsynctime. |
|---|
| 58 | |
|---|
| 59 | 9. CTL must be >1 for time steps shorter than the Lagrangian time scale |
|---|
| 60 | If CTL<0, a purely random walk simulation is done |
|---|
| 61 | |
|---|
| 62 | 10.IFINE=Reduction factor for time step used for vertical wind |
|---|
| 63 | |
|---|
| 64 | 11.IOUT determines how the output shall be made: concentration |
|---|
| 65 | (ng/m3, Bq/m3), mixing ratio (pptv), or both, or plume trajectory mode, |
|---|
| 66 | or concentration + plume trajectory mode. |
|---|
| 67 | In plume trajectory mode, output is in the form of average trajectories. |
|---|
| 68 | |
|---|
| 69 | 12.IPOUT determines whether particle positions are outputted (in addition |
|---|
| 70 | to the gridded concentrations or mixing ratios) or not. |
|---|
| 71 | 0=no output, 1 output every output interval, 2 only at end of the |
|---|
| 72 | simulation |
|---|
| 73 | |
|---|
| 74 | 13.Switch on/off subgridscale terrain parameterization (increase of |
|---|
| 75 | mixing heights due to subgridscale orographic variations) |
|---|
| 76 | |
|---|
| 77 | 14.Switch on/off the convection parameterization |
|---|
| 78 | |
|---|
| 79 | 15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES |
|---|
| 80 | must be available |
|---|
| 81 | |
|---|
| 82 | 16. If IPIN=1, a file "partposit_end" from a previous run must be available in |
|---|
| 83 | the output directory. Particle positions are read in and previous simulation |
|---|
| 84 | is continued. If IPIN=0, no particles from a previous run are used |
|---|
| 85 | |
|---|
| 86 | 17. IF IOUTPUTFOREACHRELEASE is set to 1, one output field for each location |
|---|
| 87 | in the RLEASE file is created. For backward calculation this should be |
|---|
| 88 | set to 1. For forward calculation both possibilities are applicable. |
|---|
| 89 | |
|---|
| 90 | 18. If IFLUX is set to 1, fluxes of each species through each of the output |
|---|
| 91 | boxes are calculated. Six fluxes, corresponding to northward, southward, |
|---|
| 92 | eastward, westward, upward and downward are calculated for each grid cell of |
|---|
| 93 | the output grid. The control surfaces are placed in the middle of each |
|---|
| 94 | output grid cell. If IFLUX is set to 0, no fluxes are determined. |
|---|
| 95 | |
|---|
| 96 | 19. If MDOMAINFILL is set to 1, the first box specified in file RELEASES is used |
|---|
| 97 | as the domain where domain-filling trajectory calculations are to be done. |
|---|
| 98 | Particles are initialized uniformly distributed (according to the air mass |
|---|
| 99 | distribution) in that domain at the beginning of the simulation, and are |
|---|
| 100 | created at the boundaries throughout the simulation period. |
|---|
| 101 | |
|---|
| 102 | 20. IND_SOURCE switches between different units for concentrations at the source |
|---|
| 103 | NOTE that in backward simulations the release of computational particles |
|---|
| 104 | takes place at the "receptor" and the sampling of particles at the "source". |
|---|
| 105 | 1=mass units (for bwd-runs = concentration) |
|---|
| 106 | 2=mass mixing ratio units |
|---|
| 107 | 21. IND_RECEPTOR switches between different units for concentrations at the receptor |
|---|
| 108 | 1=mass units (concentrations) |
|---|
| 109 | 2=mass mixing ratio units |
|---|
| 110 | |
|---|
| 111 | 22. MQUASILAG indicates whether particles shall be numbered consecutively (1) or |
|---|
| 112 | with their release location number (0). The first option allows tracking of |
|---|
| 113 | individual particles using the partposit output files |
|---|
| 114 | |
|---|
| 115 | 23. NESTED_OUTPUT decides whether model output shall be made also for a nested |
|---|
| 116 | output field (normally with higher resolution) |
|---|
| 117 | |
|---|
| 118 | 24. LINIT_COND determines whether, for backward runs only, the sensitivity to initial |
|---|
| 119 | conditions shall be calculated and written to output files |
|---|
| 120 | 0=no output, 1 or 2 determines in which units the initial conditions are provided. |
|---|
| 121 | |
|---|
| 122 | 25. SURF_ONLY: When set to 1, concentration/emission sensitivity is written out only |
|---|
| 123 | for the surface layer; useful for instance when only footprint emission sensitivity is needed |
|---|
| 124 | but initial conditions are needed on a full 3-D grid |
|---|
