...
 
Commits (278)
FP_ecmwf_gfortran*
*.o
*.mod
*_mod.mod
class_gribfile.mod
.DS_Store
output
src/FLEXPART_*
fort.*
CONVERT2
FLEXPART
FLEXPART_MPI
XXXXXX EMPTY LINES XXXXXXXXX
XXXXXX EMPTY LINES XXXXXXXX
YYYYMMDD HHMMSS name of the file(up to 80 characters)
20120101 000000 EA12010100 ON DISK
20120101 010000 EA12010101 ON DISK
20120101 020000 EA12010102 ON DISK
20120101 030000 EA12010103 ON DISK
20120101 040000 EA12010104 ON DISK
20120101 050000 EA12010105 ON DISK
20120101 060000 EA12010106 ON DISK
20120101 070000 EA12010107 ON DISK
20120101 080000 EA12010108 ON DISK
20120101 090000 EA12010109 ON DISK
20120101 100000 EA12010110 ON DISK
20120101 110000 EA12010111 ON DISK
20120101 120000 EA12010112 ON DISK
20120101 130000 EA12010113 ON DISK
20120101 140000 EA12010114 ON DISK
20120101 150000 EA12010115 ON DISK
20120101 160000 EA12010116 ON DISK
20120101 170000 EA12010117 ON DISK
20120101 180000 EA12010118 ON DISK
20120101 190000 EA12010119 ON DISK
20120101 200000 EA12010120 ON DISK
20120101 210000 EA12010121 ON DISK
20120101 220000 EA12010122 ON DISK
20120101 230000 EA12010123 ON DISK
# README #
This README would normally document whatever steps are necessary to get your application up and running.
### What is this repository for? ###
* This repository contains the development version of the Lagrangian model FLEXPART 9.2 migrated from svn
### How do I get set up? ###
* Summary of set up
* Configuration
* Dependencies
* Database configuration
* How to run tests
* Deployment instructions
### Contribution guidelines ###
* The version contributed should compile on a reference version of the system and compiler. The current reference is gfortran on Ubuntu 12.04 -- GNU Fortran (Ubuntu/Linaro 4.6.3-1ubuntu5) 4.6.3
* Code contribution including new features and bug fixes should be complemented with appropriate tests
* Writing tests:
An essential test consists of a set of input files and directories that allow FLEXPART to run.
A test can be accompanied by output files for verification
* Code review
* Other guidelines
### Who do I talk to? ###
* Ignacio Pisso ip@nilu.no
# README #
### What is this repository for? ###
* This repository contains the development version of the Lagrangian model FLEXPART
### How do I get set up? ###
* Configuration
Edit the makefile with the paths to libraries and include files
* Dependencies
* Jasper and grib_api or ECCodes
* NetCDF (optional)
* Compilation
```
> cd src
change the appropriate paths to libraries in the makefile
> make
```
* Deployment instructions
FLEXPART is a standalone executable
The necessary ECMWF wind fields can be obtained testing flex_ecmwf
The AVAILABLE file works with the default ERA 5 retrieved winds
In the winds are available in flex_ecmwf/work it should suffice to execute
./src/FLEXPART in the main directory
### Contribution guidelines ###
* The version contributed should compile on a reference version of the system and compiler. The current reference is gfortran 5.4 on Ubuntu 16.04
* Code contribution including new features and bug fixes should be complemented with appropriate tests
An essential test consists of a set of input files and directories that allow FLEXPART to run.
A test can be accompanied by output files for verification
* Code review
[comment]: # "### Who do I talk to? ###"
ignacio.pisso@nilu.no
FLEXPART VERSION 10.0 beta (MPI)
Description
-----------
This branch contains both the standard (serial) FLEXPART, and a parallel
version (implemented with MPI). The latter is under developement, so not
every FLEXPART option is implemented yet.
MPI related subroutines and variables are in file mpi_mod.f90.
Most of the source files are identical/shared between the serial and
parallel versions. Those that depend on the MPI module have '_mpi'
apppended to their names, e.g. 'timemanager_mpi.f90'
Installation
------------
A MPI library must be installed on the target platform, either as a
system library or compiled from source.
So far, we have tested the following freely available implementations:
mpich2 -- versions 3.0.1, 3.0.4, 3.1, 3.1.3
OpenMPI -- version 1.8.3
Based on testing so far, OpenMPI is recommended.
Compiling the parallel version (executable: FP_ecmwf_MPI) is done by
'make [-j] ecmwf-mpi'
The makefile has resolved dependencies, so 'make -j' will compile
and link in parallel.
The included makefile must be edited to match the target platform
(location of system libraries, compiler etc.).
Usage
-----
Running the parallel version with MPI is done with the "mpirun" command
(some MPI implementations may use a "mpiexec" command instead). The
simplest case is:
'mpirun -n [number] ./FP_ecmwf_MPI'
where 'number' is the number of processes to launch. Depending on the
target platform, useful options regarding process-to-processor bindings
can be specified (for performance reasons), e.g,
'mpirun --bind-to l3cache -n [number] ./FP_ecmwf_MPI'
Implementation
--------------
The current parallel model is based on distributing particles equally
among the running processes. In the code, variables like 'maxpart' and
'numpart' are complemented by variables 'maxpart_mpi' and 'numpart_mpi'
which are the run-time determined number of particles per process, i.e,
maxpart_mpi = maxpart/np, where np are the number of processes. The variable 'numpart'
is still used in the code, but redefined to mean 'number of particles
per MPI process'
The root MPI process writes concentrations to file, following a MPI
communication step where each process sends its contributions to root,
where the individual contributions are summed.
In the parallel version one can choose to set aside a process dedicated
to reading and distributing meteorological data ("windfields"). This process will
thus not participate in the calculation of trajectories. This might not be
the optimal choice when running with very few processes.
As an example, running with a total number of processes np=4 and
using one of these processes for reading windfields will normally
be faster than running with np=3 and no dedicated 'reader' process.
But it is also possible that the
program will run even faster if the 4th process is participating in
the calculation of particle trajectories instead. This will largely depend on
the problem size (total number of particles in the simulation, resolution
of grids etc) and hardware being used (disk speed/buffering, memory
bandwidth etc).
To control this
behavior, edit the parameter 'read_grp_min' in file mpi_mod.f90. This
sets the minimum number of total processes at which one will be set
aside for reading the fields. Experimentation is required to find
the optimum value. At typical NILU machines (austre.nilu.no,
dmz-proc01.nilu.no) with 24-32 cores, a value of 6-8 seems to be a
good choice.
An experimental feature, which is an extension of the functionality
described above, is to hold 3 fields in memory instead of the usual 2.
Here, the transfer of fields from the "reader" process to the "particle"
processes is done on the vacant field index, simultaneously while the
"particle" processes are calculating trajectories. To use this feature,
set 'lmp_sync=.false'. in file mpi_mod.f90 and set numwfmem=3 in file
par_mod.f90. At the moment, this method does not seem to produce faster
running code (about the same as the "2-fields" version).
Performance efficency considerations
------------------------------------
A couple of reference runs have been set up to measure performace of the
MPI version (as well as checking for errors in the implementation).
They are as follows:
Reference run 1 (REF1):
* Forward modelling (24h) of I2-131, variable number of particles
* Two release locations
* 360x720 Global grid, no nested grid
* Species file modified to include (not realistic) values for
scavenging/deposition
As the parallization is based on particles, it follows that if
FLEXPART-MPI is run with no (or just a few) particles, no performance
improvement is possible. In this case, most processing time is spent
in the 'getfields'-routine.
A) Running without dedicated reader process
----------------------------------------
Running REF1 with 100M particles on 16 processes (NILU machine 'dmz-proc04'),
a speedup close to 8 is observed (~50% efficiency).
Running REF1 with 10M particles on 8 processes (NILU machine 'dmz-proc04'),
a speedup close to 3 is observed (~40% efficiency). Running with 16
processes gives only marginal improvements (speedup ~3.5) because of the 'getfields'
bottleneck.
Running REF1 with 1M particles: Here 'getfields' consumes ~70% of the CPU
time. Running with 4 processes gives a speedup of ~1.5. Running with more
processes does not help much here.
B) Running with dedicated reader process
----------------------------------------
Running REF1 with 40M particles on 16 processes (NILU machine 'dmz-proc04'),
a speedup above 10 is observed (~63% efficiency).
:TODO: more to come...
Advice
------
From the tests referred to above, the following advice can be given:
* Do not run with too many processes.
* Do not use the parallel version when running with very few particles.
What is implemented in the MPI version
--------------------------------------
-The following should work (have been through initial testing):
* Forward runs
* OH fields
* Radioactive decay
* Particle splitting
* Dumping particle positions to file
* ECMWF data
* Wet/dry deposition
* Nested grid output
* NetCDF output
* Namelist input/output
* Domain-filling trajectory calculations
* Nested wind fields
-Implemented but untested:
* Backward runs (but not initial_cond_output.f90)
-The following will most probably not work (untested/under developement):
* Calculation/output of fluxes
-This will positively NOT work yet
* Subroutine partoutput_short (MQUASILAG = 1) will not dump particles
correctly at the moment
* Reading particle positions from file (the tools to implement this
are available in mpi_mod.f90 so it will be possible soon)
Please keep in mind that running the serial version (FP_ecmwf_gfortran)
should yield identical results as running the parallel version
(FP_ecmwf_MPI) using only one process, i.e. "mpirun -n 1 FP_ecmwf_MPI".
If not, this indicates a bug.
When running with multiple processes, statistical differences are expected
in the results.
Contact
-------
If you have questions, or wish to work with the parallel version, please
contact Espen Sollum (eso@nilu.no). Please report any errors/anomalies!
This diff is collapsed.
......@@ -11,5 +11,5 @@
*time a particle is carried in the simulation. *
* *
************************************************
1 Integer Number of age classes
1728000
1 Integer Number of age classes
1728000
......@@ -11,5 +11,5 @@
*time a particle is carried in the simulation. *
* *
************************************************
1 Integer Number of age classes
1728000
1 Integer Number of age classes
1728000
&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/",
/
......@@ -196,3 +196,4 @@
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
********************************************************************************
* *
* 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
360 NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
NUMXGRID
4. ----- 4X,I5
180 NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
NUMYGRID
5. ------.--- 4X,F10.3
1.000 GRID DISTANCE IN X DIRECTION
DXOUTLON
6. ------.--- 4X,F10.3
1.000 GRID DISTANCE IN Y DIRECTION
DYOUTLAT
10. -----.- 4X, F7.1
100.0
LEVEL 4 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
3000.0
LEVEL 4 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
50000.0
LEVEL 4 HEIGHT OF LEVEL (UPPER BOUNDARY)
********************************************************************************
* *
* 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
100.0
LEVEL 1 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
500.0
LEVEL 2 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
1000.0
LEVEL 3 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
10000.0
LEVEL 4 HEIGHT OF LEVEL (UPPER BOUNDARY)
10. -----.- 4X, F7.1
40000.0
LEVEL 5 HEIGHT OF LEVEL (UPPER BOUNDARY)
......@@ -6,25 +6,25 @@
********************************************************************************
1. ------.---- 4X,F11.4
-15.0000 GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
-2.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
35.0000 GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
1.0000 GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
OUTLATLOWER (lower boundary of the first grid cell - not its centre)
3. ----- 4X,I5
220 NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
400 NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
NUMXGRID
4. ----- 4X,I5
120 NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
200 NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
NUMYGRID
5. ------.----- 4X,F12.5
0.25000 GRID DISTANCE IN X DIRECTION
0.01000 GRID DISTANCE IN X DIRECTION
DXOUTLON
6. ------.----- 4X,F12.5
0.25000 GRID DISTANCE IN Y DIRECTION
0.01000 GRID DISTANCE IN Y DIRECTION
DYOUTLAT
......@@ -10,11 +10,11 @@
RECEPTORNAME
2. ------.---- 4X,F11.4
6.1333 GEOGRAFICAL LONGITUDE
21.1333 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
49.0833 GEOGRAFICAL LATITUDE
38.0833 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
......
&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 ,
/
......@@ -10,11 +10,11 @@
RECEPTORNAME
2. ------.---- 4X,F11.4
6.1333 GEOGRAFICAL LONGITUDE
21.1333 GEOGRAFICAL LONGITUDE
XRECEPTOR
3. ------.---- 4X,F11.4
49.0833 GEOGRAFICAL LATITUDE
38.0833 GEOGRAFICAL LATITUDE
YRECEPTOR
================================================================================
1. ---------------- 4X,A16
......
&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",
/
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.3:
---------------------------------------------------------------------------------------------------------------------------
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.
&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)
/
&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)
/
&SPECIES_PARAMS
PSPECIES="NO2", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=1.0E-05, ! 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.9E9, ! 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.6, ! Dry deposition (gases) - D
PHENRY=1.0E-02, ! Dry deposition (gases) - Henrys const.
PF0=0.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=46.0, ! molweight
POHCCONST=-0.9E9, ! 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)
/
&SPECIES_PARAMS
PSPECIES="HNO3", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=5.0E-05, ! 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.9, ! Dry deposition (gases) - D
PHENRY=1E14, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=63.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)
/
&SPECIES_PARAMS
PSPECIES="HNO2", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=1.6, ! Dry deposition (gases) - D
PHENRY=100000.0, ! Dry deposition (gases) - Henrys const.
PF0=.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=47.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)
/
&SPECIES_PARAMS
PSPECIES="H2O2", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=1.0E-04, ! 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.4, ! Dry deposition (gases) - D
PHENRY=100000.0, ! Dry deposition (gases) - Henrys const.
PF0=1.0, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=34.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)
/
&SPECIES_PARAMS
PSPECIES="HCHO", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=1,3, ! Dry deposition (gases) - D
PHENRY=6000.0, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! 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)
/
&SPECIES_PARAMS
PSPECIES="PAN", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=2.6, ! Dry deposition (gases) - D
PHENRY=3.6, ! Dry deposition (gases) - Henrys const.
PF0=0.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=121.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)
/
&SPECIES_PARAMS
PSPECIES="NH3", ! Tracer name
PDECAY=-9.9, ! Species half life
PWETA_GAS=9.9E-05, ! 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.1, ! Dry deposition (gases) - D
PHENRY=2.0E14, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=17.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)
/
&SPECIES_PARAMS
PSPECIES="SO4-aero", ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=2000.0, ! Dry deposition (particles) - rho
PDQUER=4.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E+9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=0.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)
/
&SPECIES_PARAMS
PSPECIES="NO3-aero", ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=-0.9E+9, ! Dry deposition (particles) - rho
PDQUER=4.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E+9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=2000.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)
/
&SPECIES_PARAMS
PSPECIES="I2-131", ! Tracer name
PDECAY=691200.0, ! Species half life
PWETA_GAS=8.0E-05, ! 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=2.7, ! Dry deposition (gases) - D
PHENRY=100000.0, ! Dry deposition (gases) - Henrys const.
PF0=0.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! 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)
/
&SPECIES_PARAMS
PSPECIES=""I-131", ! Tracer name
PDECAY=691200.0, ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=2500.0, ! Dry deposition (particles) - rho
PDQUER=6.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-9.9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=350.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)
/
&SPECIES_PARAMS
PSPECIES="Cs-137", ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.9, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=1900.0, ! Dry deposition (particles) - rho
PDQUER=0.4E-06, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! 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)
/
&SPECIES_PARAMS
PSPECIES="Y-91", ! Tracer name
PDECAY=5037120.0, ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=2500.0, ! Dry deposition (particles) - rho
PDQUER=6.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=350.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)
/
&SPECIES_PARAMS
PSPECIES=""Ru-106", ! Tracer name
PDECAY=31536000.0, ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=2500.0, ! Dry deposition (particles) - rho
PDQUER=6.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=350.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)
/
&SPECIES_PARAMS
PSPECIES="KR-85", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! 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)
/
&SPECIES_PARAMS
PSPECIES="Sr-90", ! 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=1.0, ! Below-cloud scavenging (particles) - Crain (crain_aero)
PCSNOW_AERO=1.0, ! Below-cloud scavenging (particles) - Csnow (csnow_aero)
PCCN_AERO=0.9, ! In-cloud scavenging (particles) - CCNeff (ccn_aero)
PIN_AERO=0.1, ! In-cloud scavenging (particles) - INeff (in_aero)
PDENSITY=2500.0, ! Dry deposition (particles) - rho
PDQUER=6.0E-07, ! Dry deposition (particles) - dquer
PDSIGMA=0.3, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! 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)
/
&SPECIES_PARAMS
PSPECIES="Xe-133", ! Tracer name
PDECAY=453168.0, ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! 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)
/
&SPECIES_PARAMS
PSPECIES="CO", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=28.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)
/
&SPECIES_PARAMS
PSPECIES="SO2", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=2.0, ! Dry deposition (gases) - D
PHENRY=100000.0, ! Dry deposition (gases) - Henrys const.
PF0=0, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=64.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)
/
&SPECIES_PARAMS
PSPECIES="AIRTRACER", ! 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=0.0, ! Dry deposition (particles) - dquer
PDSIGMA=0.0, ! Dry deposition (particles) - dsig
PDRYVEL=-9.99, ! Alternative: dry deposition velocity
PRELDIFF=-9.9, ! Dry deposition (gases) - D
PHENRY=-0.9E-9, ! Dry deposition (gases) - Henrys const.
PF0=-9, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=29.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)
/
&SPECIES_PARAMS
PSPECIES="AERO-TRACE", ! Tracer name
PDECAY=-9.9, ! Species half life