Index: options/SPECIES/SPECIES_001
===================================================================
--- options/SPECIES/SPECIES_001	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_001	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,4 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-TRACER             Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-350.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="TRACER          ",
+ PWEIGHTMOLAR=  350.000000    ,
+ /
Index: options/SPECIES/SPECIES_001.oldformat
===================================================================
--- options/SPECIES/SPECIES_001.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_001.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+TRACER             Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+350.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_002
===================================================================
--- options/SPECIES/SPECIES_002	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_002	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,9 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-O3                 Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
- 1.5               Dry deposition (gases) - D
-1.0e-02            Dry deposition (gases) - Henrys const.
-1.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 48.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="O3              ",
+ PRELDIFF=  1.50000000    ,
+ PHENRY=  1.0E-02,
+ PF0=  1.00000000    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  48.0000000    ,
+ /
Index: options/SPECIES/SPECIES_002.oldformat
===================================================================
--- options/SPECIES/SPECIES_002.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_002.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+O3                 Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+ 1.5               Dry deposition (gases) - D
+1.0e-02            Dry deposition (gases) - Henrys const.
+1.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 48.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_003
===================================================================
--- options/SPECIES/SPECIES_003	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_003	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,10 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-NO                 Tracer name
--999.9             Species half life
-8.0E-06            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
-1.2                Dry deposition (gases) - D
-3.0e-03            Dry deposition (gases) - Henrys const.
-0.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 30.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="NO              ",
+ PWETA_GAS=  8.0E-06,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  1.2    ,
+ PHENRY=  3.0E-03,
+ PF0=  0.0    ,
+ PDQUER=  0.00000000    ,
+ PWEIGHTMOLAR=  30.0000000    ,
+ /
Index: options/SPECIES/SPECIES_003.oldformat
===================================================================
--- options/SPECIES/SPECIES_003.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_003.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+NO                 Tracer name
+-999.9             Species half life
+8.0E-06            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.2                Dry deposition (gases) - D
+3.0e-03            Dry deposition (gases) - Henrys const.
+0.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 30.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_004
===================================================================
--- options/SPECIES/SPECIES_004	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_004	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,11 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-NO2                Tracer name
--999.9             Species half life
-1.0E-05            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
-1.6                Dry deposition (gases) - D
-1.0e-02            Dry deposition (gases) - Henrys const.
-0.1                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 46.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="NO2             ",
+ PWETA_GAS=  1.0E-05,
+ PWETB_GAS= 0.620000005    ,
+ PRELDIFF=  1.60000002    ,
+ PHENRY=  1.0E-02,
+ PF0= 0.1    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  46.0000000    ,
+ /
Index: options/SPECIES/SPECIES_004.oldformat
===================================================================
--- options/SPECIES/SPECIES_004.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_004.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+NO2                Tracer name
+-999.9             Species half life
+1.0E-05            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.6                Dry deposition (gases) - D
+1.0e-02            Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 46.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_005
===================================================================
--- options/SPECIES/SPECIES_005	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_005	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,11 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-HNO3               Tracer name
--999.9             Species half life
-5.0E-05            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
-1.9                Dry deposition (gases) - D
-1.0e+14            Dry deposition (gases) - Henrys const.
-0.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 63.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="HNO3            ",
+ PWETA_GAS=  5.0E-05,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  1.9    ,
+ PHENRY=  1.00000000E+14,
+ PF0=  0.00000000    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  63.0000000    ,
+ /
Index: options/SPECIES/SPECIES_005.oldformat
===================================================================
--- options/SPECIES/SPECIES_005.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_005.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+HNO3               Tracer name
+-999.9             Species half life
+5.0E-05            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.9                Dry deposition (gases) - D
+1.0e+14            Dry deposition (gases) - Henrys const.
+0.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 63.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_006
===================================================================
--- options/SPECIES/SPECIES_006	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_006	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,9 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-HNO2               Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
- 1.6               Dry deposition (gases) - D
-1.0e+05            Dry deposition (gases) - Henrys const.
-0.1                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 47.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="HNO2            ",
+ PRELDIFF=  1.6    ,
+ PHENRY=  100000.000    ,
+ PF0= 0.1    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  47.0000000    ,
+ /
Index: options/SPECIES/SPECIES_006.oldformat
===================================================================
--- options/SPECIES/SPECIES_006.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_006.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+HNO2               Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.6                Dry deposition (gases) - D
+1.0e+05            Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 47.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_007
===================================================================
--- options/SPECIES/SPECIES_007	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_007	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,12 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-H2O2               Tracer name
--999.9             Species half life
-1.0E-04            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
- 1.4               Dry deposition (gases) - D
-1.0e+05            Dry deposition (gases) - Henrys const.
-1.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 34.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="H2O2            ",
+ PDECAY= -999.900024    ,
+ PWETA_GAS=  1.0E-04,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  1.4    ,
+ PHENRY=  100000.000    ,
+ PF0=  1.00000000    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  34.0000000    ,
+ /
Index: options/SPECIES/SPECIES_007.oldformat
===================================================================
--- options/SPECIES/SPECIES_007.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_007.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+H2O2               Tracer name
+-999.9             Species half life
+1.0E-04            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.4                Dry deposition (gases) - D
+1.0e+05            Dry deposition (gases) - Henrys const.
+1.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 34.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_008
===================================================================
--- options/SPECIES/SPECIES_008	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_008	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -7,9 +7,11 @@
 SO2var             Tracer name
 -999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
- 2.0               Dry deposition (gases) - D
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+2.0                Dry deposition (gases) - D
 1.0e+05            Dry deposition (gases) - Henrys const.
 0.0                Dry deposition (gases) - f0 (reactivity)
Index: options/SPECIES/SPECIES_008.oldformat
===================================================================
--- options/SPECIES/SPECIES_008.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_008.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,60 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+SO2var             Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+2.0                Dry deposition (gases) - D
+1.0e+05            Dry deposition (gases) - Henrys const.
+0.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 64.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
+hr_start so2_area  so2_point
+ 0       0.636       0.814      0-1 local time
+ 1       0.571       0.767      1-2 local time
+ 2       0.523       0.743
+ 3       0.495       0.735
+ 4       0.507       0.751
+ 5       0.577       0.798
+ 6       0.732       0.873
+ 7       0.929       0.976
+ 8       1.015       1.037
+ 9       1.077       1.084
+10       1.153       1.126
+11       1.241       1.142
+12       1.315       1.157
+13       1.370       1.163
+14       1.422       1.171
+15       1.463       1.169
+16       1.463       1.162
+17       1.426       1.141
+18       1.326       1.117
+19       1.225       1.109
+20       1.082       1.095
+21       0.928       1.044
+22       0.804       0.959
+23       0.717       0.869      23-24 local time
+week_day so2_area  so2_point
+1        1.060       1.000      Monday
+2        1.060       1.000      Tuesday
+3        1.060       1.000      Wednesday
+4        1.060       1.000      Thursday
+5        1.060       1.000      Friday
+6        0.900       1.000      Saturday
+7        0.800       1.000      Sunday
Index: options/SPECIES/SPECIES_009
===================================================================
--- options/SPECIES/SPECIES_009	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_009	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,9 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-HCHO               Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
- 1.3               Dry deposition (gases) - D
-6.0e+03            Dry deposition (gases) - Henrys const.
-0.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 30.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="HCHO            ",
+ PRELDIFF=  1.3    ,
+ PHENRY=  6000.00000    ,
+ PF0=  0.00000000    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  30.0000000    ,
+ /
Index: options/SPECIES/SPECIES_009.oldformat
===================================================================
--- options/SPECIES/SPECIES_009.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_009.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+HCHO               Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.3                Dry deposition (gases) - D
+6.0e+03            Dry deposition (gases) - Henrys const.
+0.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 30.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_010
===================================================================
--- options/SPECIES/SPECIES_010	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_010	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,9 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-PAN                Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
- 2.6               Dry deposition (gases) - D
-3.6e+00            Dry deposition (gases) - Henrys const.
-0.1                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-121.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="PAN             ",
+ PRELDIFF=  2.6    ,
+ PHENRY=  3.6    ,
+ PF0= 0.1    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  121.000000    ,
+ /
Index: options/SPECIES/SPECIES_010.oldformat
===================================================================
--- options/SPECIES/SPECIES_010.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_010.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+PAN                Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+2.6                Dry deposition (gases) - D
+3.6e+00            Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+121.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_011
===================================================================
--- options/SPECIES/SPECIES_011	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_011	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,11 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-NH3                Tracer name
--999.9             Species half life
-9.9E-05            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
- 1.1               Dry deposition (gases) - D
-2.0e+14            Dry deposition (gases) - Henrys const.
-0.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 17.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="NH3             ",
+ PWETA_GAS=  9.9E-05,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  1.1    ,
+ PHENRY=  2.0E+14,
+ PF0=  0.00000000    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR=  17.0000000    ,
+ /
Index: options/SPECIES/SPECIES_011.oldformat
===================================================================
--- options/SPECIES/SPECIES_011.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_011.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+NH3                Tracer name
+-999.9             Species half life
+9.9E-05            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+1.1                Dry deposition (gases) - D
+2.0e+14            Dry deposition (gases) - Henrys const.
+0.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 17.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_012
===================================================================
--- options/SPECIES/SPECIES_012	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_012	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,11 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-SO4-aero           Tracer name
--999.9             Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.0E03             Dry deposition (particles) - rho
-4.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="SO4-aero        ",
+ PCRAIN_AERO=1.0    ,
+ PCSNOW_AERO=1.0    ,
+ PCCN_AERO=0.9    ,
+ PIN_AERO=0.1    ,
+ PDENSITY=2000.00000    ,
+ PDQUER=4.0E-07   ,
+ PDSIGMA=0.3    ,
+ PWEIGHTMOLAR=  -9.9    ,
+ /
Index: options/SPECIES/SPECIES_012.oldformat
===================================================================
--- options/SPECIES/SPECIES_012.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_012.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+SO4-aero           Tracer name
+-999.9             Species half life
+-999.9             Below-cloud scavenging (gases) - A (weta_gas)       
+-999.9             Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.0E03             Dry deposition (particles) - rho
+4.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_013
===================================================================
--- options/SPECIES/SPECIES_013	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_013	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,13 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-NO3-aero           Tracer name
--999.9             Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.0E03             Dry deposition (particles) - rho
-4.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="NO3-aero        ",
+ PCRAIN_AERO=  1.00000000    ,
+ PCSNOW_AERO=  1.00000000    ,
+ PCCN_AERO= 0.9    ,
+ PIN_AERO= 0.1    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY=  2000.00000    ,
+ PDQUER=  4.0E-07,
+ PDSIGMA= 0.3    ,
+ PWEIGHTMOLAR= -9.99    ,
+ /
Index: options/SPECIES/SPECIES_013.oldformat
===================================================================
--- options/SPECIES/SPECIES_013.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_013.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+NO3-aero           Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.0E03             Dry deposition (particles) - rho
+4.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_014
===================================================================
--- options/SPECIES/SPECIES_014	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_014	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,12 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-I2-131             Tracer name
-691200.0           Species half life
-8.0E-05            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi 
- 2.7               Dry deposition (gases) - D
-1.0e+05            Dry deposition (gases) - Henrys const.
-0.1                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="I2-131          ",
+ PDECAY=  691200.000    ,
+ PWETA_GAS=  8.0E-05,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  2.70    ,
+ PHENRY=  100000.000    ,
+ PF0= 0.1    ,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PWEIGHTMOLAR= -9.99    ,
+ /
Index: options/SPECIES/SPECIES_014.oldformat
===================================================================
--- options/SPECIES/SPECIES_014.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_014.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+I2-131             Tracer name
+691200.0           Species half life
+8.0E-05            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+2.7                Dry deposition (gases) - D
+1.0e+05            Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_015
===================================================================
--- options/SPECIES/SPECIES_015	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_015	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,14 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-I-131              Tracer name
-691200.0           Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.5E03             Dry deposition (particles) - rho
-6.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="I-131           ",
+ PDECAY=  691200.000    ,
+ PCRAIN_AERO=  1.00000000    ,
+ PCSNOW_AERO=  1.00000000    ,
+ PCCN_AERO= 0.9    ,
+ PIN_AERO= 0.1    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY=  2500.00000    ,
+ PDQUER=  6.0E-07,
+ PDSIGMA= 0.3    ,
+ PWEIGHTMOLAR= -9.9    ,
+ /
Index: options/SPECIES/SPECIES_015.oldformat
===================================================================
--- options/SPECIES/SPECIES_015.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_015.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+I-131              Tracer name
+691200.0           Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.5E03             Dry deposition (particles) - rho
+6.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_016
===================================================================
--- options/SPECIES/SPECIES_016	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_016	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,23 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-Cs-137             Tracer name
--999.9             Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.5E03             Dry deposition (particles) - rho
-6.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="Cs-137          ",
+ PDECAY= -999.900024    ,
+ PWETA_GAS= -9.89999962    ,
+ PWETB_GAS= -9.89999962    ,
+ PCRAIN_AERO=  1.00000000    ,
+ PCSNOW_AERO=  1.00000000    ,
+ PCCN_AERO= 0.899999976    ,
+ PIN_AERO= 0.100000001    ,
+ PRELDIFF= -9.89999962    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY=  2500.00000    ,
+ PDQUER=  6.00000021E-07,
+ PDSIGMA= 0.300000012    ,
+ PDRYVEL= -9.98999977    ,
+ PWEIGHTMOLAR= -9.98999977    ,
+ POHCCONST= -9.98999994E-09,
+ POHDCONST= -9.89999962    ,
+ POHNCONST= -9.89999962    ,
+ PSPEC_ASS= -9.00000000    ,
+ PKAO= -99.9899979    ,
+ /
Index: options/SPECIES/SPECIES_016.oldformat
===================================================================
--- options/SPECIES/SPECIES_016.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_016.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+Cs-137             Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.5E03             Dry deposition (particles) - rho
+6.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_017
===================================================================
--- options/SPECIES/SPECIES_017	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_017	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,14 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-Y-91               Tracer name
-5037120.0          Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.5E03             Dry deposition (particles) - rho
-6.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="Y-91            ",
+ PDECAY=  5037120.00    ,
+ PCRAIN_AERO=  1.00000000    ,
+ PCSNOW_AERO=  1.00000000    ,
+ PCCN_AERO= 0.9    ,
+ PIN_AERO= 0.1    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY=  2500.00000    ,
+ PDQUER=  6.0E-07,
+ PDSIGMA= 0.3    ,
+ PWEIGHTMOLAR= -9.9    ,
+ /
Index: options/SPECIES/SPECIES_017.oldformat
===================================================================
--- options/SPECIES/SPECIES_017.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_017.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+Y-91               Tracer name
+5037120.0          Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.5E03             Dry deposition (particles) - rho
+6.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_018
===================================================================
--- options/SPECIES/SPECIES_018	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_018	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,14 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-Ru-106             Tracer name
-31536000.0         Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.5E03             Dry deposition (particles) - rho
-6.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="Ru-106          ",
+ PDECAY=  31536000.0    ,
+ PCRAIN_AERO=  1.00000000    ,
+ PCSNOW_AERO=  1.00000000    ,
+ PCCN_AERO= 0.9    ,
+ PIN_AERO= 0.1    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY=  2500.00000    ,
+ PDQUER=  6.0E-07,
+ PDSIGMA= 0.3    ,
+ PWEIGHTMOLAR=  -9.9    ,
+ /
Index: options/SPECIES/SPECIES_018.oldformat
===================================================================
--- options/SPECIES/SPECIES_018.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_018.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+Ru-106             Tracer name
+31536000.0         Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.5E03             Dry deposition (particles) - rho
+6.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_019
===================================================================
--- options/SPECIES/SPECIES_019	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_019	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,5 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-Kr-85              Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
--9.99              molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="Kr-85           ",
+ PDENSITY= -9.9E+09,
+ PWEIGHTMOLAR= -9.9    ,
+ /
Index: options/SPECIES/SPECIES_019.oldformat
===================================================================
--- options/SPECIES/SPECIES_019.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_019.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+Kr-85              Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+-9.99              molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_020
===================================================================
--- options/SPECIES/SPECIES_020	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_020	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,12 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-Sr-90              Tracer name
--999.9             Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-2.5E03             Dry deposition (particles) - rho
-6.0E-7             Dry deposition (particles) - dquer
-3.0E-1             Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- -9.99             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="Sr-90           ",
+ PDECAY= -999.900024    ,
+ PCRAIN_AERO=  1.00000000    ,
+ PCSNOW_AERO=  1.00000000    ,
+ PCCN_AERO= 0.9    ,
+ PIN_AERO= 0.1    ,
+ PDENSITY=  2500.00000    ,
+ PDQUER=  6.0E-07,
+ PDSIGMA= 0.3    ,
+ PWEIGHTMOLAR=  -9.9    ,
+ /
Index: options/SPECIES/SPECIES_020.oldformat
===================================================================
--- options/SPECIES/SPECIES_020.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_020.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+Sr-90              Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+2.5E03             Dry deposition (particles) - rho
+6.0E-7             Dry deposition (particles) - dquer
+3.0E-1             Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ -9.99             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_021
===================================================================
--- options/SPECIES/SPECIES_021	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_021	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,5 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-Xe-133             Tracer name
-453168.0           Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
--9.99              molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="Xe-133          ",
+ PDECAY=  453168.000    ,
+ PWEIGHTMOLAR=  -9.9    ,
+ /
Index: options/SPECIES/SPECIES_021.oldformat
===================================================================
--- options/SPECIES/SPECIES_021.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_021.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+Xe-133             Tracer name
+453168.0           Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+-9.99              molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_022
===================================================================
--- options/SPECIES/SPECIES_022	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_022	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,4 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-CO                 Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 28.00             molweight
-??                 OH Reaction rate - C [cm^3/molecule/sec]
-??                 OH reaction rate - D [K]
-??                 OH Reaction rate - N (no unit)
--9                 number of associated specias (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="CO              ",
+ PWEIGHTMOLAR=  28.0000000    ,
+ /
Index: options/SPECIES/SPECIES_022.oldformat
===================================================================
--- options/SPECIES/SPECIES_022.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_022.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+CO                 Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 28.00             molweight
+-9.9               OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated specias (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_023
===================================================================
--- options/SPECIES/SPECIES_023	(revision d486c43c828d7d072cf249311015e51297baf0b7)
+++ options/SPECIES/SPECIES_023	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,6 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-SO2dryOH           Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
- 2.0               Dry deposition (gases) - D
-1.0e+05            Dry deposition (gases) - Henrys const.
-0.0                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 64.00             molweight
-??                 OH Reaction rate - C [cm^3/molecule/sec]
-??                 OH reaction rate - D [K]
-2.0                OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="SO2dryOH        ",
+ PRELDIFF=  2.00000000    ,
+ PHENRY=  100000.000    ,
+ PWEIGHTMOLAR=  64.0000000    ,
+ /
Index: options/SPECIES/SPECIES_023.oldformat
===================================================================
--- options/SPECIES/SPECIES_023.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_023.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+SO2dryOH           Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+2.0                Dry deposition (gases) - D
+1.0e+05            Dry deposition (gases) - Henrys const.
+0.0                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 64.00             molweight
+-9.9               OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_024
===================================================================
--- options/SPECIES/SPECIES_024	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_024	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,26 +1,23 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-AIRTRACER             Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-29.00              molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
-
+&SPECIES_PARAMS
+ PSPECIES="AIRTRACER       ",
+ PDECAY= -999.900024    ,
+ PWETA_GAS= -9.89999962    ,
+ PWETB_GAS= -9.89999962    ,
+ PCRAIN_AERO= -9.89999962    ,
+ PCSNOW_AERO= -9.89999962    ,
+ PCCN_AERO= -9.89999962    ,
+ PIN_AERO= -9.89999962    ,
+ PRELDIFF= -9.89999962    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY= -9.90000026E+09,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PDRYVEL= -9.98999977    ,
+ PWEIGHTMOLAR=  29.0000000    ,
+ POHCCONST= -9.98999994E-09,
+ POHDCONST= -9.89999962    ,
+ POHNCONST= -9.89999962    ,
+ PSPEC_ASS= -9.00000000    ,
+ PKAO= -99.9899979    ,
+ /
Index: options/SPECIES/SPECIES_024.oldformat
===================================================================
--- options/SPECIES/SPECIES_024.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_024.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+AIRTRACER             Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+29.00              molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_025
===================================================================
--- options/SPECIES/SPECIES_025	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_025	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,11 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-AERO-TRACER        Tracer name
--999.9             Species half life
-1.0                Below-cloud scavenging (gases and particles) - A
-1.0                Below-cloud scavenging (gases and particles)- B
-0.9                In-cloud scavenging (gases and particles) - Ai
-0.1                In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
-1.4E+03            Dry deposition (particles) - rho
-2.5E-07            Dry deposition (particles) - dquer
-1.25               Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
- 29.00             molweight
--9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
--9.9               OH reaction rate - D [K]
--9.9               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="AERO-TRACE      ",
+ PCRAIN_AERO=  1.0    ,
+ PCSNOW_AERO=  1.0    ,
+ PCCN_AERO= 0.9    ,
+ PIN_AERO= 0.1    ,
+ PDENSITY=  1400.00000    ,
+ PDQUER=  2.5E-07,
+ PDSIGMA=  1.25000000    ,
+ PWEIGHTMOLAR=  29.0000000    ,
+ /
Index: options/SPECIES/SPECIES_025.oldformat
===================================================================
--- options/SPECIES/SPECIES_025.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_025.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+AERO-TRACER        Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+1.0                Below-cloud scavenging (particles) - Crain (crain_aero)
+1.0                Below-cloud scavenging (particles) - Csnow (csnow_aero)
+0.9                In-cloud scavenging (particles) - CCNeff (ccn_aero)
+0.1                In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+1.4E+03            Dry deposition (particles) - rho
+2.5E-07            Dry deposition (particles) - dquer
+1.25               Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+ 29.00             molweight
+-9.99E-09          OH Reaction rate - C [cm^3/molecule/sec]
+-9.9               OH reaction rate - D [K]
+-9.9               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_026
===================================================================
--- options/SPECIES/SPECIES_026	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_026	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,8 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-CH4                Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-16.00              molweight
-9.65E-20           OH Reaction rate - C [cm^3/molecule/sec]
-1082.0             OH reaction rate - D [K]
-2.58               OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="CH4             ",
+ PDENSITY= -9.9+09,
+ PWEIGHTMOLAR=  16.0000000    ,
+ POHCCONST=  9.65E-20,
+ POHDCONST=  1082.00000    ,
+ POHNCONST=  2.58    ,
+ /
Index: options/SPECIES/SPECIES_026.oldformat
===================================================================
--- options/SPECIES/SPECIES_026.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_026.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+CH4                Tracer name
+-999.9             Species half lif
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+16.00              molweight
+9.65E-20           OH Reaction rate - C [cm^3/molecule/sec]
+1082.0             OH reaction rate - D [K]
+2.58               OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_027
===================================================================
--- options/SPECIES/SPECIES_027	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_027	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,23 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-C2H6               Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-30.00              molweight
-15.2E-18           OH Reaction rate - C [cm^3/molecule/sec]
-498.0              OH reaction rate - D [K]
-2.0                OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="C2H6            ",
+ PDECAY= -999.900024    ,
+ PWETA_GAS= -9.89999962    ,
+ PWETB_GAS= -9.89999962    ,
+ PCRAIN_AERO= -9.89999962    ,
+ PCSNOW_AERO= -9.89999962    ,
+ PCCN_AERO= -9.89999962    ,
+ PIN_AERO= -9.89999962    ,
+ PRELDIFF= -9.89999962    ,
+ PHENRY=  0.00000000    ,
+ PF0=  0.00000000    ,
+ PDENSITY= -9.90000026E+09,
+ PDQUER=  0.00000000    ,
+ PDSIGMA=  0.00000000    ,
+ PDRYVEL= -9.98999977    ,
+ PWEIGHTMOLAR=  30.0000000    ,
+ POHCCONST=  1.51999996E-17,
+ POHDCONST=  498.000000    ,
+ POHNCONST=  2.00000000    ,
+ PSPEC_ASS= -9.00000000    ,
+ PKAO= -99.9899979    ,
+ /
Index: options/SPECIES/SPECIES_027.oldformat
===================================================================
--- options/SPECIES/SPECIES_027.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_027.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+C2H6               Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+30.00              molweight
+15.2E-18           OH Reaction rate - C [cm^3/molecule/sec]
+498.0              OH reaction rate - D [K]
+2.0                OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_028
===================================================================
--- options/SPECIES/SPECIES_028	(revision 0f20c31f55c80f501e55c69aad5a2dd4081628f1)
+++ options/SPECIES/SPECIES_028	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,7 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*     This file is on the format required by version 10 (05/2015)          *
-****************************************************************************
-C3H8               Tracer name
--999.9             Species half life
--9.9               Below-cloud scavenging (gases and particles) - A
--9.9               Below-cloud scavenging (gases and particles)- B
--9.9               In-cloud scavenging (gases and particles) - Ai
--9.9               In-cloud scavenging (gases and particles) - Bi
--9.9               Dry deposition (gases) - D
-                   Dry deposition (gases) - Henrys const.
-                   Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-44.000             molweight
-15.5E-18           OH Reaction rate - C [cm^3/molecule/sec]
-61.0               OH reaction rate - D [K]
-2.0                OH Reaction rate - N (no unit)
--9                 number of associated species (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="C3H8            ",
+ PWEIGHTMOLAR=  44.0000000    ,
+ POHCCONST=  1.55-E17,
+ POHDCONST=  61.0000000    ,
+ POHNCONST=  2.00000000    ,
+ /
Index: options/SPECIES/SPECIES_028.oldformat
===================================================================
--- options/SPECIES/SPECIES_028.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_028.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*     This file is on the format required by version 10 (05/2015)          *
+****************************************************************************
+C3H8               Tracer name
+-999.9             Species half life
+-9.9               Below-cloud scavenging (gases) - A (weta_gas)
+-9.9               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+-9.9               Dry deposition (gases) - D
+                   Dry deposition (gases) - Henrys const.
+                   Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+44.000             molweight
+15.5E-18           OH Reaction rate - C [cm^3/molecule/sec]
+61.0               OH reaction rate - D [K]
+2.0                OH Reaction rate - N (no unit)
+-9                 number of associated species (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_031
===================================================================
--- options/SPECIES/SPECIES_031	(revision d4ba52840ad7900e8913b3e649483181eb493a2a)
+++ options/SPECIES/SPECIES_031	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,12 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*                                                                          *
-****************************************************************************
-PCB28              Tracer name
--999.9             Species half life
-3.5E-05            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
- 4.3               Dry deposition (gases) - D
-3.1e-00            Dry deposition (gases) - Henrys const.
-0.1                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-257.54             molweight
-1.07E-11           OH Reaction rate - C [cm^3/molecule/sec]
-1203.0             OH reaction rate - D [K]
-2.0                OH Reaction rate - N (no unit)
--9                 number of associated specias (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="PCB28           ",
+ PWETA_GAS=  3.5E-05,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  4.3    ,
+ PHENRY=  3.1    ,
+ PF0= 0.1    ,
+ PWEIGHTMOLAR=  257.54    ,
+ POHCCONST=  1.07E-11,
+ POHDCONST=  1203.00000    ,
+ POHNCONST=  2.00000000    ,
+ /
Index: options/SPECIES/SPECIES_031.oldformat
===================================================================
--- options/SPECIES/SPECIES_031.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_031.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*                                                                          *
+****************************************************************************
+PCB28              Tracer name
+-999.9             Species half life
+3.5E-05            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+4.3                Dry deposition (gases) - D
+3.1E-00            Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+257.54             molweight
+1.07E-11           OH Reaction rate - C [cm^3/molecule/sec]
+1203.0             OH reaction rate - D [K]
+2.0                OH Reaction rate - N (no unit)
+-9                 number of associated specias (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: options/SPECIES/SPECIES_034
===================================================================
--- options/SPECIES/SPECIES_034	(revision d4ba52840ad7900e8913b3e649483181eb493a2a)
+++ options/SPECIES/SPECIES_034	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -1,25 +1,12 @@
-****************************************************************************
-*                                                                          *
-*     Input file for the Lagrangian particle dispersion model FLEXPART     *
-*               Definition file of chemical species/radionuclides          *
-*                                                                          *
-****************************************************************************
-G-HCH              Tracer name
--999.9             Species half life
-3.5E-05            Below-cloud scavenging (gases and particles) - A
-0.62               Below-cloud scavenging (gases and particles)- B
-1.0                In-cloud scavenging (gases and particles) - Ai
-1.0                In-cloud scavenging (gases and particles) - Bi
-4.3                Dry deposition (gases) - D
-4.22E+02           Dry deposition (gases) - Henrys const.
-0.1                Dry deposition (gases) - f0 (reactivity)
--9.9E09            Dry deposition (particles) - rho
-                   Dry deposition (particles) - dquer
-                   Dry deposition (particles) - dsig
--9.99              Alternative: dry deposition velocity
-290.83             molweight
-6.21E-11           OH Reaction rate - C [cm^3/molecule/sec]
-1203.0             OH reaction rate - D [K]
-2.0                OH Reaction rate - N (no unit)
--9                 number of associated specias (neg. none)
--99.99             KOA - organic matter air partitioning
+&SPECIES_PARAMS
+ PSPECIES="G-HCH           ",
+ PWETA_GAS=  3.5E-05,
+ PWETB_GAS= 0.62    ,
+ PRELDIFF=  4.3    ,
+ PHENRY=  422.0    ,
+ PF0= 0.10    ,
+ PWEIGHTMOLAR=  290.83    ,
+ POHCCONST=  6.21E-11,
+ POHDCONST=  1203.00000    ,
+ POHNCONST=  2.00000000    ,
+ /
Index: options/SPECIES/SPECIES_034.oldformat
===================================================================
--- options/SPECIES/SPECIES_034.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
+++ options/SPECIES/SPECIES_034.oldformat	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -0,0 +1,27 @@
+****************************************************************************
+*                                                                          *
+*     Input file for the Lagrangian particle dispersion model FLEXPART     *
+*               Definition file of chemical species/radionuclides          *
+*                                                                          *
+****************************************************************************
+G-HCH              Tracer name
+-999.9             Species half life
+3.5E-05            Below-cloud scavenging (gases) - A (weta_gas)
+0.62               Below-cloud scavenging (gases) - B (wetb_gas)
+-9.9               Below-cloud scavenging (particles) - Crain (crain_aero)
+-9.9               Below-cloud scavenging (particles) - Csnow (csnow_aero)
+-9.9               In-cloud scavenging (particles) - CCNeff (ccn_aero)
+-9.9               In-cloud scavenging (particles) - INeff (in_aero)
+4.3                Dry deposition (gases) - D
+4.22E+02           Dry deposition (gases) - Henrys const.
+0.1                Dry deposition (gases) - f0 (reactivity)
+-9.9E09            Dry deposition (particles) - rho
+                   Dry deposition (particles) - dquer
+                   Dry deposition (particles) - dsig
+-9.99              Alternative: dry deposition velocity
+290.83             molweight
+6.21E-11           OH Reaction rate - C [cm^3/molecule/sec]
+1203.0             OH reaction rate - D [K]
+2.0                OH Reaction rate - N (no unit)
+-9                 number of associated specias (neg. none)
+-99.99             KOA - organic matter air partitioning
Index: src/com_mod.f90
===================================================================
--- src/com_mod.f90	(revision db712a8ffdd9c478928d0173162055a7c109b178)
+++ src/com_mod.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -160,7 +160,8 @@
   !real xmass(maxpoint,maxspec)
   real :: decay(maxspec)
-  real :: weta(maxspec),wetb(maxspec)
+  real :: weta_gas(maxspec),wetb_gas(maxspec)
+  real :: crain_aero(maxspec),csnow_aero(maxspec)
 ! NIK: 31.01.2013- parameters for in-cloud scavening
-  real :: weta_in(maxspec), wetb_in(maxspec)
+  real :: ccn_aero(maxspec),in_aero(maxspec)
   real :: reldiff(maxspec),henry(maxspec),f0(maxspec)
   real :: density(maxspec),dquer(maxspec),dsigma(maxspec)
@@ -197,7 +198,7 @@
 
   ! WET DEPOSITION
-  ! weta, wetb              parameters for determining below-cloud wet scavenging coefficients
-  ! weta_in, wetb_in       parameters for determining in-cloud wet scavenging coefficients
-  ! wetc_in, wetd_in       parameters for determining in-cloud wet scavenging coefficients
+  ! weta_gas, wetb_gas     parameters for below-cloud wet scavenging coefficients (gasses)
+  ! crain_aero, in_aero    parameters for in-cloud wet scavenging coefficients (aerosols)
+  ! wetc_in, wetd_in       parameters for in-cloud wet scavenging coefficients (aerosols)
 
   ! GAS DEPOSITION
@@ -216,5 +217,5 @@
   ! dquer [m]               mean diameter of particles
   ! dsigma                  dsigma=10 or dsigma=0.1 means that 68% of the
-  !                    mass are between 0.1*dquer and 10*dquer
+  !                         mass are between 0.1*dquer and 10*dquer
 
   ! fract                   mass fraction of each diameter interval
@@ -371,6 +372,6 @@
 
 !ZHG Sep 2015  
-   real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)
-   real :: clw4(0:nxmax-1,0:nymax-1,numwfmem) ! eso: =icloud_stats(:,:,4,:)
+!   real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)
+   real :: ctwc(0:nxmax-1,0:nymax-1,numwfmem) ! eso: =icloud_stats(:,:,4,:)
 
 
@@ -392,4 +393,5 @@
 
   ! pplev for the GFS version
+  ! ctwc                 total cloud water content
 
   ! 2d fields
@@ -489,5 +491,5 @@
   real,allocatable,dimension(:,:,:,:,:) :: uun, vvn, wwn, ttn, qvn, pvn,&
        & rhon, drhodzn, tthn, qvhn, clwcn, ciwcn, clwn, clwchn, ciwchn
-  real,allocatable,dimension(:,:,:,:) :: clw4n
+  real,allocatable,dimension(:,:,:,:) :: ctwcn
   integer,allocatable,dimension(:,:,:,:) :: cloudshn
   integer(kind=1),allocatable,dimension(:,:,:,:,:) :: cloudsn
@@ -802,15 +804,10 @@
     allocate(clwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
     allocate(ciwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests))
-    allocate(clw4n(0:nxmax-1,0:nymax-1,numwfmem,numbnests))
-
-!    clw4n(:,:,:,:)=0.
+    allocate(ctwcn(0:nxmax-1,0:nymax-1,numwfmem,numbnests))
+
     clwcn(:,:,:,:,:)=0.
     ciwcn(:,:,:,:,:)=0.
     clwchn(:,:,:,:,:)=0.
     ciwchn(:,:,:,:,:)=0.
-    ! clwn(:,:,:,:,:)=0.
-
-    ! cloudsn(:,:,:,:,:)=0
-    ! cloudshn(:,:,:,:)=0
     
   end subroutine com_mod_allocate_nests
Index: src/makefile
===================================================================
--- src/makefile	(revision 32b49c31198ae7800bda9ceb9125f98584d75077)
+++ src/makefile	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -66,5 +66,5 @@
 
 ## LIBRARIES
-LIBS = -lgrib_api_f90 -lgrib_api -lm -ljasper -lnetcdff -llapack  # -fopenmp
+LIBS = -lgrib_api_f90 -lgrib_api -lm -ljasper -lnetcdff # -fopenmp
 
 FFLAGS   = -I$(INCPATH1) -I$(INCPATH2) -O$(O_LEV) -g -m64 -mcmodel=medium -fconvert=little-endian -frecord-marker=4 -fmessage-length=0 -flto=jobserver -O$(O_LEV) $(FUSER) # -march=native
Index: src/mpi_mod.f90
===================================================================
--- src/mpi_mod.f90	(revision 32b49c31198ae7800bda9ceb9125f98584d75077)
+++ src/mpi_mod.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -980,5 +980,5 @@
       ! call MPI_Bcast(icloud_stats(:,:,:,li:ui),d2s1*5,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
       ! if (mp_ierr /= 0) goto 600
-      call MPI_Bcast(clw4(:,:,li:ui),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
+      call MPI_Bcast(ctwc(:,:,li:ui),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
       if (mp_ierr /= 0) goto 600
       ! call MPI_Bcast(clwc(:,:,:,li:ui),d3s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
@@ -1151,5 +1151,5 @@
       ! call MPI_Bcast(icloud_stats(:,:,:,li:ui),d2s1*5,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
       ! if (mp_ierr /= 0) goto 600
-      call MPI_Bcast(clw4n(:,:,li:ui,i),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
+      call MPI_Bcast(ctwcn(:,:,li:ui,i),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)
       if (mp_ierr /= 0) goto 600
     end if
@@ -1360,5 +1360,5 @@
         ! call MPI_Isend(icloud_stats(:,:,:,mind),d2s1*5,mp_sp,dest,tm1,&
         !      &MPI_COMM_WORLD,reqs(i),mp_ierr)
-        call MPI_Isend(clw4(:,:,mind),d2s1,mp_sp,dest,tm1,&
+        call MPI_Isend(ctwc(:,:,mind),d2s1,mp_sp,dest,tm1,&
              &MPI_COMM_WORLD,reqs(i),mp_ierr)
 
@@ -1573,5 +1573,5 @@
       ! call MPI_Irecv(icloud_stats(:,:,:,mind),d2s1*5,mp_sp,id_read,MPI_ANY_TAG,&
       !      &MPI_COMM_WORLD,reqs(j),mp_ierr)
-      call MPI_Irecv(clw4(:,:,mind),d2s1*5,mp_sp,id_read,MPI_ANY_TAG,&
+      call MPI_Irecv(ctwc(:,:,mind),d2s1*5,mp_sp,id_read,MPI_ANY_TAG,&
            &MPI_COMM_WORLD,reqs(j),mp_ierr)
       if (mp_ierr /= 0) goto 600
Index: src/netcdf_output_mod.f90
===================================================================
--- src/netcdf_output_mod.f90	(revision db712a8ffdd9c478928d0173162055a7c109b178)
+++ src/netcdf_output_mod.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -61,6 +61,6 @@
                        dx,xlon0,dy,ylat0,compoint,method,lsubgrid,lconvection,&
                        ind_source,ind_receptor,nageclass,lage,&
-                       drydep,wetdep,decay,weta,wetb, numbnests, &
-                       weta_in,wetb_in, & ! wetc_in,wetd_in, &
+                       drydep,wetdep,decay,weta_gas,wetb_gas, numbnests, &
+                       ccn_aero,in_aero, & ! wetc_in,wetd_in, &
                        reldiff,henry,f0,density,dquer,dsigma,dryvel,&
 !                       weightmolar,ohreact,spec_ass,kao,vsetaver,&
@@ -543,8 +543,8 @@
              chunksizes = dep_chunksizes))
         call nf90_err(nf90_put_att(ncid, wdsID, 'units', '1e-12 kg m-2'))
-        call nf90_err(nf90_put_att(ncid, wdsID, 'weta', weta(i)))
-        call nf90_err(nf90_put_att(ncid, wdsID, 'wetb', wetb(i)))
-        call nf90_err(nf90_put_att(ncid, wdsID, 'weta_in', weta_in(i)))
-        call nf90_err(nf90_put_att(ncid, wdsID, 'wetb_in', wetb_in(i)))
+        call nf90_err(nf90_put_att(ncid, wdsID, 'weta_gas', weta_gas(i)))
+        call nf90_err(nf90_put_att(ncid, wdsID, 'wetb_gas', wetb_gas(i)))
+        call nf90_err(nf90_put_att(ncid, wdsID, 'ccn_aero', ccn_aero(i)))
+        call nf90_err(nf90_put_att(ncid, wdsID, 'in_aero', in_aero(i)))
         ! call nf90_err(nf90_put_att(ncid, wdsID, 'wetc_in', wetc_in(i)))
         ! call nf90_err(nf90_put_att(ncid, wdsID, 'wetd_in', wetd_in(i)))
Index: src/par_mod.f90
===================================================================
--- src/par_mod.f90	(revision 32b49c31198ae7800bda9ceb9125f98584d75077)
+++ src/par_mod.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -89,5 +89,5 @@
   real,parameter :: rho_water=1000. !ZHG 2015 [kg/m3]
   !ZHG MAR2016
-  real,parameter :: incloud_ratio=1.6
+  real,parameter :: incloud_ratio=2.2
 
   ! karman                  Karman's constant
Index: src/readreleases.f90
===================================================================
--- src/readreleases.f90	(revision b0434e1968915cd3d4a8afcbf2cc275acecb9ff7)
+++ src/readreleases.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -57,7 +57,7 @@
   ! xpoint2,ypoint2     geograf. coordinates of upper right corner of release  *
   !                     area                                                   *
-  ! weta, wetb          parameters to determine the below-cloud scavenging     *
-  ! weta_in, wetb_in    parameters to determine the in-cloud scavenging        *
-  ! wetc_in, wetd_in    parameters to determine the in-cloud scavenging        *
+  ! weta_gas, wetb_gas  parameters for below-cloud scavenging (gas)            *
+  ! crain_aero, csnow_aero  parameters for below-cloud scavenging (aerosol)    *
+  ! ccn_aero, in_aero   parameters for in-cloud scavenging (aerosol)           *
   ! zpoint1,zpoint2     height range, over which release takes place           *
   ! num_min_discrete    if less, release cannot be randomized and happens at   *
@@ -363,5 +363,5 @@
   !***********************************************************
 ! NIK 15.02.2015, with new wet dep scheme either A or B parameters can be assigned a positive value which switches on wet dep
-    if (weta(i).gt.0. .or. wetb(i).gt.0.)  then
+    if (weta_gas(i).gt.0. .or. wetb_gas(i).gt.0.)  then ! :TODO: or particle
       WETDEP=.true.
       if (lroot) then
@@ -374,10 +374,10 @@
     
 ! NIK 31.01.2013 + 10.12.2013 + 15.02.2015
-    if (weta_in(i).gt.0. .or. wetb_in(i).gt.0.)  then
+    if (ccn_aero(i).gt.0. .or. in_aero(i).gt.0.)  then
       WETDEP=.true.
       if (lroot) then
         write (*,*) 'In-cloud scavenging: ON'
 !        write (*,*) 'In-cloud scavenging coefficients: ',&
-!           &weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i
+!           &ccn_aero(i),in_aero(i),i !,wetc_in(i), wetd_in(i),i
       end if
     else
Index: src/readspecies.f90
===================================================================
--- src/readspecies.f90	(revision 32b49c31198ae7800bda9ceb9125f98584d75077)
+++ src/readspecies.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -42,7 +42,7 @@
   ! decaytime(maxtable)  half time for radiological decay                      *
   ! specname(maxtable)   names of chemical species, radionuclides              *
-  ! weta, wetb           Parameters for determining below-cloud scavenging     *
-  ! weta_in              Parameter for determining in-cloud scavenging         *
-  ! wetb_in              Parameter for determining in-cloud scavenging         *
+  ! weta_gas, wetb_gas   Parameters for determining below-cloud scavenging     *
+  ! ccn_aero              Parameter for determining in-cloud scavenging         *
+  ! in_aero              Parameter for determining in-cloud scavenging         *
   ! ohcconst             OH reaction rate constant C                           *
   ! ohdconst             OH reaction rate constant D                           *
@@ -66,24 +66,24 @@
 
   character(len=16) :: pspecies
-  real :: pdecay, pweta, pwetb, preldiff, phenry, pf0, pdensity, pdquer
+  real :: pdecay, pweta_gas, pwetb_gas, preldiff, phenry, pf0, pdensity, pdquer
   real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao
-  real :: pweta_in, pwetb_in, pwetc_in, pwetd_in
+  real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero
   integer :: readerror
 
 ! declare namelist
   namelist /species_params/ &
-       pspecies, pdecay, pweta, pwetb, &
-       pweta_in, pwetb_in, pwetc_in, pwetd_in, &
+       pspecies, pdecay, pweta_gas, pwetb_gas, &
+       pcrain_aero, pcsnow_aero, pccn_aero, pin_aero, &
        preldiff, phenry, pf0, pdensity, pdquer, &
        pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao
 
-  pspecies=" "
+  pspecies="" ! read failure indicator value
   pdecay=-999.9
-  pweta=-9.9E-09
-  pwetb=0.0
-  pweta_in=-9.9E-09
-  pwetb_in=-9.9E-09
-!  pwetc_in=-9.9E-09
-!  pwetd_in=-9.9E-09
+  pweta_gas=-9.9E-09
+  pwetb_gas=0.0
+  pcrain_aero=-9.9E-09
+  pcsnow_aero=-9.9E-09
+  pccn_aero=-9.9E-09
+  pin_aero=-9.9E-09
   preldiff=-9.9
   phenry=0.0
@@ -98,5 +98,5 @@
   pspec_ass=-9
   pkao=-99.99
-  pweightmolar=-789.0 ! read failure indicator value
+  pweightmolar=-999.9
 
 ! Open the SPECIES file and read species names and properties
@@ -113,5 +113,6 @@
   close(unitspecies)
 
-  if ((pweightmolar.eq.-789.0).or.(readerror.ne.0)) then ! no namelist found
+!  if ((pweightmolar.eq.-789.0).or.(readerror.ne.0)) then ! no namelist found
+  if ((len(pspecies).eq.0).or.(readerror.ne.0)) then ! no namelist found
 
     readerror=1
@@ -127,19 +128,17 @@
     read(unitspecies,'(f18.1)',end=22) decay(pos_spec)
 !  write(*,*) decay(pos_spec)
-    read(unitspecies,'(e18.1)',end=22) weta(pos_spec)
-!  write(*,*) weta(pos_spec)
-    read(unitspecies,'(f18.2)',end=22) wetb(pos_spec)
-!  write(*,*) wetb(pos_spec)
-
+    read(unitspecies,'(e18.1)',end=22) weta_gas(pos_spec)
+!  write(*,*) weta_gas(pos_spec)
+    read(unitspecies,'(f18.2)',end=22) wetb_gas(pos_spec)
+!  write(*,*) wetb_gas(pos_spec)
+    read(unitspecies,'(e18.1)',end=22) crain_aero(pos_spec)
+!  write(*,*) crain_aero(pos_spec)
+    read(unitspecies,'(f18.2)',end=22) csnow_aero(pos_spec)
+!  write(*,*) csnow_aero(pos_spec)
 !*** NIK 31.01.2013: including in-cloud scavening parameters
-    read(unitspecies,'(e18.1)',end=22) weta_in(pos_spec)
-!  write(*,*) weta_in(pos_spec)
-    read(unitspecies,'(f18.2)',end=22) wetb_in(pos_spec)
-!  write(*,*) wetb_in(pos_spec)
-! read(unitspecies,'(f18.2)',end=22) wetc_in(pos_spec)
-!  write(*,*) wetc_in(pos_spec)
-! read(unitspecies,'(f18.2)',end=22) wetd_in(pos_spec)
-!  write(*,*) wetd_in(pos_spec)
-
+    read(unitspecies,'(e18.1)',end=22) ccn_aero(pos_spec)
+!  write(*,*) ccn_aero(pos_spec)
+    read(unitspecies,'(f18.2)',end=22) in_aero(pos_spec)
+!  write(*,*) in_aero(pos_spec)
     read(unitspecies,'(f18.1)',end=22) reldiff(pos_spec)
 !  write(*,*) reldiff(pos_spec)
@@ -171,10 +170,10 @@
     pspecies=species(pos_spec)
     pdecay=decay(pos_spec)
-    pweta=weta(pos_spec)
-    pwetb=wetb(pos_spec)
-    pweta_in=weta_in(pos_spec)
-    pwetb_in=wetb_in(pos_spec)
-!    pwetc_in=wetc_in(pos_spec)
-!    pwetd_in=wetd_in(pos_spec)
+    pweta_gas=weta_gas(pos_spec)
+    pwetb_gas=wetb_gas(pos_spec)
+    pcrain_aero=crain_aero(pos_spec)
+    pcsnow_aero=csnow_aero(pos_spec)
+    pccn_aero=ccn_aero(pos_spec)
+    pin_aero=in_aero(pos_spec)
     preldiff=reldiff(pos_spec)
     phenry=henry(pos_spec)
@@ -195,10 +194,10 @@
     species(pos_spec)=pspecies
     decay(pos_spec)=pdecay
-    weta(pos_spec)=pweta
-    wetb(pos_spec)=pwetb
-    weta_in(pos_spec)=pweta_in
-    wetb_in(pos_spec)=pwetb_in
-!    wetc_in(pos_spec)=pwetc_in
-!    wetd_in(pos_spec)=pwetd_in
+    weta_gas(pos_spec)=pweta_gas
+    wetb_gas(pos_spec)=pwetb_gas
+    crain_aero=pcrain_aero
+    csnow_aero=pcsnow_aero
+    ccn_aero(pos_spec)=pccn_aero
+    in_aero(pos_spec)=pin_aero
     reldiff(pos_spec)=preldiff
     henry(pos_spec)=phenry
@@ -222,14 +221,14 @@
 ! Check scavenging parameters given in SPECIES file
 
-  if (weta(pos_spec).gt.0.0 .or. wetb(pos_spec).gt.0.0 .or. weta_in(pos_spec).gt.0.0 .or. wetb_in(pos_spec).gt.0.0) then
+  if (weta_gas(pos_spec).gt.0.0 .or. wetb_gas(pos_spec).gt.0.0 .or. ccn_aero(pos_spec).gt.0.0 .or. in_aero(pos_spec).gt.0.0) then
 
   if (dquer(pos_spec).gt.0) then !is particle
     if (lroot) then
       write(*,'(a,f5.2)') '  Particle below-cloud scavenging parameter A &
-           &(Rain collection efficiency)  ', weta(pos_spec)
+           &(Rain collection efficiency)  ', crain_aero(pos_spec)
       write(*,'(a,f5.2)') '  Particle below-cloud scavenging parameter B &
-           &(Snow collection efficiency)  ', wetb(pos_spec)
-    end if
-    if (weta(pos_spec).gt.1.0 .or. wetb(pos_spec).gt.1.0) then
+           &(Snow collection efficiency)  ', csnow_aero(pos_spec)
+    end if
+    if (weta_gas(pos_spec).gt.1.0 .or. wetb_gas(pos_spec).gt.1.0) then
       if (lroot) then
         write(*,*) '*******************************************'
@@ -242,9 +241,9 @@
     if (lroot) then
       write(*,'(a,f5.2)') '  Particle in-cloud scavenging parameter Ai &
-           &(CCN efficiency)  ', weta_in(pos_spec)
+           &(CCN efficiency)  ', ccn_aero(pos_spec)
       write(*,'(a,f5.2)') '  Particle in-cloud scavenging parameter Bi &
-           &(IN efficiency)  ', wetb_in(pos_spec)
-    end if
-    if (weta_in(pos_spec).gt.1.0 .or. weta_in(pos_spec).lt.0.01) then
+           &(IN efficiency)  ', in_aero(pos_spec)
+    end if
+    if (ccn_aero(pos_spec).gt.1.0 .or. ccn_aero(pos_spec).lt.0.01) then
       if (lroot) then
         write(*,*) '*******************************************'
@@ -254,5 +253,5 @@
       end if
     endif
-    if (wetb_in(pos_spec).gt.1.0 .or. wetb_in(pos_spec).lt.0.01) then
+    if (in_aero(pos_spec).gt.1.0 .or. in_aero(pos_spec).lt.0.01) then
       if (lroot) then
         write(*,*) '*******************************************'
@@ -265,10 +264,10 @@
   else !is gas
     if (lroot) then
-      write(*,*) '  Gas below-cloud scavenging parameter A  ', weta(pos_spec)
-      write(*,'(a,f5.2)') '  Gas below-cloud scavenging parameter B  ', wetb(pos_spec)
+      write(*,*) '  Gas below-cloud scavenging parameter A  ', weta_gas(pos_spec)
+      write(*,'(a,f5.2)') '  Gas below-cloud scavenging parameter B  ', wetb_gas(pos_spec)
       write(*,*) ' Gas in-cloud scavenging uses default values as in Hertel et al 1995'
     end if
-    if (weta(pos_spec).gt.0.) then !if wet deposition is turned on
-      if (weta(pos_spec).gt.1E-04 .or. weta(pos_spec).lt.1E-09) then
+    if (weta_gas(pos_spec).gt.0.) then !if wet deposition is turned on
+      if (weta_gas(pos_spec).gt.1E-04 .or. weta_gas(pos_spec).lt.1E-09) then
         if (lroot) then
           write(*,*) '*******************************************'
@@ -279,6 +278,6 @@
       endif
     end if
-    if (wetb(pos_spec).gt.0.) then !if wet deposition is turned on
-      if (wetb(pos_spec).gt.0.8 .or. wetb(pos_spec).lt.0.6) then
+    if (wetb_gas(pos_spec).gt.0.) then !if wet deposition is turned on
+      if (wetb_gas(pos_spec).gt.0.8 .or. wetb_gas(pos_spec).lt.0.6) then
         if (lroot) then
           write(*,*) '*******************************************'
@@ -292,5 +291,5 @@
   endif
 
-  if (((weta(pos_spec).gt.0).or.(wetb(pos_spec).gt.0)).and.(henry(pos_spec).le.0)) then
+  if (((weta_gas(pos_spec).gt.0).or.(wetb_gas(pos_spec).gt.0)).and.(henry(pos_spec).le.0)) then
     if (dquer(pos_spec).le.0) goto 996 ! no particle, no henry set
   endif
Index: src/timemanager_mpi.f90
===================================================================
--- src/timemanager_mpi.f90	(revision 32b49c31198ae7800bda9ceb9125f98584d75077)
+++ src/timemanager_mpi.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -169,5 +169,5 @@
     if (mp_dev_mode) write(*,*) 'myid, itime: ',mp_pid,itime
     
-    if (wetdep .and. itime .ne. 0 .and. numpart .gt. 0) then
+    if (WETDEP .and. itime .ne. 0 .and. numpart .gt. 0) then
       if (verbosity.gt.0) then
         write (*,*) 'timemanager> call wetdepo'
Index: src/verttransform.f90
===================================================================
--- src/verttransform.f90	(revision db712a8ffdd9c478928d0173162055a7c109b178)
+++ src/verttransform.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -82,5 +82,5 @@
 
   integer :: ix,jy,kz,iz,n,kmin,ix1,jy1,ixp,jyp,ixm,jym,kz_inv
-  real :: f_qvsat,pressure,rh,lsp,convp,prec
+  real :: f_qvsat,pressure,rh,lsp,convp,cloudh_min,prec
   real :: ew,dz1,dz2,dz
   real :: xlon,ylat,xlonr,dzdx,dzdy
@@ -586,5 +586,6 @@
     write(*,*) 'Global ECMWF fields: using cloud water'
     clw(:,:,:,n)=0.0
-    icloud_stats(:,:,:,n)=0.0
+!    icloud_stats(:,:,:,n)=0.0
+    ctwc(:,:,n)=0.0
     clouds(:,:,:,n)=0
 ! If water/ice are read separately into clwc and ciwc, store sum in clwc
@@ -604,10 +605,13 @@
             clw(ix,jy,kz,n)=(clwc(ix,jy,kz,n)*rho(ix,jy,kz,n))*(height(kz+1)-height(kz))
             tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 
-            icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n)          ! Column cloud water [m3/m3]
-            icloud_stats(ix,jy,3,n)= min(height(kz+1),height(kz))                     ! Cloud BOT height stats      [m]
+            
+!            icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n)          ! Column cloud water [m3/m3]
+            ctwc(ix,jy,n) = ctwc(ix,jy,n)+clw(ix,jy,kz,n)
+!            icloud_stats(ix,jy,3,n)= min(height(kz+1),height(kz))                     ! Cloud BOT height stats      [m]
+            cloudh_min=min(height(kz+1),height(kz))
 !ZHG 2015 extra for testing
 !            clh(ix,jy,kz,n)=height(kz+1)-height(kz)
-            icloud_stats(ix,jy,1,n)=icloud_stats(ix,jy,1,n)+(height(kz+1)-height(kz)) ! Cloud total vertical extent [m]
-            icloud_stats(ix,jy,2,n)= max(icloud_stats(ix,jy,2,n),height(kz))          ! Cloud TOP height            [m]
+!            icloud_stats(ix,jy,1,n)=icloud_stats(ix,jy,1,n)+(height(kz+1)-height(kz)) ! Cloud total vertical extent [m]
+!            icloud_stats(ix,jy,2,n)= max(icloud_stats(ix,jy,2,n),height(kz))          ! Cloud TOP height            [m]
 !ZHG
           endif
@@ -626,5 +630,5 @@
                 clouds(ix,jy,kz,n)=2                             ! convp in-cloud
               endif                                              ! convective or large scale
-            elseif((clw(ix,jy,kz,n).le.0) .and. (icloud_stats(ix,jy,3,n).ge.height(kz)) ) then ! is below cloud
+            elseif((clw(ix,jy,kz,n).le.0) .and. (cloudh_min.ge.height(kz))) then ! is below cloud
               if (lsp.ge.convp) then
                 clouds(ix,jy,kz,n)=5                             ! lsp dominated washout
@@ -643,5 +647,5 @@
 
 ! eso: copy the relevant data to clw4 to reduce amount of communicated data for MPI
-    clw4(:,:,n) = icloud_stats(:,:,4,n)
+!    ctwc(:,:,n) = icloud_stats(:,:,4,n)
 
 !**************************************************************************
Index: src/verttransform_nests.f90
===================================================================
--- src/verttransform_nests.f90	(revision db712a8ffdd9c478928d0173162055a7c109b178)
+++ src/verttransform_nests.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -466,5 +466,5 @@
       clwn(0:nxm1,0:nym1,:,n,l)=0.0
 !    icloud_stats(0:nxm1,0:nym1,:,n)=0.0
-      clw4n(0:nxm1,0:nym1,n,l)=0.0
+      ctwcn(0:nxm1,0:nym1,n,l)=0.0
       cloudsn(0:nxm1,0:nym1,:,n,l)=0
 ! If water/ice are read separately into clwc and ciwc, store sum in clwcn
@@ -484,5 +484,5 @@
               clwn(ix,jy,kz,n,l)=(clwcn(ix,jy,kz,n,l)*rhon(ix,jy,kz,n,l))*(height(kz+1)-height(kz))
               tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 
-              clw4n(ix,jy,n,l) = clw4n(ix,jy,n,l)+clwn(ix,jy,kz,n,l)
+              ctwcn(ix,jy,n,l) = ctwcn(ix,jy,n,l)+clwn(ix,jy,kz,n,l)
 !            icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n)          ! Column cloud water [m3/m3]
 !           icloud_stats(ix,jy,3,n)= min(height(kz+1),height(kz))                     ! Cloud BOT height stats      [m]
Index: src/wetdepo.f90
===================================================================
--- src/wetdepo.f90	(revision ec7fc72438c14086c97319dbcbf75d1e0463ed5f)
+++ src/wetdepo.f90	(revision 341f4b7eb1c45f4749b96fc8dbea2dfff30140f2)
@@ -270,39 +270,46 @@
       if (clouds_v.ge.4) then !below cloud
 
-        if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then !if positive below-cloud parameters given in SPECIES file (either A or B)
+! For gas: if positive below-cloud parameters (A or B), and dquer<=0
+!******************************************************************
+        if ((dquer(ks).le.0.).and.(weta_gas(ks).gt.0..or.wetb_gas(ks).gt.0.)) then
+          !        if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then 
           blc_count=blc_count+1
-
-!GAS
-          if (dquer(ks) .le. 0.) then  !is gas
-            wetscav=weta(ks)*prec(1)**wetb(ks)
-
-!AEROSOL
-          else !is particle
+          wetscav=weta_gas(ks)*prec(1)**wetb_gas(ks)
+
+! For aerosols: if positive below-cloud parameters (Crain/Csnow or B), and dquer>0
+!*********************************************************************************
+        else if ((dquer(ks).gt.0.).and.(crain_aero(ks).gt.0..or.csnow_aero(ks).gt.0.)) then
+          blc_count=blc_count+1
+
 !NIK 17.02.2015
-! For the calculation here particle size needs to be in meter and not um as dquer is changed to in readreleases
-! for particles larger than 10 um use the largest size defined in the parameterizations (10um)
-            dquer_m=min(10.,dquer(ks))/1000000. !conversion from um to m
-            if (act_temp .ge. 273. .and. weta(ks).gt.0.)  then !Rain 
+! For the calculation here particle size needs to be in meter and not um as dquer is
+! changed in readreleases
+! For particles larger than 10 um use the largest size defined in the parameterizations (10um)
+          dquer_m=min(10.,dquer(ks))/1000000. !conversion from um to m
+
+! Rain:
+          if (act_temp .ge. 273. .and. crain_aero(ks).gt.0.)  then
+
 ! ZHG 2014 : Particle RAIN scavenging coefficient based on Laakso et al 2003, 
-! the below-cloud scavenging (rain efficienty) 
-! parameter A (=weta) from SPECIES file
-              wetscav= weta(ks)*10**(bclr(1)+ (bclr(2)*(log10(dquer_m))**(-4))+(bclr(3)*(log10(dquer_m))**(-3))+ (bclr(4)* &
-                   (log10(dquer_m))**(-2))+ (bclr(5)*(log10(dquer_m))**(-1))+ bclr(6)* (prec(1))**(0.5))
-
-            elseif (act_temp .lt. 273. .and. wetb(ks).gt.0.)  then ! Snow
+! the below-cloud scavenging (rain efficienty) parameter Crain (=crain_aero) from SPECIES file
+            wetscav=crain_aero(ks)*10**(bclr(1)+(bclr(2)*(log10(dquer_m))**(-4))+ &
+                 & (bclr(3)*(log10(dquer_m))**(-3))+ (bclr(4)*(log10(dquer_m))**(-2))+&
+                 &(bclr(5)*(log10(dquer_m))**(-1))+bclr(6)* (prec(1))**(0.5))
+
+! Snow:
+          elseif (act_temp .lt. 273. .and. csnow_aero(ks).gt.0.)  then 
 ! ZHG 2014 : Particle SNOW scavenging coefficient based on Kyro et al 2009, 
-! the below-cloud scavenging (Snow efficiency) 
-! parameter B (=wetb) from SPECIES file
-              wetscav= wetb(ks)*10**(bcls(1)+ (bcls(2)*(log10(dquer_m))**(-4))+(bcls(3)*(log10(dquer_m))**(-3))+ (bcls(4)* &
-                   (log10(dquer_m))**(-2))+ (bcls(5)*(log10(dquer_m))**(-1))+ bcls(6)* (prec(1))**(0.5))
-
-            endif
-
+! the below-cloud scavenging (Snow efficiency) parameter Csnow (=csnow_aero) from SPECIES file
+            wetscav=csnow_aero(ks)*10**(bcls(1)+(bcls(2)*(log10(dquer_m))**(-4))+&
+                 &(bcls(3)*(log10(dquer_m))**(-3))+ (bcls(4)*(log10(dquer_m))**(-2))+&
+                 &(bcls(5)*(log10(dquer_m))**(-1))+ bcls(6)* (prec(1))**(0.5))
+
+          endif
+          
 !             write(*,*) 'bl-cloud, act_temp=',act_temp, ',prec=',prec(1),',wetscav=', wetscav, ', jpart=',jpart
 
-          endif !gas or particle
-        endif ! positive below-cloud scavenging parameters given in Species file
+        endif ! gas or particle
+!      endif ! positive below-cloud scavenging parameters given in Species file
       endif !end BELOW
-
 
 !********************
@@ -310,17 +317,18 @@
 !********************
       if (clouds_v.lt.4) then ! In-cloud
-! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are given in species file
-        if (weta_in(ks).gt.0. .or. wetb_in(ks).gt.0.) then 
+! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are
+! given in species file, or if gas and positive Henry's constant
+        if ((ccn_aero(ks).gt.0. .or. in_aero(ks).gt.0.).or.(henry(ks).gt.0.and.dquer(ks).le.0)) then 
           inc_count=inc_count+1
 ! if negative coefficients (turned off) set to zero for use in equation
-          if (weta_in(ks).lt.0.) weta_in(ks)=0.
-          if (wetb_in(ks).lt.0.) wetb_in(ks)=0.
+          if (ccn_aero(ks).lt.0.) ccn_aero(ks)=0.
+          if (in_aero(ks).lt.0.) in_aero(ks)=0.
 
 !ZHG 2015 Cloud liquid & ice water (CLWC+CIWC) from ECMWF
 ! nested fields
           if (ngrid.gt.0.and.readclouds_this_nest) then
-            cl=clw4n(ix,jy,n,ngrid)*(grfraction(1)/cc)
+            cl=ctwcn(ix,jy,n,ngrid)*(grfraction(1)/cc)
           else if (ngrid.eq.0.and.readclouds) then
-            cl=clw4(ix,jy,n)*(grfraction(1)/cc)
+            cl=ctwc(ix,jy,n)*(grfraction(1)/cc)
           else                                  !parameterize cloudwater m2/m3
 !ZHG updated parameterization of cloud water to better reproduce the values coming from ECMWF
@@ -337,5 +345,5 @@
           end if
 ! ZHG: Calculate the aerosol partition based on cloud phase and Ai and Bi
-          frac_act = liq_frac*weta_in(ks) +(1-liq_frac)*wetb_in(ks)
+          frac_act = liq_frac*ccn_aero(ks) +(1-liq_frac)*in_aero(ks)
 
 !ZHG Use the activated fraction and the liqid water to calculate the washout
@@ -343,5 +351,5 @@
 ! AEROSOL
 !**************************************************
-          if (dquer(ks).gt. 0.) then ! is particle
+          if (dquer(ks).gt.0.) then ! is particle
 
             S_i= frac_act/cl
