[e200b7a] | 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 |
---|
[4fbe7a5] | 32 | 0 SURF_ONLY IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER |
---|
[e200b7a] | 33 | |
---|
| 34 | |
---|
| 35 | 1. Simulation direction, 1 for forward, -1 for backward in time |
---|
| 36 | |
---|
| 37 | 2. Beginning date and time of simulation. Must be given in format |
---|
| 38 | YYYYMMDD HHMISS, where YYYY is YEAR, MM is MONTH, DD is DAY, HH is HOUR, |
---|
| 39 | MI is MINUTE and SS is SECOND. Current version utilizes UTC. |
---|
| 40 | |
---|
| 41 | 3. Ending date and time of simulation. Same format as 3. |
---|
| 42 | |
---|
| 43 | 4. Average concentrations are calculated every SSSSS seconds. |
---|
| 44 | |
---|
| 45 | 5. The average concentrations are time averages of SSSSS seconds |
---|
| 46 | duration. If SSSSS is 0, instantaneous concentrations are outputted. |
---|
| 47 | |
---|
| 48 | 6. The concentrations are sampled every SSSSS seconds to calculate the time |
---|
| 49 | average concentration. This period must be shorter than the averaging time. |
---|
| 50 | |
---|
| 51 | 7. Time constant for particle splitting. Particles are split into two |
---|
| 52 | after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on. |
---|
| 53 | |
---|
| 54 | 8. All processes are synchronized with this time interval (lsynctime). |
---|
| 55 | Therefore, all other time constants must be multiples of this value. |
---|
| 56 | Output interval and time average of output must be at least twice lsynctime. |
---|
| 57 | |
---|
| 58 | 9. CTL must be >1 for time steps shorter than the Lagrangian time scale |
---|
| 59 | If CTL<0, a purely random walk simulation is done |
---|
| 60 | |
---|
| 61 | 10.IFINE=Reduction factor for time step used for vertical wind |
---|
| 62 | |
---|
| 63 | 11.IOUT determines how the output shall be made: concentration |
---|
| 64 | (ng/m3, Bq/m3), mixing ratio (pptv), or both, or plume trajectory mode, |
---|
| 65 | or concentration + plume trajectory mode. |
---|
| 66 | In plume trajectory mode, output is in the form of average trajectories. |
---|
| 67 | |
---|
| 68 | 12.IPOUT determines whether particle positions are outputted (in addition |
---|
| 69 | to the gridded concentrations or mixing ratios) or not. |
---|
| 70 | 0=no output, 1 output every output interval, 2 only at end of the |
---|
| 71 | simulation |
---|
| 72 | |
---|
| 73 | 13.Switch on/off subgridscale terrain parameterization (increase of |
---|
| 74 | mixing heights due to subgridscale orographic variations) |
---|
| 75 | |
---|
| 76 | 14.Switch on/off the convection parameterization |
---|
| 77 | |
---|
| 78 | 15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES |
---|
| 79 | must be available |
---|
| 80 | |
---|
| 81 | 16. If IPIN=1, a file "partposit_end" from a previous run must be available in |
---|
| 82 | the output directory. Particle positions are read in and previous simulation |
---|
| 83 | is continued. If IPIN=0, no particles from a previous run are used |
---|
| 84 | |
---|
| 85 | 17. IF IOUTPUTFOREACHRELEASE is set to 1, one output field for each location |
---|
| 86 | in the RLEASE file is created. For backward calculation this should be |
---|
| 87 | set to 1. For forward calculation both possibilities are applicable. |
---|
| 88 | |
---|
| 89 | 18. If IFLUX is set to 1, fluxes of each species through each of the output |
---|
| 90 | boxes are calculated. Six fluxes, corresponding to northward, southward, |
---|
| 91 | eastward, westward, upward and downward are calculated for each grid cell of |
---|
| 92 | the output grid. The control surfaces are placed in the middle of each |
---|
| 93 | output grid cell. If IFLUX is set to 0, no fluxes are determined. |
---|
| 94 | |
---|
| 95 | 19. If MDOMAINFILL is set to 1, the first box specified in file RELEASES is used |
---|
| 96 | as the domain where domain-filling trajectory calculations are to be done. |
---|
| 97 | Particles are initialized uniformly distributed (according to the air mass |
---|
| 98 | distribution) in that domain at the beginning of the simulation, and are |
---|
| 99 | created at the boundaries throughout the simulation period. |
---|
| 100 | |
---|
| 101 | 20. IND_SOURCE switches between different units for concentrations at the source |
---|
| 102 | NOTE that in backward simulations the release of computational particles |
---|
| 103 | takes place at the "receptor" and the sampling of particles at the "source". |
---|
| 104 | 1=mass units (for bwd-runs = concentration) |
---|
| 105 | 2=mass mixing ratio units |
---|
| 106 | 21. IND_RECEPTOR switches between different units for concentrations at the receptor |
---|
| 107 | 1=mass units (concentrations) |
---|
| 108 | 2=mass mixing ratio units |
---|
| 109 | |
---|
| 110 | 22. MQUASILAG indicates whether particles shall be numbered consecutively (1) or |
---|
| 111 | with their release location number (0). The first option allows tracking of |
---|
| 112 | individual particles using the partposit output files |
---|
| 113 | |
---|
| 114 | 23. NESTED_OUTPUT decides whether model output shall be made also for a nested |
---|
| 115 | output field (normally with higher resolution) |
---|
| 116 | |
---|
| 117 | 24. LINIT_COND determines whether, for backward runs only, the sensitivity to initial |
---|
| 118 | conditions shall be calculated and written to output files |
---|
| 119 | 0=no output, 1 or 2 determines in which units the initial conditions are provided. |
---|
[4fbe7a5] | 120 | |
---|
| 121 | 25. SURF_ONLY: When set to 1, concentration/emission sensitivity is written out only |
---|
| 122 | for the surface layer; useful for instance when only footprint emission sensitivity is needed |
---|
| 123 | but initial conditions are needed on a full 3-D grid |
---|