Skip to content
Snippets Groups Projects
Commit e978bf96 authored by Ignacio Pisso's avatar Ignacio Pisso
Browse files

add reference options dir and change default options template

parent cbceb211
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1182 additions and 1 deletion
gen_options.sh
+ 5
− 1
View file @ e978bf96
......@@ -9,9 +9,13 @@ define_options=${1:-$define_options_default}
#template is in 1 or defalt here
options_template_default=~/repos/flexpart/options
#options_template_default=~/repos/flexpart/options
options_template_default=~/repos/flex_gen_input/reference/options
options_template=${2:-$options_template_default}
# to do:
# source settings
#define options contains $options_template
source $define_options
reference/options/AGECLASSES 0 → 100644
+ 15
− 0
View file @ e978bf96
************************************************
* *
*Lagrangian particle dispersion model FLEXPART *
* Please select your options *
* *
*This file determines the ageclasses to be used*
* *
*Ages are given in seconds. The first class *
*starts at age zero and goes up to the first *
*age specified. The last age gives the maximum *
*time a particle is carried in the simulation. *
* *
************************************************
1 Integer Number of age classes
1728000
reference/options/AGECLASSES.oldformat 0 → 100644
+ 15
− 0
View file @ e978bf96
************************************************
* *
*Lagrangian particle dispersion model FLEXPART *
* Please select your options *
* *
*This file determines the ageclasses to be used*
* *
*Ages are given in seconds. The first class *
*starts at age zero and goes up to the first *
*age specified. The last age gives the maximum *
*time a particle is carried in the simulation. *
* *
************************************************
1 Integer Number of age classes
1728000
reference/options/COMMAND 0 → 100644
+ 31
− 0
View file @ e978bf96
&COMMAND
LDIRECT= 1,
IBDATE= 20170102,
IBTIME= 060000,
IEDATE= 20170102,
IETIME= 120000,
LOUTSTEP= 3600,
LOUTAVER= 3600,
LOUTSAMPLE= 900,
ITSPLIT= 99999999,
LSYNCTIME= 900,
CTL= -5.0000000,
IFINE= 4,
IOUT= 1,
IPOUT= 0,
LSUBGRID= 0,
LCONVECTION= 0,
LAGESPECTRA= 0,
IPIN= 0,
IOUTPUTFOREACHRELEASE= 0,
IFLUX= 0,
MDOMAINFILL= 0,
IND_SOURCE= 1,
IND_RECEPTOR= 1,
MQUASILAG= 0,
NESTED_OUTPUT= 0,
LINIT_COND= 0,
SURF_ONLY= 0,
CBLFLAG= 0,
OHFIELDS_PATH= "../../flexin/",
/
reference/options/COMMAND.alternative 0 → 100644
+ 123
− 0
View file @ e978bf96
********************************************************************************
* *
* Input file for the Lagrangian particle dispersion model FLEXPART *
* Please select your options *
* *
********************************************************************************
1 LDIRECT 1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION
20040720 000000 YYYYMMDD HHMISS BEGINNING DATE OF SIMULATION
20040721 120000 YYYYMMDD HHMISS ENDING DATE OF SIMULATION
10800 SSSSS OUTPUT EVERY SSSSS SECONDS
10800 SSSSS TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
900 SSSSS SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
9999999 SSSSSSS TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)
900 SSSSS SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
-5.0 CTL FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL
4 IFINE DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE
3 IOUT 1 CONC. (RESID. TIME FOR BACKWARD RUNS) OUTPUT,2 MIX. RATIO OUTPUT,3 BOTH,4 PLUME TRAJECT.,5=1+4
0 IPOUT PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END
1 LSUBGRID SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO
1 LCONVECTION CONVECTION: 1 YES, 0 NO
0 LAGESPECTRA AGE SPECTRA: 1 YES, 0 NO
0 IPIN CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO
0 IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
0 IFLUX CALCULATE FLUXES: 1 YES, 0 NO
0 MDOMAINFILL DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO
1 IND_SOURCE 1=MASS UNIT , 2=MASS MIXING RATIO UNIT
1 IND_RECEPTOR 1=MASS UNIT , 2=MASS MIXING RATIO UNIT
0 MQUASILAG QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO
0 NESTED_OUTPUT SHALL NESTED OUTPUT BE USED? YES, 0 NO
2 LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
0 SURF_ONLY IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER
1. Simulation direction, 1 for forward, -1 for backward in time
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. Time constant for particle splitting. Particles are split into two
after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on.
8. 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.
9. CTL must be >1 for time steps shorter than the Lagrangian time scale
If CTL<0, a purely random walk simulation is done
10.IFINE=Reduction factor for time step used for vertical wind
11.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.
12.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
13.Switch on/off subgridscale terrain parameterization (increase of
mixing heights due to subgridscale orographic variations)
14.Switch on/off the convection parameterization
15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES
must be available
16. 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
17. 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.
18. 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.
19. 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.
20. 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
21. IND_RECEPTOR switches between different units for concentrations at the receptor
1=mass units (concentrations)
2=mass mixing ratio units
22. 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
23. NESTED_OUTPUT decides whether model output shall be made also for a nested
output field (normally with higher resolution)
24. LINIT_COND determines whether, for backward runs only, the sensitivity to initial
conditions shall be calculated and written to output files
0=no output, 1 or 2 determines in which units the initial conditions are provided.
25. SURF_ONLY: When set to 1, concentration/emission sensitivity is written out only
for the surface layer; useful for instance when only footprint emission sensitivity is needed
but initial conditions are needed on a full 3-D grid
reference/options/COMMAND.namelist 0 → 100644
+ 31
− 0
View file @ e978bf96
&COMMAND
LDIRECT= 1,
IBDATE= 20140911,
IBTIME= 050000,
IEDATE= 20140911,
IETIME= 150000,
LOUTSTEP= 3600,
LOUTAVER= 3600,
LOUTSAMPLE= 900,
ITSPLIT= 99997200,
LSYNCTIME= 900,
CTL= -5.0000000 ,
IFINE= 4,
IOUT= 1,
IPOUT= 0,
LSUBGRID= 1,
LCONVECTION= 1,
LAGESPECTRA= 0,
IPIN= 0,
IOUTPUTFOREACHRELEASE= 0,
IFLUX= 0,
MDOMAINFILL= 0,
IND_SOURCE= 1,
IND_RECEPTOR= 1,
MQUASILAG= 0,
NESTED_OUTPUT= 1,
LINIT_COND= 0,
SURF_ONLY= 0,
CBLFLAG= 0,
OHFIELDS_PATH= "../../flexin/",
/
reference/options/COMMAND.oldformat 0 → 100644
+ 199
− 0
View file @ e978bf96
********************************************************************************
* *
* 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
20140911 050000
YYYYMMDD HHMISS BEGINNING DATE OF SIMULATION
3. ________ ______ 3X, I8, 1X, I6
20140911 150000
YYYYMMDD HHMISS ENDING DATE OF SIMULATION
4. _____ 3X, I5
3600
SSSSS OUTPUT EVERY SSSSS SECONDS
5. _____ 3X, I5
3600
SSSSS TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
6. _____ 3X, I5
900
SSSSS SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
7. _________ 3X, I9
999999999
SSSSSSSSS TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)
8. _____ 3X, I5
900
SSSSS SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
9. ---.-- 4X, F6.4
-5.0
CTL FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL
10. --- 4X, I3
4
IFINE DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE
11. - 4X, I1
1
IOUT 1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4
12. - 4X, I1
0
IPOUT PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END
13. _ 4X, I1
1
LSUBGRID SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO
14. _ 4X, I1
1
LCONVECTION CONVECTION: 1 YES, 0 NO
15. _ 4X, I1
0
LAGESPECTRA AGE SPECTRA: 1 YES, 0 NO
16. _ 4X, I1
0
IPIN CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO
17. _
0 4X,I1
IOFR IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
18. _ 4X, I1
0
IFLUX CALCULATE FLUXES: 1 YES, 0 NO
19. _ 4X, I1
0
MDOMAINFILL DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT. O3 TRACER
20. _ 4X, I1
1
IND_SOURCE 1=MASS UNIT , 2=MASS MIXING RATIO UNIT
21. _ 4X, I1
1
IND_RECEPTOR 1=MASS UNIT , 2=MASS MIXING RATIO UNIT
22. _ 4X, I1
0
MQUASILAG QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO
23. _ 4X, I1
1
NESTED_OUTPUT SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO
24. _ 4X, I1
0
LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
25. _ 4X, I1
0
SURF_ONLY IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER
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. Time constant for particle splitting. Particles are split into two
after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on.
8. 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.
9. CTL must be >1 for time steps shorter than the Lagrangian time scale
If CTL<0, a purely random walk simulation is done
10.IFINE=Reduction factor for time step used for vertical wind
11.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.
12.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
13.Switch on/off subgridscale terrain parameterization (increase of
mixing heights due to subgridscale orographic variations)
14.Switch on/off the convection parameterization
15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES
must be available
16. 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
17. 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.
18. 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.
19. 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.
20. 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
21. IND_RECEPTOR switches between different units for concentrations at the receptor
1=mass units (concentrations)
2=mass mixing ratio units
22. 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
23. NESTED_OUTPUT decides whether model output shall be made also for a nested
output field (normally with higher resolution)
24. LINIT_COND determines whether, for backward runs only, the sensitivity to initial
conditions shall be calculated and written to output files
0=no output, 1 or 2 determines in which units the initial conditions are provided.
25. SURF_ONLY: When set to 1, concentration/emission sensitivity is written out only
for the surface layer; useful for instance when only footprint emission sensitivity is needed
but initial conditions are needed on a full 3-D grid
reference/options/COMMAND.reference 0 → 100644
+ 199
− 0
View file @ e978bf96
********************************************************************************
* *
* 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
20110310 000000
YYYYMMDD HHMISS BEGINNING DATE OF SIMULATION
3. ________ ______ 3X, I8, 1X, I6
20110310 120000
YYYYMMDD HHMISS ENDING DATE OF SIMULATION
4. _____ 3X, I5
10800
SSSSS OUTPUT EVERY SSSSS SECONDS
5. _____ 3X, I5
10800
SSSSS TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
6. _____ 3X, I5
900
SSSSS SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
7. _________ 3X, I9
999999999
SSSSSSSSS TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)
8. _____ 3X, I5
900
SSSSS SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
9. ---.-- 4X, F6.4
-5.0
CTL FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL
10. --- 4X, I3
4
IFINE DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE
11. - 4X, I1
3
IOUT 1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4
12. - 4X, I1
0
IPOUT PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END
13. _ 4X, I1
1
LSUBGRID SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO
14. _ 4X, I1
1
LCONVECTION CONVECTION: 1 YES, 0 NO
15. _ 4X, I1
0
LAGESPECTRA AGE SPECTRA: 1 YES, 0 NO
16. _ 4X, I1
0
IPIN CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO
17. _
0 4X,I1
IOFR IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
18. _ 4X, I1
0
IFLUX CALCULATE FLUXES: 1 YES, 0 NO
19. _ 4X, I1
0
MDOMAINFILL DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT. O3 TRACER
20. _ 4X, I1
1
IND_SOURCE 1=MASS UNIT , 2=MASS MIXING RATIO UNIT
21. _ 4X, I1
1
IND_RECEPTOR 1=MASS UNIT , 2=MASS MIXING RATIO UNIT
22. _ 4X, I1
0
MQUASILAG QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO
23. _ 4X, I1
0
NESTED_OUTPUT SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO
24. _ 4X, I1
2
LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
25. _ 4X, I1
0
SURF_ONLY IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER
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. Time constant for particle splitting. Particles are split into two
after SSSSS seconds, 2xSSSSS seconds, 4xSSSSS seconds, and so on.
8. 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.
9. CTL must be >1 for time steps shorter than the Lagrangian time scale
If CTL<0, a purely random walk simulation is done
10.IFINE=Reduction factor for time step used for vertical wind
11.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.
12.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
13.Switch on/off subgridscale terrain parameterization (increase of
mixing heights due to subgridscale orographic variations)
14.Switch on/off the convection parameterization
15.Switch on/off the calculation of age spectra: if yes, the file AGECLASSES
must be available
16. 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
17. 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.
18. 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.
19. 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.
20. 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
21. IND_RECEPTOR switches between different units for concentrations at the receptor
1=mass units (concentrations)
2=mass mixing ratio units
22. 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
23. NESTED_OUTPUT decides whether model output shall be made also for a nested
output field (normally with higher resolution)
24. LINIT_COND determines whether, for backward runs only, the sensitivity to initial
conditions shall be calculated and written to output files
0=no output, 1 or 2 determines in which units the initial conditions are provided.
25. SURF_ONLY: When set to 1, concentration/emission sensitivity is written out only
for the surface layer; useful for instance when only footprint emission sensitivity is needed
but initial conditions are needed on a full 3-D grid
reference/options/IGBP_int1.dat 0 → 100644
+ 0
− 0
View file @ e978bf96
File added
reference/options/OUTGRID 0 → 100644
+ 156
− 0
View file @ e978bf96
********************************************************************************
* *
* Input file for the Lagrangian particle dispersion model FLEXPART *
* Please specify your output grid *
* *
********************************************************************************
1. ------.---- 4X,F11.4
-179.0000 GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
OUTLONLEFT (left boundary of the first grid cell - not its centre)
2. ------.---- 4X,F11.4
-90.0000 GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
OUTLATLOWER (lower boundary of the first grid cell - not its centre)
3. ----- 4X,I5
720 NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
NUMXGRID
4. ----- 4X,I5
360 NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
NUMYGRID
5. ------.--- 4X,F10.3
0.50 GRID DISTANCE IN X DIRECTION
DXOUTLON
6. ------.--- 4X,F10.3
0.50 GRID DISTANCE IN Y DIRECTION
DYOUTLAT
10. -----.- 4X, F7.1
1000.0
LEVEL 1 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
2000.0
LEVEL 2 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
3000.0
LEVEL 3 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
4000.0
LEVEL 4 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
5000.0
LEVEL 5 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
6000.0
LEVEL 6 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
7000.0
LEVEL 7 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
8000.0
LEVEL 8 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
9000.0
LEVEL 9 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
10000.0
LEVEL10 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
11000.0
LEVEL11 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
12000.0
LEVEL12 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
13000.0
LEVEL13 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
14000.0
LEVEL14 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
15000.0
LEVEL15 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
16000.0
LEVEL16 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
17000.0
LEVEL17 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
18000.0
LEVEL18 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
19000.0
LEVEL19 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
20000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
21000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
22000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
23000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
24000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
25000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
26000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
27000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
28000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
29000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
30000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
40000.0
LEVEL20 HEIGHT OF LEVEL (UPPER BOUNDARY)
reference/options/OUTGRID_NEST 0 → 100644
+ 30
− 0
View file @ e978bf96
********************************************************************************
* *
* Input file for the Lagrangian particle dispersion model FLEXPART *
* Please specify your output grid *
* *
********************************************************************************
1. ------.---- 4X,F11.4
-74.0000 GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
OUTLONLEFT (left boundary of the first grid cell - not its centre)
2. ------.---- 4X,F11.4
37.0000 GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
OUTLATLOWER (lower boundary of the first grid cell - not its centre)
3. ----- 4X,I5
1000 NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
NUMXGRID
4. ----- 4X,I5
200 NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
NUMYGRID
5. ------.----- 4X,F12.5
0.01000 GRID DISTANCE IN X DIRECTION
DXOUTLON
6. ------.----- 4X,F12.5
0.01000 GRID DISTANCE IN Y DIRECTION
DYOUTLAT
reference/options/RECEPTORS 0 → 100644
+ 10
− 0
View file @ e978bf96
&RECEPTORS
RECEPTOR="receptor 1 ",
LON= 0.0000000 ,
LAT= 0.0000000 ,
/
&RECEPTORS
RECEPTOR="receptor 2 ",
LON= 1.00000000 ,
LAT= 0.00000000 ,
/
reference/options/RECEPTORS.namelist 0 → 100644
+ 30
− 0
View file @ e978bf96
&RECEPTORS
RECEPTOR="F15 ",
LON= 21.1333008 ,
LAT= 38.0833015 ,
/
&RECEPTORS
RECEPTOR="NL01 ",
LON= 5.78329992 ,
LAT= 50.9166985 ,
/
&RECEPTORS
RECEPTOR="B05 ",
LON= 4.34999990 ,
LAT= 50.7999992 ,
/
&RECEPTORS
RECEPTOR="D27 ",
LON= 11.1332998 ,
LAT= 52.9667015 ,
/
&RECEPTORS
RECEPTOR="D08 ",
LON= 8.69999981 ,
LAT= 53.8666992 ,
/
&RECEPTORS
RECEPTOR="DK05 ",
LON= 8.13329983 ,
LAT= 56.0000000 ,
/
reference/options/RECEPTORS.oldformat 0 → 100644
+ 79
− 0
View file @ e978bf96
********************************************************************************
* *
* Input file for the Lagrangian particle dispersion model FLEXPART *
* Please specify your receptor points *
* For the receptor points, ground level concentrations are calculated *
* *
********************************************************************************
1. ---------------- 4X,A16
F15 NAME OF RECEPTOR POINT
RECEPTORNAME
2. ------.---- 4X,F11.4
21.1333 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
38.0833 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
NL01 NAME OF RECEPTOR POINT
RECEPTORNAME
2. ------.---- 4X,F11.4
5.7833 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
50.9167 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
B05 NAME OF RECEPTOR POINT
RECEPTORNAME
2. ------.---- 4X,F11.4
4.3500 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
50.8000 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
D27 NAME OF RECEPTOR POINT
RECEPTORNAME
2. ------.---- 4X,F11.4
11.1333 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
52.9667 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
D08 NAME OF RECEPTOR POINT
RECEPTORNAME
2. ------.---- 4X,F11.4
8.7000 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
53.8667 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
DK05 NAME OF RECEPTOR POINT
RECEPTORNAME
2. ------.---- 4X,F11.4
8.1333 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
56.0000 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
reference/options/RELEASES 0 → 100644
+ 20
− 0
View file @ e978bf96
&RELEASES_CTRL
NSPEC= 1,
SPECNUM_REL= 24,
/
&RELEASE
IDATE1= 20170102,
ITIME1= 090000,
IDATE2= 20170102,
ITIME2= 090000,
LON1= 0.000 ,
LON2= 0.000 ,
LAT1= 0.000 ,
LAT2= 0.000 ,
Z1= 50.000 ,
Z2= 50.000 ,
ZKIND= 1,
MASS= 1.0000E8 ,
PARTS= 10000
COMMENT="TEST1",
/
reference/options/RELEASES.alternative 0 → 100644
+ 40
− 0
View file @ e978bf96
*************************************************************************
* *
* *
* *
* Input file for the Lagrangian particle dispersion model FLEXPART *
* Please select your options *
* *
* *
* *
*************************************************************************
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1 Total number of species emitted
24 Index of species in file SPECIES
=========================================================================
20011028 150007
20011028 150046
9.40480
48.5060
9.50670
48.5158
2
6933.60
6950.40
20000
1.00000
FLIGHT_11242
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
20011028 150047
20011028 150107
9.30380
48.5158
9.40480
48.5906
2
6833.50
6950.40
20000
1.00000
FLIGHT_11185
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
reference/options/RELEASES.reference 0 → 100644
+ 90
− 0
View file @ e978bf96
*************************************************************************
* *
* *
* *
* Input file for the Lagrangian particle dispersion model FLEXPART *
* Please select your options *
* *
* *
* *
*************************************************************************
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1
___ i3 Total number of species emitted
24
___ i3 Index of species in file SPECIES
=========================================================================
20040720 0
________ ______ i8,1x,i6 Beginning date and time of release
20040720 120000
________ ______ i8,1x,i6 Ending date and time of release
9.4048
____.____ f9.4 Longitude [DEG] of lower left corner
48.5060
____.____ f9.4 Latitude [DEG] of lower left corner
9.5067
____.____ f9.4 Longitude [DEG] of upper right corner
48.5158
____.____ f9.4 Latitude [DEG] of upper right corner
2
_________ i9 1 for m above ground, 2 for m above sea level
6933.60
_____.___ f10.3 Lower z-level (in m agl or m asl)
6950.40
_____.___ f10.3 Upper z-level (in m agl or m asl)
20000
_________ i9 Total number of particles to be released
1.0000E00
_.____E__ e9.4 Total mass emitted
RELEASE_TEST1
________________________________________ character*40 comment
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
20040720 0
________ ______ i8,1x,i6 Beginning date and time of release
20040720 120000
________ ______ i8,1x,i6 Ending date and time of release
9.3038
____.____ f9.4 Longitude [DEG] of lower left corner
48.5158
____.____ f9.4 Latitude [DEG] of lower left corner
9.4048
____.____ f9.4 Longitude [DEG] of upper right corner
48.5906
____.____ f9.4 Latitude [DEG] of upper right corner
2
_________ i9 1 for m above ground, 2 for m above sea level
6833.50
_____.___ f10.3 Lower z-level (in m agl or m asl)
6950.40
_____.___ f10.3 Upper z-level (in m agl or m asl)
20000
_________ i9 Total number of particles to be released
1.0000E00
_.____E__ e9.4 Total mass emitted
RELEASE_TEST2
________________________________________ character*40 comment
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
reference/options/SPECIES/SPECIES.README 0 → 100644
+ 67
− 0
View file @ e978bf96
Specifies all physico-chemical properties for the species.
To generate an overview table compile the specoverview.f90 program
gfortran specoverview.f90 -o specoverview
./specoverview
UPDATE FOR FLEXPART VERSION 10.0 DEV:
---------------------------------------------------------------------------------------------------------------------------
Please note that the format of the SPECIES files has changed from the previous version (9.2) of FLEXPART and it is important
to update all SPECIES files accordingly. The use of SPECIES files from older FLEXPART versions may lead to run-time
errors or erroneous results.
---------------------------------------------------------------------------------------------------------------------------
Half life due to radioactive or chemical decay and can be switched off by specifying a negative value.
Molecular weight (molweight of the species, needed for mixing ratio output)
---------------------------------------------------------------------------------------------------------------------------
WET DEPOSITION
Below-cloud scavenging
- Gases:
A and B are the parameters for below-cloud scavenging defined by Eq. scavenging coefficient = A x precipitation-rate^B
If you migrate from FLEXPART version 9.2 and below this is equal to A and B as it was.
- Particles:
A "rain collection efficiency"
B "snow collection efficiency" are both defined by Equations 46-47.
Suggested values for particle scavenging are A=1.0,B=1.0.
Using values A=0.2 and B=0.2 will give results close to the old scavenging scheme of FLEXPART version 9.2 and below.
In-cloud scavenging
- Gases:
Is switched on if Ai or Bi are given greater than 0
If you migrate from FLEXPART version 9.2 and below this is equal to Ai and Bi should be set to 1
- Particles:
Ai ("Cloud Condensation Nuclei (CCN) efficiency")
Bi ("Ice Nuclei (IN) efficiency") are both from Eq. 51.
Suggested values for particle scavenging are Ai=0.9 and Bi=0.1.
Using values Ai=0.9 and Bi=0.9 should give results cto the old scavenging scheme of FLEXPART version 9.2 and below.
---------------------------------------------------------------------------------------------------------------------------
DRY DEPOSITION
- Gases:
D is the diffusivity of water vapor and
D_i is the diffusivity of the species,
H is the effective Henry’s constant,
f0 varies between 0 and 1 and gives the reactivity of a species relative to that of ozone. For non-reactive species f0 is 0, for slightly reactive it is 0.1 and for highly reactive it is 1.
Dry deposition of gases is switched off by negative D.
- Particles:
rho specifies the density of the substance,
dquer its mean diameter
dsig the measure of variation.
If no detailed information for deposition velocity calculation is available, a constant deposition velocity vd [cm s-1] can be used.
Dry deposition of particles is switched off by negative rho.
---------------------------------------------------------------------------------------------------------------------------
OH REACTION:
can be turned on by giving parameter C [cm^3/molecule/sec]], D [K] and N [no unit] positive values, defined by Eq. 39.
OH reaction is switched off by negative C.
---------------------------------------------------------------------------------------------------------------------------
Optionally emission variation information
can be added at the end of the file, defined as following: Since FLEXPART version 6.0, emission factors can be defined that change the temporal variation of particle releases. This is useful, for instance,
to simulate the typical daily and weekly cycle of anthropogenic emissions. The emission factors are given in the file of the
corresponding species SPECIES_nnn, where nnn is the species number defined in file RELEASES. If no emission variation
information is given, emission rates for species nnn are taken as constant. Release rates can vary with the hour of the day
and with the day of the week, according to the local time at the release location. Emission factors must be 1 on average. 24
hourly as well as 7 daily values must be specified. Furthermore, different disaggregation factors must be given for area sources
and for point sources. FLEXPART distinguishes between the two using the lower altitude of the release box: area sources are
assumed to start below 0.5 m above the ground, whereas point sources are assumed to be higher. Please note that when this
option is used, it is not so easy to determine the maximum number of particles present at a particular time of the model run. It
might then be necessary to increase the parameter maxpart to a higher value than what would otherwise be needed.
reference/options/SPECIES/SPECIES_002 0 → 100644
+ 21
− 0
View file @ e978bf96
&SPECIES_PARAMS
PSPECIES="Ozone", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=-0.9E-9, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=-9.9, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=-9.9, ! Dry deposition (particles) - dquer
PDSIGMA=-9.9, ! Dry deposition (particles) - dsig
PDRYVEL=-9.9, ! Alternative: dry deposition velocity
PRELDIFF=1.6, ! Dry deposition (gases) - D
PHENRY=1.0E-2, ! Dry deposition (gases) - Henrys const.
PF0=1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=48.0, ! molweight
POHCCONST=-0.9E-9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9, ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
reference/options/SPECIES/SPECIES_003 0 → 100644
+ 21
− 0
View file @ e978bf96
&SPECIES_PARAMS
PSPECIES="NO", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=8.0E-6, ! Below-cloud scavenging (gases) - A (weta_gas)
PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas)
PCRAIN_AERO=-9.9, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=-9.9, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=-9.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=-9.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=1.2, ! Dry deposition (gases) - D
PHENRY=3.0E-03, ! Dry deposition (gases) - Henrys const.
PF0=0.0, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=30.0, ! molweight
POHCCONST=-0.9E+9, ! OH Reaction rate - C [cm^3/molecule/sec]
POHDCONST=-9.9 ! OH Reaction rate - D [K]
POHNCONST=2.0, ! OH Reaction rate - N (no unit)
/
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment