********************************************************************************
*                                                                              *
*      Input file for the Lagrangian particle dispersion model FLEXPART        *
*                           Please select your options                         *
*                                                                              *
********************************************************************************

1. __                3X, I2
    1       
   LDIRECT           1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION

2. ________ ______   3X, I8, 1X, I6
   20040101 000000
   YYYYMMDD HHMISS   BEGINNING DATE OF SIMULATION

3. ________ ______   3X, I8, 1X, I6
   20040201 000000 
   YYYYMMDD HHMISS   ENDING DATE OF SIMULATION

4. _____             3X, I5
 2678400
   SSSSS             GRID OUTPUT EVERY SSSSS SECONDS

5. _____             3X, I5
 2678400
   SSSSS             TIME AVERAGE OF GRID OUTPUT (IN SSSSS SECONDS)

6. _____             3X, I5
     900 
   SSSSS             SAMPLING RATE OF GRID OUTPUT (IN SSSSS SECONDS)

7. _____             3X, I5
       1 
                     Number of UNCERTAINTY CLASSES OF GRID OUTPUT

8. _____             3X, I5
   86400 
   SSSSS             RECEPTOR OUTPUT EVERY SSSSS SECONDS

9. _____             3X, I5
   86400
   SSSSS             TIME AVERAGE OF RECEPTOR OUTPUT (IN SSSSS SECONDS)

10._____             3X, I5
     900 
   SSSSS             SAMPLING RATE OF RECEPTOR OUTPUT (IN SSSSS SECONDS)

11._____             3X, I5
       9 
                     Number of UNCERTAINTY CLASSES OF RECEPTOR OUTPUT

12._________         3X, I9
   999999999
   SSSSSSSSS         TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)

13._____             3X, I5
     900 
   SSSSS             SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)

14. ---.--           4X, F6.4
     -5.0
    CTL              FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL

15. ---              4X, I3
      4
    IFINE            DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE

16. -                4X, I1
    2  
    IOUT             1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4

17. -                4X, I1
    2  
    IPOUT            PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END

18. _                4X, I1
    1
    LSUBGRID         SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO

19. _                4X, I1
    1
    LCONVECTION      CONVECTION: 1 YES, 0 NO

20. _                4X, I1
    0
    LAGESPECTRA      AGE SPECTRA: 1 YES, 0 NO

21. _                4X, I1
     0
    IPIN             CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO

22. _                
    0                4X,I1 
    IOFR             IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO

23. _                4X, I1
    0
    IFLUX            CALCULATE FLUXES: 1 YES, 0 NO

24. _                4X, I1
    1
    MDOMAINFILL      DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT. O3 TRACER

25. _                4X, I1
    1
    IND_SOURCE       1=MASS UNIT , 2=MASS MIXING RATIO UNIT 

26. _                4X, I1
    2
    IND_RECEPTOR     1=MASS UNIT , 2=MASS MIXING RATIO UNIT 

27. _                4X, I1
    0
    MQUASILAG        QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO

28. _                4X, I1
    0
    NESTED_OUTPUT    SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO

29. __________       4X, I10 
   2678400
    SSSSSSSSSS       Interval for mass shift (0: no mass shift)

30. _                A
/mnt/storage/spc134/FLEXPART80CTM/input/URMEL/emissions/
    EM_PATH          PATH TO EMISSION FIELDS, IF EMPTY EMISSION SWITCHED OFF, (include / at the end)

31. _                A
/mnt/storage/spc134/FLEXPART80CTM/input/URMEL/regions/
    REG_PATH          PATH TO REGION MASKS (include / at the end)

32. _                A
/mnt/storage/spc134/FLEXPART80CTM/input/URMEL/OH/
    OH_PATH          PATH TO OH fields (include / at the end)
    
33._                 A
/mnt/storage/spc134/FLEXPART80CTM/input/URMEL/CLR/
    CLR_PATH         PATH TO add. chemical loss rate fields (include / at the end)

1. Simulation direction, 1 for forward, -1 for backward in time 
	(consult Seibert and Frank, 2004 for backward runs)

2. Beginning date and time of simulation. Must be given in format
   YYYYMMDD HHMISS, where YYYY is YEAR, MM is MONTH, DD is DAY, HH is HOUR,
   MI is MINUTE and SS is SECOND. Current version utilizes UTC.

3. Ending date and time of simulation. Same format as 3.

4. Average concentrations are calculated every SSSSS seconds.

5. The average concentrations are time averages of SSSSS seconds
   duration. If SSSSS is 0, instantaneous concentrations are outputted.

6. The concentrations are sampled every SSSSS seconds to calculate the time
   average concentration. This period must be shorter than the averaging time.

7. Number of uncertainty classes used for grid output. Has to be equal to, a multiple of, or 
   a divisor of the number of uncertainty classes for receptor output. 

8. Average receptor concentrations are calculated every SSSSS seconds.

9. The average receptor concentrations are time averages of SSSSS seconds
   duration. If SSSSS is 0, instantaneous concentrations are outputted.

10. The receptor concentrations are sampled every SSSSS seconds to calculate the time
   average concentration. This period must be shorter than the averaging time.

11. Number of uncertainty classes used for receptor output. Has to be equal to, a multiple of, or 
   a divisor of the number of uncertainty classes for grid output. 

12. Time constant for particle splitting. Particles are split into two
   after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on.

13. All processes are synchronized with this time interval (lsynctime).
   Therefore, all other time constants must be multiples of this value.
   Output interval and time average of output must be at least twice lsynctime.

14. CTL must be >1 for time steps shorter than the  Lagrangian time scale
   If CTL<0, a purely random walk simulation is done

15.IFINE=Reduction factor for time step used for vertical wind

16.IOUT determines how the output shall be made: concentration
   (ng/m3, Bq/m3), mixing ratio (pptv), or both, or plume trajectory mode,
   or concentration + plume trajectory mode.
   In plume trajectory mode, output is in the form of average trajectories.

17.IPOUT determines whether particle positions are outputted (in addition
   to the gridded concentrations or mixing ratios) or not.
   0=no output, 1 output every output interval, 2 only at end of the
   simulation

18.Switch on/off subgridscale terrain parameterization (increase of
   mixing heights due to subgridscale orographic variations)

19.Switch on/off the convection parameterization

20.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES
   must be available

21. If IPIN=1, a file "partposit_end" from a previous run must be available in
    the output directory. Particle positions are read in and previous simulation
    is continued. If IPIN=0, no particles from a previous run are used

22. IF IOUTPUTFOREACHRELEASE is set to 1, one output field for each location
    in the RLEASE file is created. For backward calculation this should be
    set to 1. For forward calculation both possibilities are applicable.

23. If IFLUX is set to 1, fluxes of each species through each of the output
    boxes are calculated. Six fluxes, corresponding to northward, southward,
    eastward, westward, upward and downward are calculated for each grid cell of
    the output grid. The control surfaces are placed in the middle of each
    output grid cell. If IFLUX is set to 0, no fluxes are determined.

24. If MDOMAINFILL is set to 1, the first box specified in file RELEASES is used
    as the domain where domain-filling trajectory calculations are to be done.
    Particles are initialized uniformly distributed (according to the air mass
    distribution) in that domain at the beginning of the simulation, and are
    created at the boundaries throughout the simulation period.

25. IND_SOURCE switches between different units for concentrations at the source
    NOTE that in backward simulations the release of computational particles 
    takes place at the "receptor" and the sampling of particles at the "source".
          1=mass units (for bwd-runs = concentration)
          2=mass mixing ratio units 
26. IND_RECEPTOR switches between different units for concentrations at the receptor
          1=mass units (concentrations)
          2=mass mixing ratio units 

27. MQUASILAG indicates whether particles shall be numbered consecutively (1) or
    with their release location number (0). The first option allows tracking of
    individual particles using the partposit output files

28. NESTED_OUTPUT decides whether model output shall be made also for a nested
    output field (normally with higher resolution)

29. PATH TO EMISSION FIELDS, IF EMPTY EMISSION SWITCHED OFF, (include "/" at the end of path name)

30. PATH TO REGION MASKS, (include "/" at the end of path name)

31. PATH TO OH fields, (include "/" at the end of path name) 
