Changeset 341f4b7 in flexpart.git
- Timestamp:
- Apr 20, 2016, 1:39:51 PM (8 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 61df8d9
- Parents:
- 32b49c3
- Files:
-
- 30 added
- 41 edited
-
options/SPECIES/SPECIES_001 (modified) (1 diff)
-
options/SPECIES/SPECIES_001.oldformat (added)
-
options/SPECIES/SPECIES_002 (modified) (1 diff)
-
options/SPECIES/SPECIES_002.oldformat (added)
-
options/SPECIES/SPECIES_003 (modified) (1 diff)
-
options/SPECIES/SPECIES_003.oldformat (added)
-
options/SPECIES/SPECIES_004 (modified) (1 diff)
-
options/SPECIES/SPECIES_004.oldformat (added)
-
options/SPECIES/SPECIES_005 (modified) (1 diff)
-
options/SPECIES/SPECIES_005.oldformat (added)
-
options/SPECIES/SPECIES_006 (modified) (1 diff)
-
options/SPECIES/SPECIES_006.oldformat (added)
-
options/SPECIES/SPECIES_007 (modified) (1 diff)
-
options/SPECIES/SPECIES_007.oldformat (added)
-
options/SPECIES/SPECIES_008 (modified) (1 diff)
-
options/SPECIES/SPECIES_008.oldformat (added)
-
options/SPECIES/SPECIES_009 (modified) (1 diff)
-
options/SPECIES/SPECIES_009.oldformat (added)
-
options/SPECIES/SPECIES_010 (modified) (1 diff)
-
options/SPECIES/SPECIES_010.oldformat (added)
-
options/SPECIES/SPECIES_011 (modified) (1 diff)
-
options/SPECIES/SPECIES_011.oldformat (added)
-
options/SPECIES/SPECIES_012 (modified) (1 diff)
-
options/SPECIES/SPECIES_012.oldformat (added)
-
options/SPECIES/SPECIES_013 (modified) (1 diff)
-
options/SPECIES/SPECIES_013.oldformat (added)
-
options/SPECIES/SPECIES_014 (modified) (1 diff)
-
options/SPECIES/SPECIES_014.oldformat (added)
-
options/SPECIES/SPECIES_015 (modified) (1 diff)
-
options/SPECIES/SPECIES_015.oldformat (added)
-
options/SPECIES/SPECIES_016 (modified) (1 diff)
-
options/SPECIES/SPECIES_016.oldformat (added)
-
options/SPECIES/SPECIES_017 (modified) (1 diff)
-
options/SPECIES/SPECIES_017.oldformat (added)
-
options/SPECIES/SPECIES_018 (modified) (1 diff)
-
options/SPECIES/SPECIES_018.oldformat (added)
-
options/SPECIES/SPECIES_019 (modified) (1 diff)
-
options/SPECIES/SPECIES_019.oldformat (added)
-
options/SPECIES/SPECIES_020 (modified) (1 diff)
-
options/SPECIES/SPECIES_020.oldformat (added)
-
options/SPECIES/SPECIES_021 (modified) (1 diff)
-
options/SPECIES/SPECIES_021.oldformat (added)
-
options/SPECIES/SPECIES_022 (modified) (1 diff)
-
options/SPECIES/SPECIES_022.oldformat (added)
-
options/SPECIES/SPECIES_023 (modified) (1 diff)
-
options/SPECIES/SPECIES_023.oldformat (added)
-
options/SPECIES/SPECIES_024 (modified) (1 diff)
-
options/SPECIES/SPECIES_024.oldformat (added)
-
options/SPECIES/SPECIES_025 (modified) (1 diff)
-
options/SPECIES/SPECIES_025.oldformat (added)
-
options/SPECIES/SPECIES_026 (modified) (1 diff)
-
options/SPECIES/SPECIES_026.oldformat (added)
-
options/SPECIES/SPECIES_027 (modified) (1 diff)
-
options/SPECIES/SPECIES_027.oldformat (added)
-
options/SPECIES/SPECIES_028 (modified) (1 diff)
-
options/SPECIES/SPECIES_028.oldformat (added)
-
options/SPECIES/SPECIES_031 (modified) (1 diff)
-
options/SPECIES/SPECIES_031.oldformat (added)
-
options/SPECIES/SPECIES_034 (modified) (1 diff)
-
options/SPECIES/SPECIES_034.oldformat (added)
-
src/com_mod.f90 (modified) (7 diffs)
-
src/makefile (modified) (1 diff)
-
src/mpi_mod.f90 (modified) (4 diffs)
-
src/netcdf_output_mod.f90 (modified) (2 diffs)
-
src/par_mod.f90 (modified) (1 diff)
-
src/readreleases.f90 (modified) (3 diffs)
-
src/readspecies.f90 (modified) (13 diffs)
-
src/timemanager_mpi.f90 (modified) (1 diff)
-
src/verttransform.f90 (modified) (5 diffs)
-
src/verttransform_nests.f90 (modified) (2 diffs)
-
src/wetdepo.f90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
options/SPECIES/SPECIES_001
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 TRACER Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 350.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="TRACER ", 3 PWEIGHTMOLAR= 350.000000 , 4 / -
options/SPECIES/SPECIES_002
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 O3 Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 1.5 Dry deposition (gases) - D 14 1.0e-02 Dry deposition (gases) - Henrys const. 15 1.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 48.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="O3 ", 3 PRELDIFF= 1.50000000 , 4 PHENRY= 1.0E-02, 5 PF0= 1.00000000 , 6 PDQUER= 0.00000000 , 7 PDSIGMA= 0.00000000 , 8 PWEIGHTMOLAR= 48.0000000 , 9 / -
options/SPECIES/SPECIES_003
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 NO Tracer name 8 -999.9 Species half life 9 8.0E-06 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 1.2 Dry deposition (gases) - D 14 3.0e-03 Dry deposition (gases) - Henrys const. 15 0.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 30.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="NO ", 3 PWETA_GAS= 8.0E-06, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.2 , 6 PHENRY= 3.0E-03, 7 PF0= 0.0 , 8 PDQUER= 0.00000000 , 9 PWEIGHTMOLAR= 30.0000000 , 10 / -
options/SPECIES/SPECIES_004
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 NO2 Tracer name 8 -999.9 Species half life 9 1.0E-05 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 1.6 Dry deposition (gases) - D 14 1.0e-02 Dry deposition (gases) - Henrys const. 15 0.1 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 46.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="NO2 ", 3 PWETA_GAS= 1.0E-05, 4 PWETB_GAS= 0.620000005 , 5 PRELDIFF= 1.60000002 , 6 PHENRY= 1.0E-02, 7 PF0= 0.1 , 8 PDQUER= 0.00000000 , 9 PDSIGMA= 0.00000000 , 10 PWEIGHTMOLAR= 46.0000000 , 11 / -
options/SPECIES/SPECIES_005
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 HNO3 Tracer name 8 -999.9 Species half life 9 5.0E-05 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 1.9 Dry deposition (gases) - D 14 1.0e+14 Dry deposition (gases) - Henrys const. 15 0.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 63.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="HNO3 ", 3 PWETA_GAS= 5.0E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.9 , 6 PHENRY= 1.00000000E+14, 7 PF0= 0.00000000 , 8 PDQUER= 0.00000000 , 9 PDSIGMA= 0.00000000 , 10 PWEIGHTMOLAR= 63.0000000 , 11 / -
options/SPECIES/SPECIES_006
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 HNO2 Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 1.6 Dry deposition (gases) - D 14 1.0e+05 Dry deposition (gases) - Henrys const. 15 0.1 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 47.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="HNO2 ", 3 PRELDIFF= 1.6 , 4 PHENRY= 100000.000 , 5 PF0= 0.1 , 6 PDQUER= 0.00000000 , 7 PDSIGMA= 0.00000000 , 8 PWEIGHTMOLAR= 47.0000000 , 9 / -
options/SPECIES/SPECIES_007
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 H2O2 Tracer name 8 -999.9 Species half life 9 1.0E-04 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 1.4 Dry deposition (gases) - D 14 1.0e+05 Dry deposition (gases) - Henrys const. 15 1.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 34.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="H2O2 ", 3 PDECAY= -999.900024 , 4 PWETA_GAS= 1.0E-04, 5 PWETB_GAS= 0.62 , 6 PRELDIFF= 1.4 , 7 PHENRY= 100000.000 , 8 PF0= 1.00000000 , 9 PDQUER= 0.00000000 , 10 PDSIGMA= 0.00000000 , 11 PWEIGHTMOLAR= 34.0000000 , 12 / -
options/SPECIES/SPECIES_008
r0f20c31 r341f4b7 7 7 SO2var Tracer name 8 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 2.0 Dry deposition (gases) - D 9 -9.9 Below-cloud scavenging (gases) - A (weta_gas) 10 -9.9 Below-cloud scavenging (gases) - B (wetb_gas) 11 -9.9 Below-cloud scavenging (particles) - Crain (crain_aero) 12 -9.9 Below-cloud scavenging (particles) - Csnow (csnow_aero) 13 -9.9 In-cloud scavenging (particles) - CCNeff (ccn_aero) 14 -9.9 In-cloud scavenging (particles) - INeff (in_aero) 15 2.0 Dry deposition (gases) - D 14 16 1.0e+05 Dry deposition (gases) - Henrys const. 15 17 0.0 Dry deposition (gases) - f0 (reactivity) -
options/SPECIES/SPECIES_009
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 HCHO Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 1.3 Dry deposition (gases) - D 14 6.0e+03 Dry deposition (gases) - Henrys const. 15 0.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 30.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="HCHO ", 3 PRELDIFF= 1.3 , 4 PHENRY= 6000.00000 , 5 PF0= 0.00000000 , 6 PDQUER= 0.00000000 , 7 PDSIGMA= 0.00000000 , 8 PWEIGHTMOLAR= 30.0000000 , 9 / -
options/SPECIES/SPECIES_010
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 PAN Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 2.6 Dry deposition (gases) - D 14 3.6e+00 Dry deposition (gases) - Henrys const. 15 0.1 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 121.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="PAN ", 3 PRELDIFF= 2.6 , 4 PHENRY= 3.6 , 5 PF0= 0.1 , 6 PDQUER= 0.00000000 , 7 PDSIGMA= 0.00000000 , 8 PWEIGHTMOLAR= 121.000000 , 9 / -
options/SPECIES/SPECIES_011
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 NH3 Tracer name 8 -999.9 Species half life 9 9.9E-05 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 1.1 Dry deposition (gases) - D 14 2.0e+14 Dry deposition (gases) - Henrys const. 15 0.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 17.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="NH3 ", 3 PWETA_GAS= 9.9E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 1.1 , 6 PHENRY= 2.0E+14, 7 PF0= 0.00000000 , 8 PDQUER= 0.00000000 , 9 PDSIGMA= 0.00000000 , 10 PWEIGHTMOLAR= 17.0000000 , 11 / -
options/SPECIES/SPECIES_012
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 SO4-aero Tracer name 8 -999.9 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.0E03 Dry deposition (particles) - rho 17 4.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="SO4-aero ", 3 PCRAIN_AERO=1.0 , 4 PCSNOW_AERO=1.0 , 5 PCCN_AERO=0.9 , 6 PIN_AERO=0.1 , 7 PDENSITY=2000.00000 , 8 PDQUER=4.0E-07 , 9 PDSIGMA=0.3 , 10 PWEIGHTMOLAR= -9.9 , 11 / -
options/SPECIES/SPECIES_013
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 NO3-aero Tracer name 8 -999.9 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.0E03 Dry deposition (particles) - rho 17 4.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="NO3-aero ", 3 PCRAIN_AERO= 1.00000000 , 4 PCSNOW_AERO= 1.00000000 , 5 PCCN_AERO= 0.9 , 6 PIN_AERO= 0.1 , 7 PHENRY= 0.00000000 , 8 PF0= 0.00000000 , 9 PDENSITY= 2000.00000 , 10 PDQUER= 4.0E-07, 11 PDSIGMA= 0.3 , 12 PWEIGHTMOLAR= -9.99 , 13 / -
options/SPECIES/SPECIES_014
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 I2-131 Tracer name 8 691200.0 Species half life 9 8.0E-05 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 2.7 Dry deposition (gases) - D 14 1.0e+05 Dry deposition (gases) - Henrys const. 15 0.1 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="I2-131 ", 3 PDECAY= 691200.000 , 4 PWETA_GAS= 8.0E-05, 5 PWETB_GAS= 0.62 , 6 PRELDIFF= 2.70 , 7 PHENRY= 100000.000 , 8 PF0= 0.1 , 9 PDQUER= 0.00000000 , 10 PDSIGMA= 0.00000000 , 11 PWEIGHTMOLAR= -9.99 , 12 / -
options/SPECIES/SPECIES_015
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 I-131 Tracer name 8 691200.0 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.5E03 Dry deposition (particles) - rho 17 6.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="I-131 ", 3 PDECAY= 691200.000 , 4 PCRAIN_AERO= 1.00000000 , 5 PCSNOW_AERO= 1.00000000 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PHENRY= 0.00000000 , 9 PF0= 0.00000000 , 10 PDENSITY= 2500.00000 , 11 PDQUER= 6.0E-07, 12 PDSIGMA= 0.3 , 13 PWEIGHTMOLAR= -9.9 , 14 / -
options/SPECIES/SPECIES_016
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 Cs-137 Tracer name 8 -999.9 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.5E03 Dry deposition (particles) - rho 17 6.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="Cs-137 ", 3 PDECAY= -999.900024 , 4 PWETA_GAS= -9.89999962 , 5 PWETB_GAS= -9.89999962 , 6 PCRAIN_AERO= 1.00000000 , 7 PCSNOW_AERO= 1.00000000 , 8 PCCN_AERO= 0.899999976 , 9 PIN_AERO= 0.100000001 , 10 PRELDIFF= -9.89999962 , 11 PHENRY= 0.00000000 , 12 PF0= 0.00000000 , 13 PDENSITY= 2500.00000 , 14 PDQUER= 6.00000021E-07, 15 PDSIGMA= 0.300000012 , 16 PDRYVEL= -9.98999977 , 17 PWEIGHTMOLAR= -9.98999977 , 18 POHCCONST= -9.98999994E-09, 19 POHDCONST= -9.89999962 , 20 POHNCONST= -9.89999962 , 21 PSPEC_ASS= -9.00000000 , 22 PKAO= -99.9899979 , 23 / -
options/SPECIES/SPECIES_017
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 Y-91 Tracer name 8 5037120.0 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.5E03 Dry deposition (particles) - rho 17 6.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="Y-91 ", 3 PDECAY= 5037120.00 , 4 PCRAIN_AERO= 1.00000000 , 5 PCSNOW_AERO= 1.00000000 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PHENRY= 0.00000000 , 9 PF0= 0.00000000 , 10 PDENSITY= 2500.00000 , 11 PDQUER= 6.0E-07, 12 PDSIGMA= 0.3 , 13 PWEIGHTMOLAR= -9.9 , 14 / -
options/SPECIES/SPECIES_018
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 Ru-106 Tracer name 8 31536000.0 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.5E03 Dry deposition (particles) - rho 17 6.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="Ru-106 ", 3 PDECAY= 31536000.0 , 4 PCRAIN_AERO= 1.00000000 , 5 PCSNOW_AERO= 1.00000000 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PHENRY= 0.00000000 , 9 PF0= 0.00000000 , 10 PDENSITY= 2500.00000 , 11 PDQUER= 6.0E-07, 12 PDSIGMA= 0.3 , 13 PWEIGHTMOLAR= -9.9 , 14 / -
options/SPECIES/SPECIES_019
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 Kr-85 Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="Kr-85 ", 3 PDENSITY= -9.9E+09, 4 PWEIGHTMOLAR= -9.9 , 5 / -
options/SPECIES/SPECIES_020
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 Sr-90 Tracer name 8 -999.9 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 2.5E03 Dry deposition (particles) - rho 17 6.0E-7 Dry deposition (particles) - dquer 18 3.0E-1 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="Sr-90 ", 3 PDECAY= -999.900024 , 4 PCRAIN_AERO= 1.00000000 , 5 PCSNOW_AERO= 1.00000000 , 6 PCCN_AERO= 0.9 , 7 PIN_AERO= 0.1 , 8 PDENSITY= 2500.00000 , 9 PDQUER= 6.0E-07, 10 PDSIGMA= 0.3 , 11 PWEIGHTMOLAR= -9.9 , 12 / -
options/SPECIES/SPECIES_021
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 Xe-133 Tracer name 8 453168.0 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 -9.99 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="Xe-133 ", 3 PDECAY= 453168.000 , 4 PWEIGHTMOLAR= -9.9 , 5 / -
options/SPECIES/SPECIES_022
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 CO Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 28.00 molweight 21 ?? OH Reaction rate - C [cm^3/molecule/sec] 22 ?? OH reaction rate - D [K] 23 ?? OH Reaction rate - N (no unit) 24 -9 number of associated specias (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="CO ", 3 PWEIGHTMOLAR= 28.0000000 , 4 / -
options/SPECIES/SPECIES_023
rd486c43 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 SO2dryOH Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 2.0 Dry deposition (gases) - D 14 1.0e+05 Dry deposition (gases) - Henrys const. 15 0.0 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 64.00 molweight 21 ?? OH Reaction rate - C [cm^3/molecule/sec] 22 ?? OH reaction rate - D [K] 23 2.0 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="SO2dryOH ", 3 PRELDIFF= 2.00000000 , 4 PHENRY= 100000.000 , 5 PWEIGHTMOLAR= 64.0000000 , 6 / -
options/SPECIES/SPECIES_024
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 AIRTRACER Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 29.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 26 1 &SPECIES_PARAMS 2 PSPECIES="AIRTRACER ", 3 PDECAY= -999.900024 , 4 PWETA_GAS= -9.89999962 , 5 PWETB_GAS= -9.89999962 , 6 PCRAIN_AERO= -9.89999962 , 7 PCSNOW_AERO= -9.89999962 , 8 PCCN_AERO= -9.89999962 , 9 PIN_AERO= -9.89999962 , 10 PRELDIFF= -9.89999962 , 11 PHENRY= 0.00000000 , 12 PF0= 0.00000000 , 13 PDENSITY= -9.90000026E+09, 14 PDQUER= 0.00000000 , 15 PDSIGMA= 0.00000000 , 16 PDRYVEL= -9.98999977 , 17 PWEIGHTMOLAR= 29.0000000 , 18 POHCCONST= -9.98999994E-09, 19 POHDCONST= -9.89999962 , 20 POHNCONST= -9.89999962 , 21 PSPEC_ASS= -9.00000000 , 22 PKAO= -99.9899979 , 23 / -
options/SPECIES/SPECIES_025
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 AERO-TRACER Tracer name 8 -999.9 Species half life 9 1.0 Below-cloud scavenging (gases and particles) - A 10 1.0 Below-cloud scavenging (gases and particles)- B 11 0.9 In-cloud scavenging (gases and particles) - Ai 12 0.1 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 1.4E+03 Dry deposition (particles) - rho 17 2.5E-07 Dry deposition (particles) - dquer 18 1.25 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 29.00 molweight 21 -9.99E-09 OH Reaction rate - C [cm^3/molecule/sec] 22 -9.9 OH reaction rate - D [K] 23 -9.9 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="AERO-TRACE ", 3 PCRAIN_AERO= 1.0 , 4 PCSNOW_AERO= 1.0 , 5 PCCN_AERO= 0.9 , 6 PIN_AERO= 0.1 , 7 PDENSITY= 1400.00000 , 8 PDQUER= 2.5E-07, 9 PDSIGMA= 1.25000000 , 10 PWEIGHTMOLAR= 29.0000000 , 11 / -
options/SPECIES/SPECIES_026
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 CH4 Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 16.00 molweight 21 9.65E-20 OH Reaction rate - C [cm^3/molecule/sec] 22 1082.0 OH reaction rate - D [K] 23 2.58 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="CH4 ", 3 PDENSITY= -9.9+09, 4 PWEIGHTMOLAR= 16.0000000 , 5 POHCCONST= 9.65E-20, 6 POHDCONST= 1082.00000 , 7 POHNCONST= 2.58 , 8 / -
options/SPECIES/SPECIES_027
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 C2H6 Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 30.00 molweight 21 15.2E-18 OH Reaction rate - C [cm^3/molecule/sec] 22 498.0 OH reaction rate - D [K] 23 2.0 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="C2H6 ", 3 PDECAY= -999.900024 , 4 PWETA_GAS= -9.89999962 , 5 PWETB_GAS= -9.89999962 , 6 PCRAIN_AERO= -9.89999962 , 7 PCSNOW_AERO= -9.89999962 , 8 PCCN_AERO= -9.89999962 , 9 PIN_AERO= -9.89999962 , 10 PRELDIFF= -9.89999962 , 11 PHENRY= 0.00000000 , 12 PF0= 0.00000000 , 13 PDENSITY= -9.90000026E+09, 14 PDQUER= 0.00000000 , 15 PDSIGMA= 0.00000000 , 16 PDRYVEL= -9.98999977 , 17 PWEIGHTMOLAR= 30.0000000 , 18 POHCCONST= 1.51999996E-17, 19 POHDCONST= 498.000000 , 20 POHNCONST= 2.00000000 , 21 PSPEC_ASS= -9.00000000 , 22 PKAO= -99.9899979 , 23 / -
options/SPECIES/SPECIES_028
r0f20c31 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * This file is on the format required by version 10 (05/2015) * 6 **************************************************************************** 7 C3H8 Tracer name 8 -999.9 Species half life 9 -9.9 Below-cloud scavenging (gases and particles) - A 10 -9.9 Below-cloud scavenging (gases and particles)- B 11 -9.9 In-cloud scavenging (gases and particles) - Ai 12 -9.9 In-cloud scavenging (gases and particles) - Bi 13 -9.9 Dry deposition (gases) - D 14 Dry deposition (gases) - Henrys const. 15 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 44.000 molweight 21 15.5E-18 OH Reaction rate - C [cm^3/molecule/sec] 22 61.0 OH reaction rate - D [K] 23 2.0 OH Reaction rate - N (no unit) 24 -9 number of associated species (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="C3H8 ", 3 PWEIGHTMOLAR= 44.0000000 , 4 POHCCONST= 1.55-E17, 5 POHDCONST= 61.0000000 , 6 POHNCONST= 2.00000000 , 7 / -
options/SPECIES/SPECIES_031
rd4ba528 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * * 6 **************************************************************************** 7 PCB28 Tracer name 8 -999.9 Species half life 9 3.5E-05 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 4.3 Dry deposition (gases) - D 14 3.1e-00 Dry deposition (gases) - Henrys const. 15 0.1 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 257.54 molweight 21 1.07E-11 OH Reaction rate - C [cm^3/molecule/sec] 22 1203.0 OH reaction rate - D [K] 23 2.0 OH Reaction rate - N (no unit) 24 -9 number of associated specias (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="PCB28 ", 3 PWETA_GAS= 3.5E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 4.3 , 6 PHENRY= 3.1 , 7 PF0= 0.1 , 8 PWEIGHTMOLAR= 257.54 , 9 POHCCONST= 1.07E-11, 10 POHDCONST= 1203.00000 , 11 POHNCONST= 2.00000000 , 12 / -
options/SPECIES/SPECIES_034
rd4ba528 r341f4b7 1 **************************************************************************** 2 * * 3 * Input file for the Lagrangian particle dispersion model FLEXPART * 4 * Definition file of chemical species/radionuclides * 5 * * 6 **************************************************************************** 7 G-HCH Tracer name 8 -999.9 Species half life 9 3.5E-05 Below-cloud scavenging (gases and particles) - A 10 0.62 Below-cloud scavenging (gases and particles)- B 11 1.0 In-cloud scavenging (gases and particles) - Ai 12 1.0 In-cloud scavenging (gases and particles) - Bi 13 4.3 Dry deposition (gases) - D 14 4.22E+02 Dry deposition (gases) - Henrys const. 15 0.1 Dry deposition (gases) - f0 (reactivity) 16 -9.9E09 Dry deposition (particles) - rho 17 Dry deposition (particles) - dquer 18 Dry deposition (particles) - dsig 19 -9.99 Alternative: dry deposition velocity 20 290.83 molweight 21 6.21E-11 OH Reaction rate - C [cm^3/molecule/sec] 22 1203.0 OH reaction rate - D [K] 23 2.0 OH Reaction rate - N (no unit) 24 -9 number of associated specias (neg. none) 25 -99.99 KOA - organic matter air partitioning 1 &SPECIES_PARAMS 2 PSPECIES="G-HCH ", 3 PWETA_GAS= 3.5E-05, 4 PWETB_GAS= 0.62 , 5 PRELDIFF= 4.3 , 6 PHENRY= 422.0 , 7 PF0= 0.10 , 8 PWEIGHTMOLAR= 290.83 , 9 POHCCONST= 6.21E-11, 10 POHDCONST= 1203.00000 , 11 POHNCONST= 2.00000000 , 12 / -
src/com_mod.f90
rdb712a8 r341f4b7 160 160 !real xmass(maxpoint,maxspec) 161 161 real :: decay(maxspec) 162 real :: weta(maxspec),wetb(maxspec) 162 real :: weta_gas(maxspec),wetb_gas(maxspec) 163 real :: crain_aero(maxspec),csnow_aero(maxspec) 163 164 ! NIK: 31.01.2013- parameters for in-cloud scavening 164 real :: weta_in(maxspec), wetb_in(maxspec)165 real :: ccn_aero(maxspec),in_aero(maxspec) 165 166 real :: reldiff(maxspec),henry(maxspec),f0(maxspec) 166 167 real :: density(maxspec),dquer(maxspec),dsigma(maxspec) … … 197 198 198 199 ! WET DEPOSITION 199 ! weta , wetb parameters for determining below-cloud wet scavenging coefficients200 ! weta_in, wetb_in parameters for determining in-cloud wet scavenging coefficients201 ! wetc_in, wetd_in parameters for determining in-cloud wet scavenging coefficients200 ! weta_gas, wetb_gas parameters for below-cloud wet scavenging coefficients (gasses) 201 ! crain_aero, in_aero parameters for in-cloud wet scavenging coefficients (aerosols) 202 ! wetc_in, wetd_in parameters for in-cloud wet scavenging coefficients (aerosols) 202 203 203 204 ! GAS DEPOSITION … … 216 217 ! dquer [m] mean diameter of particles 217 218 ! dsigma dsigma=10 or dsigma=0.1 means that 68% of the 218 ! mass are between 0.1*dquer and 10*dquer219 ! mass are between 0.1*dquer and 10*dquer 219 220 220 221 ! fract mass fraction of each diameter interval … … 371 372 372 373 !ZHG Sep 2015 373 real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem)374 real :: c lw4(0:nxmax-1,0:nymax-1,numwfmem) ! eso: =icloud_stats(:,:,4,:)374 ! real :: icloud_stats(0:nxmax-1,0:nymax-1,5,numwfmem) 375 real :: ctwc(0:nxmax-1,0:nymax-1,numwfmem) ! eso: =icloud_stats(:,:,4,:) 375 376 376 377 … … 392 393 393 394 ! pplev for the GFS version 395 ! ctwc total cloud water content 394 396 395 397 ! 2d fields … … 489 491 real,allocatable,dimension(:,:,:,:,:) :: uun, vvn, wwn, ttn, qvn, pvn,& 490 492 & rhon, drhodzn, tthn, qvhn, clwcn, ciwcn, clwn, clwchn, ciwchn 491 real,allocatable,dimension(:,:,:,:) :: c lw4n493 real,allocatable,dimension(:,:,:,:) :: ctwcn 492 494 integer,allocatable,dimension(:,:,:,:) :: cloudshn 493 495 integer(kind=1),allocatable,dimension(:,:,:,:,:) :: cloudsn … … 802 804 allocate(clwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests)) 803 805 allocate(ciwchn(0:nxmaxn-1,0:nymaxn-1,nuvzmax,numwfmem,numbnests)) 804 allocate(clw4n(0:nxmax-1,0:nymax-1,numwfmem,numbnests)) 805 806 ! clw4n(:,:,:,:)=0. 806 allocate(ctwcn(0:nxmax-1,0:nymax-1,numwfmem,numbnests)) 807 807 808 clwcn(:,:,:,:,:)=0. 808 809 ciwcn(:,:,:,:,:)=0. 809 810 clwchn(:,:,:,:,:)=0. 810 811 ciwchn(:,:,:,:,:)=0. 811 ! clwn(:,:,:,:,:)=0.812 813 ! cloudsn(:,:,:,:,:)=0814 ! cloudshn(:,:,:,:)=0815 812 816 813 end subroutine com_mod_allocate_nests -
src/makefile
r32b49c3 r341f4b7 66 66 67 67 ## LIBRARIES 68 LIBS = -lgrib_api_f90 -lgrib_api -lm -ljasper -lnetcdff -llapack# -fopenmp68 LIBS = -lgrib_api_f90 -lgrib_api -lm -ljasper -lnetcdff # -fopenmp 69 69 70 70 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 -
src/mpi_mod.f90
r32b49c3 r341f4b7 980 980 ! call MPI_Bcast(icloud_stats(:,:,:,li:ui),d2s1*5,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr) 981 981 ! if (mp_ierr /= 0) goto 600 982 call MPI_Bcast(c lw4(:,:,li:ui),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)982 call MPI_Bcast(ctwc(:,:,li:ui),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr) 983 983 if (mp_ierr /= 0) goto 600 984 984 ! call MPI_Bcast(clwc(:,:,:,li:ui),d3s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr) … … 1151 1151 ! call MPI_Bcast(icloud_stats(:,:,:,li:ui),d2s1*5,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr) 1152 1152 ! if (mp_ierr /= 0) goto 600 1153 call MPI_Bcast(c lw4n(:,:,li:ui,i),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr)1153 call MPI_Bcast(ctwcn(:,:,li:ui,i),d2s1,mp_sp,id_read,MPI_COMM_WORLD,mp_ierr) 1154 1154 if (mp_ierr /= 0) goto 600 1155 1155 end if … … 1360 1360 ! call MPI_Isend(icloud_stats(:,:,:,mind),d2s1*5,mp_sp,dest,tm1,& 1361 1361 ! &MPI_COMM_WORLD,reqs(i),mp_ierr) 1362 call MPI_Isend(c lw4(:,:,mind),d2s1,mp_sp,dest,tm1,&1362 call MPI_Isend(ctwc(:,:,mind),d2s1,mp_sp,dest,tm1,& 1363 1363 &MPI_COMM_WORLD,reqs(i),mp_ierr) 1364 1364 … … 1573 1573 ! call MPI_Irecv(icloud_stats(:,:,:,mind),d2s1*5,mp_sp,id_read,MPI_ANY_TAG,& 1574 1574 ! &MPI_COMM_WORLD,reqs(j),mp_ierr) 1575 call MPI_Irecv(c lw4(:,:,mind),d2s1*5,mp_sp,id_read,MPI_ANY_TAG,&1575 call MPI_Irecv(ctwc(:,:,mind),d2s1*5,mp_sp,id_read,MPI_ANY_TAG,& 1576 1576 &MPI_COMM_WORLD,reqs(j),mp_ierr) 1577 1577 if (mp_ierr /= 0) goto 600 -
src/netcdf_output_mod.f90
rdb712a8 r341f4b7 61 61 dx,xlon0,dy,ylat0,compoint,method,lsubgrid,lconvection,& 62 62 ind_source,ind_receptor,nageclass,lage,& 63 drydep,wetdep,decay,weta ,wetb, numbnests, &64 weta_in,wetb_in, & ! wetc_in,wetd_in, &63 drydep,wetdep,decay,weta_gas,wetb_gas, numbnests, & 64 ccn_aero,in_aero, & ! wetc_in,wetd_in, & 65 65 reldiff,henry,f0,density,dquer,dsigma,dryvel,& 66 66 ! weightmolar,ohreact,spec_ass,kao,vsetaver,& … … 543 543 chunksizes = dep_chunksizes)) 544 544 call nf90_err(nf90_put_att(ncid, wdsID, 'units', '1e-12 kg m-2')) 545 call nf90_err(nf90_put_att(ncid, wdsID, 'weta ', weta(i)))546 call nf90_err(nf90_put_att(ncid, wdsID, 'wetb ', wetb(i)))547 call nf90_err(nf90_put_att(ncid, wdsID, ' weta_in', weta_in(i)))548 call nf90_err(nf90_put_att(ncid, wdsID, ' wetb_in', wetb_in(i)))545 call nf90_err(nf90_put_att(ncid, wdsID, 'weta_gas', weta_gas(i))) 546 call nf90_err(nf90_put_att(ncid, wdsID, 'wetb_gas', wetb_gas(i))) 547 call nf90_err(nf90_put_att(ncid, wdsID, 'ccn_aero', ccn_aero(i))) 548 call nf90_err(nf90_put_att(ncid, wdsID, 'in_aero', in_aero(i))) 549 549 ! call nf90_err(nf90_put_att(ncid, wdsID, 'wetc_in', wetc_in(i))) 550 550 ! call nf90_err(nf90_put_att(ncid, wdsID, 'wetd_in', wetd_in(i))) -
src/par_mod.f90
r32b49c3 r341f4b7 89 89 real,parameter :: rho_water=1000. !ZHG 2015 [kg/m3] 90 90 !ZHG MAR2016 91 real,parameter :: incloud_ratio= 1.691 real,parameter :: incloud_ratio=2.2 92 92 93 93 ! karman Karman's constant -
src/readreleases.f90
rb0434e1 r341f4b7 57 57 ! xpoint2,ypoint2 geograf. coordinates of upper right corner of release * 58 58 ! area * 59 ! weta , wetb parameters to determine the below-cloud scavenging*60 ! weta_in, wetb_in parameters to determine the in-cloud scavenging*61 ! wetc_in, wetd_in parameters to determine the in-cloud scavenging*59 ! weta_gas, wetb_gas parameters for below-cloud scavenging (gas) * 60 ! crain_aero, csnow_aero parameters for below-cloud scavenging (aerosol) * 61 ! ccn_aero, in_aero parameters for in-cloud scavenging (aerosol) * 62 62 ! zpoint1,zpoint2 height range, over which release takes place * 63 63 ! num_min_discrete if less, release cannot be randomized and happens at * … … 363 363 !*********************************************************** 364 364 ! 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 365 if (weta (i).gt.0. .or. wetb(i).gt.0.) then365 if (weta_gas(i).gt.0. .or. wetb_gas(i).gt.0.) then ! :TODO: or particle 366 366 WETDEP=.true. 367 367 if (lroot) then … … 374 374 375 375 ! NIK 31.01.2013 + 10.12.2013 + 15.02.2015 376 if ( weta_in(i).gt.0. .or. wetb_in(i).gt.0.) then376 if (ccn_aero(i).gt.0. .or. in_aero(i).gt.0.) then 377 377 WETDEP=.true. 378 378 if (lroot) then 379 379 write (*,*) 'In-cloud scavenging: ON' 380 380 ! write (*,*) 'In-cloud scavenging coefficients: ',& 381 ! & weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i381 ! &ccn_aero(i),in_aero(i),i !,wetc_in(i), wetd_in(i),i 382 382 end if 383 383 else -
src/readspecies.f90
r32b49c3 r341f4b7 42 42 ! decaytime(maxtable) half time for radiological decay * 43 43 ! specname(maxtable) names of chemical species, radionuclides * 44 ! weta , wetbParameters for determining below-cloud scavenging *45 ! weta_inParameter for determining in-cloud scavenging *46 ! wetb_inParameter for determining in-cloud scavenging *44 ! weta_gas, wetb_gas Parameters for determining below-cloud scavenging * 45 ! ccn_aero Parameter for determining in-cloud scavenging * 46 ! in_aero Parameter for determining in-cloud scavenging * 47 47 ! ohcconst OH reaction rate constant C * 48 48 ! ohdconst OH reaction rate constant D * … … 66 66 67 67 character(len=16) :: pspecies 68 real :: pdecay, pweta , pwetb, preldiff, phenry, pf0, pdensity, pdquer68 real :: pdecay, pweta_gas, pwetb_gas, preldiff, phenry, pf0, pdensity, pdquer 69 69 real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao 70 real :: p weta_in, pwetb_in, pwetc_in, pwetd_in70 real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero 71 71 integer :: readerror 72 72 73 73 ! declare namelist 74 74 namelist /species_params/ & 75 pspecies, pdecay, pweta , pwetb, &76 p weta_in, pwetb_in, pwetc_in, pwetd_in, &75 pspecies, pdecay, pweta_gas, pwetb_gas, & 76 pcrain_aero, pcsnow_aero, pccn_aero, pin_aero, & 77 77 preldiff, phenry, pf0, pdensity, pdquer, & 78 78 pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao 79 79 80 pspecies=" "80 pspecies="" ! read failure indicator value 81 81 pdecay=-999.9 82 pweta =-9.9E-0983 pwetb =0.084 p weta_in=-9.9E-0985 p wetb_in=-9.9E-0986 ! pwetc_in=-9.9E-0987 ! pwetd_in=-9.9E-0982 pweta_gas=-9.9E-09 83 pwetb_gas=0.0 84 pcrain_aero=-9.9E-09 85 pcsnow_aero=-9.9E-09 86 pccn_aero=-9.9E-09 87 pin_aero=-9.9E-09 88 88 preldiff=-9.9 89 89 phenry=0.0 … … 98 98 pspec_ass=-9 99 99 pkao=-99.99 100 pweightmolar=- 789.0 ! read failure indicator value100 pweightmolar=-999.9 101 101 102 102 ! Open the SPECIES file and read species names and properties … … 113 113 close(unitspecies) 114 114 115 if ((pweightmolar.eq.-789.0).or.(readerror.ne.0)) then ! no namelist found 115 ! if ((pweightmolar.eq.-789.0).or.(readerror.ne.0)) then ! no namelist found 116 if ((len(pspecies).eq.0).or.(readerror.ne.0)) then ! no namelist found 116 117 117 118 readerror=1 … … 127 128 read(unitspecies,'(f18.1)',end=22) decay(pos_spec) 128 129 ! write(*,*) decay(pos_spec) 129 read(unitspecies,'(e18.1)',end=22) weta(pos_spec) 130 ! write(*,*) weta(pos_spec) 131 read(unitspecies,'(f18.2)',end=22) wetb(pos_spec) 132 ! write(*,*) wetb(pos_spec) 133 130 read(unitspecies,'(e18.1)',end=22) weta_gas(pos_spec) 131 ! write(*,*) weta_gas(pos_spec) 132 read(unitspecies,'(f18.2)',end=22) wetb_gas(pos_spec) 133 ! write(*,*) wetb_gas(pos_spec) 134 read(unitspecies,'(e18.1)',end=22) crain_aero(pos_spec) 135 ! write(*,*) crain_aero(pos_spec) 136 read(unitspecies,'(f18.2)',end=22) csnow_aero(pos_spec) 137 ! write(*,*) csnow_aero(pos_spec) 134 138 !*** NIK 31.01.2013: including in-cloud scavening parameters 135 read(unitspecies,'(e18.1)',end=22) weta_in(pos_spec) 136 ! write(*,*) weta_in(pos_spec) 137 read(unitspecies,'(f18.2)',end=22) wetb_in(pos_spec) 138 ! write(*,*) wetb_in(pos_spec) 139 ! read(unitspecies,'(f18.2)',end=22) wetc_in(pos_spec) 140 ! write(*,*) wetc_in(pos_spec) 141 ! read(unitspecies,'(f18.2)',end=22) wetd_in(pos_spec) 142 ! write(*,*) wetd_in(pos_spec) 143 139 read(unitspecies,'(e18.1)',end=22) ccn_aero(pos_spec) 140 ! write(*,*) ccn_aero(pos_spec) 141 read(unitspecies,'(f18.2)',end=22) in_aero(pos_spec) 142 ! write(*,*) in_aero(pos_spec) 144 143 read(unitspecies,'(f18.1)',end=22) reldiff(pos_spec) 145 144 ! write(*,*) reldiff(pos_spec) … … 171 170 pspecies=species(pos_spec) 172 171 pdecay=decay(pos_spec) 173 pweta =weta(pos_spec)174 pwetb =wetb(pos_spec)175 p weta_in=weta_in(pos_spec)176 p wetb_in=wetb_in(pos_spec)177 ! pwetc_in=wetc_in(pos_spec)178 ! pwetd_in=wetd_in(pos_spec)172 pweta_gas=weta_gas(pos_spec) 173 pwetb_gas=wetb_gas(pos_spec) 174 pcrain_aero=crain_aero(pos_spec) 175 pcsnow_aero=csnow_aero(pos_spec) 176 pccn_aero=ccn_aero(pos_spec) 177 pin_aero=in_aero(pos_spec) 179 178 preldiff=reldiff(pos_spec) 180 179 phenry=henry(pos_spec) … … 195 194 species(pos_spec)=pspecies 196 195 decay(pos_spec)=pdecay 197 weta (pos_spec)=pweta198 wetb (pos_spec)=pwetb199 weta_in(pos_spec)=pweta_in200 wetb_in(pos_spec)=pwetb_in201 ! wetc_in(pos_spec)=pwetc_in 202 ! wetd_in(pos_spec)=pwetd_in 196 weta_gas(pos_spec)=pweta_gas 197 wetb_gas(pos_spec)=pwetb_gas 198 crain_aero=pcrain_aero 199 csnow_aero=pcsnow_aero 200 ccn_aero(pos_spec)=pccn_aero 201 in_aero(pos_spec)=pin_aero 203 202 reldiff(pos_spec)=preldiff 204 203 henry(pos_spec)=phenry … … 222 221 ! Check scavenging parameters given in SPECIES file 223 222 224 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) then223 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 225 224 226 225 if (dquer(pos_spec).gt.0) then !is particle 227 226 if (lroot) then 228 227 write(*,'(a,f5.2)') ' Particle below-cloud scavenging parameter A & 229 &(Rain collection efficiency) ', weta(pos_spec)228 &(Rain collection efficiency) ', crain_aero(pos_spec) 230 229 write(*,'(a,f5.2)') ' Particle below-cloud scavenging parameter B & 231 &(Snow collection efficiency) ', wetb(pos_spec)232 end if 233 if (weta (pos_spec).gt.1.0 .or. wetb(pos_spec).gt.1.0) then230 &(Snow collection efficiency) ', csnow_aero(pos_spec) 231 end if 232 if (weta_gas(pos_spec).gt.1.0 .or. wetb_gas(pos_spec).gt.1.0) then 234 233 if (lroot) then 235 234 write(*,*) '*******************************************' … … 242 241 if (lroot) then 243 242 write(*,'(a,f5.2)') ' Particle in-cloud scavenging parameter Ai & 244 &(CCN efficiency) ', weta_in(pos_spec)243 &(CCN efficiency) ', ccn_aero(pos_spec) 245 244 write(*,'(a,f5.2)') ' Particle in-cloud scavenging parameter Bi & 246 &(IN efficiency) ', wetb_in(pos_spec)247 end if 248 if ( weta_in(pos_spec).gt.1.0 .or. weta_in(pos_spec).lt.0.01) then245 &(IN efficiency) ', in_aero(pos_spec) 246 end if 247 if (ccn_aero(pos_spec).gt.1.0 .or. ccn_aero(pos_spec).lt.0.01) then 249 248 if (lroot) then 250 249 write(*,*) '*******************************************' … … 254 253 end if 255 254 endif 256 if ( wetb_in(pos_spec).gt.1.0 .or. wetb_in(pos_spec).lt.0.01) then255 if (in_aero(pos_spec).gt.1.0 .or. in_aero(pos_spec).lt.0.01) then 257 256 if (lroot) then 258 257 write(*,*) '*******************************************' … … 265 264 else !is gas 266 265 if (lroot) then 267 write(*,*) ' Gas below-cloud scavenging parameter A ', weta (pos_spec)268 write(*,'(a,f5.2)') ' Gas below-cloud scavenging parameter B ', wetb (pos_spec)266 write(*,*) ' Gas below-cloud scavenging parameter A ', weta_gas(pos_spec) 267 write(*,'(a,f5.2)') ' Gas below-cloud scavenging parameter B ', wetb_gas(pos_spec) 269 268 write(*,*) ' Gas in-cloud scavenging uses default values as in Hertel et al 1995' 270 269 end if 271 if (weta (pos_spec).gt.0.) then !if wet deposition is turned on272 if (weta (pos_spec).gt.1E-04 .or. weta(pos_spec).lt.1E-09) then270 if (weta_gas(pos_spec).gt.0.) then !if wet deposition is turned on 271 if (weta_gas(pos_spec).gt.1E-04 .or. weta_gas(pos_spec).lt.1E-09) then 273 272 if (lroot) then 274 273 write(*,*) '*******************************************' … … 279 278 endif 280 279 end if 281 if (wetb (pos_spec).gt.0.) then !if wet deposition is turned on282 if (wetb (pos_spec).gt.0.8 .or. wetb(pos_spec).lt.0.6) then280 if (wetb_gas(pos_spec).gt.0.) then !if wet deposition is turned on 281 if (wetb_gas(pos_spec).gt.0.8 .or. wetb_gas(pos_spec).lt.0.6) then 283 282 if (lroot) then 284 283 write(*,*) '*******************************************' … … 292 291 endif 293 292 294 if (((weta (pos_spec).gt.0).or.(wetb(pos_spec).gt.0)).and.(henry(pos_spec).le.0)) then293 if (((weta_gas(pos_spec).gt.0).or.(wetb_gas(pos_spec).gt.0)).and.(henry(pos_spec).le.0)) then 295 294 if (dquer(pos_spec).le.0) goto 996 ! no particle, no henry set 296 295 endif -
src/timemanager_mpi.f90
r32b49c3 r341f4b7 169 169 if (mp_dev_mode) write(*,*) 'myid, itime: ',mp_pid,itime 170 170 171 if ( wetdep.and. itime .ne. 0 .and. numpart .gt. 0) then171 if (WETDEP .and. itime .ne. 0 .and. numpart .gt. 0) then 172 172 if (verbosity.gt.0) then 173 173 write (*,*) 'timemanager> call wetdepo' -
src/verttransform.f90
rdb712a8 r341f4b7 82 82 83 83 integer :: ix,jy,kz,iz,n,kmin,ix1,jy1,ixp,jyp,ixm,jym,kz_inv 84 real :: f_qvsat,pressure,rh,lsp,convp, prec84 real :: f_qvsat,pressure,rh,lsp,convp,cloudh_min,prec 85 85 real :: ew,dz1,dz2,dz 86 86 real :: xlon,ylat,xlonr,dzdx,dzdy … … 586 586 write(*,*) 'Global ECMWF fields: using cloud water' 587 587 clw(:,:,:,n)=0.0 588 icloud_stats(:,:,:,n)=0.0 588 ! icloud_stats(:,:,:,n)=0.0 589 ctwc(:,:,n)=0.0 589 590 clouds(:,:,:,n)=0 590 591 ! If water/ice are read separately into clwc and ciwc, store sum in clwc … … 604 605 clw(ix,jy,kz,n)=(clwc(ix,jy,kz,n)*rho(ix,jy,kz,n))*(height(kz+1)-height(kz)) 605 606 tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 606 icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n) ! Column cloud water [m3/m3] 607 icloud_stats(ix,jy,3,n)= min(height(kz+1),height(kz)) ! Cloud BOT height stats [m] 607 608 ! icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n) ! Column cloud water [m3/m3] 609 ctwc(ix,jy,n) = ctwc(ix,jy,n)+clw(ix,jy,kz,n) 610 ! icloud_stats(ix,jy,3,n)= min(height(kz+1),height(kz)) ! Cloud BOT height stats [m] 611 cloudh_min=min(height(kz+1),height(kz)) 608 612 !ZHG 2015 extra for testing 609 613 ! clh(ix,jy,kz,n)=height(kz+1)-height(kz) 610 icloud_stats(ix,jy,1,n)=icloud_stats(ix,jy,1,n)+(height(kz+1)-height(kz)) ! Cloud total vertical extent [m]611 icloud_stats(ix,jy,2,n)= max(icloud_stats(ix,jy,2,n),height(kz)) ! Cloud TOP height [m]614 ! icloud_stats(ix,jy,1,n)=icloud_stats(ix,jy,1,n)+(height(kz+1)-height(kz)) ! Cloud total vertical extent [m] 615 ! icloud_stats(ix,jy,2,n)= max(icloud_stats(ix,jy,2,n),height(kz)) ! Cloud TOP height [m] 612 616 !ZHG 613 617 endif … … 626 630 clouds(ix,jy,kz,n)=2 ! convp in-cloud 627 631 endif ! convective or large scale 628 elseif((clw(ix,jy,kz,n).le.0) .and. ( icloud_stats(ix,jy,3,n).ge.height(kz))) then ! is below cloud632 elseif((clw(ix,jy,kz,n).le.0) .and. (cloudh_min.ge.height(kz))) then ! is below cloud 629 633 if (lsp.ge.convp) then 630 634 clouds(ix,jy,kz,n)=5 ! lsp dominated washout … … 643 647 644 648 ! eso: copy the relevant data to clw4 to reduce amount of communicated data for MPI 645 clw4(:,:,n) = icloud_stats(:,:,4,n)649 ! ctwc(:,:,n) = icloud_stats(:,:,4,n) 646 650 647 651 !************************************************************************** -
src/verttransform_nests.f90
rdb712a8 r341f4b7 466 466 clwn(0:nxm1,0:nym1,:,n,l)=0.0 467 467 ! icloud_stats(0:nxm1,0:nym1,:,n)=0.0 468 c lw4n(0:nxm1,0:nym1,n,l)=0.0468 ctwcn(0:nxm1,0:nym1,n,l)=0.0 469 469 cloudsn(0:nxm1,0:nym1,:,n,l)=0 470 470 ! If water/ice are read separately into clwc and ciwc, store sum in clwcn … … 484 484 clwn(ix,jy,kz,n,l)=(clwcn(ix,jy,kz,n,l)*rhon(ix,jy,kz,n,l))*(height(kz+1)-height(kz)) 485 485 tot_cloud_h=tot_cloud_h+(height(kz+1)-height(kz)) 486 c lw4n(ix,jy,n,l) = clw4n(ix,jy,n,l)+clwn(ix,jy,kz,n,l)486 ctwcn(ix,jy,n,l) = ctwcn(ix,jy,n,l)+clwn(ix,jy,kz,n,l) 487 487 ! icloud_stats(ix,jy,4,n)= icloud_stats(ix,jy,4,n)+clw(ix,jy,kz,n) ! Column cloud water [m3/m3] 488 488 ! icloud_stats(ix,jy,3,n)= min(height(kz+1),height(kz)) ! Cloud BOT height stats [m] -
src/wetdepo.f90
rec7fc72 r341f4b7 270 270 if (clouds_v.ge.4) then !below cloud 271 271 272 if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then !if positive below-cloud parameters given in SPECIES file (either A or B) 272 ! For gas: if positive below-cloud parameters (A or B), and dquer<=0 273 !****************************************************************** 274 if ((dquer(ks).le.0.).and.(weta_gas(ks).gt.0..or.wetb_gas(ks).gt.0.)) then 275 ! if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then 273 276 blc_count=blc_count+1 274 275 !GAS 276 if (dquer(ks) .le. 0.) then !is gas 277 wetscav=weta(ks)*prec(1)**wetb(ks) 278 279 !AEROSOL 280 else !is particle 277 wetscav=weta_gas(ks)*prec(1)**wetb_gas(ks) 278 279 ! For aerosols: if positive below-cloud parameters (Crain/Csnow or B), and dquer>0 280 !********************************************************************************* 281 else if ((dquer(ks).gt.0.).and.(crain_aero(ks).gt.0..or.csnow_aero(ks).gt.0.)) then 282 blc_count=blc_count+1 283 281 284 !NIK 17.02.2015 282 ! For the calculation here particle size needs to be in meter and not um as dquer is changed to in readreleases 283 ! for particles larger than 10 um use the largest size defined in the parameterizations (10um) 284 dquer_m=min(10.,dquer(ks))/1000000. !conversion from um to m 285 if (act_temp .ge. 273. .and. weta(ks).gt.0.) then !Rain 285 ! For the calculation here particle size needs to be in meter and not um as dquer is 286 ! changed in readreleases 287 ! For particles larger than 10 um use the largest size defined in the parameterizations (10um) 288 dquer_m=min(10.,dquer(ks))/1000000. !conversion from um to m 289 290 ! Rain: 291 if (act_temp .ge. 273. .and. crain_aero(ks).gt.0.) then 292 286 293 ! ZHG 2014 : Particle RAIN scavenging coefficient based on Laakso et al 2003, 287 ! the below-cloud scavenging (rain efficienty) 288 ! parameter A (=weta) from SPECIES file 289 wetscav= weta(ks)*10**(bclr(1)+ (bclr(2)*(log10(dquer_m))**(-4))+(bclr(3)*(log10(dquer_m))**(-3))+ (bclr(4)* & 290 (log10(dquer_m))**(-2))+ (bclr(5)*(log10(dquer_m))**(-1))+ bclr(6)* (prec(1))**(0.5)) 291 292 elseif (act_temp .lt. 273. .and. wetb(ks).gt.0.) then ! Snow 294 ! the below-cloud scavenging (rain efficienty) parameter Crain (=crain_aero) from SPECIES file 295 wetscav=crain_aero(ks)*10**(bclr(1)+(bclr(2)*(log10(dquer_m))**(-4))+ & 296 & (bclr(3)*(log10(dquer_m))**(-3))+ (bclr(4)*(log10(dquer_m))**(-2))+& 297 &(bclr(5)*(log10(dquer_m))**(-1))+bclr(6)* (prec(1))**(0.5)) 298 299 ! Snow: 300 elseif (act_temp .lt. 273. .and. csnow_aero(ks).gt.0.) then 293 301 ! ZHG 2014 : Particle SNOW scavenging coefficient based on Kyro et al 2009, 294 ! the below-cloud scavenging (Snow efficiency) 295 ! parameter B (=wetb) from SPECIES file 296 wetscav= wetb(ks)*10**(bcls(1)+ (bcls(2)*(log10(dquer_m))**(-4))+(bcls(3)*(log10(dquer_m))**(-3))+ (bcls(4)*&297 (log10(dquer_m))**(-2))+(bcls(5)*(log10(dquer_m))**(-1))+ bcls(6)* (prec(1))**(0.5))298 299 300 302 ! the below-cloud scavenging (Snow efficiency) parameter Csnow (=csnow_aero) from SPECIES file 303 wetscav=csnow_aero(ks)*10**(bcls(1)+(bcls(2)*(log10(dquer_m))**(-4))+& 304 &(bcls(3)*(log10(dquer_m))**(-3))+ (bcls(4)*(log10(dquer_m))**(-2))+& 305 &(bcls(5)*(log10(dquer_m))**(-1))+ bcls(6)* (prec(1))**(0.5)) 306 307 endif 308 301 309 ! write(*,*) 'bl-cloud, act_temp=',act_temp, ',prec=',prec(1),',wetscav=', wetscav, ', jpart=',jpart 302 310 303 endif !gas or particle304 311 endif ! gas or particle 312 ! endif ! positive below-cloud scavenging parameters given in Species file 305 313 endif !end BELOW 306 307 314 308 315 !******************** … … 310 317 !******************** 311 318 if (clouds_v.lt.4) then ! In-cloud 312 ! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are given in species file 313 if (weta_in(ks).gt.0. .or. wetb_in(ks).gt.0.) then 319 ! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are 320 ! given in species file, or if gas and positive Henry's constant 321 if ((ccn_aero(ks).gt.0. .or. in_aero(ks).gt.0.).or.(henry(ks).gt.0.and.dquer(ks).le.0)) then 314 322 inc_count=inc_count+1 315 323 ! if negative coefficients (turned off) set to zero for use in equation 316 if ( weta_in(ks).lt.0.) weta_in(ks)=0.317 if ( wetb_in(ks).lt.0.) wetb_in(ks)=0.324 if (ccn_aero(ks).lt.0.) ccn_aero(ks)=0. 325 if (in_aero(ks).lt.0.) in_aero(ks)=0. 318 326 319 327 !ZHG 2015 Cloud liquid & ice water (CLWC+CIWC) from ECMWF 320 328 ! nested fields 321 329 if (ngrid.gt.0.and.readclouds_this_nest) then 322 cl=c lw4n(ix,jy,n,ngrid)*(grfraction(1)/cc)330 cl=ctwcn(ix,jy,n,ngrid)*(grfraction(1)/cc) 323 331 else if (ngrid.eq.0.and.readclouds) then 324 cl=c lw4(ix,jy,n)*(grfraction(1)/cc)332 cl=ctwc(ix,jy,n)*(grfraction(1)/cc) 325 333 else !parameterize cloudwater m2/m3 326 334 !ZHG updated parameterization of cloud water to better reproduce the values coming from ECMWF … … 337 345 end if 338 346 ! ZHG: Calculate the aerosol partition based on cloud phase and Ai and Bi 339 frac_act = liq_frac* weta_in(ks) +(1-liq_frac)*wetb_in(ks)347 frac_act = liq_frac*ccn_aero(ks) +(1-liq_frac)*in_aero(ks) 340 348 341 349 !ZHG Use the activated fraction and the liqid water to calculate the washout … … 343 351 ! AEROSOL 344 352 !************************************************** 345 if (dquer(ks).gt. 353 if (dquer(ks).gt.0.) then ! is particle 346 354 347 355 S_i= frac_act/cl
Note: See TracChangeset
for help on using the changeset viewer.
