Skip to content
Snippets Groups Projects
Commit 10d41309 authored by Lucie Bakels's avatar Lucie Bakels
Browse files

Added GeoSphere SPECIES by Marie D. Mulder and Christian Maurer

parent bb1732de
No related branches found
No related tags found
No related merge requests found
Showing
with 349 additions and 9 deletions
...@@ -119,7 +119,7 @@ And for each release: ...@@ -119,7 +119,7 @@ And for each release:
**Note:** the RELEASES file is no longer necessary when using [IPIN](configuration.md#ipin)=3, giving full control to the user to decide where and when particles of different species are being released (see [User-defined initial conditions](configuration.md#ic)). **Note:** the RELEASES file is no longer necessary when using [IPIN](configuration.md#ipin)=3, giving full control to the user to decide where and when particles of different species are being released (see [User-defined initial conditions](configuration.md#ic)).
### <a name="species"></a>SPECIES ### <a name="species"></a>SPECIES
The subdirectory options/SPECIES/ needs to contain one or more files named SPECIES_nnn. For each species nnn listed in the header section of the RELEASES file, such a SPECIES_nnn file must exist. The parameters in the SPECIES_nnn file, contained in the namelist &SPECIES_PARAMS, set the species name and define the physicochemical properties of the species; they are described in Table 10. These are important for simulating radioactive or chemical decay, wet deposition (scavenging) for gases and aerosols, dry deposition for gases and aerosols, particle settling, and chemical reaction with the OH radical. Some parameters are only necessary for gas tracers and some are only necessary for aerosol tracers; thus, a namelist does not need to contain all parameters for both gases and particles. Optionally, since FLEXPART version 6.0, information about temporal emission variations can be added at the end of the file. The subdirectory options/SPECIES/ needs to contain one or more files named SPECIES_nnn. In options/SPECIES/ templates for several species are given. These come with no warranty and will have to be renamed to SPECIES_nnn where nnn are three digits. For each species nnn listed in the header section of the RELEASES file, such a SPECIES_nnn file must exist. The parameters in the SPECIES_nnn file, contained in the namelist &SPECIES_PARAMS, set the species name and define the physicochemical properties of the species; they are described in Table 10. These are important for simulating radioactive or chemical decay, wet deposition (scavenging) for gases and aerosols, dry deposition for gases and aerosols, particle settling, and chemical reaction with the OH radical. Some parameters are only necessary for gas tracers and some are only necessary for aerosol tracers; thus, a namelist does not need to contain all parameters for both gases and particles. Optionally, since FLEXPART version 6.0, information about temporal emission variations can be added at the end of the file.
The following specifies the parameters associated with each physicochemical process simulated. The following specifies the parameters associated with each physicochemical process simulated.
......
...@@ -14,6 +14,9 @@ ...@@ -14,6 +14,9 @@
LOUTAVER= 3600, ! Interval of output averaging (s) LOUTAVER= 3600, ! Interval of output averaging (s)
LOUTSAMPLE= 900, ! Interval of output sampling (s), higher stat. accuracy with shorter intervals LOUTSAMPLE= 900, ! Interval of output sampling (s), higher stat. accuracy with shorter intervals
LOUTRESTART= 86400, ! Interval of writing restart files (s), switched off when set to -1 LOUTRESTART= 86400, ! Interval of writing restart files (s), switched off when set to -1
LRECOUTSTEP= 3600, ! Interval of model output at receptors (s)
LRECOUTAVER= 3600, ! Interval of receptor output averaging (s)
LRECOUTSAMPLE= 1200, ! Interval of receptor output sampling (s)
LSYNCTIME= 900, ! All processes are synchronized to this time interval (s) LSYNCTIME= 900, ! All processes are synchronized to this time interval (s)
CTL= -5.0000000, ! CTL>1, ABL time step = (Lagrangian timescale (TL))/CTL, uses LSYNCTIME if CTL<0 CTL= -5.0000000, ! CTL>1, ABL time step = (Lagrangian timescale (TL))/CTL, uses LSYNCTIME if CTL<0
IFINE= 4, ! Reduction for time step in vertical transport, used only if CTL>1 IFINE= 4, ! Reduction for time step in vertical transport, used only if CTL>1
...@@ -41,6 +44,5 @@ ...@@ -41,6 +44,5 @@
MAXTHREADGRID= 1, ! Set maximum number of threads for doing grid computations. Recommended to set this no higher than 16. High numbers create more overhead and a larger memory footprint, 1=no parallelisation on grid. MAXTHREADGRID= 1, ! Set maximum number of threads for doing grid computations. Recommended to set this no higher than 16. High numbers create more overhead and a larger memory footprint, 1=no parallelisation on grid.
MAXFILESIZE= 10000, ! Maximum output of each partoutput NetCDF-4 file in Mb before a new one is created MAXFILESIZE= 10000, ! Maximum output of each partoutput NetCDF-4 file in Mb before a new one is created
LOGVERTINTERP= 0, ! Flag to set all vertical interpolation to logarithmic instead of linear LOGVERTINTERP= 0, ! Flag to set all vertical interpolation to logarithmic instead of linear
BCSCHEME= 1, ! Below-cloud scheme; [1] Laakso & Kyro [2] Wang et al [3] modified Wang et al
LCMOUTPUT= 0, ! Switch for the Linear Chemistry Module; [0] off [1] on LCMOUTPUT= 0, ! Switch for the Linear Chemistry Module; [0] off [1] on
/ /
File moved
...@@ -11,6 +11,6 @@ ...@@ -11,6 +11,6 @@
!******************************************************************************* !*******************************************************************************
&REAGENT_PARAMS &REAGENT_PARAMS
PREAGENT="OH", PREAGENT="OH",
PREAG_PATH="/xnilu_wrk/users/rlt/NILU/FLEXDATA/REGAME/flexpart_ctm/oh_fields/", PREAG_PATH="./oh_fields/",
PHOURLY=1, PHOURLY=1,
/ /
&RECEPTORS &RECEPTORS
RECEPTOR="receptor 1", RECEPTOR="BRW ",
LON = 0.0000000, LAT= 71.3199997 ,
LAT = 0.0000000, LON= -156.610001 ,
ALT= 16.0000000 ,
TIME= 2459032.0000000000 ,
/ /
&RECEPTORS &RECEPTORS
RECEPTOR="receptor 2", RECEPTOR="CRV ",
LON = 1.00000000, LAT= 64.9899979 ,
LAT = 0.00000000, LON= -147.600006 ,
ALT= 32.0000000 ,
TIME= 2459032.0000000000 ,
/ /
&RECEPTORS
RECEPTOR="LEF ",
LAT= 45.9500008 ,
LON= -90.2699966 ,
ALT= 396.000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="MLO ",
LAT= 19.5400009 ,
LON= -155.580002 ,
ALT= 300.000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="SCT ",
LAT= 33.4099998 ,
LON= -81.8300018 ,
ALT= 305.000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="WGC ",
LAT= 38.2599983 ,
LON= -121.489998 ,
ALT= 484.000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="MHD ",
LAT= 53.3300018 ,
LON= -9.98999977 ,
ALT= 11.0000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="THD ",
LAT= 41.0499992 ,
LON= -124.150002 ,
ALT= 15.0000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="RPB ",
LAT= 13.1700001 ,
LON= -59.4300003 ,
ALT= 13.0000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="SMO ",
LAT= -14.2500000 ,
LON= -170.559998 ,
ALT= 35.0000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="CGO ",
LAT= -40.6800003 ,
LON= 144.690002 ,
ALT= 80.0000000 ,
TIME= 2459032.0000000000 ,
/
&RECEPTORS
RECEPTOR="BRW ",
LAT= 71.3199997 ,
LON= -156.610001 ,
ALT= 16.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="CRV ",
LAT= 64.9899979 ,
LON= -147.600006 ,
ALT= 32.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="LEF ",
LAT= 45.9500008 ,
LON= -90.2699966 ,
ALT= 396.000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="MLO ",
LAT= 19.5400009 ,
LON= -155.580002 ,
ALT= 300.000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="SCT ",
LAT= 33.4099998 ,
LON= -81.8300018 ,
ALT= 305.000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="WGC ",
LAT= 38.2599983 ,
LON= -121.489998 ,
ALT= 484.000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="MHD ",
LAT= 53.3300018 ,
LON= -9.98999977 ,
ALT= 11.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="THD ",
LAT= 41.0499992 ,
LON= -124.150002 ,
ALT= 15.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="RPB ",
LAT= 13.1700001 ,
LON= -59.4300003 ,
ALT= 13.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="SMO ",
LAT= -14.2500000 ,
LON= -170.559998 ,
ALT= 35.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="CGO ",
LAT= -40.6800003 ,
LON= 144.690002 ,
ALT= 80.0000000 ,
TIME= 2459032.0416666665 ,
/
&RECEPTORS
RECEPTOR="BRW ",
LAT= 71.3199997 ,
LON= -156.610001 ,
ALT= 16.0000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="CRV ",
LAT= 64.9899979 ,
LON= -147.600006 ,
ALT= 32.0000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="LEF ",
LAT= 45.9500008 ,
LON= -90.2699966 ,
ALT= 396.000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="MLO ",
LAT= 19.5400009 ,
LON= -155.580002 ,
ALT= 300.000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="SCT ",
LAT= 33.4099998 ,
LON= -81.8300018 ,
ALT= 305.000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="WGC ",
LAT= 38.2599983 ,
LON= -121.489998 ,
ALT= 484.000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="MHD ",
LAT= 53.3300018 ,
LON= -9.98999977 ,
ALT= 11.0000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="THD ",
LAT= 41.0499992 ,
LON= -124.150002 ,
ALT= 15.0000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="RPB ",
LAT= 13.1700001 ,
LON= -59.4300003 ,
ALT= 13.0000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="SMO ",
LAT= -14.2500000 ,
LON= -170.559998 ,
ALT= 35.0000000 ,
TIME= 2459032.0833333335 ,
/
&RECEPTORS
RECEPTOR="CGO ",
LAT= -40.6800003 ,
LON= 144.690002 ,
ALT= 80.0000000 ,
TIME= 2459032.0833333335 ,
/
File moved
&SPECIES_PARAMS
PSPECIES="Ba-140", ! Tracer name: BARIUM 140
PDECAY=1101721.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [Van Leuven et al. (2023)]
PDIA=0.65E-06, ! Dry deposition (particles) - dquer [Van Leuven et al. (2023)]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [Van Leuven et al. (2023)]
/
&SPECIES_PARAMS
PSPECIES="CH-3-I-131", ! Tracer name: METHYL IODIDE
PDECAY=693382.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
PWETA_GAS=8.0E-07, ! Below-cloud scavenging (gases) - A (weta_gas) [Asman (1995), Paesler-Sauer (2003)]
PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) [Asman (1995)]
PDRYVEL=0.1E-3, ! Alternative: dry deposition velocity [Mueller and Proehl (1993), average of grass and soil]
/
&SPECIES_PARAMS
PSPECIES="Cs-137", ! Tracer name: CESIUM 137
PDECAY=-9.9, ! Species half life (948886704 sec) [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [Van Leuven et al. (2023)]
PDIA=0.65E-06, ! Dry deposition (particles) - dquer [Van Leuven et al. (2023)]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [Van Leuven et al. (2023)]
PWEIGHTMOLAR=-9.9, ! molweight
/
&SPECIES_PARAMS
PSPECIES="I-131", ! Tracer name: PARTICLE-BOUND IODINE 131
PDECAY=693382.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [Van Leuven et al. (2023)]
PDIA=5.0E-07, ! Dry deposition (particles) - dquer [Baklanov and Sorensen (2001)]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [Baklanov and Sorensen (2001)]
/
&SPECIES_PARAMS
PSPECIES="I2-131", ! Tracer name: ELEMENTAL IODINE 131
PDECAY=693382.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
PWETA_GAS=8.0E-05, ! Below-cloud scavenging (gases) - A (weta_gas) [Asman (1995), Paesler-Sauer (2003)]
PWETB_GAS=0.62, ! Below-cloud scavenging (gases) - B (wetb_gas) [Asman (1995)]
PRELDIFF=2.7, ! Dry deposition (gases) - D
PHENRY=100000.0, ! Dry deposition (gases) - Henrys const.
PF0=0.1, ! Dry deposition (gases) - f0 (reactivity)
PWEIGHTMOLAR=-9.9, ! molweight
/
&SPECIES_PARAMS
PSPECIES="Kr-85", ! Tracer name: KRYPTON 85
PDECAY=-9.9, ! Species half life (338669519 sec) [https://www.internetchemie.info/chemische-elemente/isotope.php]
PWEIGHTMOLAR=-9.9, ! molweight
/
&SPECIES_PARAMS
PSPECIES="Particulate-Radionuclide", ! Tracer name: e.g. Y-92, Ru-106, I-131, Cs-137
PDECAY=-9.9, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [arbitrary value]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [arbitrary value]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [arbitrary value]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [arbitrary value]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [arbitrary value]
PDIA=0.65E-06, ! Dry deposition (particles) - dquer [arbitrary value]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [arbitrary value]
/
&SPECIES_PARAMS
PSPECIES="Ru-106", ! Tracer name: RUTHENIUM 106
PDECAY=-9.9, ! Species half life (32125075 sec) [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [Van Leuven et al. (2023)]
PDIA=0.65E-06, ! Dry deposition (particles) - dquer [Van Leuven et al. (2023)]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [Van Leuven et al. (2023)]
/
&SPECIES_PARAMS
PSPECIES="Sr-90", ! Tracer name: STRONTIUM 90
PDECAY=-9.9, ! Species half life (911800368 sec) [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [Van Leuven et al. (2023)]
PDIA=0.65E-06, ! Dry deposition (particles) - dquer [Van Leuven et al. (2023)]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [Van Leuven et al. (2023)]
/
&SPECIES_PARAMS
PSPECIES="Xe-131m", ! Tracer name: XENON 131 metastable
PDECAY=1023322.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
/
&SPECIES_PARAMS
PSPECIES="Xe-133", ! Tracer name: XENON 133
PDECAY=453388, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
/
&SPECIES_PARAMS
PSPECIES="Xe-133m", ! Tracer name: XENON 133 metastable
PDECAY=189918.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
/
&SPECIES_PARAMS
PSPECIES="Xe-135", ! Tracer name: XENON 135
PDECAY=32911.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
/
&SPECIES_PARAMS
PSPECIES="Y-91", ! Tracer name: YTTRIUM 91
PDECAY=5055782.0, ! Species half life [https://www.internetchemie.info/chemische-elemente/isotope.php]
PCRAIN_AERO=3.6, ! Below-cloud scavenging (particles) - Crain (crain_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCSNOW_AERO=1.4, ! Below-cloud scavenging (particles) - Csnow (csnow_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PCCN_AERO=1.8, ! In-cloud scavenging (particles) - CCNeff (ccn_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PIN_AERO=1.6, ! In-cloud scavenging (particles) - INeff (in_aero) [Grythe et al. (2017), Van Leuven et al. (2023)]
PDENSITY=1900.0, ! Dry deposition (particles) - rho [Van Leuven et al. (2023)]
PDIA=0.65E-06, ! Dry deposition (particles) - dquer [Van Leuven et al. (2023)]
PDSIGMA=3.0, ! Dry deposition (particles) - dsig [Van Leuven et al. (2023)]
/
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment