icon-art guide - User contributions [en]

Stratospheric Ozone Chemistry and Polar Stratospheric Clouds

2023-09-07T08:37:06Z

Luca R: Syntaxhighlighting

- under construction! -

In this example a simulation with only stratospheric (simplified) ozone chemistry is performed, also polar stratospheric clouds are taken into account. The tutorial teaches you...
* the implementation of stratospheric specific chemistry (here on thee example of ozone)
* applying linearized ozone chemistry (LINOZ) in a simulation
* the implementation of polar stratospheric clouds (PSCs)
No emission data will be included in this simulation.

== Configuration case ==
=== Stratospheric ozone with linearized ozone chemistry (LINOZ) ===
On the one hand the depicted case is dealing with simulating stratospheric ozone which is calculated with a linearized ozone chemistry (LINOZ). Here a linearized version of the differential equation for production or depletion like <math>\frac{\mathrm{d}c_i}{\mathrm{d}t}=P_i-\frac{c_i}{\tau _i}</math> is used. After calculating the Taylor expansion of first order the computed ozone concentrations are then anomalies from temperature and ozone column climatologies. This LINOZ chemistry is applied in heights of 10km and higher. Below that height the lifetime of ozone is set to constant 28 days so a constant climatology is applied. Without including LINOZ at this point, a complete ozone chemistry must be included in the simulation what would result in higher computational effort.
=== Polar Stratospheric Clouds (PSCs) ===
Further, the simulation of polar stratospheric clouds (PSC), are also considered. Polar stratospheric clouds are forming under very cold conditions in the antarctic region during the south hemispheric winter. These cold conditions are reached due to the underlying orographic conditions in Antarctica and a resulting very strong and stable polar vortex. In May or June the temperature is normally low enough to form polar stratospheric clouds. Then the reaction

<chem>ClNO3 + HCl -> Cl2 + HNO3</chem>

is effektively producing high <chem>HNO3</chem> concentrations in the stratosphere resulting in a formation of <chem>HNO3*H20</chem> phases, identifying as polar stratospheric clouds. The educts are existing due to the regular <chem>ClOx</chem> family reactions (see e.g. [https://onlinelibrary.wiley.com/doi/10.1002/ciuz.200700418#pane-pcw-references Dameris et al., 2007]). Later ion the year, when more light reaches Antarctica, the <chem>Cl2</chem> product can then be photolysed and finally, due to the <chem>ClOx</chem> cycle, Ozone can be depleted:

<chem>Cl2 + hv -> 2Cl</chem>

<chem>Cl + O3 -> ClO + O2</chem>

That's also why the ozone hole has formed: The more chlorofluorocarbons are emitted, the more <chem>ClNO3</chem> as well as <chem>HCl</chem> are produced and the more ozone can be depleted.
In this simulation, PSCs are considered to calculate the stratospheric ozone concentration but not to simulate the clouds itself.

== Setting up the Runscript ==
Let's start with the runscript that has to be prepared. Please note that the in the following explained parts have to be printed in one runscript-file with the naming designation "xyz.run". Here it is named <code>exp.testsuite.strat_ozone_linoz_psc.run</code> but of course you can call it differently as well. The runscript can be stored under the following path in your icon directory: <code>/icon-kit/run</code>

If you've already walked through the example of the [[Simplified Chemistry]], you can use nearly the same runscript. Note that you have to change the paths from Part 1, the timing settings from Part 2, the output variables from Part 3, the emission settings as well as the path of the chemtracer-xml-file in the ART-settings from Part 4 and finally the timing in the job settings from Part 5. Details can be found below.

Inside of that, first check in part 1 that all your paths to your directories are correct, probably they have to be adjusted. Note that "hp8526" is a name of a specific account here, make sure to double check especially these lines. Abbreviations used here are the following:
*CENTER: Your organization
*EXPNAME: name of your ICON-Simulation
*OUTDIR: Directory where the simulation output will be stored
*ARTFOLDER: Directory where the ICON-ART code is stored
*INDIR: Directory where the necessary Input data are stored
*EXP:
*lart: For ICON-ART Simulation that has to be switched to <code>.True.</code>.
<div class="toccolours mw-collapsible mw-collapsed">
Part 1: Runscript Directory Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
#!/bin/bash
CENTER=IMK
workspace=/hkfs/work/workspace/scratch/hp8526-strat_ozone_linoz_psc
basedir=${workspace}/icon-kit-testsuite
icon_data_poolFolder=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT/AMIP/amip_input
EXPNAME=liftime_tracer_test
OUTDIR=${workspace}/output/${EXPNAME}
ICONFOLDER=/home/hk-project-iconart/hp8526/icon-kit
ARTFOLDER=${ICONFOLDER}/externals/art
INDIR=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT
EXP=ECHAM_AMIP_LIFETIME
lart=.True.

FILETYPE=4
COMPILER=intel
restart=.False.
read_restart_namelists=.False.

# Remove folder from OUTDIR for postprocessing output
OUTDIR_PREFIX=${workspace}

# Create output directory and go to this directory

if [ ! -d $OUTDIR ]; then
mkdir -p $OUTDIR
fi

cd $OUTDIR

#input for global domain
ln -sf ${INDIR}/../INPUT/GRID/icon_grid_0014_R02B05_G.nc iconR2B05_DOM01.nc
ln -sf ${INDIR}/../INPUT/EXTPAR/icon_extpar_0014_R02B05_G.nc extpar_iconR2B05_DOM01.nc
ln -sf /hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/output/0014_R02B05/uc1_ifs_t1279_grb2_remap_rev832_0014_R02B05_2018010100.nc ifs2icon_R2B05_DOM01.nc

ln -sf $ICONFOLDER/data/rrtmg_lw.nc rrtmg_lw.nc
ln -sf $ICONFOLDER/data/ECHAM6_CldOptProps.nc ECHAM6_CldOptProps.nc

ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/mozart_coord.nc ${OUTDIR}/mozart_coord.nc
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Linoz2004Br.dat ${OUTDIR}/Linoz2004Br.dat
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Simnoy2002.dat ${OUTDIR}/Simnoy2002.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-aer.dat ${OUTDIR}/FJX_scat-aer.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j.dat ${OUTDIR}/FJX_j2j.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-cld.dat ${OUTDIR}/FJX_scat-cld.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-ssa.dat ${OUTDIR}/FJX_scat-ssa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-UMa.dat ${OUTDIR}/FJX_scat-UMa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec.dat ${OUTDIR}/FJX_spec.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_std.dat ${OUTDIR}/atmos_std.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_h2och4.dat ${OUTDIR}/atmos_h2och4.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j_extended.dat ${OUTDIR}/FJX_j2j_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended_lyman.dat ${OUTDIR}/FJX_spec_extended_lyman.dat

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# the namelist filename
atmo_namelist=NAMELIST_${EXPNAME}
#
#-----------------------------------------------------------------------------
</syntaxhighlight>
</div>

Additionally in the next lines of code you set the timing. In this simulation we simulate one week (01 October 2012 until 07 October 2012). To really catch all day times and so the time dependent solar radiation, the output interval is set to 10 hours to calculate <chem>O3</chem> to every time of the day. Because of the Photolysis dependency of <chem>O3</chem> this is particularly important.
<div class="toccolours mw-collapsible mw-collapsed">
Part 2: Runscript Timing Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
! run_nml: general switches ----------
&# model timing

start_date=${start_date:="2012-10-01T00:00:00Z"}
end_date=${end_date:="2012-10-08T00:00:00Z"}
output_start=${start_date:="2012-10-01T00:00:00Z"}
output_end=${end_date:="2012-10-08T00:00:00Z"}
output_interval="PT10H"
modelTimeStep="PT6M"
leadtime="P8H"
checkpoint_interval="P30D"
</syntaxhighlight>
</div>

Further, all the namelist parameters (from the regular ICON model without ART-extension) have to be set. For a regular ICON-ART-Simulation the following settings are recommended - if not stated differently. For a detailed description, check out the ICON Documentation ([https://code.mpimet.mpg.de/attachments/download/19568/ICON_tutorial_2019.pdf Drill et. al. (2019)]).

Since we make use of LINOZ, in the radiation namelist section <code>&radiation_nml</code> the namelist parameter <code>irad_o3</code> has to be set to 10:

<div class="toccolours mw-collapsible mw-collapsed">
Part 3: Runscript ICON-Parameter and -Namelist Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
# model parameters
model_equations=3 # equation system
# 1=hydrost. atm. T
# 1=hydrost. atm. theta dp
# 3=non-hydrost. atm.,
# 0=shallow water model
# -1=hydrost. ocean
#-----------------------------------------------------------------------------
# the grid parameters
declare -a atmo_dyn_grids=("iconR2B04_DOM01.nc" "iconR2B05_DOM02.nc" "iconR2B06_DOM03.nc")
# "iconR2B08_DOM03.nc"
#atmo_dyn_grids="iconR2B07_DOM02.nc","iconR2B08_DOM03.nc"
#atmo_rad_grids="iconR2B06_DOM01.nc"
#-----------------------------------------------------------------------------

# create ICON master namelist
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf

no_of_models=1

cat > icon_master.namelist << EOF
&master_nml
lRestart = .false.
/
&master_time_control_nml
experimentStartDate = "$start_date"
experimentStopDate = "$end_date"
forecastLeadTime = "$leadtime"
checkpointTimeIntval = "$checkpoint_interval"
/
&master_model_nml
model_type=1
model_name="ATMO"
model_namelist_filename="$atmo_namelist"
model_min_rank=1
model_max_rank=65536
model_inc_rank=1
/
EOF

#-----------------------------------------------------------------------------
#

#-----------------------------------------------------------------------------
#
# write ICON namelist parameters
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf
#
# ------------------------
# reconstrcuct the grid parameters in namelist form
dynamics_grid_filename=""
for gridfile in ${atmo_dyn_grids}; do
dynamics_grid_filename="${dynamics_grid_filename} '${gridfile}',"
done
radiation_grid_filename=""
for gridfile in ${atmo_rad_grids}; do
radiation_grid_filename="${radiation_grid_filename} '${gridfile}',"
done

# ------------------------

cat > ${atmo_namelist} << EOF
&parallel_nml
nproma = 8 ! optimal setting 8 for CRAY; use 16 or 24 for IBM
p_test_run = .false.
l_test_openmp = .false.
l_log_checks = .false.
num_io_procs = 0 ! up to one PE per output stream is possible
itype_comm = 1
iorder_sendrecv = 3 ! best value for CRAY (slightly faster than option 1)
/
&grid_nml
dynamics_grid_filename = 'iconR2B05_DOM01.nc'
dynamics_parent_grid_id = 0
!radiation_grid_filename = ${radiation_grid_filename}
lredgrid_phys = .false.
lfeedback = .true.
ifeedback_type = 2
/
&initicon_nml
lconsistency_checks = .false.
init_mode = 2 ! operation mode 2: IFS
zpbl1 = 500.
zpbl2 = 1000.
! l_sst_in = .true.
/
&run_nml
num_lev = 90
lvert_nest = .true. ! use vertical nesting if a nest is active
! nsteps = ${nsteps} ! 50 ! 1200 ! 7200 !
! dtime = ${dtime} ! timestep in seconds
modelTimeStep = "${modelTimeStep}"
ldynamics = .TRUE. ! dynamics
ltransport = .true.
iforcing = 3 ! NWP forcing
ltestcase = .false. ! false: run with real data
msg_level = 7 ! print maximum wind speeds every 5 time steps
ltimer = .true. ! set .TRUE. for timer output
timers_level = 10 ! can be increased up to 10 for detailed timer output
output = "nml"
lart = ${lart}
/
&nwp_phy_nml
inwp_gscp = 1
inwp_convection = 1
inwp_radiation = 1
inwp_cldcover = 1
inwp_turb = 1
inwp_satad = 1
inwp_sso = 1
inwp_gwd = 1
inwp_surface = 1
icapdcycl = 3 ! apply CAPE modification to improve diurnalcycle over tropical land (optimizes NWP scores)
latm_above_top = .false., .true. ! the second entry refers to the nested domain (if present)
efdt_min_raylfric = 7200.
itype_z0 = 2
icpl_aero_conv = 1
icpl_aero_gscp = 1
! resolution-dependent settings - please choose the appropriate one
dt_rad = 2160.
dt_conv = 720.
dt_sso = 1440.
dt_gwd = 1440.
/
&nwp_tuning_nml
tune_zceff_min = 0.075 ! ** default value to be used for R3B7; use 0.05 for R2B6 in order to get similar temperature biases in upper troposphere **
itune_albedo = 1 ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara in order to reduce cold bias
/
&turbdiff_nml
tkhmin = 0.75 ! new default since rev. 16527
tkmmin = 0.75 ! ""
pat_len = 100.
c_diff = 0.2
rat_sea = 8.5 ! ** new since r20191: 8.5 for R3B7, 8.0 for R2B6 in order to get similar temperature biases in the tropics **
ltkesso = .true.
frcsmot = 0.2 ! these 2 switches together apply vertical smoothing of the TKE source terms
imode_frcsmot = 2 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects:
itype_sher = 3
ltkeshs = .true.
a_hshr = 2.0
/
&lnd_nml
ntiles = 3 !!! 1 for assimilation cycle and forecast
nlev_snow = 3 !!! 1 for assimilation cycle and forecast
lmulti_snow = .true. !!! .false. for assimilation cycle and forecast
itype_heatcond = 2
idiag_snowfrac = 2
lsnowtile = .false. !! later on .true. if GRIB encoding issues are solved
lseaice = .true.
llake = .false.
itype_lndtbl = 3 ! minimizes moist/cold bias in lower tropical troposphere
itype_root = 2
/
&radiation_nml
irad_o3 = 10
irad_aero = 6
albedo_type = 2 ! Modis albedo
vmr_co2 = 390.e-06 ! values representative for 2012
vmr_ch4 = 1800.e-09
vmr_n2o = 322.0e-09
vmr_o2 = 0.20946
vmr_cfc11 = 240.e-12
vmr_cfc12 = 532.e-12
/
&nonhydrostatic_nml
iadv_rhotheta = 2
ivctype = 2
itime_scheme = 4
exner_expol = 0.333
vwind_offctr = 0.2
damp_height = 50000.
rayleigh_coeff = 0.10
lhdiff_rcf = .true.
divdamp_order = 24 ! for data assimilation runs, '2' provides extra-strong filtering of gravity waves
divdamp_type = 32 !!! optional: 2 for assimilation cycle if very strong gravity-wave filtering is desired
divdamp_fac = 0.004
l_open_ubc = .false.
igradp_method = 3
l_zdiffu_t = .true.
thslp_zdiffu = 0.02
thhgtd_zdiffu = 125.
htop_moist_proc= 22500.
hbot_qvsubstep = 22500. ! use 19000. with R3B7
/
&sleve_nml
min_lay_thckn = 20.
max_lay_thckn = 400. ! maximum layer thickness below htop_thcknlimit
htop_thcknlimit = 14000. ! this implies that the upcoming COSMO-EU nest will have 60 levels
top_height = 75000.
stretch_fac = 0.9
decay_scale_1 = 4000.
decay_scale_2 = 2500.
decay_exp = 1.2
flat_height = 16000.
/
&dynamics_nml
iequations = 3
idiv_method = 1
divavg_cntrwgt = 0.50
lcoriolis = .TRUE.
/
&transport_nml
! qv, qc, qi, qr, qs
itype_vlimit = 1,1,1,1,1
ivadv_tracer = 3, 3, 3, 3, 3
itype_hlimit = 3, 4, 4, 4 , 4
ihadv_tracer = 52, 2,2,2,2
iadv_tke = 0
/
&diffusion_nml
hdiff_order = 5
itype_vn_diffu = 1
itype_t_diffu = 2
hdiff_efdt_ratio = 24.0
hdiff_smag_fac = 0.025
lhdiff_vn = .TRUE.
lhdiff_temp = .TRUE.
/
&interpol_nml
nudge_zone_width = 8
lsq_high_ord = 3
l_intp_c2l = .true.
l_mono_c2l = .true.
support_baryctr_intp = .false.
/
&extpar_nml
itopo = 1
n_iter_smooth_topo = 1
heightdiff_threshold = 3000.
/
&io_nml
itype_pres_msl = 4 ! IFS method with bug fix for self-consistency between SLP and geopotential
itype_rh = 1 ! RH w.r.t. water
/
&output_nml
filetype = ${FILETYPE} ! output format: 2=GRIB2, 4=NETCDFv2
dom = -1 ! write all domains
output_start = "${output_start}"
output_end = "${output_end}"
output_interval = "${output_interval}"
steps_per_file = 1
include_last = .TRUE.
output_filename = 'icon-art-${EXPNAME}-chem' ! file name base
ml_varlist = 'temp','pres','group:ART_CHEMISTRY','u','v'
output_grid = .TRUE.
remap = 1
reg_lon_def = -180.,1,180.
reg_lat_def = -90.,1,90.
/
</syntaxhighlight>
</div>
Please note in the last namelist section "output_nml" that you can set all output variables that you need to postprocess your data later. All assigned variables here will be written in the output netCDF-files as well. To learn more about post processing your data, check out a later chapter of this article or the [[Postprocessing]] article.

Now, we're getting to the ICON-ART settings. To enable chemistry in an ICON-ART Simulation inn general, the switch <code>lart_chem</code> has to be set to <code>.TRUE.</code>. With <code>lart_diag_out</code> output of the diagnostic fields can be enabled. Due to setting <code>lart_chem=.TRUE.</code> either <code>lart_chemtracer</code> or <code>lart_mecca</code> has to be set to <code>.TRUE.</code>. Because we want to perform a simulation with simplified chemistry, we have to switch on <code>lart_chemtracer</code>. If this namelist parameter is set to <code>.TRUE.</code>, also <code>cart_chemtracer_xml</code> has to be fulfilled. Here you enter the path of your xml-file which describes the tracers occurring and their properties in the simulation, which is only <chem>O3</chem> in this case. How to create this xml-file is explained in the next chapter. Since PSCs shall be considered, the <code>lart_psc</code> switch has to be set to <code>.TRUE.</code>. Please also note the extra information that have to be given for <chem>O3</chem> in the xml-file that PSCs are considered correctly (see next chapter).
An example configuration for this part is shown in the following:
<div class="toccolours mw-collapsible mw-collapsed">
Part 4: Runscript ICON-ART Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
&art_nml
lart_chem = .TRUE.
lart_diag_out = .TRUE.
lart_aerosol = .FALSE.
lart_mecca = .FALSE.
lart_chemtracer = .TRUE.
lart_psc = .TRUE.
cart_chemtracer_xml = '/home/hk-project-iconart/hp8526/icon-kit/externals/art/runctrl_examples/xml_ctrl/chemtracer_lifetime_test.xml'
cart_input_folder = '${OUTDIR}'
cart_io_suffix = '0014'
iart_init_gas = 0
/
EOF
</syntaxhighlight>
</div>

Depending on the used HPC-System, some parameter concerning the running job like maximum running time and used nodes can be set. For this case study the following settings can be copied. Note that this is valid for the HoreKa HPC system and that it can differ to other systems.
<div class="toccolours mw-collapsible mw-collapsed">
Part 5: Runscript job Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
cp ${ICONFOLDER}/bin/icon ./icon.exe

cat > job_ICON << ENDFILE
#!/bin/bash -x
#SBATCH --nodes=4
#SBATCH --time=06:00:00
#SBATCH --ntasks-per-node=76
#SBATCH --partition=cpuonly
#SBATCH -A hk-project-iconart
###SBATCH --constraint=LSDF

module load compiler/intel/2022.0.2 mpi/openmpi/4.0 lib/netcdf/4.9_serial lib/hdf5/1.12_serial lib/netcdf-fortran/4.5_serial lib/eccodes/2.25.0 numlib/mkl/2022.0.2

mpirun --bind-to core --map-by core --report-bindings ./icon.exe

ENDFILE

chmod +x job_ICON
sbatch job_ICON
</syntaxhighlight>
</div>

To conclude and to double check, in the following box the complete runscript is shown once again.
<div class="toccolours mw-collapsible mw-collapsed">
Complete example configuration of the runscript
<syntaxhighlight line lang=bash class="mw-collapsible-content">
#!/bin/bash
CENTER=IMK
workspace=/hkfs/work/workspace/scratch/hp8526-strat_ozone_linoz_psc
basedir=${workspace}/icon-kit-testsuite
icon_data_poolFolder=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT/AMIP/amip_input
EXPNAME=liftime_tracer_test
OUTDIR=${workspace}/output/${EXPNAME}
ICONFOLDER=/home/hk-project-iconart/hp8526/icon-kit
ARTFOLDER=${ICONFOLDER}/externals/art
INDIR=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT
EXP=ECHAM_AMIP_LIFETIME
lart=.True.

FILETYPE=4
COMPILER=intel
restart=.False.
read_restart_namelists=.False.

# Remove folder from OUTDIR for postprocessing output
OUTDIR_PREFIX=${workspace}

# Create output directory and go to this directory

if [ ! -d $OUTDIR ]; then
mkdir -p $OUTDIR
fi

cd $OUTDIR

#input for global domain
ln -sf ${INDIR}/../INPUT/GRID/icon_grid_0014_R02B05_G.nc iconR2B05_DOM01.nc
ln -sf ${INDIR}/../INPUT/EXTPAR/icon_extpar_0014_R02B05_G.nc extpar_iconR2B05_DOM01.nc
ln -sf /hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/output/0014_R02B05/uc1_ifs_t1279_grb2_remap_rev832_0014_R02B05_2018010100.nc ifs2icon_R2B05_DOM01.nc

ln -sf $ICONFOLDER/data/rrtmg_lw.nc rrtmg_lw.nc
ln -sf $ICONFOLDER/data/ECHAM6_CldOptProps.nc ECHAM6_CldOptProps.nc

ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/mozart_coord.nc ${OUTDIR}/mozart_coord.nc
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Linoz2004Br.dat ${OUTDIR}/Linoz2004Br.dat
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Simnoy2002.dat ${OUTDIR}/Simnoy2002.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-aer.dat ${OUTDIR}/FJX_scat-aer.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j.dat ${OUTDIR}/FJX_j2j.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-cld.dat ${OUTDIR}/FJX_scat-cld.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-ssa.dat ${OUTDIR}/FJX_scat-ssa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-UMa.dat ${OUTDIR}/FJX_scat-UMa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec.dat ${OUTDIR}/FJX_spec.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_std.dat ${OUTDIR}/atmos_std.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_h2och4.dat ${OUTDIR}/atmos_h2och4.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j_extended.dat ${OUTDIR}/FJX_j2j_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended_lyman.dat ${OUTDIR}/FJX_spec_extended_lyman.dat

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# the namelist filename
atmo_namelist=NAMELIST_${EXPNAME}
#
#-----------------------------------------------------------------------------
! run_nml: general switches ----------
&# model timing

start_date=${start_date:="2012-10-01T00:00:00Z"}
end_date=${end_date:="2012-10-08T00:00:00Z"}
output_start=${start_date:="2012-10-01T00:00:00Z"}
output_end=${end_date:="2012-10-08T00:00:00Z"}
output_interval="PT10H"
modelTimeStep="PT6M"
leadtime="P8H"
checkpoint_interval="P30D"

# model parameters
model_equations=3 # equation system
# 1=hydrost. atm. T
# 1=hydrost. atm. theta dp
# 3=non-hydrost. atm.,
# 0=shallow water model
# -1=hydrost. ocean
#-----------------------------------------------------------------------------
# the grid parameters
declare -a atmo_dyn_grids=("iconR2B04_DOM01.nc" "iconR2B05_DOM02.nc" "iconR2B06_DOM03.nc")
# "iconR2B08_DOM03.nc"
#atmo_dyn_grids="iconR2B07_DOM02.nc","iconR2B08_DOM03.nc"
#atmo_rad_grids="iconR2B06_DOM01.nc"
#-----------------------------------------------------------------------------

# create ICON master namelist
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf

no_of_models=1

cat > icon_master.namelist << EOF
&master_nml
lRestart = .false.
/
&master_time_control_nml
experimentStartDate = "$start_date"
experimentStopDate = "$end_date"
forecastLeadTime = "$leadtime"
checkpointTimeIntval = "$checkpoint_interval"
/
&master_model_nml
model_type=1
model_name="ATMO"
model_namelist_filename="$atmo_namelist"
model_min_rank=1
model_max_rank=65536
model_inc_rank=1
/
EOF

#-----------------------------------------------------------------------------
#

#-----------------------------------------------------------------------------
#
# write ICON namelist parameters
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf
#
# ------------------------
# reconstrcuct the grid parameters in namelist form
dynamics_grid_filename=""
for gridfile in ${atmo_dyn_grids}; do
dynamics_grid_filename="${dynamics_grid_filename} '${gridfile}',"
done
radiation_grid_filename=""
for gridfile in ${atmo_rad_grids}; do
radiation_grid_filename="${radiation_grid_filename} '${gridfile}',"
done

# ------------------------

cat > ${atmo_namelist} << EOF
&parallel_nml
nproma = 8 ! optimal setting 8 for CRAY; use 16 or 24 for IBM
p_test_run = .false.
l_test_openmp = .false.
l_log_checks = .false.
num_io_procs = 0 ! up to one PE per output stream is possible
itype_comm = 1
iorder_sendrecv = 3 ! best value for CRAY (slightly faster than option 1)
/
&grid_nml
dynamics_grid_filename = 'iconR2B05_DOM01.nc'
dynamics_parent_grid_id = 0
!radiation_grid_filename = ${radiation_grid_filename}
lredgrid_phys = .false.
lfeedback = .true.
ifeedback_type = 2
/
&initicon_nml
lconsistency_checks = .false.
init_mode = 2 ! operation mode 2: IFS
zpbl1 = 500.
zpbl2 = 1000.
! l_sst_in = .true.
/
&run_nml
num_lev = 90
lvert_nest = .true. ! use vertical nesting if a nest is active
! nsteps = ${nsteps} ! 50 ! 1200 ! 7200 !
! dtime = ${dtime} ! timestep in seconds
modelTimeStep = "${modelTimeStep}"
ldynamics = .TRUE. ! dynamics
ltransport = .true.
iforcing = 3 ! NWP forcing
ltestcase = .false. ! false: run with real data
msg_level = 7 ! print maximum wind speeds every 5 time steps
ltimer = .true. ! set .TRUE. for timer output
timers_level = 10 ! can be increased up to 10 for detailed timer output
output = "nml"
lart = ${lart}
/
&nwp_phy_nml
inwp_gscp = 1
inwp_convection = 1
inwp_radiation = 1
inwp_cldcover = 1
inwp_turb = 1
inwp_satad = 1
inwp_sso = 1
inwp_gwd = 1
inwp_surface = 1
icapdcycl = 3 ! apply CAPE modification to improve diurnalcycle over tropical land (optimizes NWP scores)
latm_above_top = .false., .true. ! the second entry refers to the nested domain (if present)
efdt_min_raylfric = 7200.
itype_z0 = 2
icpl_aero_conv = 1
icpl_aero_gscp = 1
! resolution-dependent settings - please choose the appropriate one
dt_rad = 2160.
dt_conv = 720.
dt_sso = 1440.
dt_gwd = 1440.
/
&nwp_tuning_nml
tune_zceff_min = 0.075 ! ** default value to be used for R3B7; use 0.05 for R2B6 in order to get similar temperature biases in upper troposphere **
itune_albedo = 1 ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara in order to reduce cold bias
/
&turbdiff_nml
tkhmin = 0.75 ! new default since rev. 16527
tkmmin = 0.75 ! ""
pat_len = 100.
c_diff = 0.2
rat_sea = 8.5 ! ** new since r20191: 8.5 for R3B7, 8.0 for R2B6 in order to get similar temperature biases in the tropics **
ltkesso = .true.
frcsmot = 0.2 ! these 2 switches together apply vertical smoothing of the TKE source terms
imode_frcsmot = 2 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects:
itype_sher = 3
ltkeshs = .true.
a_hshr = 2.0
/
&lnd_nml
ntiles = 3 !!! 1 for assimilation cycle and forecast
nlev_snow = 3 !!! 1 for assimilation cycle and forecast
lmulti_snow = .true. !!! .false. for assimilation cycle and forecast
itype_heatcond = 2
idiag_snowfrac = 2
lsnowtile = .false. !! later on .true. if GRIB encoding issues are solved
lseaice = .true.
llake = .false.
itype_lndtbl = 3 ! minimizes moist/cold bias in lower tropical troposphere
itype_root = 2
/
&radiation_nml
irad_o3 = 10
irad_aero = 6
albedo_type = 2 ! Modis albedo
vmr_co2 = 390.e-06 ! values representative for 2012
vmr_ch4 = 1800.e-09
vmr_n2o = 322.0e-09
vmr_o2 = 0.20946
vmr_cfc11 = 240.e-12
vmr_cfc12 = 532.e-12
/
&nonhydrostatic_nml
iadv_rhotheta = 2
ivctype = 2
itime_scheme = 4
exner_expol = 0.333
vwind_offctr = 0.2
damp_height = 50000.
rayleigh_coeff = 0.10
lhdiff_rcf = .true.
divdamp_order = 24 ! for data assimilation runs, '2' provides extra-strong filtering of gravity waves
divdamp_type = 32 !!! optional: 2 for assimilation cycle if very strong gravity-wave filtering is desired
divdamp_fac = 0.004
l_open_ubc = .false.
igradp_method = 3
l_zdiffu_t = .true.
thslp_zdiffu = 0.02
thhgtd_zdiffu = 125.
htop_moist_proc= 22500.
hbot_qvsubstep = 22500. ! use 19000. with R3B7
/
&sleve_nml
min_lay_thckn = 20.
max_lay_thckn = 400. ! maximum layer thickness below htop_thcknlimit
htop_thcknlimit = 14000. ! this implies that the upcoming COSMO-EU nest will have 60 levels
top_height = 75000.
stretch_fac = 0.9
decay_scale_1 = 4000.
decay_scale_2 = 2500.
decay_exp = 1.2
flat_height = 16000.
/
&dynamics_nml
iequations = 3
idiv_method = 1
divavg_cntrwgt = 0.50
lcoriolis = .TRUE.
/
&transport_nml
! qv, qc, qi, qr, qs
itype_vlimit = 1,1,1,1,1
ivadv_tracer = 3, 3, 3, 3, 3
itype_hlimit = 3, 4, 4, 4 , 4
ihadv_tracer = 52, 2,2,2,2
iadv_tke = 0
/
&diffusion_nml
hdiff_order = 5
itype_vn_diffu = 1
itype_t_diffu = 2
hdiff_efdt_ratio = 24.0
hdiff_smag_fac = 0.025
lhdiff_vn = .TRUE.
lhdiff_temp = .TRUE.
/
&interpol_nml
nudge_zone_width = 8
lsq_high_ord = 3
l_intp_c2l = .true.
l_mono_c2l = .true.
support_baryctr_intp = .false.
/
&extpar_nml
itopo = 1
n_iter_smooth_topo = 1
heightdiff_threshold = 3000.
/
&io_nml
itype_pres_msl = 4 ! IFS method with bug fix for self-consistency between SLP and geopotential
itype_rh = 1 ! RH w.r.t. water
/
&output_nml
filetype = ${FILETYPE} ! output format: 2=GRIB2, 4=NETCDFv2
dom = -1 ! write all domains
output_start = "${output_start}"
output_end = "${output_end}"
output_interval = "${output_interval}"
steps_per_file = 1
include_last = .TRUE.
output_filename = 'icon-art-${EXPNAME}-chem' ! file name base
ml_varlist = 'temp','pres','group:ART_CHEMISTRY','u','v'
output_grid = .TRUE.
remap = 1
reg_lon_def = -180.,1,180.
reg_lat_def = -90.,1,90.
/

&art_nml
lart_chem = .TRUE.
lart_diag_out = .TRUE.
lart_aerosol = .FALSE.
lart_mecca = .FALSE.
lart_chemtracer = .TRUE.
lart_psc = .TRUE.
cart_chemtracer_xml = '/home/hk-project-iconart/hp8526/icon-kit/externals/art/runctrl_examples/xml_ctrl/chemtracer_ozone_psc.xml'
cart_input_folder = '${OUTDIR}'
cart_io_suffix = '0014'
iart_init_gas = 0
/
EOF

cp ${ICONFOLDER}/bin/icon ./icon.exe

cat > job_ICON << ENDFILE
#!/bin/bash -x
#SBATCH --nodes=4
#SBATCH --time=06:00:00
#SBATCH --ntasks-per-node=76
#SBATCH --partition=cpuonly
#SBATCH -A hk-project-iconart
###SBATCH --constraint=LSDF

module load compiler/intel/2022.0.2 mpi/openmpi/4.0 lib/netcdf/4.9_serial lib/hdf5/1.12_serial lib/netcdf-fortran/4.5_serial lib/eccodes/2.25.0 numlib/mkl/2022.0.2

mpirun --bind-to core --map-by core --report-bindings ./icon.exe

ENDFILE

chmod +x job_ICON
sbatch job_ICON
</syntaxhighlight>
</div>

== Setting up the xml-file ==
An xml-file describes the chemical components of the simulation which means that all trace gases or aerosols and their properties, that are relevant for the simulation, are listed here. Since we perform a simulation with simplified and linearized ozone chemistry, we need the matching chemtracer-xml-file where only <chem>O3</chem> is mentioned. Because of no specific other settings (e.g. emissions), this is the only xml-file needed here.

To prepare the xml-file we can select the matching species from the general previously generated xml-file <code>standard_chemtracer.xml</code>.

<div class="toccolours mw-collapsible mw-collapsed">
Standard Chemtracer-xml-file for simulations without emission data
<syntaxhighlight line lang=xml class="mw-collapsible-content">
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE tracers SYSTEM "tracers.dtd">

<tracers>
<chemtracer id="TRCH4" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.604E-2</mol_weight>
<?source_lifetime Hayman et al., ACP, 2014 ?>
<lifetime type="real">286977600</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRCO;2.*TRH2O</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">CH4</init_name>
</chemtracer>
<chemtracer id="TRC2H6" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">3.006E-2</mol_weight>
<?source_lifetime Hodnebrog et al, Atmos. Sci. Lett., 2018 ?>
<lifetime type="real">5011200</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRC3H8" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.40956E-2</mol_weight>
<?source_lifetime Hodnebrog et al, Atmos. Sci. Lett., 2018 ?>
<lifetime type="real">1123200</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">0.736*TRCH3COCH3</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRC5H8" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.812E-2</mol_weight>
<?source_lifetime Weimer (2015), p. 16 ?>
<lifetime type="real">8640</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCH3COCH3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">5.808E-2</mol_weight>
<?source_lifetime Weimer (2015), p. 9 ?>
<lifetime type="real">1728000</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCO" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">2.801E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">5184000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRCO2</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">CO</init_name>
</chemtracer>
<chemtracer id="TRCO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.401E-2</mol_weight>
<?source_lifetime Houghton et al., IPCC, Cambridge University Press, 2001 ?>
<lifetime type="real">3153600000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRCO</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRH2O" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.802E-2</mol_weight>
<?source_lifetime Just a placeholder (not used in the code) ?>
<lifetime type="real">2592000000</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">H2O</init_name>
</chemtracer>
<chemtracer id="TRCHBr3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">2.527E-1</mol_weight>
<?source_lifetime Rieger et al., GMD, 2015 ?>
<lifetime type="real">2073600</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCH2Br2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.738E-1</mol_weight>
<?source_lifetime Rieger et al., GMD, 2015 ?>
<lifetime type="real">10627200</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRNH3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.70E-2</mol_weight>
<?source_lifetime Pinder et al., GRL, 2008?>
<lifetime type="real">86400</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRNO2</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">NH3</init_name>
</chemtracer>
<chemtracer id="TRNO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.601E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">259200</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRHNO3</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">NO2</init_name>
</chemtracer>
<chemtracer id="TRHNO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.30E-2</mol_weight>
<?source_lifetime Day et al., ACP, 2008?>
<lifetime type="real">21600</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">HNO3</init_name>
</chemtracer>
<chemtracer id="TRSO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.40E-2</mol_weight>
<?source_lifetime Von Glasow, Chemical Geology, 2009 ?>
<lifetime type="real">1209600</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRH2SO4</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">SO2</init_name>
</chemtracer>
<chemtracer id="TROCS" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.01E-2</mol_weight>
<?source_lifetime Ullwer (2017) ?>
<lifetime type="real">504576000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRSO2</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">OCS</init_name>
</chemtracer>
<chemtracer id="TRDMS" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.21E-2</mol_weight>
<?source_lifetime Ullwer (2017) ?>
<lifetime type="real">216000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">0.993*TRSO2;0.007*TROCS</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">DMS</init_name>
</chemtracer>
<chemtracer id="TRH2SO4" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">9.80E-2</mol_weight>
<?source_lifetime Fiedler et al., ACP, 2005?>
<lifetime type="real">1800</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<dataset01 type="char" cvar_name="H2SO4" rbottom_height="10000" rupper_height="40000">CCMI-ETH_MPIC1.1</dataset01>
</chemtracer>
<chemtracer id="TRHCL" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">3.60E-2</mol_weight>
<lifetime type="real">259200</lifetime>

<c_solve type="char">param</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TR_cold">
<unit type="char">none</unit>
<mol_weight type="real">1.000E+0</mol_weight>
<lifetime type="real">432000</lifetime>

<c_solve type="char">cold</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.800E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">1576800</lifetime>

<c_solve type="char">linoz</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">O3</init_name>
</chemtracer>
<chemtracer id="TRN2O" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.401E-2</mol_weight>
<?source_lifetime Wypych, 2017, Atlas of Material Damage?>
<lifetime type="real">4730400000</lifetime>

<c_solve type="char">simnoy</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">N2O</init_name>
</chemtracer>
<chemtracer id="TRNOy" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.301E-1</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">259200</lifetime>

<c_solve type="char">simnoy</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRAGE" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">none</unit>
<mol_weight type="real">1.000E+0</mol_weight>
<?source_lifetime Just a placeholder, not used in the code?>
<lifetime type="real">25920000</lifetime>

<c_solve type="char">passive</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
</tracers>
</syntaxhighlight>
</div>

After identifying the needed chemtracers (in our case only <chem>O3</chem> (in xml-file: <code>TRO3</code>) we copy the needed code lines into our new xml-file <code>chemtracer_ozone_psc.xml</code> needed for the simulation. Additionally we add the three namelist parameter <code>polarchem</code>, <code>Thet</code> and <code>lt_het</code> as noted below. Also double check that the <code>c_solve</code> parameter is set to <code>linoz</code> to really mark that this Ozone tracer has to be computed with LINOZ chemistry and that the <code>transport</code> parameter is not set to <code>stdaero</code> but to <code>hadv52aero</code> which improves the accuracy in several cases.
<div class="toccolours mw-collapsible mw-collapsed">
Chemtracer-xml-file for stratospheric ozone simulation (Example configuration)
<syntaxhighlight line lang=xml class="mw-collapsible-content">
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE tracers SYSTEM "tracers.dtd">

<tracers>
<chemtracer id="TRO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.800E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">1576800</lifetime>

<c_solve type="char">linoz</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">hadv52aero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">O3</init_name>
<polarchem type="char">on</polarchem>
<Thet type="real">195</Thet>
<lt_het type="real">864000</lt_het>
</chemtracer>
</tracers>
</syntaxhighlight>
</div>

== Running the simulation ==
Double check all filled in paths and namelist - especially the ART-namelists. If every namelist parameter in the runscript is filled in correctly, the runscript has to be saved. Afterwards by typing
./exp.testsuite.strat_ozone_linoz_psc.run
a job can be submitted to the respective HPC-System. Type the terminal command
squeue
to view a list of your submitted and currently running and jobs.
By changing in the output directory (which is according to our runscript <code>/hkfs/work/workspace/scratch/hp8526-liftime_tracer_test</code> you can check the slurm file for possible errors and run times after your job has been run through.

In the output directory you can also find all output data for postprocessing in netCDF format.

Simulating a Point Source

2023-09-07T08:32:39Z

Luca R: Syntaxhighlighting

___ Work in Progress ___

In this Tutorial, the steps to simulate a Volcanic Eruption via a point source are given.
<gallery mode="slideshow">
File:Pointsource.gif | Pointsource simulation of the Raikoke eruption in June 2019.
</gallery>
== Setting up ICON-ART ==

The first thing to do is having a working installation of ICON-art. To check if your icon version has been built correctly, you can check if
your_icon_folder/bin/icon.exe
exists.

== Setting up Directories ==

Now a directory structure has to be set up. Usually the following directories are used:
* '''Working Directory''': Here the files that are needed for an individual run are saved. This usually includes the runscript and the relevant .xml files.
* '''Icon Code directory''': This is where the Icon code is stored.
* '''External Data directory''': Here external files which are needed for a run are stored. For example, to parametrize the optical properties of clouds a files like ECHAM6_CldOptProps.nc is used. These files of course can be switched out for others, however in most applications the same ones are used. A list of the used files is given in the Runscript, which then creates a link of these files in the Output directory.
* '''Output directory''': This is where the new Simulation data will be stored. Since most of the time large amounts of data are produced, this is stored in the work or scratch partitions on most HPC Systems. The Namelists produced by the runfile are also stored here.

== Setting up the Runscript ==

The function of the runscript is to set up the directory structure,link the relevant files and create the namelists.

== Setting up the .xml's ==
Here the .xml
pntSRC.xml
is set up. It contains all the necessary information to describe the emission of the here defined aerosols into the atmosphere.
The following information is contained:
* where is the Pointsource
* when is it emitting
* what substances are emitted
* how much of each substance is emitted
* what size are the emitted substances (median and standard deviation)

This can be adapted as needed.

In this example, the Raikoke eruption on 21. of June 2019 is Simulated, as can be seen in the following .xml:

<syntaxhighlight line lang=xml class="mw-collapsible-content">

<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE tracers SYSTEM "sources_selTrnsp.dtd">

<sources>
<pntSrc id="Raikoke-SO2">
<lon type="real">153.24</lon>
<lat type="real">49.29</lat>
<substance type="char">TRSO2</substance>
<source_strength type="real">46300.0</source_strength>
<height type="real">14000</height>
<height_bot type="real">600</height_bot>
<unit type="char">kg s-1</unit>
<startTime type="char">2019-06-21T13:00:00</startTime>
<endTime type="char">2019-06-21T22:00:00</endTime>
</pntSrc>
<pntSrc id="Raikoke-ashacc">
<lon type="real">153.24</lon>
<lat type="real">49.29</lat>
<substance type="char">ash_insol_acc</substance>
<dg3_emiss type="real">0.8E-6</dg3_emiss>
<sigma_emiss type="real">1.4</sigma_emiss>
<source_strength type="real">19700.0</source_strength>
<height type="real">14000</height>
<height_bot type="real">600</height_bot>
<unit type="char">kg s-1</unit>
<startTime type="char">2019-06-21T13:00:00</startTime>
<endTime type="char">2019-06-21T22:00:00</endTime>
</pntSrc>
<pntSrc id="Raikoke-ashcoa">
<lon type="real">153.24</lon>
<lat type="real">49.29</lat>
<substance type="char">ash_insol_coa</substance>
<dg3_emiss type="real">2.98E-6</dg3_emiss>
<sigma_emiss type="real">1.4</sigma_emiss>
<source_strength type="real">19700.0</source_strength>
<height type="real">14000</height>
<height_bot type="real">600</height_bot>
<unit type="char">kg s-1</unit>
<startTime type="char">2019-06-21T13:00:00</startTime>
<endTime type="char">2019-06-21T22:00:00</endTime>
</pntSrc>
<pntSrc id="Raikoke-ashgiant">
<lon type="real">153.24</lon>
<lat type="real">49.29</lat>
<substance type="char">ash_giant</substance>
<dg3_emiss type="real">11.35E-6</dg3_emiss>
<sigma_emiss type="real">1.4</sigma_emiss>
<source_strength type="real">19700.0</source_strength>
<height type="real">14000</height>
<height_bot type="real">600</height_bot>
<unit type="char">kg s-1</unit>
<startTime type="char">2019-06-21T13:00:00</startTime>
<endTime type="char">2019-06-21T22:00:00</endTime>
</pntSrc>
</sources>
</syntaxhighlight>

== Running the Simulation ==

== Inspecting the Output ==

== Making a Plot of the Simulation Data ==

Simplified Chemistry

2023-09-07T08:31:19Z

Luca R: Syntaxhighlighting

- work in progress -

In this first example it is shown how to perform a simulation of with simplified chemistry in ICON-ART. This tutorial teaches you...
* the basics of setting up an ICON runscript with ICON-ART settings
* the use of the most simple ICON-ART namelist parameter
*the implementation of the desired chemical species in a simulation by setting up a chemtracer xml-data for simplified chemistry simulations
* the implementation of emission data in a simulation
Emission data will be applied on the simulation as well.

== Configuration case ==
The depicted case is about simulating the tropospheric hydroxyl radical (OH), one of the most important oxidants of the atmosphere. It's main source in the lower troposphere is the photolysis of ozone and its consequent reaction of an excited oxygen atom with the surrounding water vapor:

<chem>O3 + hv -> O2 + O(^1D)</chem>
<chem>O(^1D) + H2O -> 2OH</chem>

Additionally the excited Oxygen atom reacts further with Nitrogen and Oxygen:

<chem>O3 + hv -> O2 + O(^1D)</chem>
<chem>O(^1D) + H2O -> 2OH</chem>

The main sink of OH in the Troposphere is methane and carbon monooxide:

<chem>OH + CH4 -> H2O + CH3 -> ... -> CO + HO2</chem>
<chem>OH + CO -> H + CO2</chem>
<chem>OH + CO -> HOCO</chem>

Now the OH concentrations are calculated with the respective kinetic and photolysis constants, based on chemical kinetic laws:
<math>[\ce{OH}]=\frac{\mathrm{2[\ce{O(^1D)}]}k_{\ce{H2O}}[\ce{H2O}]}{k_{\ce{CH4}}[\ce{CH4}]+(k_{\ce{CO,1}}+k_{\ce{CO,2}})[\ce{CO}]}</math>
with <math>[\ce{O(^1D)}]=\frac{J_{\ce{O3}}[\ce{O3}]}{k_{\ce{O2}}[\ce{O2}]+k_{\ce{N2}}[\ce{N2}]+k_{\ce{H2O}}[\ce{H2O}]}</math>.

Additionally emission data of the main sinks of OH are implemented. Since the simulation is performed on a R2B05-grid the following emission data are the most suitable ones for the respective trace gases:
*<chem>CH4</chem>: anthropogenic (EDGAR-432 monthly), biomass-burning (GFED3), biogenic (MEGAN-MACC)
*<chem>CO</chem>: anthropogenic (EDGAR-432 monthly)
For more information on recommended emission data see the abstract, dealing with [[Input|Emission Data]].
Since emission data are relatively large, they can also just be left out.

== Setting up the Runscript ==
Let's start with the runscript that has to be prepared. Please note that the in the following explained parts have to be printed in one runscript-file with the naming designation "xyz.run". Here it is named <code>exp.testsuite.ohsim_simple_icon.run</code> but of course you can call it differently as well. The runscript can be stored under the following path in your icon directory: <code>/icon-kit/run</code>

Inside of that, first check in part 1 that all your paths to your directories are correct, probably they have to be adjusted. Note that "hp8526" is a name of a specific account here, make sure to double check especially these lines. Abbreviations used here are the following:
*CENTER: Your organization
*EXPNAME: name of your ICON-Simulation
*OUTDIR: Directory where the simulation output will be stored
*ARTFOLDER: Directory where the ICON-ART code is stored
*INDIR: Directory where the necessary Input data are stored
*EXP:
*lart: For ICON-ART Simulation that has to be switched to <code>.True.</code>.

<div class="toccolours mw-collapsible mw-collapsed">
Part 1: Runscript Directory Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
#!/bin/bash
CENTER=IMK
workspace=/hkfs/work/workspace/scratch/hp8526-ws_icon_oh
basedir=${workspace}/icon-kit-testsuite
icon_data_poolFolder=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT/AMIP/amip_input
EXPNAME=ohsim_icon_simple_atom1
OUTDIR=${workspace}/output/${EXPNAME}
ICONFOLDER=/home/hk-project-iconart/hp8526/icon-kit
ARTFOLDER=${ICONFOLDER}/externals/art
INDIR=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT
EXP=ECHAM_AMIP_LIFETIME
lart=.True.

FILETYPE=4
COMPILER=intel
restart=.False.
read_restart_namelists=.False.

# Remove folder from OUTDIR for postprocessing output
OUTDIR_PREFIX=${workspace}

# Create output directory and go to this directory

if [ ! -d $OUTDIR ]; then
mkdir -p $OUTDIR
fi

cd $OUTDIR

#input for global domain
ln -sf ${INDIR}/../INPUT/GRID/icon_grid_0014_R02B05_G.nc iconR2B05_DOM01.nc
ln -sf ${INDIR}/../INPUT/EXTPAR/icon_extpar_0014_R02B05_G.nc extpar_iconR2B05_DOM01.nc
ln -sf /hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/output/0014_R02B05/uc1_ifs_t1279_grb2_remap_rev832_0014_R02B05_2018010100.nc ifs2icon_R2B05_DOM01.nc

ln -sf $ICONFOLDER/data/rrtmg_lw.nc rrtmg_lw.nc
ln -sf $ICONFOLDER/data/ECHAM6_CldOptProps.nc ECHAM6_CldOptProps.nc

ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/mozart_coord.nc ${OUTDIR}/mozart_coord.nc
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Linoz2004Br.dat ${OUTDIR}/Linoz2004Br.dat
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Simnoy2002.dat ${OUTDIR}/Simnoy2002.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-aer.dat ${OUTDIR}/FJX_scat-aer.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j.dat ${OUTDIR}/FJX_j2j.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-cld.dat ${OUTDIR}/FJX_scat-cld.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-ssa.dat ${OUTDIR}/FJX_scat-ssa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-UMa.dat ${OUTDIR}/FJX_scat-UMa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec.dat ${OUTDIR}/FJX_spec.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_std.dat ${OUTDIR}/atmos_std.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_h2och4.dat ${OUTDIR}/atmos_h2och4.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j_extended.dat ${OUTDIR}/FJX_j2j_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended_lyman.dat ${OUTDIR}/FJX_spec_extended_lyman.dat

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# the namelist filename
atmo_namelist=NAMELIST_${EXPNAME}
#
#-----------------------------------------------------------------------------
</syntaxhighlight>
</div>

Additionally in the next lines of code you set the timing. In this simulation we only simulate a few days. Because OH is dependent from the solar radiation, the output interval is set to 10 hours to calculate OH to every time of the day.
<div class="toccolours mw-collapsible mw-collapsed">
Part 2: Runscript Timing Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
! run_nml: general switches ----------
&# model timing

start_date=${start_date:="2016-07-29T00:00:00Z"}
end_date=${end_date:="2016-08-23T00:00:00Z"}
output_start=${start_date:="2016-07-29T00:00:00Z"}
output_end=${end_date:="2016-08-23T00:00:00Z"}
output_interval="PT10H"
modelTimeStep="PT6M"
leadtime="P10D"
checkpoint_interval="P30D"
</syntaxhighlight>
</div>

Further, all the namelist parameters (from the regular ICON model without ART-extension) have to be set. For a regular ICON-ART-Simulation the following settings are recommended - if not stated differently. For a detailed description, check out the ICON Documentation ([https://code.mpimet.mpg.de/attachments/download/19568/ICON_tutorial_2019.pdf Drill et. al. (2019)]).

Since we make use of LINOZ, in the radiation namelist section <code>&radiation_nml</code> the namelist parameter <code>irad_o3</code> has to be set to 10:

<div class="toccolours mw-collapsible mw-collapsed">
Part 3: Runscript ICON-Parameter and -Namelist Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
# model parameters
model_equations=3 # equation system
# 1=hydrost. atm. T
# 1=hydrost. atm. theta dp
# 3=non-hydrost. atm.,
# 0=shallow water model
# -1=hydrost. ocean
#-----------------------------------------------------------------------------
# the grid parameters
declare -a atmo_dyn_grids=("iconR2B04_DOM01.nc" "iconR2B05_DOM02.nc" "iconR2B06_DOM03.nc")
# "iconR2B08_DOM03.nc"
#atmo_dyn_grids="iconR2B07_DOM02.nc","iconR2B08_DOM03.nc"
#atmo_rad_grids="iconR2B06_DOM01.nc"
#-----------------------------------------------------------------------------

# create ICON master namelist
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf

no_of_models=1

cat > icon_master.namelist << EOF
&master_nml
lRestart = .false.
/
&master_time_control_nml
experimentStartDate = "$start_date"
experimentStopDate = "$end_date"
forecastLeadTime = "$leadtime"
checkpointTimeIntval = "$checkpoint_interval"
/
&master_model_nml
model_type=1
model_name="ATMO"
model_namelist_filename="$atmo_namelist"
model_min_rank=1
model_max_rank=65536
model_inc_rank=1
/
EOF

#-----------------------------------------------------------------------------
#

#-----------------------------------------------------------------------------
#
# write ICON namelist parameters
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf
#
# ------------------------
# reconstrcuct the grid parameters in namelist form
dynamics_grid_filename=""
for gridfile in ${atmo_dyn_grids}; do
dynamics_grid_filename="${dynamics_grid_filename} '${gridfile}',"
done
radiation_grid_filename=""
for gridfile in ${atmo_rad_grids}; do
radiation_grid_filename="${radiation_grid_filename} '${gridfile}',"
done

# ------------------------

cat > ${atmo_namelist} << EOF
&parallel_nml
nproma = 8 ! optimal setting 8 for CRAY; use 16 or 24 for IBM
p_test_run = .false.
l_test_openmp = .false.
l_log_checks = .false.
num_io_procs = 0 ! up to one PE per output stream is possible
itype_comm = 1
iorder_sendrecv = 3 ! best value for CRAY (slightly faster than option 1)
/
&grid_nml
dynamics_grid_filename = 'iconR2B05_DOM01.nc'
dynamics_parent_grid_id = 0
!radiation_grid_filename = ${radiation_grid_filename}
lredgrid_phys = .false.
lfeedback = .true.
ifeedback_type = 2
/
&initicon_nml
lconsistency_checks = .false.
init_mode = 2 ! operation mode 2: IFS
zpbl1 = 500.
zpbl2 = 1000.
! l_sst_in = .true.
/
&run_nml
num_lev = 90
lvert_nest = .true. ! use vertical nesting if a nest is active
! nsteps = ${nsteps} ! 50 ! 1200 ! 7200 !
! dtime = ${dtime} ! timestep in seconds
modelTimeStep = "${modelTimeStep}"
ldynamics = .TRUE. ! dynamics
ltransport = .true.
iforcing = 3 ! NWP forcing
ltestcase = .false. ! false: run with real data
msg_level = 7 ! print maximum wind speeds every 5 time steps
ltimer = .true. ! set .TRUE. for timer output
timers_level = 10 ! can be increased up to 10 for detailed timer output
output = "nml"
lart = ${lart}
/
&nwp_phy_nml
inwp_gscp = 1
inwp_convection = 1
inwp_radiation = 1
inwp_cldcover = 1
inwp_turb = 1
inwp_satad = 1
inwp_sso = 1
inwp_gwd = 1
inwp_surface = 1
icapdcycl = 3 ! apply CAPE modification to improve diurnalcycle over tropical land (optimizes NWP scores)
latm_above_top = .false., .true. ! the second entry refers to the nested domain (if present)
efdt_min_raylfric = 7200.
itype_z0 = 2
icpl_aero_conv = 1
icpl_aero_gscp = 1
! resolution-dependent settings - please choose the appropriate one
dt_rad = 2160.
dt_conv = 720.
dt_sso = 1440.
dt_gwd = 1440.
/
&nwp_tuning_nml
tune_zceff_min = 0.075 ! ** default value to be used for R3B7; use 0.05 for R2B6 in order to get similar temperature biases in upper troposphere **
itune_albedo = 1 ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara in order to reduce cold bias
/
&turbdiff_nml
tkhmin = 0.75 ! new default since rev. 16527
tkmmin = 0.75 ! ""
pat_len = 100.
c_diff = 0.2
rat_sea = 8.5 ! ** new since r20191: 8.5 for R3B7, 8.0 for R2B6 in order to get similar temperature biases in the tropics **
ltkesso = .true.
frcsmot = 0.2 ! these 2 switches together apply vertical smoothing of the TKE source terms
imode_frcsmot = 2 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects:
itype_sher = 3
ltkeshs = .true.
a_hshr = 2.0
/
&lnd_nml
ntiles = 3 !!! 1 for assimilation cycle and forecast
nlev_snow = 3 !!! 1 for assimilation cycle and forecast
lmulti_snow = .true. !!! .false. for assimilation cycle and forecast
itype_heatcond = 2
idiag_snowfrac = 2
lsnowtile = .false. !! later on .true. if GRIB encoding issues are solved
lseaice = .true.
llake = .false.
itype_lndtbl = 3 ! minimizes moist/cold bias in lower tropical troposphere
itype_root = 2
/
&radiation_nml
irad_o3 = 10
irad_aero = 6
albedo_type = 2 ! Modis albedo
vmr_co2 = 390.e-06 ! values representative for 2012
vmr_ch4 = 1800.e-09
vmr_n2o = 322.0e-09
vmr_o2 = 0.20946
vmr_cfc11 = 240.e-12
vmr_cfc12 = 532.e-12
/
&nonhydrostatic_nml
iadv_rhotheta = 2
ivctype = 2
itime_scheme = 4
exner_expol = 0.333
vwind_offctr = 0.2
damp_height = 50000.
rayleigh_coeff = 0.10
lhdiff_rcf = .true.
divdamp_order = 24 ! for data assimilation runs, '2' provides extra-strong filtering of gravity waves
divdamp_type = 32 !!! optional: 2 for assimilation cycle if very strong gravity-wave filtering is desired
divdamp_fac = 0.004
l_open_ubc = .false.
igradp_method = 3
l_zdiffu_t = .true.
thslp_zdiffu = 0.02
thhgtd_zdiffu = 125.
htop_moist_proc= 22500.
hbot_qvsubstep = 22500. ! use 19000. with R3B7
/
&sleve_nml
min_lay_thckn = 20.
max_lay_thckn = 400. ! maximum layer thickness below htop_thcknlimit
htop_thcknlimit = 14000. ! this implies that the upcoming COSMO-EU nest will have 60 levels
top_height = 75000.
stretch_fac = 0.9
decay_scale_1 = 4000.
decay_scale_2 = 2500.
decay_exp = 1.2
flat_height = 16000.
/
&dynamics_nml
iequations = 3
idiv_method = 1
divavg_cntrwgt = 0.50
lcoriolis = .TRUE.
/
&transport_nml
! qv, qc, qi, qr, qs
itype_vlimit = 1,1,1,1,1
ivadv_tracer = 3, 3, 3, 3, 3
itype_hlimit = 3, 4, 4, 4 , 4
ihadv_tracer = 52, 2,2,2,2
iadv_tke = 0
/
&diffusion_nml
hdiff_order = 5
itype_vn_diffu = 1
itype_t_diffu = 2
hdiff_efdt_ratio = 24.0
hdiff_smag_fac = 0.025
lhdiff_vn = .TRUE.
lhdiff_temp = .TRUE.
/
&interpol_nml
nudge_zone_width = 8
lsq_high_ord = 3
l_intp_c2l = .true.
l_mono_c2l = .true.
support_baryctr_intp = .false.
/
&extpar_nml
itopo = 1
n_iter_smooth_topo = 1
heightdiff_threshold = 3000.
/
&io_nml
itype_pres_msl = 4 ! IFS method with bug fix for self-consistency between SLP and geopotential
itype_rh = 1 ! RH w.r.t. water
/
&output_nml
filetype = ${FILETYPE} ! output format: 2=GRIB2, 4=NETCDFv2
dom = -1 ! write all domains
output_start = "${output_start}"
output_end = "${output_end}"
output_interval = "${output_interval}"
steps_per_file = 1
include_last = .TRUE.
output_filename = 'icon-art-${EXPNAME}-chem' ! file name base
ml_varlist = 'temp','pres','group:ART_CHEMISTRY','u','v','OH_Nconc'
output_grid = .TRUE.
remap = 1
reg_lon_def = -180.,1,180.
reg_lat_def = -90.,1,90.
/
</syntaxhighlight>
</div>
Please note in the last namelist section "output_nml" that you can set all output variables that you need to postprocess your data later. All assigned variables here will be written in the output netCDF-files as well. To learn more about post processing your data, check out a later chapter of this article or the [[Postprocessing]] article.

Now, we're getting to the ICON-ART settings. To enable chemistry in an ICON-ART Simulation inn general, the switch <code>lart_chem</code> has to be set to <code>.TRUE.</code>. With <code>lart_diag_out</code> output of the diagnostic fields can be enabled. Due to setting <code>lart_chem=.TRUE.</code> either <code>lart_chemtracer</code> or <code>lart_mecca</code> have to be set to <code>.TRUE.</code>. Because we want to perform a simulation with simplified chemistry, we have to switch on <code>lart_chemtracer</code>. If this namelist parameter is set to <code>.TRUE.</code>, also <code>cart_chemtracer_xml</code> has to be fulfilled. Here you enter the path of your xml-file which describes the tracers occurring and their properties in the simulation. How to create this xml-file is explained in the next chapter. Because this xml-file will contain information about included emission data of certain chemical species, in <code>cart_emiss_xml_file</code> the path of a second emission-xml-file has to be set.
An example configuration for this part is shown in the following:
<div class="toccolours mw-collapsible mw-collapsed">
Part 4: Runscript ICON-ART Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
&art_nml
lart_chem = .TRUE.
lart_diag_out = .TRUE.
lart_aerosol = .FALSE.
lart_mecca = .FALSE.
lart_chemtracer = .TRUE.
cart_emiss_xml_file = '${ARTFOLDER}/runctrl_examples/emiss_ctrl/emissions_R2B05_0014_cs.xml'
cart_chemtracer_xml = '/home/hk-project-iconart/hp8526/icon-kit/externals/art/runctrl_examples/xml_ctrl/chemtracer_reimus.xml'
cart_input_folder = '${OUTDIR}'
cart_io_suffix = '0014'
iart_init_gas = 0
/
EOF
</syntaxhighlight>
</div>
Please note that there are also several other namelist parameter you can select from (see [[Namelist]] article) but to perform our we're done for this part.

Depending on the used HPC-System, some parameter concerning the running job like maximum running time and used nodes can be set. For this case study the following settings can be copied. Note that this is valid for the HoreKa HPC system and that it can differ to other systems.
<div class="toccolours mw-collapsible mw-collapsed">
Part 5: Runscript job Settings (Example configuration)
<syntaxhighlight line lang=bash class="mw-collapsible-content">
cp ${ICONFOLDER}/bin/icon ./icon.exe

cat > job_ICON << ENDFILE
#!/bin/bash -x
#SBATCH --nodes=4
#SBATCH --time=12:00:00
#SBATCH --ntasks-per-node=76
#SBATCH --partition=cpuonly
#SBATCH -A hk-project-iconart
###SBATCH --constraint=LSDF

module load compiler/intel/2022.0.2 mpi/openmpi/4.0 lib/netcdf/4.9_serial lib/hdf5/1.12_serial lib/netcdf-fortran/4.5_serial lib/eccodes/2.25.0 numlib/mkl/2022.0.2

mpirun --bind-to core --map-by core --report-bindings ./icon.exe

ENDFILE

chmod +x job_ICON
sbatch job_ICON
</syntaxhighlight>
</div>

To conclude and to double check, in the following box the complete runscript is shown once again.
<div class="toccolours mw-collapsible mw-collapsed">
Complete example configuration of the runscript
<syntaxhighlight line lang=bash class="mw-collapsible-content">
#!/bin/bash
CENTER=IMK
workspace=/hkfs/work/workspace/scratch/hp8526-ws_icon_oh
basedir=${workspace}/icon-kit-testsuite
icon_data_poolFolder=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT/AMIP/amip_input
EXPNAME=ohsim_icon_mecca_atom1
OUTDIR=${workspace}/output/${EXPNAME}
ICONFOLDER=/home/hk-project-iconart/hp8526/icon-kit
ARTFOLDER=${ICONFOLDER}/externals/art
INDIR=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT
EXP=ECHAM_AMIP_LIFETIME
lart=.True.

FILETYPE=4
COMPILER=intel
restart=.False.
read_restart_namelists=.False.

# Remove folder from OUTDIR for postprocessing output
OUTDIR_PREFIX=${workspace}

# Create output directory and go to this directory

if [ ! -d $OUTDIR ]; then
mkdir -p $OUTDIR
fi

cd $OUTDIR

#input for global domain
ln -sf ${INDIR}/../INPUT/GRID/icon_grid_0014_R02B05_G.nc iconR2B05_DOM01.nc
ln -sf ${INDIR}/../INPUT/EXTPAR/icon_extpar_0014_R02B05_G.nc extpar_iconR2B05_DOM01.nc
ln -sf /hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/output/0014_R02B05/uc1_ifs_t1279_grb2_remap_rev832_0014_R02B05_2018010100.nc ifs2icon_R2B05_DOM01.nc

ln -sf $ICONFOLDER/data/rrtmg_lw.nc rrtmg_lw.nc
ln -sf $ICONFOLDER/data/ECHAM6_CldOptProps.nc ECHAM6_CldOptProps.nc

ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/mozart_coord.nc ${OUTDIR}/mozart_coord.nc
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Linoz2004Br.dat ${OUTDIR}/Linoz2004Br.dat
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Simnoy2002.dat ${OUTDIR}/Simnoy2002.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-aer.dat ${OUTDIR}/FJX_scat-aer.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j.dat ${OUTDIR}/FJX_j2j.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-cld.dat ${OUTDIR}/FJX_scat-cld.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-ssa.dat ${OUTDIR}/FJX_scat-ssa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-UMa.dat ${OUTDIR}/FJX_scat-UMa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec.dat ${OUTDIR}/FJX_spec.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_std.dat ${OUTDIR}/atmos_std.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_h2och4.dat ${OUTDIR}/atmos_h2och4.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j_extended.dat ${OUTDIR}/FJX_j2j_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended_lyman.dat ${OUTDIR}/FJX_spec_extended_lyman.dat

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# the namelist filename
atmo_namelist=NAMELIST_${EXPNAME}
#
#-----------------------------------------------------------------------------
# global timing
ndays_restart=366
dt_restart=`expr ${ndays_restart} \* 86400`
#
#-----------------------------------------------------------------------------
# model timing

start_date=${start_date:="2016-07-29T00:00:00Z"}
end_date=${end_date:="2016-08-23T00:00:00Z"}
output_start=${start_date:="2016-07-29T00:00:00Z"}
output_end=${end_date:="2016-08-23T00:00:00Z"}
output_interval="PT10H"
modelTimeStep="PT6M"
leadtime="P10D"
checkpoint_interval="P30D"

#
#-----------------------------------------------------------------------------
# model parameters
model_equations=3 # equation system
# 1=hydrost. atm. T
# 1=hydrost. atm. theta dp
# 3=non-hydrost. atm.,
# 0=shallow water model
# -1=hydrost. ocean
#-----------------------------------------------------------------------------
# the grid parameters
declare -a atmo_dyn_grids=("iconR2B04_DOM01.nc" "iconR2B05_DOM02.nc" "iconR2B06_DOM03.nc")
# "iconR2B08_DOM03.nc"
#atmo_dyn_grids="iconR2B07_DOM02.nc","iconR2B08_DOM03.nc"
#atmo_rad_grids="iconR2B06_DOM01.nc"
#-----------------------------------------------------------------------------

# create ICON master namelist
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf

no_of_models=1

cat > icon_master.namelist << EOF
&master_nml
lRestart = .false.
/
&master_time_control_nml
experimentStartDate = "$start_date"
experimentStopDate = "$end_date"
forecastLeadTime = "$leadtime"
checkpointTimeIntval = "$checkpoint_interval"
/
&master_model_nml
model_type=1
model_name="ATMO"
model_namelist_filename="$atmo_namelist"
model_min_rank=1
model_max_rank=65536
model_inc_rank=1
/
EOF

#-----------------------------------------------------------------------------
#

#-----------------------------------------------------------------------------
#
# write ICON namelist parameters
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf
#
# ------------------------
# reconstrcuct the grid parameters in namelist form
dynamics_grid_filename=""
for gridfile in ${atmo_dyn_grids}; do
dynamics_grid_filename="${dynamics_grid_filename} '${gridfile}',"
done
radiation_grid_filename=""
for gridfile in ${atmo_rad_grids}; do
radiation_grid_filename="${radiation_grid_filename} '${gridfile}',"
done

# ------------------------

cat > ${atmo_namelist} << EOF
&parallel_nml
nproma = 8 ! optimal setting 8 for CRAY; use 16 or 24 for IBM
p_test_run = .false.
l_test_openmp = .false.
l_log_checks = .false.
num_io_procs = 0 ! up to one PE per output stream is possible
itype_comm = 1
iorder_sendrecv = 3 ! best value for CRAY (slightly faster than option 1)
/
&grid_nml
dynamics_grid_filename = 'iconR2B05_DOM01.nc'
dynamics_parent_grid_id = 0
!radiation_grid_filename = ${radiation_grid_filename}
lredgrid_phys = .false.
lfeedback = .true.
ifeedback_type = 2
/
&initicon_nml
lconsistency_checks = .false.
init_mode = 2 ! operation mode 2: IFS
zpbl1 = 500.
zpbl2 = 1000.
! l_sst_in = .true.
/
&run_nml
num_lev = 90
lvert_nest = .true. ! use vertical nesting if a nest is active
! nsteps = ${nsteps} ! 50 ! 1200 ! 7200 !
! dtime = ${dtime} ! timestep in seconds
modelTimeStep = "${modelTimeStep}"
ldynamics = .TRUE. ! dynamics
ltransport = .true.
iforcing = 3 ! NWP forcing
ltestcase = .false. ! false: run with real data
msg_level = 7 ! print maximum wind speeds every 5 time steps
ltimer = .true. ! set .TRUE. for timer output
timers_level = 10 ! can be increased up to 10 for detailed timer output
output = "nml"
lart = ${lart}
/
&nwp_phy_nml
inwp_gscp = 1
inwp_convection = 1
inwp_radiation = 1
inwp_cldcover = 1
inwp_turb = 1
inwp_satad = 1
inwp_sso = 1
inwp_gwd = 1
inwp_surface = 1
icapdcycl = 3 ! apply CAPE modification to improve diurnalcycle over tropical land (optimizes NWP scores)
latm_above_top = .false., .true. ! the second entry refers to the nested domain (if present)
efdt_min_raylfric = 7200.
itype_z0 = 2
icpl_aero_conv = 1
icpl_aero_gscp = 1
! resolution-dependent settings - please choose the appropriate one
dt_rad = 2160.
dt_conv = 720.
dt_sso = 1440.
dt_gwd = 1440.
/
&nwp_tuning_nml
tune_zceff_min = 0.075 ! ** default value to be used for R3B7; use 0.05 for R2B6 in order to get similar temperature biases in upper troposphere **
itune_albedo = 1 ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara in order to reduce cold bias
/
&turbdiff_nml
tkhmin = 0.75 ! new default since rev. 16527
tkmmin = 0.75 ! ""
pat_len = 100.
c_diff = 0.2
rat_sea = 8.5 ! ** new since r20191: 8.5 for R3B7, 8.0 for R2B6 in order to get similar temperature biases in the tropics **
ltkesso = .true.
frcsmot = 0.2 ! these 2 switches together apply vertical smoothing of the TKE source terms
imode_frcsmot = 2 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects:
itype_sher = 3
ltkeshs = .true.
a_hshr = 2.0
/
&lnd_nml
ntiles = 3 !!! 1 for assimilation cycle and forecast
nlev_snow = 3 !!! 1 for assimilation cycle and forecast
lmulti_snow = .true. !!! .false. for assimilation cycle and forecast
itype_heatcond = 2
idiag_snowfrac = 2
lsnowtile = .false. !! later on .true. if GRIB encoding issues are solved
lseaice = .true.
llake = .false.
itype_lndtbl = 3 ! minimizes moist/cold bias in lower tropical troposphere
itype_root = 2
/
&radiation_nml
irad_o3 = 10
irad_aero = 6
albedo_type = 2 ! Modis albedo
vmr_co2 = 390.e-06 ! values representative for 2012
vmr_ch4 = 1800.e-09
vmr_n2o = 322.0e-09
vmr_o2 = 0.20946
vmr_cfc11 = 240.e-12
vmr_cfc12 = 532.e-12
/
&nonhydrostatic_nml
iadv_rhotheta = 2
ivctype = 2
itime_scheme = 4
exner_expol = 0.333
vwind_offctr = 0.2
damp_height = 50000.
rayleigh_coeff = 0.10
lhdiff_rcf = .true.
divdamp_order = 24 ! for data assimilation runs, '2' provides extra-strong filtering of gravity waves
divdamp_type = 32 !!! optional: 2 for assimilation cycle if very strong gravity-wave filtering is desired
divdamp_fac = 0.004
l_open_ubc = .false.
igradp_method = 3
l_zdiffu_t = .true.
thslp_zdiffu = 0.02
thhgtd_zdiffu = 125.
htop_moist_proc= 22500.
hbot_qvsubstep = 22500. ! use 19000. with R3B7
/
&sleve_nml
min_lay_thckn = 20.
max_lay_thckn = 400. ! maximum layer thickness below htop_thcknlimit
htop_thcknlimit = 14000. ! this implies that the upcoming COSMO-EU nest will have 60 levels
top_height = 75000.
stretch_fac = 0.9
decay_scale_1 = 4000.
decay_scale_2 = 2500.
decay_exp = 1.2
flat_height = 16000.
/
&dynamics_nml
iequations = 3
idiv_method = 1
divavg_cntrwgt = 0.50
lcoriolis = .TRUE.
/
&transport_nml
! qv, qc, qi, qr, qs
itype_vlimit = 1,1,1,1,1
ivadv_tracer = 3, 3, 3, 3, 3
itype_hlimit = 3, 4, 4, 4 , 4
ihadv_tracer = 52, 2,2,2,2
iadv_tke = 0
/
&diffusion_nml
hdiff_order = 5
itype_vn_diffu = 1
itype_t_diffu = 2
hdiff_efdt_ratio = 24.0
hdiff_smag_fac = 0.025
lhdiff_vn = .TRUE.
lhdiff_temp = .TRUE.
/
&interpol_nml
nudge_zone_width = 8
lsq_high_ord = 3
l_intp_c2l = .true.
l_mono_c2l = .true.
support_baryctr_intp = .false.
/
&extpar_nml
itopo = 1
n_iter_smooth_topo = 1
heightdiff_threshold = 3000.
/
&io_nml
itype_pres_msl = 4 ! IFS method with bug fix for self-consistency between SLP and geopotential
itype_rh = 1 ! RH w.r.t. water
/
&output_nml
filetype = ${FILETYPE} ! output format: 2=GRIB2, 4=NETCDFv2
dom = -1 ! write all domains
output_start = "${output_start}"
output_end = "${output_end}"
output_interval = "${output_interval}"
steps_per_file = 1
include_last = .TRUE.
output_filename = 'icon-art-${EXPNAME}-chem' ! file name base
ml_varlist = 'z_ifc','temp','pres','group:ART_CHEMISTRY','u','v','OH_Nconc'
output_grid = .TRUE.
remap = 1
reg_lon_def = -180.,1,180.
reg_lat_def = -90.,1,90.
/
&art_nml
lart_diag_out = .TRUE.
lart_aerosol = .FALSE.
lart_chem = .TRUE.
lart_mecca = .TRUE.
lart_chemtracer = .TRUE.
cart_emiss_xml_file = '${ARTFOLDER}/runctrl_examples/emiss_ctrl/emissions_R2B05_0014_cs.xml'
cart_mecca_xml = '${ARTFOLDER}/runctrl_examples/xml_ctrl/tracers_oh_reimus.xml'
cart_chemtracer_xml = '/home/hk-project-iconart/hp8526/icon-kit/externals/art/runctrl_examples/xml_ctrl/chemtracer_reimus.xml'
cart_input_folder = '${OUTDIR}'
cart_io_suffix = '0014'
iart_init_gas = 0
/
EOF

cp ${ICONFOLDER}/bin/icon ./icon.exe

cat > job_ICON << ENDFILE
#!/bin/bash -x
#SBATCH --nodes=4
#SBATCH --time=12:00:00
#SBATCH --ntasks-per-node=76
#SBATCH --partition=cpuonly
#SBATCH -A hk-project-iconart
###SBATCH --constraint=LSDF

module load compiler/intel/2022.0.2 mpi/openmpi/4.0 lib/netcdf/4.9_serial lib/hdf5/1.12_serial lib/netcdf-fortran/4.5_serial lib/eccodes/2.25.0 numlib/mkl/2022.0.2

mpirun --bind-to core --map-by core --report-bindings ./icon.exe

ENDFILE

chmod +x job_ICON
sbatch job_ICON
</syntaxhighlight>
</div>

== Setting up the xml-files ==
An xml-file describes the chemical components of the simulation which means that all trace gases or aerosols and their properties that are relevant for the simulation are listed here. Since we perform a simulation with simplified ICON-ART chemistry we need the matching chemtracer-xml-file. Additionally we need to create an emission-xml-file since emission data of certain chemical species is included in our simulation.

=== Chemtracer-xml-file ===
This file contains all the necessary information to describe the chemical mechanism with the respective important chemical species.
The following information are given per chemical species:
* properties like mol weight, units or lifetime of the tracer
* main sink
* main reaction product after reacting with the main sink after a given lifetime
* partly information about included emission data (ANT=anthropogenic emission, BIO=biogenic emission, BBE=biomass burning emission, more information about available emission data, see [[Emission Data]])
<div class="toccolours mw-collapsible mw-collapsed">
Chemtracer-xml-file for OH chemistry (Example configuration)
<syntaxhighlight line lang=xml class="mw-collapsible-content">
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE tracers SYSTEM "tracers.dtd">

<tracers>
<chemtracer id="TRCH4" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">1.604E-2</mol_weight>
<?source_lifetime Hayman et al., ACP, 2017 ?>
<lifetime type="real">286977600</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<c_solve type="char">OH</c_solve>
<products type="char">TRCO</products>
<emiss_ANT type="char" inum_levs="1">CH4_ANT_EDGAR432-monthly</emiss_ANT>
<emiss_BBE type="char" inum_levs="1">CH4_BBE_GFED3</emiss_BBE>
<emiss_BIO type="char" inum_levs="1">CH4_BIO_MEGAN-MACC</emiss_BIO>
<init_mode type="int">0</init_mode>
<init_name type="char">CH4</init_name>
</chemtracer>
<chemtracer id="TRC2H6" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">3.006E-2</mol_weight>
<?source_lifetime Hodnebrog et al, Atmos. Sci. Lett., 2018 ?>
<lifetime type="real">5011200</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<c_solve type="char">OH</c_solve>
<init_mode type="int">0</init_mode>
<init_name type="char">C2H6</init_name>
</chemtracer>
<chemtracer id="TRC3H8" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">4.40956E-2</mol_weight>
<?source_lifetime Hodnebrog et al, Atmos. Sci. Lett., 2018 ?>
<lifetime type="real">1123200</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<c_solve type="char">OH</c_solve>
<products type="char">0.736*TRCH3COCH3</products>
<init_mode type="int">0</init_mode>
<init_name type="char">C3H8</init_name>
</chemtracer>
<chemtracer id="TRC5H8" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">6.812E-2</mol_weight>
<?source_lifetime Weimer (2015), p. 16 ?>
<lifetime type="real">8640</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
</chemtracer>
<chemtracer id="TRCH3COCH3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">5.808E-2</mol_weight>
<?source_lifetime Weimer (2015), p. 9 ?>
<lifetime type="real">1728000</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
</chemtracer>
<chemtracer id="TRCO" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">2.801E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">5184000</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<c_solve type="char">OH</c_solve>
<products type="char">TRCO2</products>
<emiss_ANT type="char" inum_levs="1">CO_ANT_EDGAR432-monthly</emiss_ANT>
<init_mode type="int">0</init_mode>
<init_name type="char">CO</init_name>
</chemtracer>
<chemtracer id="TRCO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">4.401E-2</mol_weight>
<?source_lifetime Houghton et al., IPCC, Cambridge University Press, 2001 ?>
<lifetime type="real">3153600000</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">lt</c_solve>
<products type="char">TRCO</products>
</chemtracer>
<chemtracer id="TRNH3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">1.70E-2</mol_weight>
<?source_lifetime Pinder et al., GRL, 2008?>
<lifetime type="real">86400</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
<products type="char">TRNO2</products>
</chemtracer>
<chemtracer id="TRNO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">4.601E-2</mol_weight>
<?source_lifetime Wypych, 2017, Atlas of Material Damage?>
<lifetime type="real">4730400000</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
<products type="char">TRHNO3</products>
</chemtracer>
<chemtracer id="TRSO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">6.40E-2</mol_weight>
<?source_lifetime Von Glasow, Chemical Geology, 2009 ?>
<lifetime type="real">1209600</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
<products type="char">TRH2SO4</products>
</chemtracer>
<chemtracer id="TROCS" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">6.01E-2</mol_weight>
<?source_lifetime Ullwer (2017) ?>
<lifetime type="real">504576000</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
<products type="char">TRSO2</products>
</chemtracer>
<chemtracer id="TRDMS" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">6.21E-2</mol_weight>
<?source_lifetime Ullwer (2017) ?>
<lifetime type="real">216000</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">OH</c_solve>
<products type="char">0.993*TRSO2;0.007*TROCS</products>
</chemtracer>
<chemtracer id="TRHNO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">6.30E-2</mol_weight>
<?source_lifetime Day et al., ACP, 2008?>
<lifetime type="real">21600</lifetime>
<transport type="char"> hadv52aero </transport>
<init_name type="char"> HNO3 </init_name>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">lt</c_solve>
</chemtracer>
<chemtracer id="TRH2SO4" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">9.80E-2</mol_weight>
<?source_lifetime Fiedler et al., ACP, 2005?>
<lifetime type="real">1800</lifetime>
<transport type="char"> hadv52aero </transport>
<unit type="char">mol mol-1</unit>
<init_mode type="int">0</init_mode>
<c_solve type="char">lt</c_solve>
</chemtracer>
<chemtracer id="TRO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<mol_weight type="real">4.800E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">1555200</lifetime>
<transport type="char"> hadv52aero </transport>
<init_mode type="int"> 0 </init_mode>
<init_name type="char">O3</init_name>
<unit type="char">mol mol-1</unit>
<c_solve type="char">linoz</c_solve>
</chemtracer>
</tracers>
</syntaxhighlight>
</div>

=== Emission-xml-file ===
To link the information from the chemtracer-xml-file that we want to include emission data of a specific chemical species with the specific information about the respective emission, ICON-ART is working with, we need a second xml-file: The emission-xml-file. It is normally a standard grid-dependant xml-file you can use where nothing has to be changed. Since this Simulation is performed on a R2B05 grid, we use the <code>emissions_R2B05_0014_cs.xml</code>. Have a look at the article about [[Emission Data]] to check out all the available emission-xml-files for different resolutions.

== Running the simulation ==
Double check all filled in paths and namelist - especially the ART-namelists. If every namelist parameter in the runscript is filled in correctly, the runscript has to be saved. Afterwards by typing
./exp.testsuite.ohsim_simple_icon.run
a job can be submitted to the respective HPC-System. Type the terminal command
squeue
to view a list of your submitted and currently running and jobs.
By changing in the output directory (which is according to our runscript <code>/hkfs/work/workspace/scratch/hp8526-ws_icon_oh/output/ohsim_icon_simple_atom1</code> you can check the slurm file for possible errors and run times after your job has been run through.

In the output directory you can also find all output data for postprocessing in netCDF format.

Postprocessing

2023-09-07T08:16:53Z

Luca R: Syntaxhighlighting

== Output Checks with SAMOA ==

SAMOA performs a sanity check on all model outputs that can be read by CDO. It checks if a variable lies in-between a predefined range and if the minimum and maximum value of each variable are the same. For this purpose CDO version 1.6.2rc3 is required currently (see https://code.zmaw.de/projects/cdo).

For more information about the usage please refer to the README-file within the SAMOA package. You can get a copy of the SAMOA script by writing an e-mail to the contact person of the ART code (see http://icon-art.imk-tro.kit.edu). SAMOA is licensed under the GNU GENERAL PUBLIC LICENSE Version 3.

As SAMOA is primarily developed for the usage with COSMO-ART and COSMO-CLM, you have to do a minor change before using it. The latest version of SAMOA has a list for the usage of SAMOA with ICON-ART output included but not loaded automatically. This list is called samoa_list_icon-art. You have to replace the default (COSMO) list that is used by SAMOA by editing samoa.sh:

Search for the following lines:

<syntaxhighlight lang=bash># Path to the list with variables (is overwritten when -l specified)
# Assumed to be on the same path as script

path_list=$SCRIPTPATH/list </syntaxhighlight>
Change the name of the list to:

<syntaxhighlight lang=bash>path_list=$SCRIPTPATH/samoa_list_icon-art </syntaxhighlight>
Now you may use SAMOA with the ICON-ART output file out.nc with the following command:

<syntaxhighlight lang=bash>./samoa.sh out.nc</syntaxhighlight>
For all options see:

<syntaxhighlight lang=bash>./samoa.sh --help</syntaxhighlight>

== Data Processing ==
Some ways to handle and modify the outputs of ICON-ART are described here.

=== CDO ===
[https://code.mpimet.mpg.de/projects/cdo CDO] is a command line tool to interact with netCDF files, enabling to list the variables and do calculations efficiently.
Different examples for its usage can be found in the [https://code.mpimet.mpg.de/projects/cdo/wiki/FAQ CDO FAQ]
=== ICONTOOLS ===
[https://wiki.c2sm.ethz.ch/pub/MODELS/ICONDwdIconTools/doc_icontools.pdf ICONTOOLS] is a command line tool used to produce the required external parameters data of a simulation and is especially good at remapping data from one grid to another.

== Visualisation ==
There are many ways to visualize data produced by ICON-ART. In general, there are two possibilities: The output may exist on the ICON grid or it may exist on an interpolated longitude/latitude grid. This can be chosen by adaptions of the output namelist. Although it comes along with a loss in information, sometimes it is only possible to use interpolated output. By this, the visualization is much easier to handle.

In the following sections some tools are introduced which can be used to visualize ICON output. Note, that '''only NETCDF''' is supported by ICON-ART so far. With the tool Ncview it is very easy to have a quick look into the interpolated model output. With ParaView and Met3D a nice-looking three-dimensional visualization can be created. Python in the recent years has become the standard tool for 2D visualisation and its capabilities far exceed those of Ncview.

=== Ncview ===

"Ncview is a visual browser for netCDF format files. Typically you would use ncview to get a quick and easy, push-button look at your netCDF files. You can view simple movies of the data, view along various dimensions, take a look at the actual data values, change color maps, invert the data, etc." (http://meteora.ucsd.edu/~pierce/ncview_home_page.html)

Please note that Ncview does only work for latitude-longitude grid data and cannot be used for RXXBXX-style ICON grids without remapping.

=== Met3d ===

Met3D is a visualization software that can be used on a HPC system to visualize ICON as well as ICON-ART data in 3D. An online documentation can be found [https://met3d.wavestoweather.de/met-3d.html here]. Met3D needs remapped ICON fields to Latitude-Longitude . Furthermore, it is important to have the pressure (pres - variable) as model output, as Met3D transforms model levels in pressure levels for plotting.

=== ParaView ===

"ParaView is an open-source, multi-platform data analysis and visualization application. ParaView users can quickly build visualizations to analyze their data using qualitative and quantitative techniques. The data exploration can be done interactively in 3D or programmatically using ParaView’s batch processing capabilities.

ParaView was developed to analyze extremely large datasets using distributed memory computing resources. It can be run on supercomputers to analyze datasets of exascale size as well as on laptops for smaller data." (http://www.paraview.org/)

Paraview can be used for both Latitude-Longitude grids as well as RXXBXX-style ICON grids.

=== Python ===

On the [https://www.python.org/about/ official Website] Python describes itself as

<pre>Python is powerful... and fast;
plays well with others;
runs everywhere;
is friendly & easy to learn;
is Open.</pre>
Using Python is a simple but effective way to display ICON-ART model output data. There is a large number of Packages available to help with Visualisation, the most useful Packages for visualising ICON-ART data are given in this [[#tab:pythonpackages|Table]]

<div id="tab:pythonpackages">

{| class="wikitable" style="text-align:left;"
|+ Helpful Python Packages and their primary usage for Visualisation
|align="center"| numpy
|align="center"| predefined Mathematical functions
|-
|align="center"| matplotlib
|align="center"| Plotting framework
|-
|align="center"| xarray
|align="center"| reading in and processing netcdf datasets
|-
|align="center"| cartopy
|align="center"| Include country borders in plots
|-
|align="center"| ...
|align="center"|
|}

Python can be used for both Latitude-Longitude grids as well as RXXBXX-style ICON grids.
</div>

Output

2023-09-07T08:15:20Z

Luca R: Syntaxhighlighting

__TOC__

== General Remarks ==

In principle, output of ICON-ART variables works the same way as for ICON variables. As described in , the following five quantities of the output have to be specified:

* The time interval between two model outputs.
* The name of the output file.
* The name of the variable(s) and/or variable group(s).
* The type of vertical output grid.
* The type of horizontal output grid.

For the best results it is recommended to use NETCDF output on the icosahedral grid which ICON-ART is using. However in some applications remapping the grid to a latitude-longitude grid may be required, which can be set via the <code>remap</code> option. A corresponding output namelist for sea salt on model levels can be seen here:

<syntaxhighlight lang=bash line>
NAMELIST EXAMPLE
&output_nml
filetype = 4 ! output format: 2=GRIB2, 4=NETCDFv2
dom = 1 ! write output for domain 1
output_start = 'JJJJ-MM-DDTHH:MM:SSZ' !put date in
output_end = 'JJJJ-MM-DDTHH:MM:SSZ' !put date in
output_interval = 'PT1H' ! see https://en.wikipedia.org/wiki/ISO_8601
steps_per_file = 1 ! max. num. of time steps within one file
mode = 1 ! 1: forecast mode (relative t-axis)
include_last = .TRUE. ! include the last time step
output_filename = 'INSERTFILENAME' ! file name base
ml_varlist = 'seasa','seasb','seasc',
'seasa0','seasb0','seasc0'
remap = 1 ! output is transferred to lat long grid
reg_lon_def = -180.,0.5,179.5 !start, incr., end, in deg.
reg_lat_def = 90.,-0.5, -90. !start, incr., end, in deg.
</syntaxhighlight>

There is an option to obtain all diagnostic Variables of a certain Group without having to specifying all of them. For example, you may use the group ART_DIAGNOSTICS.

To include a group of Variables in the output file change the namelist variable ml_varlist from the example above to the following:

<pre>ml_varlist = 'group:ART_AEROSOL'</pre>
The output variables that are associated to this group will be written to the output file. You can check the groups of output variables in this [[#OutputTable | Table]] .

== Aerosol Naming Conventions ==
The following table contains an overview of the possible output variables.

There are several ways to choose the Naming of the output variables, depending on your application

* '''Externally mixed Aerosols:''' : The Tracers for Dust, Seasalt, Ash and Soot are combined with the three Possible modes a, b and c, which correspond to the different size bins of the particles <div id="OutputTable"></div>

::{| class="wikitable" style="text-align:left;"
|+ Externally mixed Tracer modes
!|| dust||seasalt ||ash ||soot
|-
!a
| dusta ||seasa ||asha ||soota
|-
!b
| dustb ||seasb ||ashb ||-
|-
!c
| dustc||seasc ||ashc ||-
|-
|}

* '''Internally mixed Aerosols ([[AERODYN]]):''' Here a tracer is defined in a different way, with the goal being to have a more flexible framework for various applications. In this framework modes are created in a different way, as illustrated int the table below:

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ AERODYN mode configurations
! || aitken || accumulation || coarse || giant
|-
! soluble
| sol_aitken || sol_acc || sol_coarse || sol_giant
|-
! insoluble
|insol_aitken || insol_acc || insol_coarse || insol_giant
|-
! mixed
| mixed_aitken ||mixed_acc || mixed_coarse || mixed_giant
|-
|}

The modes are then combined with The Tracername to obtain the name of the Variable using <code>varname = 'Tracer' + '_' + 'mode from Table'</code>.
Example : <code> dust_insol_acc </code>

* '''Monodisperse Aerosols'''

== Available Output Variables ==
The following Table contains an overview over the diagnostic Icon-ART Variables. Expressions in Brackets are Placeholders which can be used to construct the name of the actual variables:

[aeronet wavelength] => [340, 380, 440, 500, 550, 675, 870, 1020, 1064]

[ceilo_wavelength] => [355,532,1064]

[pollen] => [ALNU,BETU,...] to be defined in diagnostics.xml

{| class="wikitable"
|+ Output Overview
|-
!
! varname
! groups
! unit
! descripition
! namelist switch
! required xml
|-
!rowspan="19" style="vertical-align:top;" | Aerosols
| diam_[mode]
| ART_DIAGNOSTICS
| m
| with AERODYN : aerosol diameter
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| diam_[aerosol][mode]
| ART_DIAGNOSTICS
| m
| WITHOUT AERODYN: aerosol diameter
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| aod_[aerosol]_[aeronet wavelength]nm
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| Layer-1
| [AEROSOL] optical depth
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| bsc_[arosol]_[ceilo_wavelength]nm
| ART_DIAGNOSTICS
| m-1 sr-1
| [AEROSOL] backscatter
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| ceil_[arosol]_[ceilo_wavelength]nm
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| m-1 sr-1
| [AEROSOL] Attenuated Backscatter Ceilometer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| sat_[arosol]_[ceilo_wavelength]nm
| ART_DIAGNOSTICS
| m-1 sr-1
| [AEROSOL] Attenuated Backscatter Satellite
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| AOD_550_so4_sol
| ART_DIAGNOSTICS
| layer-1
| SO4 sol Optical Depth
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| AOD_550_ash_insol
| ART_DIAGNOSTICS
| layer-1
| Ash insol Optical Depth
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| AOD_550_ash_mixed
| ART_DIAGNOSTICS
| layer-1
| Ash mixed Optical Depth
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| AOD_550_ash_giant
| ART_DIAGNOSTICS
| layer-1
| Ash giant Optical Depth
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| ustar_thres
| ART_ROUTINE_DIAG
| m s-1
| threshold friction velocity for dust emission'
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| ustar
| ART_ROUTINE_DIAG
| m s-1
| Friction velocity
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| acc_drydepo_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2
| accumulated dry deposition of tracer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| acc_sedim_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2
| accumulated sedimentation of tracer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| acc_wetdepo_gscp_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2
| accumulated wet deposition by grid scale precipitation of tracer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| acc_wetdepo_con_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2
| accumulated wet deposition by convective precipitation of tracer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| acc_wetdepo_rrsfc_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2
| accumulated wet deposition of tracer if precipitation reaches surface
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| emiss_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2 s-1
| emission of tracer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
| acc_emiss_[tracer]
| ART_DIAGNOSTICS , ART_ROUTINE_DIAG
| tracer-unit m-2
| accumulated emission of tracer
| lart_aerosol=True and lart_diag_out=True
| requires cart_diagnostics_xml file
|-
!rowspan="9" style="vertical-align:top;"| Pollen
| [pollen]rprec
| ART_ROUTINE_DIAG
| m-2
| precipitation reservoir of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]reso
| ART_ROUTINE_DIAG
| m-2
| Pollen reservoir (previous timestep) of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]ress
| ART_ROUTINE_DIAG
| m-2
| Pollen reservoir (daily sum) of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]sdes
| ART_ROUTINE_DIAG
| -
| State of pollen season of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]ctsum
| ART_ROUTINE_DIAG
| K
| Cumulated weighted 2m temperature sum of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]saisn
| ART_ROUTINE_DIAG
| days
| Number of days since start of pollen season of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]saisl
| ART_ROUTINE_DIAG
| days
| length of pollen season of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
| [pollen]saisa
| ART_ROUTINE_DIAG
| days
| Number of days since the start of pollen season of [pollen]. if present day is out of the season: length of current season
| REQUIRES diagnostics.xml
|
|-
| [pollen]fe
| ART_ROUTINE_DIAG
| m-2 s-1
| Emission flux of [pollen]
| iart_pollen>0
| REQUIRES diagnostics.xml
|-
!rowspan="18" style="vertical-align:top;" | Chemistry
| reac_rates
| ART_DIAGNOSTICS
| s-1
| MECCA reaction rates
| lart_mecca=True, lart_diag_out=True
|
|-
| art_o3
| kg/kg
| Ozone mass mixing ratio
| lart_chem =True, lart_diag_out=True
|
|
|-
| OH_Nconc
| # / cm3
| OH number concentration
| lart_chem =TRUE
|
|
|-
| photo
| -
| s-1
| photolysis rates
| lart_chem=TRUE, lart_mecca=TRUE
|
|-
| art_full_chemistry_o3_col
| -
| DU
| Ozone column
| lart_chem=TRUE, lart_mecca=TRUE
|
|-
| sts_liqsur
| ART_DIAGNOSTICS
| cm2 cm-3
| liquid area density of STS
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| cgaml
| ART_DIAGNOSTICS
| -
| STS uptake coefficient of the reaction
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| dens_ice
| ART_DIAGNOSTICS
| m-3
| number density of ice particles
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| radius_ice
| ART_DIAGNOSTICS
| m
| radius of ice particles
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| radius_STS
| ART_DIAGNOSTICS
| m
| radius of STS particles
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| dens_NAT
| ART_DIAGNOSTICS
| m-3
| number density of NAT particles
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| radius_NAT
| ART_DIAGNOSTICS
| m
| radius of NAT particles
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| HNO3_Nconc_s
| ART_DIAGNOSTICS
| cm-3
| number concentration of HNO3 in NAT
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| HNO3_Nconc_l
| ART_DIAGNOSTICS
| cm-3
| number concentration of HNO3 in STS
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| ice_vmr_Marti
| ART_DIAGNOSTICS
| mol mol-1
| volume mixing ratio of solid water by Marti and Mauersberger
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| NAT_sedi_rel_difference
| ART_DIAGNOSTICS
| -
| relative difference of NAT mass bef and aft sedi (aft - bef) * 2 / (aft + bef)
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| NAT_sedi_vel
| ART_DIAGNOSTICS
| m s-1
| sedimentation velocity of NAT particles
| lart_chem=TRUE , lart_psc=TRUE
|
|-
| art_so2_col
| ART_DIAGNOSTICS
| DU
| SO2 column
| lat_chem=TRUE , lart_chemtracer=TRUE
|
|-
!rowspan="3" style="vertical-align:top;"| Radioactive Tracer Diagnostics
| wet deposition of xml defined tracer
| ART_DIAGNOSTICS, ART_ROUTINE_DIAG
| Bq m-2
| wet deposition of xml defined tracer
| lart_aerosol=True and iart_radioact=1
|
|-
| dry deposition of xml defined tracer
| ART_DIAGNOSTICS, ART_ROUTINE_DIAG
| Bq m-2
| dry deposition of xml defined tracer
| lart_aerosol=True and iart_radioact=1
|
|-
| Averaged air concentration of xml defined traer
| ART_DIAGNOSTICS, ART_ROUTINE_DIAG
| Bq m-3
| Averaged air concentration of xml defined traer
| lart_aerosol=True and iart_radioact=1
|
|-
!rowspan="4" style="vertical-align:top;"| FPLUME Output
| plume_height
| ART_FPLUME
| m
| plume height
| iart_fplume/=0
|
|-
| plume_MFR
| ART_FPLUME
| kg s-1
| plume MFR
| iart_fplume/=0
|
|-
| MER_transport
| ART_FPLUME
| kg s-1
| Amount of very fine ash for transport
| iart_fplume/=0
|
|-
| solution_with
| ART_FPLUME
| -
| FPlume off, Mastin, or FPlume
| iart_fplume/=0
|
|}

MECCA-based (full) Chemistry

2023-09-07T08:11:01Z

Luca R: Syntaxhighlighting

- work in progress -

In this configuration example a simulation with (full) MECCA-based chemistry is performed. If we talk about MECCA-based chemistry we mean a full gas phase chemistry that can be applied additionally to the existing standard parametrized chemistry from ICON-ART (explained in the article [[Atmospheric Chemistry]], example: [[Simplified Chemistry|click here]]). MECCA uses the Kinetic PreProcessor (KPP) to convert the chemical equations into differential equations that can be used then e.g., in ICON-ART. MECCA already contains a comprehensive chemical mech-
anism, but it also allows to create an own mechanism by changing equations or only picking a subset of equations.
To perform this example, a complete reaction mechanism is created and transferred to ICON-ART.

This article teaches you...
*the implementation of (full) MECCA-based chemistry in ICON-ART
*the creation of a chemical mechanism and the selection of the respective desired chemical species and their reactions in MECCA
*the implementation of not yet in MECCA implemented reactions in your mechanism
*the creation of the to your mechanism belonging Mecca-xml data to link MECCA calculations with ICON-ART

== Introductioin ==
The MECCA(=Module Efficiently Calculating the Chemistry of the Atmosphere) based chemistry describes a full gas phase chemistry that can be applied as an extension to the parametrized [[Atmospheric Chemistry|Simplified Chemistry]] (see above). MECCA based chemistry is generally more exact in the concentration values but the overall runtime is longer compared to purely simplified chemistry simulations. MECCA itself is originally a submodule of the CAABA box model where an air parcel is described as a box and outgoing from this model all exchange processes in- and outward of the box are calculated. As MECCA is part of this model, it contains a wide collection of the most important reactions, including Ozone-, Methane-, HOx-, NOx-, Carbonhydrogen-, Halogene- and Sulfur chemistry. MECCA is available in a [http://www.geosci-model-dev.net/4/373/2011/gmd-4-373-2011-supplement.zip supplement], available to download for free and containing all auxiliaries to perform MECCA-simulations.

In this case we are going to have a look at the implementation of the fully available chemistry, meaning all in default mode available reactions will be implemented. This could also be adjusted easily if necessary and will be showed later. If we talk about "adding" MECCA chemistry to ICON-ART, it means that it is additionally calculated to the in any way calculated parametrized simplified chemistry in ICON-ART.

(Note: Adding MECCA-chemistry to ICON-ART also means that some namelist parameters will be overwritten in the runscript which means that some extra options (e.g. LINOZ-chemistry) are only available by setting them manually in the ICON-ART code.)

[In this configuration case a regular simulation with a MECCA chemistry implementation that will be applied in most cases is explained. The implementation and link to ICON-ART works with an xml-file that has to created with help of the [http://www.geosci-model-dev.net/4/373/2011/gmd-4-373-2011-supplement.zip caaba3.0-supplement]. The created xml-file has to be included in the runscript as well. To get a better overview about the upcoming steps you can also check out the MECCA chemistry part in the [[Atmospheric Chemistry]] article.

== Preparing the MECCA-xml-file ==

=== Step 0: Download and open the caaba3.0-supplement ===
Once you have downloaded the [http://www.geosci-model-dev.net/4/373/2011/gmd-4-373-2011-supplement.zip caaba3.0-supplement], you can open it in your preferred directory browser or terminal.

=== Step 1: Setting up the chemical mechanism ===
If you browse the <code>mecca</code> directory you can check out the <code>gas.eqn</code> file. Just use the text editor if you use the directory browser or use the command
vi gas.eqn
in the terminal. Once you have opened it, you can see all available reactions implemented in MECCA with their respective reaction codes which represents the full MECCA mechanism. In this configuration case we will implement a MECCA chemistry with all available reactions. That's why the <code>gas.eqn</code> can be left as it is.

==== What to do if not all reactions are wanted ====
First, never edit the <code>gas.eqn</code> itself! Better copy and rename it for your respective chemical mechanism that you want to create, e.g. <code>gas_Mechanism1.eqn</code>. Afterwards you can open it and delete all reactions that are not wanted (for Terminal users: to activate writing inside vi, use command <code>:w</code>, for closing and saving: <code>:wq</code>).

==== What to do if you want to edit existing reactions of the gas.eqn or add new reactions to your mechanism ====
If a reaction of the <code>gas.eqn</code> is only similar to that one that you want to implement or a specific reaction is not implemented in the <code>gas.eqn</code>, you can make use of replacement-files. Inside the <code>mecca</code> directory, select the <code>rpl</code> directory in which you can copy and rename the <code>example.rpl</code> first of all. If you open it (again with the text editor or with the <code>vi</code>-command you can edit and add your reactions, depending on your scientific goal.
If you want to edit a reaction, use the <code>#REPLACE</code> command as well as the number of the respective reaction (e.g. <G4110>) in the first line, then write <code><a></code> (<code><nowiki><b></nowiki></code>, <code><c></code>,...) for your first (second, third,...) reaction which belongs to the same reaction number. Below some examples are shown.
* <u>Editing one reaction:</u> Reaction <code>G4110</code> originally looks like this inside the <code>gas.eqn</code>:
<pre>
<G4110> CO + OH = H + CO2 : {%StTrG} (1.57E-13+cair*3.54E-33){§1.15}; {&1628}
</pre>

Inside the replacement file you can change it for example like this:
<pre>
#REPLACE <G4110>
<a> CO + OH = HO2 + CO2 : {%StTrG} 1.57E-13 + cair*3.54E-33 {&1628}
#ENDREPLACE
</pre>
* <u>Split one reaction into two (or more):</u> If you want to split a reaction into several subreactions (e.g. reaction <code>G4101</code>). Originally <code>G4110</code> looks like this inside the <code>gas.eqn</code>:
<G4101> CH4 + OH {+O2}= CH3O2 + H2O : {%StTrG} 1.85E-20{§1.2}*EXP(2.82*log(temp)-987./temp); {&1627}
The splitting can then be written as:
#REPLACE <G4101>
<a> CH4 + OH = CH3 + H2O : {%StTrG} 1.85E-20*EXP(2.82*log(temp)-987./temp); {&&1627}
<nowiki><b></nowiki> CH3 + O2 = CH3O2 : {%StTrG} 1E-999; {&&}
#ENDREPLACE
* <u>Adding new reaction:</u> If you want to add a new reaction, define a not yet existing reaction number and write your new reaction:
#REPLACE <>
<G9876JD> XYZ + OH = RO2 + H2O : {%StG} 1.57E-13; {&&}
#ENDREPLACE
Note:
* The added reaction only works for your specific project, it is not implemented in MECCA then and it can't be applied from the <code>gas.eqn</code> afterwards.
* if you try to edit a reaction with a not existing reaction number, an error message will be given

If you're finished, don't forget to save everything (for Terminal users: to activate writing inside vi, use command <code>:w</code>, for closing and saving: <code>:wq</code>)!

=== Step 2: Setting up the batch file ===
Inside the <code>mecca</code> directory of the supplement, check out the <code>example.bat</code> file in the <code>batch</code> directory. The batch file is some kind of runscript of the respective MECCA-mechanism to create the belonging Fortran Code and the ICON-Code afterwards. Copy and rename it according to your project. After opening it afterwards you will see the following code lines:
<div class="toccolours mw-collapsible mw-collapsed">
Example batch file (<code>example.bat</code>)
<syntaxhighlight lang=bash line class="mw-collapsible-content">
# -*- Shell-script -*-

# The shell variables defined here will be used by xmecca
# when it is run in batch mode (i.e. not interactive).

set apn = 2 # number of aerosol phases [0...99, default=0]
set gaseqnfile = gas.eqn
set rplfile = # no replacements
set wanted = "Tr && (G || (Aa && Mbl)) && \!I && \!Hg"
set mcfct = n # Monte-Carlo factor?
set diagtracfile = # diagnostic tracers?
set rxnrates = n # calculate accumulated reaction rates?
set tagdbl = n # tagging, doubling, both, none ??
set kppoption = k # k=kpp, 4=kp4, q=quit
set integr = rosenbrock_posdef # integrator
set vlen = 256 # only for kp4 and integr=rosenbrock_vec
set decomp = n # remove indirect indexing
# kp4: 0/1/2/3/q; kpp: y/n/q
set deltmpkp4 = y # delete temporary kp4 files
set latex = n # latex list of reactions
set graphviz = n # graphviz plots?
set deltmp = n # delete temporary xmecca files?
</syntaxhighlight>
</div>
The most parameters can be left as they are. Important is to rename the gas-equation-file and the replacement-file (if you have created one?).
Here is a list added with a short explanation of each parameter:
* <b>apn:</b>
* <b>gaseqnfile:</b>
* <b>rplfile:</b>
* <b>wanted:</b>
* <b>mcfct:</b>
* <b>diagtracfile:</b>
* <b>rxnrates:</b>
* <b>tagdbl:</b>
* <b>kppoption:</b>
* <b>integr:</b>
* <b>vlen:</b>
* <b>decomp:</b>
* <b>deltmpkp4:</b>
* <b>latex:</b>
* <b>graphviz:</b>
* <b>deltmp:</b>

If you're finished, don't forget to save everything (for Terminal users: to activate writing inside vi, use command <code>:w</code>, for closing and saving: <code>:wq</code>)!

=== Step 3: Execute the Mecca script ===
The first most important preparations are done now. Let's execute the <code>xmecca</code> script by typing
./xmecca
inside the <code>mecca</code> directory and follow the steps. Select your previously created batch file and run it. Now the Fortran files of the selected mechanism are created.

=== Step4: Transfer MECCA mechanism to ICON ===
Since the created Fortran code is only located inside Mecca and not in ICON-ART so far, a transfer has to be carried out. A script that performs this transfer can be obtained via
git clone https://gitlab.dkrz.de/art/mecca preproc.git
Now a new directory <code>Mecca_preproc</code> has been generated where the script <code>create_icon_code4.sh</code> can be found. By executing
./create_icon_code4.sh -h
the following view will be provided:
usage: ./create_icon_code4.sh [-m <mecca_home>] [-i <icon_home>] [-o <file_name>] [-h]
-m <mecca_home> : path to Mecca home directory, e.g. '~/caaba_3.0/mecca'
-i <icon_home> : path to icon home directory, e.g. '~/icon-kit/src/ICON-ART'
-o <file_name> : name of the XML file for ICON-ART, e.g. 'tracers_full.xml'
-h : display help

Here all paths to the Mecca- and ICON home directories have to be provided (type <code>pwd</code> in terminal to find out the path of your current directory) as well as a name for the XML-file that is going to be linked in the runscript later.

Now the Mecca-XML-file is generated and can be found in ICON in the directory <code>/icon home>/runctrl examples/xml ctrl</code>.

== Involving Mecca chemistry in a ICON-ART simulation ==
Now, in the respective runscript the namelist parameter <code>lart_mecca</code> has be set to <code>.TRUE</code> and for <code>cart_mecca_xml</code> the path to the Mecca file can be provided.
'''Important before executing the simulation:''' As a final step, the ICON code has to be recompiled with the command
./config/dkrz/levante.intel --enable-art --enable-ecrad
and executed afterwards with the command
./make -j8

Afterwards the respective simulation can be executed, e.g. as explained in [[Lifetime Tracer Simulation]].
Double check all filled in paths and namelist - especially the ART-namelists. If every namelist parameter in the runscript is filled in correctly, the runscript has to be saved. Afterwards by typing
./exp.testsuite.meccasim_test.run
a job can be submitted to the respective HPC-System. Type the terminal command
squeue
to view a list of your submitted and currently running and jobs.
By changing in the output directory you can check the slurm file for possible errors and run times after your job has been run through. In the output directory you can also find all output data for postprocessing in netCDF format.

Lifetime Tracer Simulation

2023-09-07T08:07:02Z

Luca R: set up syntaxhiglight

- under construction -

In this example it is shown how to simulate a lifetime driven tracer with simplified chemistry in ICON-ART. This tutorial teaches you...
* the setup of the runscript
* the correct setup of the xml-file for such tracers.
* the structure of lifetime tracers in xml-files
* editing output variables
No emission data will be included in this simulation.

== Configuration case ==
The depicted case is dealing with biogenic very short-lived species (VSLS) with a very short lifetime, more accurately Bromoform (<chem>CHBr3</chem>) and Dibromomethane (<chem>CH2Br2</chem>). Like of most VSLSs the major source of <chem>CHBr3</chem> and <chem>CH2Br2</chem> is the ocean which leads too a large gradient with increasing height in the concentration of these tracers. Bromoform is mainly depleted by photolysis in the Troposphere whereas Dibromomethanes main loss is due to Hydroxylradicals (OH). To assess the ability to simulate the transport of VSLS from the surface to the lower stratosphere, this case study uses an idealized chemical tracer approach.

The simulation is modeling the 01 October 2012 and is initialized with data from the ECMWF Integrated Forecast System (IFS) and includes boundary conditions and chemical lifetimes from the WMO Ozone assessment 2010. The boundary conditions and the chemical lifetimes are recalculated in a sort of rate at which the substances are depleted from the atmosphere with help of the implicit solution of the balance equation:
<math>\frac{\partial(\bar{\rho)\hat{\Psi_{l}}/{\partial t} = -\nabla\cdot (\hat{v}\bar{\rho}\hat{\Psi_{l}})- \nabla\cdot\bar{(\rho v''\Psi_{g,l}''}+P_l-L_l+E_l</math>

Here <math>\rho</math> is the density of air, <math>\hat{\Psi_{l}}</math> the barycentric-averaged mass mixing ratio, <math>\nabla\cdot (\hat{v}\bar{\rho}\hat{\Psi_{l}})</math> indicates the flux divergence that includes the horizontal and vertical advection of the gaseous compound l and <math>\nabla\cdot\bar{(\rho v''\Psi_{l}'')}</math> indicates the change due to turbulent fluxes. Further <math>P_l</math> describes the production rate due to chemical reactions, <math>L_l</math> the respective loss rate and emissions are noted with <math>E_l</math>. Everything is related to the respective compound l.

In this example no emission data is used.

== Setting up the Runscript ==

Let's start with the runscript that has to be prepared. Please note that the in the following explained parts have to be printed in one runscript-file with the naming designation "xyz.run". Here it is named <code>exp.testsuite.lifetime_tracer_test.run</code> but of course you can call it differently as well. The runscript can be stored under the following path in your icon directory: <code>/icon-kit/run</code>

If you've already walked through the example of the [[Simplified Chemistry]], you can use nearly the same runscript. Note that you have to change the paths from Part 1, the timing settings from Part 2, the output variables from Part 3, the emission settings as well as the path of the chemtracer-xml-file in the ART-settings from Part 4 and finally the timing in the job settings from Part 5. Details can be found below.

Inside of that, first check in part 1 that all your paths to your directories are correct, probably they have to be adjusted. Note that "hp8526" is a name of a specific account here, make sure to double check especially these lines. Abbreviations used here are the following:
*CENTER: Your organization
*EXPNAME: name of your ICON-Simulation
*OUTDIR: Directory where the simulation output will be stored
*ARTFOLDER: Directory where the ICON-ART code is stored
*INDIR: Directory where the necessary Input data are stored
*EXP:
*lart: For ICON-ART Simulation that has to be switched to <code>.True.</code>.

<div class="toccolours mw-collapsible mw-collapsed">
Part 1: Runscript Directory Settings (Example configuration)
<syntaxhighlight lang=bash line class="mw-collapsible-content">
#!/bin/bash
CENTER=IMK
workspace=/hkfs/work/workspace/scratch/hp8526-liftime_tracer_test
basedir=${workspace}/icon-kit-testsuite
icon_data_poolFolder=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT/AMIP/amip_input
EXPNAME=liftime_tracer_test
OUTDIR=${workspace}/output/${EXPNAME}
ICONFOLDER=/home/hk-project-iconart/hp8526/icon-kit
ARTFOLDER=${ICONFOLDER}/externals/art
INDIR=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT
EXP=ECHAM_AMIP_LIFETIME
lart=.True.

FILETYPE=4
COMPILER=intel
restart=.False.
read_restart_namelists=.False.

# Remove folder from OUTDIR for postprocessing output
OUTDIR_PREFIX=${workspace}

# Create output directory and go to this directory

if [ ! -d $OUTDIR ]; then
mkdir -p $OUTDIR
fi

cd $OUTDIR

#input for global domain
ln -sf ${INDIR}/../INPUT/GRID/icon_grid_0014_R02B05_G.nc iconR2B05_DOM01.nc
ln -sf ${INDIR}/../INPUT/EXTPAR/icon_extpar_0014_R02B05_G.nc extpar_iconR2B05_DOM01.nc
ln -sf /hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/output/0014_R02B05/uc1_ifs_t1279_grb2_remap_rev832_0014_R02B05_2018010100.nc ifs2icon_R2B05_DOM01.nc

ln -sf $ICONFOLDER/data/rrtmg_lw.nc rrtmg_lw.nc
ln -sf $ICONFOLDER/data/ECHAM6_CldOptProps.nc ECHAM6_CldOptProps.nc

ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/mozart_coord.nc ${OUTDIR}/mozart_coord.nc
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Linoz2004Br.dat ${OUTDIR}/Linoz2004Br.dat
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Simnoy2002.dat ${OUTDIR}/Simnoy2002.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-aer.dat ${OUTDIR}/FJX_scat-aer.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j.dat ${OUTDIR}/FJX_j2j.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-cld.dat ${OUTDIR}/FJX_scat-cld.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-ssa.dat ${OUTDIR}/FJX_scat-ssa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-UMa.dat ${OUTDIR}/FJX_scat-UMa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec.dat ${OUTDIR}/FJX_spec.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_std.dat ${OUTDIR}/atmos_std.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_h2och4.dat ${OUTDIR}/atmos_h2och4.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j_extended.dat ${OUTDIR}/FJX_j2j_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended_lyman.dat ${OUTDIR}/FJX_spec_extended_lyman.dat

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# the namelist filename
atmo_namelist=NAMELIST_${EXPNAME}
#
#-----------------------------------------------------------------------------
</syntaxhighlight>
</div>

Additionally in the next lines of code you set the timing. In this simulation we only simulate one day (01 October 2012). To calculate <chem>CHBr3</chem> and <chem>CH2Br2</chem> to every time of the day the output interval is set to 1 hour.
<div class="toccolours mw-collapsible mw-collapsed">
Part 2: Runscript Timing Settings (Example configuration)
<syntaxhighlight lang=bash line class="mw-collapsible-content">
! run_nml: general switches ----------
&# model timing

start_date=${start_date:="2012-10-01T00:00:00Z"}
end_date=${end_date:="2012-10-02T00:00:00Z"}
output_start=${start_date:="2012-10-01T00:00:00Z"}
output_end=${end_date:="2012-10-02T00:00:00Z"}
output_interval="PT1H"
modelTimeStep="PT6M"
leadtime="P8H"
checkpoint_interval="P30D"
</syntaxhihlight>
</div>

Further, all the namelist parameters (from the regular ICON model without ART-extension) have to be set. For a regular ICON-ART-Simulation the following settings are recommended - if not stated differently. For a detailed description, check out the ICON Documentation ([https://code.mpimet.mpg.de/attachments/download/19568/ICON_tutorial_2019.pdf Drill et. al. (2019)]).

<div class="toccolours mw-collapsible mw-collapsed">
Part 3: Runscript ICON-Parameter and -Namelist Settings (Example configuration)
<syntaxhighlight lang=bash line class="mw-collapsible-content">
# model parameters
model_equations=3 # equation system
# 1=hydrost. atm. T
# 1=hydrost. atm. theta dp
# 3=non-hydrost. atm.,
# 0=shallow water model
# -1=hydrost. ocean
#-----------------------------------------------------------------------------
# the grid parameters
declare -a atmo_dyn_grids=("iconR2B04_DOM01.nc" "iconR2B05_DOM02.nc" "iconR2B06_DOM03.nc")
# "iconR2B08_DOM03.nc"
#atmo_dyn_grids="iconR2B07_DOM02.nc","iconR2B08_DOM03.nc"
#atmo_rad_grids="iconR2B06_DOM01.nc"
#-----------------------------------------------------------------------------

# create ICON master namelist
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf

no_of_models=1

cat > icon_master.namelist << EOF
&master_nml
lRestart = .false.
/
&master_time_control_nml
experimentStartDate = "$start_date"
experimentStopDate = "$end_date"
forecastLeadTime = "$leadtime"
checkpointTimeIntval = "$checkpoint_interval"
/
&master_model_nml
model_type=1
model_name="ATMO"
model_namelist_filename="$atmo_namelist"
model_min_rank=1
model_max_rank=65536
model_inc_rank=1
/
EOF

#-----------------------------------------------------------------------------
#

#-----------------------------------------------------------------------------
#
# write ICON namelist parameters
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf
#
# ------------------------
# reconstrcuct the grid parameters in namelist form
dynamics_grid_filename=""
for gridfile in ${atmo_dyn_grids}; do
dynamics_grid_filename="${dynamics_grid_filename} '${gridfile}',"
done
radiation_grid_filename=""
for gridfile in ${atmo_rad_grids}; do
radiation_grid_filename="${radiation_grid_filename} '${gridfile}',"
done

# ------------------------

cat > ${atmo_namelist} << EOF
&parallel_nml
nproma = 8 ! optimal setting 8 for CRAY; use 16 or 24 for IBM
p_test_run = .false.
l_test_openmp = .false.
l_log_checks = .false.
num_io_procs = 0 ! up to one PE per output stream is possible
itype_comm = 1
iorder_sendrecv = 3 ! best value for CRAY (slightly faster than option 1)
/
&grid_nml
dynamics_grid_filename = 'iconR2B05_DOM01.nc'
dynamics_parent_grid_id = 0
!radiation_grid_filename = ${radiation_grid_filename}
lredgrid_phys = .false.
lfeedback = .true.
ifeedback_type = 2
/
&initicon_nml
lconsistency_checks = .false.
init_mode = 2 ! operation mode 2: IFS
zpbl1 = 500.
zpbl2 = 1000.
! l_sst_in = .true.
/
&run_nml
num_lev = 90
lvert_nest = .true. ! use vertical nesting if a nest is active
! nsteps = ${nsteps} ! 50 ! 1200 ! 7200 !
! dtime = ${dtime} ! timestep in seconds
modelTimeStep = "${modelTimeStep}"
ldynamics = .TRUE. ! dynamics
ltransport = .true.
iforcing = 3 ! NWP forcing
ltestcase = .false. ! false: run with real data
msg_level = 7 ! print maximum wind speeds every 5 time steps
ltimer = .true. ! set .TRUE. for timer output
timers_level = 10 ! can be increased up to 10 for detailed timer output
output = "nml"
lart = ${lart}
/
&nwp_phy_nml
inwp_gscp = 1
inwp_convection = 1
inwp_radiation = 1
inwp_cldcover = 1
inwp_turb = 1
inwp_satad = 1
inwp_sso = 1
inwp_gwd = 1
inwp_surface = 1
icapdcycl = 3 ! apply CAPE modification to improve diurnalcycle over tropical land (optimizes NWP scores)
latm_above_top = .false., .true. ! the second entry refers to the nested domain (if present)
efdt_min_raylfric = 7200.
itype_z0 = 2
icpl_aero_conv = 1
icpl_aero_gscp = 1
! resolution-dependent settings - please choose the appropriate one
dt_rad = 2160.
dt_conv = 720.
dt_sso = 1440.
dt_gwd = 1440.
/
&nwp_tuning_nml
tune_zceff_min = 0.075 ! ** default value to be used for R3B7; use 0.05 for R2B6 in order to get similar temperature biases in upper troposphere **
itune_albedo = 1 ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara in order to reduce cold bias
/
&turbdiff_nml
tkhmin = 0.75 ! new default since rev. 16527
tkmmin = 0.75 ! ""
pat_len = 100.
c_diff = 0.2
rat_sea = 8.5 ! ** new since r20191: 8.5 for R3B7, 8.0 for R2B6 in order to get similar temperature biases in the tropics **
ltkesso = .true.
frcsmot = 0.2 ! these 2 switches together apply vertical smoothing of the TKE source terms
imode_frcsmot = 2 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects:
itype_sher = 3
ltkeshs = .true.
a_hshr = 2.0
/
&lnd_nml
ntiles = 3 !!! 1 for assimilation cycle and forecast
nlev_snow = 3 !!! 1 for assimilation cycle and forecast
lmulti_snow = .true. !!! .false. for assimilation cycle and forecast
itype_heatcond = 2
idiag_snowfrac = 2
lsnowtile = .false. !! later on .true. if GRIB encoding issues are solved
lseaice = .true.
llake = .false.
itype_lndtbl = 3 ! minimizes moist/cold bias in lower tropical troposphere
itype_root = 2
/
&radiation_nml
irad_o3 = 10
irad_aero = 6
albedo_type = 2 ! Modis albedo
vmr_co2 = 390.e-06 ! values representative for 2012
vmr_ch4 = 1800.e-09
vmr_n2o = 322.0e-09
vmr_o2 = 0.20946
vmr_cfc11 = 240.e-12
vmr_cfc12 = 532.e-12
/
&nonhydrostatic_nml
iadv_rhotheta = 2
ivctype = 2
itime_scheme = 4
exner_expol = 0.333
vwind_offctr = 0.2
damp_height = 50000.
rayleigh_coeff = 0.10
lhdiff_rcf = .true.
divdamp_order = 24 ! for data assimilation runs, '2' provides extra-strong filtering of gravity waves
divdamp_type = 32 !!! optional: 2 for assimilation cycle if very strong gravity-wave filtering is desired
divdamp_fac = 0.004
l_open_ubc = .false.
igradp_method = 3
l_zdiffu_t = .true.
thslp_zdiffu = 0.02
thhgtd_zdiffu = 125.
htop_moist_proc= 22500.
hbot_qvsubstep = 22500. ! use 19000. with R3B7
/
&sleve_nml
min_lay_thckn = 20.
max_lay_thckn = 400. ! maximum layer thickness below htop_thcknlimit
htop_thcknlimit = 14000. ! this implies that the upcoming COSMO-EU nest will have 60 levels
top_height = 75000.
stretch_fac = 0.9
decay_scale_1 = 4000.
decay_scale_2 = 2500.
decay_exp = 1.2
flat_height = 16000.
/
&dynamics_nml
iequations = 3
idiv_method = 1
divavg_cntrwgt = 0.50
lcoriolis = .TRUE.
/
&transport_nml
! qv, qc, qi, qr, qs
itype_vlimit = 1,1,1,1,1
ivadv_tracer = 3, 3, 3, 3, 3
itype_hlimit = 3, 4, 4, 4 , 4
ihadv_tracer = 52, 2,2,2,2
iadv_tke = 0
/
&diffusion_nml
hdiff_order = 5
itype_vn_diffu = 1
itype_t_diffu = 2
hdiff_efdt_ratio = 24.0
hdiff_smag_fac = 0.025
lhdiff_vn = .TRUE.
lhdiff_temp = .TRUE.
/
&interpol_nml
nudge_zone_width = 8
lsq_high_ord = 3
l_intp_c2l = .true.
l_mono_c2l = .true.
support_baryctr_intp = .false.
/
&extpar_nml
itopo = 1
n_iter_smooth_topo = 1
heightdiff_threshold = 3000.
/
&io_nml
itype_pres_msl = 4 ! IFS method with bug fix for self-consistency between SLP and geopotential
itype_rh = 1 ! RH w.r.t. water
/
&output_nml
filetype = ${FILETYPE} ! output format: 2=GRIB2, 4=NETCDFv2
dom = -1 ! write all domains
output_start = "${output_start}"
output_end = "${output_end}"
output_interval = "${output_interval}"
steps_per_file = 1
include_last = .TRUE.
output_filename = 'icon-art-${EXPNAME}-chem' ! file name base
ml_varlist = 'temp','pres','group:ART_CHEMISTRY','u','v'
output_grid = .TRUE.
remap = 1
reg_lon_def = -180.,1,180.
reg_lat_def = -90.,1,90.
/
</syntaxhighlight>
</div>
Please note in the last namelist section "output_nml" that you can set all output variables that you need to postprocess your data later. All assigned variables here will be written in the output netCDF-files as well. To learn more about post processing your data, check out a later chapter of this article or the [[Postprocessing]] article.

Now, we're getting to the ICON-ART settings. To enable chemistry in an ICON-ART Simulation in general, the switch <code>lart_chem</code> has to be set to <code>.TRUE.</code>. With <code>lart_diag_out</code> output of the diagnostic fields can be enabled. Due to setting <code>lart_chem=.TRUE.</code> either <code>lart_chemtracer</code> or <code>lart_mecca</code> has to be set to <code>.TRUE.</code>. Because we want to perform a simulation with simplified chemistry, we have to switch on <code>lart_chemtracer</code>. If this namelist parameter is set to <code>.TRUE.</code>, also <code>cart_chemtracer_xml</code> has to be fulfilled. Here you enter the path of your xml-file which describes the tracers occurring and their properties in the simulation. How to create this xml-file is explained in the next chapter.
An example configuration for this part is shown in the following:
<div class="toccolours mw-collapsible mw-collapsed">
Part 4: Runscript ICON-ART Settings (Example configuration)
<syntaxhighlight lang=bash line class="mw-collapsible-content">
&art_nml
lart_chem = .TRUE.
lart_diag_out = .TRUE.
lart_aerosol = .FALSE.
lart_mecca = .FALSE.
lart_chemtracer = .TRUE.
cart_chemtracer_xml = '/home/hk-project-iconart/hp8526/icon-kit/externals/art/runctrl_examples/xml_ctrl/chemtracer_lifetime_test.xml'
cart_input_folder = '${OUTDIR}'
cart_io_suffix = '0014'
iart_init_gas = 0
/
EOF
</syntaxhighlight>
</div>
Please note that there are also several other namelist parameter you can select from (see [[Namelist]] article) but to perform our case study we're done for the ART setting at this point.

Depending on the used HPC-System, some parameter concerning the running job like maximum running time and used nodes can be set. For this case study the following settings can be copied. Note that this is valid for the HoreKa HPC system and that it can differ to other systems.
<div class="toccolours mw-collapsible mw-collapsed">
Part 5: Runscript job Settings (Example configuration)
<syntaxhighlight lang=bash line class="mw-collapsible-content">
cp ${ICONFOLDER}/bin/icon ./icon.exe

cat > job_ICON << ENDFILE
#!/bin/bash -x
#SBATCH --nodes=2
#SBATCH --time=06:00:00
#SBATCH --ntasks-per-node=76
#SBATCH --partition=cpuonly
#SBATCH -A hk-project-iconart
###SBATCH --constraint=LSDF

module load compiler/intel/2022.0.2 mpi/openmpi/4.0 lib/netcdf/4.9_serial lib/hdf5/1.12_serial lib/netcdf-fortran/4.5_serial lib/eccodes/2.25.0 numlib/mkl/2022.0.2

mpirun --bind-to core --map-by core --report-bindings ./icon.exe

ENDFILE

chmod +x job_ICON
sbatch job_ICON
</syntaxhighlight>
</div>

To conclude and to double check, in the following box the complete runscript is shown once again.

<div class="toccolours mw-collapsible mw-collapsed">
Complete example configuration of the runscript
<syntaxhighlight lang=bash line class="mw-collapsible-content">
#!/bin/bash
CENTER=IMK
workspace=/hkfs/work/workspace/scratch/hp8526-liftime_tracer_test
basedir=${workspace}/icon-kit-testsuite
icon_data_poolFolder=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT/AMIP/amip_input
EXPNAME=liftime_tracer_test
OUTDIR=${workspace}/output/${EXPNAME}
ICONFOLDER=/home/hk-project-iconart/hp8526/icon-kit
ARTFOLDER=${ICONFOLDER}/externals/art
INDIR=/hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/INPUT
EXP=ECHAM_AMIP_LIFETIME
lart=.True.

FILETYPE=4
COMPILER=intel
restart=.False.
read_restart_namelists=.False.

# Remove folder from OUTDIR for postprocessing output
OUTDIR_PREFIX=${workspace}

# Create output directory and go to this directory

if [ ! -d $OUTDIR ]; then
mkdir -p $OUTDIR
fi

cd $OUTDIR

#input for global domain
ln -sf ${INDIR}/../INPUT/GRID/icon_grid_0014_R02B05_G.nc iconR2B05_DOM01.nc
ln -sf ${INDIR}/../INPUT/EXTPAR/icon_extpar_0014_R02B05_G.nc extpar_iconR2B05_DOM01.nc
ln -sf /hkfs/work/workspace/scratch/fb4738-dwd_ozone/icon/output/0014_R02B05/uc1_ifs_t1279_grb2_remap_rev832_0014_R02B05_2018010100.nc ifs2icon_R2B05_DOM01.nc

ln -sf $ICONFOLDER/data/rrtmg_lw.nc rrtmg_lw.nc
ln -sf $ICONFOLDER/data/ECHAM6_CldOptProps.nc ECHAM6_CldOptProps.nc

ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/mozart_coord.nc ${OUTDIR}/mozart_coord.nc
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Linoz2004Br.dat ${OUTDIR}/Linoz2004Br.dat
ln -sf ${ARTFOLDER}/runctrl_examples/init_ctrl/Simnoy2002.dat ${OUTDIR}/Simnoy2002.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-aer.dat ${OUTDIR}/FJX_scat-aer.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j.dat ${OUTDIR}/FJX_j2j.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-cld.dat ${OUTDIR}/FJX_scat-cld.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-ssa.dat ${OUTDIR}/FJX_scat-ssa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_scat-UMa.dat ${OUTDIR}/FJX_scat-UMa.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec.dat ${OUTDIR}/FJX_spec.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_std.dat ${OUTDIR}/atmos_std.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/atmos_h2och4.dat ${OUTDIR}/atmos_h2och4.dat

ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_j2j_extended.dat ${OUTDIR}/FJX_j2j_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended.dat ${OUTDIR}/FJX_spec_extended.dat
ln -sf ${ARTFOLDER}/runctrl_examples/photo_ctrl/FJX_spec_extended_lyman.dat ${OUTDIR}/FJX_spec_extended_lyman.dat

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# this if condition is necessary because otherwise
# a new link in ${INDIR}/${EXP}/emiss_minimal is generated
# linking to itself
if [ ! -L ${OUTDIR}/emissions ]; then
ln -sd /home/hk-project-iconart/hp8526/emissions ${OUTDIR}/emissions
fi

# the namelist filename
atmo_namelist=NAMELIST_${EXPNAME}
#
#-----------------------------------------------------------------------------
! run_nml: general switches ----------
&# model timing

start_date=${start_date:="2012-10-01T00:00:00Z"}
end_date=${end_date:="2012-10-02T00:00:00Z"}
output_start=${start_date:="2012-10-01T00:00:00Z"}
output_end=${end_date:="2012-10-02T00:00:00Z"}
output_interval="PT1H"
modelTimeStep="PT6M"
leadtime="P8H"
checkpoint_interval="P30D"

# model parameters
model_equations=3 # equation system
# 1=hydrost. atm. T
# 1=hydrost. atm. theta dp
# 3=non-hydrost. atm.,
# 0=shallow water model
# -1=hydrost. ocean
#-----------------------------------------------------------------------------
# the grid parameters
declare -a atmo_dyn_grids=("iconR2B04_DOM01.nc" "iconR2B05_DOM02.nc" "iconR2B06_DOM03.nc")
# "iconR2B08_DOM03.nc"
#atmo_dyn_grids="iconR2B07_DOM02.nc","iconR2B08_DOM03.nc"
#atmo_rad_grids="iconR2B06_DOM01.nc"
#-----------------------------------------------------------------------------

# create ICON master namelist
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf

no_of_models=1

cat > icon_master.namelist << EOF
&master_nml
lRestart = .false.
/
&master_time_control_nml
experimentStartDate = "$start_date"
experimentStopDate = "$end_date"
forecastLeadTime = "$leadtime"
checkpointTimeIntval = "$checkpoint_interval"
/
&master_model_nml
model_type=1
model_name="ATMO"
model_namelist_filename="$atmo_namelist"
model_min_rank=1
model_max_rank=65536
model_inc_rank=1
/
EOF

#-----------------------------------------------------------------------------
#

#-----------------------------------------------------------------------------
#
# write ICON namelist parameters
# ------------------------
# For a complete list see Namelist_overview and Namelist_overview.pdf
#
# ------------------------
# reconstrcuct the grid parameters in namelist form
dynamics_grid_filename=""
for gridfile in ${atmo_dyn_grids}; do
dynamics_grid_filename="${dynamics_grid_filename} '${gridfile}',"
done
radiation_grid_filename=""
for gridfile in ${atmo_rad_grids}; do
radiation_grid_filename="${radiation_grid_filename} '${gridfile}',"
done

# ------------------------

cat > ${atmo_namelist} << EOF
&parallel_nml
nproma = 8 ! optimal setting 8 for CRAY; use 16 or 24 for IBM
p_test_run = .false.
l_test_openmp = .false.
l_log_checks = .false.
num_io_procs = 0 ! up to one PE per output stream is possible
itype_comm = 1
iorder_sendrecv = 3 ! best value for CRAY (slightly faster than option 1)
/
&grid_nml
dynamics_grid_filename = 'iconR2B05_DOM01.nc'
dynamics_parent_grid_id = 0
!radiation_grid_filename = ${radiation_grid_filename}
lredgrid_phys = .false.
lfeedback = .true.
ifeedback_type = 2
/
&initicon_nml
lconsistency_checks = .false.
init_mode = 2 ! operation mode 2: IFS
zpbl1 = 500.
zpbl2 = 1000.
! l_sst_in = .true.
/
&run_nml
num_lev = 90
lvert_nest = .true. ! use vertical nesting if a nest is active
! nsteps = ${nsteps} ! 50 ! 1200 ! 7200 !
! dtime = ${dtime} ! timestep in seconds
modelTimeStep = "${modelTimeStep}"
ldynamics = .TRUE. ! dynamics
ltransport = .true.
iforcing = 3 ! NWP forcing
ltestcase = .false. ! false: run with real data
msg_level = 7 ! print maximum wind speeds every 5 time steps
ltimer = .true. ! set .TRUE. for timer output
timers_level = 10 ! can be increased up to 10 for detailed timer output
output = "nml"
lart = ${lart}
/
&nwp_phy_nml
inwp_gscp = 1
inwp_convection = 1
inwp_radiation = 1
inwp_cldcover = 1
inwp_turb = 1
inwp_satad = 1
inwp_sso = 1
inwp_gwd = 1
inwp_surface = 1
icapdcycl = 3 ! apply CAPE modification to improve diurnalcycle over tropical land (optimizes NWP scores)
latm_above_top = .false., .true. ! the second entry refers to the nested domain (if present)
efdt_min_raylfric = 7200.
itype_z0 = 2
icpl_aero_conv = 1
icpl_aero_gscp = 1
! resolution-dependent settings - please choose the appropriate one
dt_rad = 2160.
dt_conv = 720.
dt_sso = 1440.
dt_gwd = 1440.
/
&nwp_tuning_nml
tune_zceff_min = 0.075 ! ** default value to be used for R3B7; use 0.05 for R2B6 in order to get similar temperature biases in upper troposphere **
itune_albedo = 1 ! somewhat reduced albedo (w.r.t. MODIS data) over Sahara in order to reduce cold bias
/
&turbdiff_nml
tkhmin = 0.75 ! new default since rev. 16527
tkmmin = 0.75 ! ""
pat_len = 100.
c_diff = 0.2
rat_sea = 8.5 ! ** new since r20191: 8.5 for R3B7, 8.0 for R2B6 in order to get similar temperature biases in the tropics **
ltkesso = .true.
frcsmot = 0.2 ! these 2 switches together apply vertical smoothing of the TKE source terms
imode_frcsmot = 2 ! in the tropics (only), which reduces the moist bias in the tropical lower troposphere
! use horizontal shear production terms with 1/SQRT(Ri) scaling to prevent unwanted side effects:
itype_sher = 3
ltkeshs = .true.
a_hshr = 2.0
/
&lnd_nml
ntiles = 3 !!! 1 for assimilation cycle and forecast
nlev_snow = 3 !!! 1 for assimilation cycle and forecast
lmulti_snow = .true. !!! .false. for assimilation cycle and forecast
itype_heatcond = 2
idiag_snowfrac = 2
lsnowtile = .false. !! later on .true. if GRIB encoding issues are solved
lseaice = .true.
llake = .false.
itype_lndtbl = 3 ! minimizes moist/cold bias in lower tropical troposphere
itype_root = 2
/
&radiation_nml
irad_o3 = 10
irad_aero = 6
albedo_type = 2 ! Modis albedo
vmr_co2 = 390.e-06 ! values representative for 2012
vmr_ch4 = 1800.e-09
vmr_n2o = 322.0e-09
vmr_o2 = 0.20946
vmr_cfc11 = 240.e-12
vmr_cfc12 = 532.e-12
/
&nonhydrostatic_nml
iadv_rhotheta = 2
ivctype = 2
itime_scheme = 4
exner_expol = 0.333
vwind_offctr = 0.2
damp_height = 50000.
rayleigh_coeff = 0.10
lhdiff_rcf = .true.
divdamp_order = 24 ! for data assimilation runs, '2' provides extra-strong filtering of gravity waves
divdamp_type = 32 !!! optional: 2 for assimilation cycle if very strong gravity-wave filtering is desired
divdamp_fac = 0.004
l_open_ubc = .false.
igradp_method = 3
l_zdiffu_t = .true.
thslp_zdiffu = 0.02
thhgtd_zdiffu = 125.
htop_moist_proc= 22500.
hbot_qvsubstep = 22500. ! use 19000. with R3B7
/
&sleve_nml
min_lay_thckn = 20.
max_lay_thckn = 400. ! maximum layer thickness below htop_thcknlimit
htop_thcknlimit = 14000. ! this implies that the upcoming COSMO-EU nest will have 60 levels
top_height = 75000.
stretch_fac = 0.9
decay_scale_1 = 4000.
decay_scale_2 = 2500.
decay_exp = 1.2
flat_height = 16000.
/
&dynamics_nml
iequations = 3
idiv_method = 1
divavg_cntrwgt = 0.50
lcoriolis = .TRUE.
/
&transport_nml
! qv, qc, qi, qr, qs
itype_vlimit = 1,1,1,1,1
ivadv_tracer = 3, 3, 3, 3, 3
itype_hlimit = 3, 4, 4, 4 , 4
ihadv_tracer = 52, 2,2,2,2
iadv_tke = 0
/
&diffusion_nml
hdiff_order = 5
itype_vn_diffu = 1
itype_t_diffu = 2
hdiff_efdt_ratio = 24.0
hdiff_smag_fac = 0.025
lhdiff_vn = .TRUE.
lhdiff_temp = .TRUE.
/
&interpol_nml
nudge_zone_width = 8
lsq_high_ord = 3
l_intp_c2l = .true.
l_mono_c2l = .true.
support_baryctr_intp = .false.
/
&extpar_nml
itopo = 1
n_iter_smooth_topo = 1
heightdiff_threshold = 3000.
/
&io_nml
itype_pres_msl = 4 ! IFS method with bug fix for self-consistency between SLP and geopotential
itype_rh = 1 ! RH w.r.t. water
/
&output_nml
filetype = ${FILETYPE} ! output format: 2=GRIB2, 4=NETCDFv2
dom = -1 ! write all domains
output_start = "${output_start}"
output_end = "${output_end}"
output_interval = "${output_interval}"
steps_per_file = 1
include_last = .TRUE.
output_filename = 'icon-art-${EXPNAME}-chem' ! file name base
ml_varlist = 'temp','pres','group:ART_CHEMISTRY','u','v'
output_grid = .TRUE.
remap = 1
reg_lon_def = -180.,1,180.
reg_lat_def = -90.,1,90.
/

&art_nml
lart_chem = .TRUE.
lart_diag_out = .TRUE.
lart_aerosol = .FALSE.
lart_mecca = .FALSE.
lart_chemtracer = .TRUE.
cart_chemtracer_xml = '/home/hk-project-iconart/hp8526/icon-kit/externals/art/runctrl_examples/xml_ctrl/chemtracer_lifetime_test.xml'
cart_input_folder = '${OUTDIR}'
cart_io_suffix = '0014'
iart_init_gas = 0
/
EOF

cp ${ICONFOLDER}/bin/icon ./icon.exe

cat > job_ICON << ENDFILE
#!/bin/bash -x
#SBATCH --nodes=2
#SBATCH --time=06:00:00
#SBATCH --ntasks-per-node=76
#SBATCH --partition=cpuonly
#SBATCH -A hk-project-iconart
###SBATCH --constraint=LSDF

module load compiler/intel/2022.0.2 mpi/openmpi/4.0 lib/netcdf/4.9_serial lib/hdf5/1.12_serial lib/netcdf-fortran/4.5_serial lib/eccodes/2.25.0 numlib/mkl/2022.0.2

mpirun --bind-to core --map-by core --report-bindings ./icon.exe

ENDFILE

chmod +x job_ICON
sbatch job_ICON
</syntaxhighlight>
</div>

== Setting up the xml-file ==
An xml-file describes the chemical components of the simulation which means that all trace gases or aerosols and their properties that are relevant for the simulation are listed here. Since we perform a simulation with simplified ICON-ART chemistry of lifetime dependent tracers we need the matching chemtracer-xml-file where only these tracers have to mentioned. Then, the concentrations will be calculated according to their lifetimes with no relation to other tracers. Because of no specific other settings (e.g. emissions), this is the only xml-file needed here.

To prepare the xml-file we can select the matching species from the general previously generated xml-file <code>standard_chemtracer.xml</code>.

<div class="toccolours mw-collapsible mw-collapsed">
Standard Chemtracer-xml-file for simulations without emission data
<syntaxhighlight lang=xml line class="mw-collapsible-content">
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE tracers SYSTEM "tracers.dtd">

<tracers>
<chemtracer id="TRCH4" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.604E-2</mol_weight>
<?source_lifetime Hayman et al., ACP, 2014 ?>
<lifetime type="real">286977600</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRCO;2.*TRH2O</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">CH4</init_name>
</chemtracer>
<chemtracer id="TRC2H6" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">3.006E-2</mol_weight>
<?source_lifetime Hodnebrog et al, Atmos. Sci. Lett., 2018 ?>
<lifetime type="real">5011200</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRC3H8" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.40956E-2</mol_weight>
<?source_lifetime Hodnebrog et al, Atmos. Sci. Lett., 2018 ?>
<lifetime type="real">1123200</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">0.736*TRCH3COCH3</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRC5H8" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.812E-2</mol_weight>
<?source_lifetime Weimer (2015), p. 16 ?>
<lifetime type="real">8640</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCH3COCH3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">5.808E-2</mol_weight>
<?source_lifetime Weimer (2015), p. 9 ?>
<lifetime type="real">1728000</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCO" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">2.801E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">5184000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRCO2</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">CO</init_name>
</chemtracer>
<chemtracer id="TRCO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.401E-2</mol_weight>
<?source_lifetime Houghton et al., IPCC, Cambridge University Press, 2001 ?>
<lifetime type="real">3153600000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRCO</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRH2O" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.802E-2</mol_weight>
<?source_lifetime Just a placeholder (not used in the code) ?>
<lifetime type="real">2592000000</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">H2O</init_name>
</chemtracer>
<chemtracer id="TRCHBr3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">2.527E-1</mol_weight>
<?source_lifetime Rieger et al., GMD, 2015 ?>
<lifetime type="real">2073600</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCH2Br2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.738E-1</mol_weight>
<?source_lifetime Rieger et al., GMD, 2015 ?>
<lifetime type="real">10627200</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRNH3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.70E-2</mol_weight>
<?source_lifetime Pinder et al., GRL, 2008?>
<lifetime type="real">86400</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRNO2</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">NH3</init_name>
</chemtracer>
<chemtracer id="TRNO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.601E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">259200</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRHNO3</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">NO2</init_name>
</chemtracer>
<chemtracer id="TRHNO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.30E-2</mol_weight>
<?source_lifetime Day et al., ACP, 2008?>
<lifetime type="real">21600</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">HNO3</init_name>
</chemtracer>
<chemtracer id="TRSO2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.40E-2</mol_weight>
<?source_lifetime Von Glasow, Chemical Geology, 2009 ?>
<lifetime type="real">1209600</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRH2SO4</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">SO2</init_name>
</chemtracer>
<chemtracer id="TROCS" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.01E-2</mol_weight>
<?source_lifetime Ullwer (2017) ?>
<lifetime type="real">504576000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">TRSO2</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">OCS</init_name>
</chemtracer>
<chemtracer id="TRDMS" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.21E-2</mol_weight>
<?source_lifetime Ullwer (2017) ?>
<lifetime type="real">216000</lifetime>

<c_solve type="char">lt</c_solve>
<products type="char">0.993*TRSO2;0.007*TROCS</products>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">DMS</init_name>
</chemtracer>
<chemtracer id="TRH2SO4" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">9.80E-2</mol_weight>
<?source_lifetime Fiedler et al., ACP, 2005?>
<lifetime type="real">1800</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<dataset01 type="char" cvar_name="H2SO4" rbottom_height="10000" rupper_height="40000">CCMI-ETH_MPIC1.1</dataset01>
</chemtracer>
<chemtracer id="TRHCL" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">3.60E-2</mol_weight>
<lifetime type="real">259200</lifetime>

<c_solve type="char">param</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TR_cold">
<unit type="char">none</unit>
<mol_weight type="real">1.000E+0</mol_weight>
<lifetime type="real">432000</lifetime>

<c_solve type="char">cold</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRO3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.800E-2</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">1576800</lifetime>

<c_solve type="char">linoz</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">O3</init_name>
</chemtracer>
<chemtracer id="TRN2O" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">4.401E-2</mol_weight>
<?source_lifetime Wypych, 2017, Atlas of Material Damage?>
<lifetime type="real">4730400000</lifetime>

<c_solve type="char">simnoy</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
<init_mode type="int">0</init_mode>
<init_name type="char">N2O</init_name>
</chemtracer>
<chemtracer id="TRNOy" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">6.301E-1</mol_weight>
<?source_lifetime Ehhalt et al., IPCC, 2001, Chapter 4 ?>
<lifetime type="real">259200</lifetime>

<c_solve type="char">simnoy</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRAGE" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">none</unit>
<mol_weight type="real">1.000E+0</mol_weight>
<?source_lifetime Just a placeholder, not used in the code?>
<lifetime type="real">25920000</lifetime>

<c_solve type="char">passive</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
</tracers>
</syntaxhighlight>
</div>

After identifying the needed chemtracers (in our case <chem>CHBr3</chem> (in xml-file: <code>TRCHBr3</code>) and <chem>CH2Br2</chem> (in xml-file: <code>TRCH2Br2</code>)) we copy the needed code lines into our new xml-file <code>chemtracer_lifetime_chbr3_ch2br2.xml</code> needed for the simulation.
<div class="toccolours mw-collapsible mw-collapsed">
Chemtracer-xml-file for lifetime tracer simulation (Example configuration)
<syntaxhighlight lang=xml line class="mw-collapsible-content">
<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE tracers SYSTEM "tracers.dtd">

<tracers>
<chemtracer id="TRCHBr3" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">2.527E-1</mol_weight>
<?source_lifetime Rieger et al., GMD, 2015 ?>
<lifetime type="real">2073600</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
<chemtracer id="TRCH2Br2" full="FALSE" chemtr="TRUE">
<tag001 type="char">chemtr</tag001>
<unit type="char">mol mol-1</unit>
<mol_weight type="real">1.738E-1</mol_weight>
<?source_lifetime Rieger et al., GMD, 2015 ?>
<lifetime type="real">10627200</lifetime>

<c_solve type="char">lt</c_solve>

<lfeedback type='int'>0</lfeedback>
<transport type="char">stdaero</transport>
</chemtracer>
</tracers>
</syntaxhighlight>
</div>

== Running the simulation ==
Double check all filled in paths and namelist - especially the ART-namelists. If every namelist parameter in the runscript is filled in correctly, the runscript has to be saved. Afterwards by typing
./exp.testsuite.lifetime_tracer_test.run
a job can be submitted to the respective HPC-System. Type the terminal command
squeue
to view a list of your submitted and currently running and jobs.
By changing in the output directory (which is according to our runscript <code>/hkfs/work/workspace/scratch/hp8526-liftime_tracer_test</code> you can check the slurm file for possible errors and run times after your job has been run through.

In the output directory you can also find all output data for postprocessing in netCDF format.

Input

2023-09-07T06:57:41Z

Luca R: added syntax Highlighting for ruscripts

== Requirements for a Simulation ==

To run a Simulation with ICON-ART there are three main points to consider: Setting the namelist parameters, preparing the xmls, and optionally prepare additional Input data.

== Namelist Inputs ==

To run a simulation with ICON-ART installed the first thing to prepare is the runfile. It is usually best to start with an existing and working runfile and adapt parts as needed. The runfile contains the majority of parameters to run a simulation, like the length, timesteps and grids used, as well as a variety of options for methods and parametrisations used in the model. The runfile contents are then split up in several namelists that the model reads at the beginning of a simulation. An overview of the Namelist Parameters can be found in [[Namelist]].

To enable ART in an ICON simulation, the switch <code>lart = .TRUE.</code> has to be set in the section <code>&run_nml</code>. This is the global on/off switch for ICON-ART. This is how this could look like in the context of a runfile:
<syntaxhighlight lang=bash line>
! run_nml: general switches ----------
&run_nml
ltestcase = .FALSE.
num_lev = 50
ltransport = .TRUE.
.............
lart = .TRUE.
</syntaxhighlight>
Generally the first letters of a namelist switch refer to its type, the "l" in "lart" stands for logical, meaning it has to be either True or False. Here are some examples.

<div id="tab:vartypes">

{| class="wikitable" style="text-align:left;"
|+ Some namelist switches and their data types.
! namelist switch
! type
|-
| <code>lart</code>
| logical
|-
| <code>cart_aerosol_xml</code>
| character
|-
| <code>iart_init_aero </code>
| integer
|}

</div>
The namelist <code>&art_nml</code> is used for general options of the ART simulation. To run a certain kind of simulation the according switch has to be set to <code>.TRUE.</code>. For Example to include a point source the switch <code>cart_aerosol_xml</code> has to be set to <code>.TRUE.</code>.

== XML Inputs ==

This enables the inclusion of a .xml file containing additional information like location and strength of the point source. The table [[#tab:art_nml-params|below]] contains the most important <code>&art_nml</code> namelist parameters and additional namelist parameters required if they are set to <code>.TRUE.</code>.

<div id="tab:art_nml-params">

{| class="wikitable" style="text-align:left;"
|+ XML files and their namelist dependencies
! XML File
! Description
! Namelist parameter dependency
! Default
! Further Information
|---
| <code>cart_chemtracer_xml</code>
| Switch for simple OH chemistry
| <code>lart_chemtracer</code>
| .FALSE.
|[[#Chemistry Tracers|Chemistry Tracers]]
|-
| <code>cart_mecca_xml</code>
| Switch for kpp chemistry
| <code>lart_mecca</code>
| .FALSE.
|[[Atmospheric Chemistry]]
|-
| <code>cart_pntSrc_xml</code>
| Enables creation of point sources emitting given Aerosols at a given rate
| <code>lart_pntSrc</code>
| .FALSE.
|[[#Point Source|Point Source]]
|-
| <code>cart_aerosol_xml</code>
| Main switch for the treatment of atmospheric aerosols
| <code>lart_aerosol</code>
| .FALSE.
|[[#Aerosol Tracers|Aerosol Tracers]]
|-
| <code>cart_modes_xml</code>
| Main switch for the treatment of atmospheric aerosols
| <code>lart_aerosol</code>
| .FALSE.
|[[#Aerosol Modes|Aerosol Modes]]
|-
| <code>cart_diagnostics_xml</code>
| Enables diagnostic output fields
| <code>lart_diag_out</code>
| .FALSE.
| -
|-
| <code>cart_emiss_xml_file</code>
| XML File for emission metadata
| -
| -
|[[#Aerosol Emission|Aerosol Emission]]
|-
| <code>cart_ext_data_xml</code>
| XML File for metadata of datasets prescribing tracers
| -
| -
| -
|-
| <code> cart_coag_xml </code>
| XML File containing additional information about coagulation
| -
| -
| -
|}

</div>

{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Aerosol Tracers">Aerosol Tracers</span>
|-
! tag
! type
! options
! description
|-
| htop_proc
| real
| in m
| top height for processes
|-
| initc
| character
| file
| initialize from Input file
|-
| inucl
| integer
| 0 (off), 1 (on)
| H2SO4 nucleation for so4 tracer (default=1); 1 for so4_sol_ait, 0 for other so4 tracer)
|-
| label
| character
| e.g., dusta
| allows to name tracers individually
|-
| latbc
| character
| file
| read data for LBC
|-
| lfeedback
| integer
| 0 (off), 1 (on)
| child -> parent feedback in nested simulations (default=0)
|-
| '''mode'''
| character
| insol_acc, mixed_acc,..
| indicates in which modes the tracer occurs
|-
| '''mol_weight'''
| real
| in kg/mol
| value for molar weight
|-
| '''moment'''
| integer
| 0, 3
| zeroth (number) or third (mass) moment
|-
| '''rho'''
| real
| in g/m3
| density of tracer, not needed for zeroth moment
|-
| '''sol'''
| integer
| 0 (no), 1 (yes)
| indicates whether the tracer is soluble or not
|-
| '''transport'''
| character
| stdaero, stdchem, ..., off
| choice of transport template
|-
| '''unit'''
| character
| e.g., mug kg-1, kg-1
| unit of tracer
|-
| colspan="4"| '''bold''' letters indicate which tags are always required.
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Aerosol Modes">Aerosol Modes</span>
|-
! tag
! type
! options
! description
|-
| condensation
| integer
| 0 (off), 1 (on)
| condensation of H2SO4 on this mode?
|-
| '''d_gn'''
| real
| in m
| value for the initial median diameter of the number distribution
|-
| dissfac_mean
| real
|
| dissociation factor (needed with ikoehler=1)
|-
| '''icoag'''
| integer
| 0 (off), 1 (on)
| mode involved in coagulation? If 1 for any mode, then provide coagulate.xml
|-
| ikoehler
| integer
| 0 (off), 1 (on)
| Activation via Köhler theory (warm clouds), needs dissfac_mean tag
|-
| '''kind'''
| character
| 1mom or 2mom
| 1-moment or 2-moment description of distribution
|-
| '''sigma_g'''
| real
|
| standard deviation of the distribution
|-
| shift2larger
| character
| e.g., sol_acc
| Name of larger mode to be shifted to, when diameter threshold (shift_diam) exceeded
|-
| shift2mixed
| character
| e.g., mixed_acc
| Name of mixed mode to be shifted to, when soluble mass threshold of 5% exceeded
|-
| shift_diam
| real
| in m
| diameter threshold for shift2larger
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Aerosol Emission">Aerosol Emission</span>
|-
! tag
! type
! options
! description
|-
| '''nmodes'''
| integer
|
| number of emission modes
|-
| '''d_g0_*'''
| real
|
| median diameter of number distribution of mode * (e.g., d_g0_1, d_g0_2, d_g0_3)
|-
| '''d_g3_*'''
| real
|
| median diameter of mass distribution of mode * (e.g., d_g3_1, d_g3_2, d_g3_3)
|-
| '''rho'''
| real
| in kg/m3
| particle density (same for all modes)
|-
| '''sigma_g_*'''
| real
|
| standard deviation of mode *
|-
| substance
| character
| ash, dust, na, cl, soot
| emitted substance
|-
| colspan="4"| '''bold''' letters indicate which tags are always required. <br \> routine options: volc, volc fplume, dust, biomass burn, seas smith, seas monahan, seas martensson, seas mode1, seas mode2, seas mode3
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Point Source">Point Source</span>
|-
! tag
! type
! options
! description
|-
| dg3_emiss
| real
| in m
| median diameter of aerosol mass distribution
|-
| emiss_profile
| character
|
| anti-derivative of emission profile
|-
| endTime
| character
|
| end time of emission (default=9999-12-31T00:00:00)
|-
| '''height'''
| real
| in m
| emission height
|-
| height_bot
| real
| in m
| bottom height
|-
| '''lat'''
| real
| in degree
| latitude
|-
| '''lon'''
| real
| in degree
| longitude
|-
| sigma_emiss
| real
|
| standard deviation of aerosol distribution
|-
| startTime
| character
|
| start time of emission (default=1582-10-15T00:00:00)
|-
| '''source_strength'''
| real
|
| emission source strength
|-
| '''substance'''
| character
| e.g., TRSO2
| substance nme from tracer xml
|-
| '''unit'''
| character
| e.g., kg s-1
| unit of source strength
|-
| colspan="4"| '''bold''' letters indicate which tags are always required.
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Chemistry Tracers">Chemistry Tracers</span>
|-
! tag
! type
! options
! description
|-
| '''c_solve'''
| character
| param, lt, cold, OH, linoz, simnoy, passive
| solving mechanism/strategy
|-
| emissions
|
| anthropogenic, biogenic, biomassBurning
| usage see in tracers_chemtracer_amip.xml (*)
|-
| htop_proc
| real
| in m
| top height for processes
|-
| iconv
| integer
| 0 (off), 1 (on)
| transport by convection (default=1)
|-
| initc
| character
| file
| initialize from Input file
|-
| init_mode
| integer
| 0 (off), 1 (on)
| initialize tracer
|-
| init_name
| character
|
| name of tracer in initialization file
|-
| iturb
| integer
| 0 (off), 1 (on)
| transport by turbulence (default=1)
|-
| latbc
| character
| file
| read data for LBC
|-
| lfeedback
| integer
| 0 (off), 1 (on)
| child -> parent feedback in nested simulations (default=0)
|-
| lifetime
| real
| in s
| value for lifetime
|-
| '''mol_weight'''
| real
| in kg/mol
| value for molar weight
|-
| products
| character
| name of tracer
| name of resulting tracer after depletion
|-
| tag001,...
| character
|
| name of tag to be added to tracer name
|-
| '''transport'''
| character
| stdaero, stdchem, ..., off
| choice of transport template
|-
| '''unit'''
| character
| e.g., mol mol-1
| unit of tracer
|-
| colspan="4"| '''bold''' letters indicate which tags are always required. <br\> (*) in icon-kit/externals/art/runctrl_examples/xml_ctrl/
|}

The reason for the use of those additional .xml files is that the ART variables they contain (sea salt, mineral dust etc.) might be different for every run which differs from the Icon Variables (Temperature, Pressure, etc.) which usually don’t change between runs. .xml files are readable for both humans and machines, which makes them easy to tweak and integrate. An Example for the contents of an .xml file adding ash particles can be seen below.
<syntaxhighlight lang=xml line>
<modes>
<aerosol id="asha">
<kind type="char">2mom</kind>
<d_gn type="real">1.190E-6</d_gn>
<sigma_g type="real">1.410E+0</sigma_g>
<rho type="real">2.600E+3</rho>
</aerosol>
</modes>
</syntaxhighlight>

Here is an example for a more complex .xml using [[AERODYN]]:
This is an example for the number and mass concentration of dust in a tracer .xml.
<syntaxhighlight lang=xml line>
<tracers>
<aerosol id="nmb">
<moment type="int">0</moment>
<mode type="char">insol_acc,insol_coa</mode>
<unit type="char">kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="dust">
<moment type="int">3</moment>
<mode type="char">insol_acc,insol_coa</mode>
<sol type="real">1.0</sol>
<mol_weight type="real">50.00E-3</mol_weight>
<rho type="real">2.650E3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
</tracers>
</syntaxhighlight>

Further .xml examples can be fount in <code>/your_ART_Directory/runctrl_examples/xml_ctrl. </code>

== Input Data ==

Depending on the type of simulation there might be additional input files required. These are essential files that are not contained in classical ICON initialisation data. For example for the emission of mineral dust there has to be information about the soil types supplied to the model. The additional input files should be renamed to a netcdf file and follow the naming convention shown in fig [[#input-names|1.1]]. Please note that the XXX has to be replaced by one of the indices mentionend in Table [[#tab:input-init|input-init]] and [[#tab:input-emissions|input-emissions]].

[[File:input-names.png|800px|none|alt=ICON-ART input file naming convention.|ICON-ART input file naming convention.]]

<div id="tab:input-init">

{| class="wikitable" style="text-align:left;"
|+ Additional input files for the initialisation
! Species
! Namelist switch
! Options
! XXX
|-
| Gas
| <code>iart_init_gas</code>
| 0 (cold start), 5 (from file)
| IGX
|-
| Aerosol
| <code>iart_init_aero</code>
| 0 (cold start), 5 (from file)
| IAE
|}

</div>
<div id="tab:input-emissions">

=== Obtaining Input Data ===

The 2 ways of obtaining input data are to generate it yourself or download it.

[[How to generate input data]]

[[How to download input data]]

=== Emission Data ===
In every ICON-ART, there is the possibility to add additional input data like emission data that correspond with different sources. A quick overview about them can be found below.
{| class="wikitable" style="text-align:left;"
|+ Additional input files for the emissions
! Type
! Data
! XXX
|-
| Point souces
| XML-file
| -
|-
| Sea salt
| no extra data necessary
| -
|-
| Mineral Dust
| Soil Type Data
| ART_STY
|-
| Biogenic VOCs
| Emissions or Vegetatiom
| ART_STY
|-
| Athropogenic emissions
| Emission data sets
| ART_BIO ART_ANT
|-
| Biomass burning
| Satellite data
| ART_BCF
|}

==== Chemical Tracer ====
Emission data can be obtained from several sources. The following table should give an overview about which emission data are available for a corresponding tracer. To find out when to use which emission data type we recommend respective further reading.
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Emission Data Sources">Emission Data Sources</span>
|-
! Tracer !! Emission Type !! Emission !! Resolution
|-
| C2H6 || ANT || EDGA || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || BBE || GFED || R2B04_ECHAM
|-
| || || GFED3 || R2B06_0024
|-
| || || || R2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || BIO || MEGA || R2B04_ECHAM
|-
| || || MEGAN-MACC || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0024
|-
| || || || R2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| C3H8 || ANT || EDGA || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR432-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BBE || GFED || R2B04_ECHAM
|-
| || || GFED3 || R2B06_0024
|-
| || || || R2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || BIO || MEGA || R2B04_ECHAM
|-
| || || MEGAN-MACC || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0024
|-
| || || || R2B07_nest
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| C5H8 || BBE || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BIO || MEGAN-MACC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CFCl3 || ANT || GEIA || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH3CN || BBE || GFED.1s_Akagi_daily || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || GFED.1s_Akagi_monthly || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || GFED.1s_Stockwell_daily || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || GFED.1s_Stockwell_monthly || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH3COCH3 || ANT || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || || R3B07_0026
|-
| || || POET || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GICC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GUESS-ES || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || POET || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || RETRO || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BIO || MEGAN-MACC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MEGANv2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH3I || BIO || Bell || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || Ziska || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH4 || ANT || EDGA || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR432-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || || R3B08_0049
|-
| || || EDGARv4.2 || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || || EDGARv4.3.1 || R2B04_0012
|-
| || || || R2B05_0014
|-
| CHBr3 || BIO || Liang || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || Ordonez || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || Ziska || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CO || ANT || EDGAR || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR432-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || EDGARv4.2 || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || POET || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GICC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GUESS-ES || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || RETRO || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BIO || MEGAN || R2B04_ECHAM
|-
| || || MEGAN-MACC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MEGANv2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CO2 || ANT || EDGARv4.2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0024
|-
| || || || R2B06_0031_nest
|-
| || || || R2B07_0018
|-
| || || || R2B07_nest
|-
| || || || R3B07_0022
|-
| || || || B3B07_0026
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| DMS || BBE || GFED3 || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| N2O || ANT || EDGARv4.2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| NH3 || ANT || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MACCity_agriculture || R2B06_0016
|-
| || || MACCity_agric_waste || R2B06_0016
|-
| || || MACCity_energy || R2B06_0016
|-
| || || MACCity_industrial || R2B06_0016
|-
| || || MACCity_residential || R2B06_0016
|-
| || || MACCity_transport || R2B06_0016
|-
| SF6 || ANT || EDGARLevin || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| SO2 || ANT || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_nest
|-
| || || || R3B07_0022
|-
| || || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|}

==== Remapping Emissions ====
A comprehensible manual can be found [https://gitlab.dkrz.de/art/kit-wiki/-/wikis/uploads/44c05db0bff5da516d7812292de3dff8/MECCA_Emissions01.pdf here]. The document was provided by M. Weimer (June 2019).
This document provides an overview of the workflow to be done in order to remap a set of emission data onto your own ICON grid. The raw emission data can be taken from emission inventories such as Edgar, MACCity, etc. (see above).
The desired files can be copied to an own directory where they serve as input for the remap procedure described in the manual.
Additional remarks:
* The mentioned workflow was initially designed for FH2. Should be tested on other machines as well
* The automatic addition of emission tags to the tracer.xml (add_emissions_to_tracer_xml.py) is very sensitive to tracer names and emission species. In doubt, add emission tags manually.
* In any case, double check if all emission tags have been assigned to the correct tracer

== Running a Limited Area Meteorology (LAM) Simulation ==

=== General ===

Here are some notes on setting up an ICON-ART LAM simulation. Theses settings are important if you use initial data and boundary data from different sources. It is preferable to use data from the same source to be consistent. However, in certain situations this is not possible due to limitations of the model (e.g. initialization routines).
=== Required data for LAM domain ===

Grid of LAM domain external parameters of LAM domain external parameters containing soil parameters (only necessary for dust simulations) initial data (ICON-ART or IFS)

=== Required data for LAM boundaries ===

Auxiliary grid (grid containing boundary area of the LAM domain, generated during remapping process with ICONtools) forcing data for the boundaries

===Initialization ===

There are two different possible methods to read in the dust during initialization. You can either pass a file containing meteorological variables and a second file containing dust data. The vertical levels may differ between these two files and the dust must be delivered as ART_IAE file. The corresponding namelist setting in<code> &art_nml </code> is <code> iart_init_aero=5 </code>

The other possibility is to pass all variables required for the initialization in a single file. The vertical levels must all be consistent and the corresponding namelist setting in <code> &art_nml</code> is <code> iart_init_aero=0 </code>. Furthermore you have to add file in the tracer xml file.

=== Boundary Data ===

The boundary data can only be passed to the model as one single file per time step. The vertical levels for all time steps must be the same. Otherwise an error occurs. If you use data from a different source than the one used for initialization, it is crucial to decouple the reading of the boundary data from the reading of initial data. During the start of the simulation it is possible to read the first boundary data from the initial data when using ICON-ART data. To prevent this and to read the boundary data from a separate file during initialization, set <code>init_latbc_from_fg = .FALSE. </code>in <code> &limarea_nml </code>. Additionally you have to add file in the tracer xml file.

------
[[#Aerosol Tracers]]

== Creating A Nested ICON-Grid ==
There are four steps to create a grid. The steps have to be run separately as they are dependent on each other.
=== Graph Generation ===
The first step is creating the graph. Ensure that you specify R and B for the finest nest! I.e. if you plan a global R2B6 grid with a R2B7 Nest, you have to set R=2 and B=7.
An example namelist in a runscript looks like this:
<syntaxhighlight lang=bash line>
cat > NAMELIST_GRAPH << EOF
&graph_ini
nroot = ${R}
grid_levels = ${B}
/
EOF

echo global_graph_generator null > $commandFile
job_submit -c p -p 1 -t 60 -m 32000 ${run_commmand}
</Syntaxhighlight>
You might have to increase the allocated memory (-m 64000) if the process crashes (-p 1 is maximum because of lacking parallelization)
=== Grid generation ===
The second step is the generation of all (global) grids. I.e., if you choose R=2 and B=7, you get global grids for R2B1, R2B2, R2B3, R2B4, R2B5, R2B6 and R2B7. Even if you want R2B7 to be your nest, you have to do this step down to R2B7! This means, that the values of R and B for step 1 and 2 must not differ.
As far as I could figure out, the spring dynamics optimization is the one to choose. Therefore, you should choose @itype_optimize = 4@.
An example namelist looks like this:
<syntaxhighlight lang=bash line>
cat > NAMELIST_GRID << EOF
&grid_ini
nroot = ${R}
grid_levels = ${B}
/
&grid_options
itype_optimize = 4
/
EOF

echo global_grid_generator null > $commandFile
job_submit -c p -p 1 -t 600 -m 32000 ${run_commmand}
</syntaxhighlight>
If you want to generate a finer grid (e.g. R2B10) you might have to increase the allocated memory (-m 256000).
=== Modify the filenames ===
The spring-dynamics-optimized files carry this information within their filename. In order to continue, the names have to be changed to the standard names of grids. This can be done within a script as shown in the following. (@maxlev_optim@ is a parameter, that specifies the maximum level to which optimizations are applied. This is set in the previous step within the grid_options namelist. As the default is 100, there is usually no need to change this. You just have to set the variable @maxlev_optim@ within the script for the copying):
<syntaxhighlight lang=bash line>
level=1;
while [[ $level -le $maxlev_optim ]] ; do
cp iconR${R}B0${level}-grid_spr0.90.nc iconR${R}B0${level}-grid.nc
((level=$level+1))
done
</syntaxhighlight>
=== Nested grid creation ===
As a last step, you have to specify the nests. In the following example, three nests are added to a global R2B6 grid. Therefore, @start_lev = 6@ and @n_dom = 4@. As these nests are subsequent, the @parent_id@ of each nest is the one of the domain with one rank higher. I.e. R2B6 has the ID 1, therefore the @parent_id@ of R2B7 is 1. R2B8 has the ID of the R2B7 as @parent_id@ and therefore 2. The different domains are seperated by commas in the namelist. The global domain does of course not show up (you produced the global grid files in step 2).
In this example, the further namelist variables mean the following:
@l_circ@ gives the nests a circular instead of an rectangular shape.
@l_plot@ provides output which can be used to plot the grids with GMT scripts.
@radius, center_lon, center_lat@ define the location of the nests.
With @lsep_gridref_info = .true.@ the grid information is stored within an additional grid description file. This needs then to be specified within ICON!
<syntaxhighlight lang=bash line>
cat > NAMELIST_GRIDREF << EOF
&gridref_ini
grid_root = 2
start_lev = 6
n_dom = 4
parent_id = 1,2,3
l_circ = .true.
l_plot = .true.
radius = 20.,12.,12.
center_lon = 10.,5.,5.
center_lat = 40.,47.5,47.5
bdy_indexing_depth = 14
lsep_gridref_info = .false.
/
EOF

echo global_grid_refine null > $commandFile
job_submit -c p -p 64 -t 60 -m 16000 ${run_commmand}
</syntaxhighlight>
If you want to generate a finer grid (e.g. R2B10) you might have to increase the allocated memory (-m 64000).

Tutorial Examples

2023-09-07T06:07:40Z

Luca R: moved training slides to own article

- under construction! -

In this article several examples from different topics and application areas in ICON-ART are presented, concerning the setup of a simulation as well as post processing and plotting. Just click on the highlighted text to get to the page of the respecting example. For every example a configuration is explained to try out typical cases in ICON-ART yourselves. Thereby, these pages can be viewed as a Tutorial in ICON-ART.

== Overview of Simulation topics ==

=== Atmospheric Chemistry Simulations ===
==== [[Simplified Chemistry]] ====
* here a simulation of the Hydroxylradical (OH) is performed with simplified chemistry
* this example teaches you...
** the basics of setting up an ICON runscript with ICON-ART settings
** the use of the most simple ICON-ART namelist parameter
** the implementation of the desired chemical species in a simulation by setting up a chemtracer xml-data for simplified chemistry simulations
** the implementation of emission data in a simulation
* click [[Simplified Chemistry|here]] to see more details

==== [[MECCA-based (full) Chemistry]] ====
* here a simulation with (full) MECCA-based chemistry is performed, for that a complete reaction mechanism is created and transferred to ICON-ART
* this example teaches you...
** the implementation of (full) MECCA-based chemistry in ICON-ART
** the creation of a chemical mechanism and the selection of the respective desired chemical species and their reactions in MECCA
** the implementation of not yet in MECCA implemented reactions in your mechanism
** the creation of the to your mechanism belonging Mecca-xml data to link MECCA calculations with ICON-ART
* click [[MECCA-based (full) Chemistry|here]] to see more details

==== [[Stratospheric Ozone Chemistry and Polar Stratospheric Clouds]] ====
* here a simulation with only stratospheric (simplified) ozone chemistry is performed, also polar stratospheric clouds are taken into account
* this example teaches you...
** the implementation of stratospheric specific chemistry (here on thee example of ozone)
** applying linearized ozone chemistry (LINOZ) in a simulation
** the implementation of polar stratospheric clouds (PSCs)
* click [[Stratospheric Ozone Chemistry and Polar Stratospheric Clouds|here]] to see more details

==== [[Lifetime Tracer Simulation]] ====
* here it is shown how to simulate a lifetime driven tracer with simplified chemistry in ICON-ART. Here the very short-lived substances Bromoform (<chem>CHBr3</chem>) and Dibromomethane (<chem>CH2Br2</chem>) are simulated.
* this example teaches you...
** the correct setup of the xml-file for such tracers.
** the structure and correct usage of namelist parameter settings of lifetime tracers in xml-files
** editing output variables
* click [[Lifetime Tracer Simulation|here]] to see more details

---------

=== Aerosol Simulations ===
==== [[Simulating a Point Source]] ====
* In this example a Point Source is installed at the location of the volcano Raikoke and a constant emission of SO2 and Ash is emitted into the atmosphere.
* this example teaches you...
** How to setup and run your first Simulation
** The setup of the directory structure when doing an ICON-ART Simulation
** How to modify a runscript to suit your needs
** How to modify .xml data to set up a Simulation
** A possible approach to visualise the output data

== Overview of Post processing examples ==

=== Chemical tracer distributions ===

=== Aerosol Point Sources ===

Training Course

2023-09-07T06:07:14Z

Luca R: created Page

== Slides ICON-Training ==

[[File:ICON-ART-TC2023.pdf|left]]

TODO's

2023-09-06T20:34:06Z

Luca R: Created page with "== Currently Being Worked on == * Unify Formatting of all Articles == To Do's == * Insert references using Cite in the current articles * Install Math extension * Write Po..."

== Currently Being Worked on ==

* Unify Formatting of all Articles

== To Do's ==

* Insert references using Cite in the current articles
* Install Math extension
* Write Post-Processing article
* Write Code, Folder and Data structure article
* Programming Article
* Input Data Generation Article
* ISORROPIA article
* Article Explaining Tet suites
* ICON ART Simulation Articles
** Dust
** Pollen
** Seasalt
** ...

Input

2023-08-31T08:46:45Z

Luca R: made emission data sources Table collapsible

== Requirements for a Simulation ==

To run a Simulation with ICON-ART there are three main points to consider: Setting the namelist parameters, preparing the xmls, and optionally prepare additional Input data.

== Namelist Inputs ==

To run a simulation with ICON-ART installed the first thing to prepare is the runfile. It is usually best to start with an existing and working runfile and adapt parts as needed. The runfile contains the majority of parameters to run a simulation, like the length, timesteps and grids used, as well as a variety of options for methods and parametrisations used in the model. The runfile contents are then split up in several namelists that the model reads at the beginning of a simulation. An overview of the Namelist Parameters can be found in [[Namelist]].

To enable ART in an ICON simulation, the switch <code>lart = .TRUE.</code> has to be set in the section <code>&run_nml</code>. This is the global on/off switch for ICON-ART. This is how this could look like in the context of a runfile:

<pre>! run_nml: general switches ----------
&run_nml
ltestcase = .FALSE.
num_lev = 50
ltransport = .TRUE.
.............

\textcolor{red}{lart = .TRUE.}</pre>
Generally the first letters of a namelist switch refer to its type, the "l" in "lart" stands for logical, meaning it has to be either True or False. Here are some examples.

<div id="tab:vartypes">

{| class="wikitable" style="text-align:left;"
|+ Some namelist switches and their data types.
! namelist switch
! type
|-
| <code>lart</code>
| logical
|-
| <code>cart_aerosol_xml</code>
| character
|-
| <code>iart_init_aero </code>
| integer
|}

</div>
The namelist <code>&art_nml</code> is used for general options of the ART simulation. To run a certain kind of simulation the according switch has to be set to <code>.TRUE.</code>. For Example to include a point source the switch <code>cart_aerosol_xml</code> has to be set to <code>.TRUE.</code>.

== XML Inputs ==

This enables the inclusion of a .xml file containing additional information like location and strength of the point source. The table [[#tab:art_nml-params|below]] contains the most important <code>&art_nml</code> namelist parameters and additional namelist parameters required if they are set to <code>.TRUE.</code>.

<div id="tab:art_nml-params">

{| class="wikitable" style="text-align:left;"
|+ XML files and their namelist dependencies
! XML File
! Description
! Namelist parameter dependency
! Default
! Further Information
|---
| <code>cart_chemtracer_xml</code>
| Switch for simple OH chemistry
| <code>lart_chemtracer</code>
| .FALSE.
|[[#Chemistry Tracers|Chemistry Tracers]]
|-
| <code>cart_mecca_xml</code>
| Switch for kpp chemistry
| <code>lart_mecca</code>
| .FALSE.
|[[Atmospheric Chemistry]]
|-
| <code>cart_pntSrc_xml</code>
| Enables creation of point sources emitting given Aerosols at a given rate
| <code>lart_pntSrc</code>
| .FALSE.
|[[#Point Source|Point Source]]
|-
| <code>cart_aerosol_xml</code>
| Main switch for the treatment of atmospheric aerosols
| <code>lart_aerosol</code>
| .FALSE.
|[[#Aerosol Tracers|Aerosol Tracers]]
|-
| <code>cart_modes_xml</code>
| Main switch for the treatment of atmospheric aerosols
| <code>lart_aerosol</code>
| .FALSE.
|[[#Aerosol Modes|Aerosol Modes]]
|-
| <code>cart_diagnostics_xml</code>
| Enables diagnostic output fields
| <code>lart_diag_out</code>
| .FALSE.
| -
|-
| <code>cart_emiss_xml_file</code>
| XML File for emission metadata
| -
| -
|[[#Aerosol Emission|Aerosol Emission]]
|-
| <code>cart_ext_data_xml</code>
| XML File for metadata of datasets prescribing tracers
| -
| -
| -
|-
| <code> cart_coag_xml </code>
| XML File containing additional information about coagulation
| -
| -
| -
|}

</div>

{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Aerosol Tracers">Aerosol Tracers</span>
|-
! tag
! type
! options
! description
|-
| htop_proc
| real
| in m
| top height for processes
|-
| initc
| character
| file
| initialize from Input file
|-
| inucl
| integer
| 0 (off), 1 (on)
| H2SO4 nucleation for so4 tracer (default=1); 1 for so4_sol_ait, 0 for other so4 tracer)
|-
| label
| character
| e.g., dusta
| allows to name tracers individually
|-
| latbc
| character
| file
| read data for LBC
|-
| lfeedback
| integer
| 0 (off), 1 (on)
| child -> parent feedback in nested simulations (default=0)
|-
| '''mode'''
| character
| insol_acc, mixed_acc,..
| indicates in which modes the tracer occurs
|-
| '''mol_weight'''
| real
| in kg/mol
| value for molar weight
|-
| '''moment'''
| integer
| 0, 3
| zeroth (number) or third (mass) moment
|-
| '''rho'''
| real
| in g/m3
| density of tracer, not needed for zeroth moment
|-
| '''sol'''
| integer
| 0 (no), 1 (yes)
| indicates whether the tracer is soluble or not
|-
| '''transport'''
| character
| stdaero, stdchem, ..., off
| choice of transport template
|-
| '''unit'''
| character
| e.g., mug kg-1, kg-1
| unit of tracer
|-
| colspan="4"| '''bold''' letters indicate which tags are always required.
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Aerosol Modes">Aerosol Modes</span>
|-
! tag
! type
! options
! description
|-
| condensation
| integer
| 0 (off), 1 (on)
| condensation of H2SO4 on this mode?
|-
| '''d_gn'''
| real
| in m
| value for the initial median diameter of the number distribution
|-
| dissfac_mean
| real
|
| dissociation factor (needed with ikoehler=1)
|-
| '''icoag'''
| integer
| 0 (off), 1 (on)
| mode involved in coagulation? If 1 for any mode, then provide coagulate.xml
|-
| ikoehler
| integer
| 0 (off), 1 (on)
| Activation via Köhler theory (warm clouds), needs dissfac_mean tag
|-
| '''kind'''
| character
| 1mom or 2mom
| 1-moment or 2-moment description of distribution
|-
| '''sigma_g'''
| real
|
| standard deviation of the distribution
|-
| shift2larger
| character
| e.g., sol_acc
| Name of larger mode to be shifted to, when diameter threshold (shift_diam) exceeded
|-
| shift2mixed
| character
| e.g., mixed_acc
| Name of mixed mode to be shifted to, when soluble mass threshold of 5% exceeded
|-
| shift_diam
| real
| in m
| diameter threshold for shift2larger
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Aerosol Emission">Aerosol Emission</span>
|-
! tag
! type
! options
! description
|-
| '''nmodes'''
| integer
|
| number of emission modes
|-
| '''d_g0_*'''
| real
|
| median diameter of number distribution of mode * (e.g., d_g0_1, d_g0_2, d_g0_3)
|-
| '''d_g3_*'''
| real
|
| median diameter of mass distribution of mode * (e.g., d_g3_1, d_g3_2, d_g3_3)
|-
| '''rho'''
| real
| in kg/m3
| particle density (same for all modes)
|-
| '''sigma_g_*'''
| real
|
| standard deviation of mode *
|-
| substance
| character
| ash, dust, na, cl, soot
| emitted substance
|-
| colspan="4"| '''bold''' letters indicate which tags are always required. <br \> routine options: volc, volc fplume, dust, biomass burn, seas smith, seas monahan, seas martensson, seas mode1, seas mode2, seas mode3
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Point Source">Point Source</span>
|-
! tag
! type
! options
! description
|-
| dg3_emiss
| real
| in m
| median diameter of aerosol mass distribution
|-
| emiss_profile
| character
|
| anti-derivative of emission profile
|-
| endTime
| character
|
| end time of emission (default=9999-12-31T00:00:00)
|-
| '''height'''
| real
| in m
| emission height
|-
| height_bot
| real
| in m
| bottom height
|-
| '''lat'''
| real
| in degree
| latitude
|-
| '''lon'''
| real
| in degree
| longitude
|-
| sigma_emiss
| real
|
| standard deviation of aerosol distribution
|-
| startTime
| character
|
| start time of emission (default=1582-10-15T00:00:00)
|-
| '''source_strength'''
| real
|
| emission source strength
|-
| '''substance'''
| character
| e.g., TRSO2
| substance nme from tracer xml
|-
| '''unit'''
| character
| e.g., kg s-1
| unit of source strength
|-
| colspan="4"| '''bold''' letters indicate which tags are always required.
|}
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Chemistry Tracers">Chemistry Tracers</span>
|-
! tag
! type
! options
! description
|-
| '''c_solve'''
| character
| param, lt, cold, OH, linoz, simnoy, passive
| solving mechanism/strategy
|-
| emissions
|
| anthropogenic, biogenic, biomassBurning
| usage see in tracers_chemtracer_amip.xml (*)
|-
| htop_proc
| real
| in m
| top height for processes
|-
| iconv
| integer
| 0 (off), 1 (on)
| transport by convection (default=1)
|-
| initc
| character
| file
| initialize from Input file
|-
| init_mode
| integer
| 0 (off), 1 (on)
| initialize tracer
|-
| init_name
| character
|
| name of tracer in initialization file
|-
| iturb
| integer
| 0 (off), 1 (on)
| transport by turbulence (default=1)
|-
| latbc
| character
| file
| read data for LBC
|-
| lfeedback
| integer
| 0 (off), 1 (on)
| child -> parent feedback in nested simulations (default=0)
|-
| lifetime
| real
| in s
| value for lifetime
|-
| '''mol_weight'''
| real
| in kg/mol
| value for molar weight
|-
| products
| character
| name of tracer
| name of resulting tracer after depletion
|-
| tag001,...
| character
|
| name of tag to be added to tracer name
|-
| '''transport'''
| character
| stdaero, stdchem, ..., off
| choice of transport template
|-
| '''unit'''
| character
| e.g., mol mol-1
| unit of tracer
|-
| colspan="4"| '''bold''' letters indicate which tags are always required. <br\> (*) in icon-kit/externals/art/runctrl_examples/xml_ctrl/
|}

The reason for the use of those additional .xml files is that the ART variables they contain (sea salt, mineral dust etc.) might be different for every run which differs from the Icon Variables (Temperature, Pressure, etc.) which usually don’t change between runs. .xml files are readable for both humans and machines, which makes them easy to tweak and integrate. An Example for the contents of an .xml file adding ash particles can be seen below.
<syntaxhighlight lang=xml line>
<modes>
<aerosol id="asha">
<kind type="char">2mom</kind>
<d_gn type="real">1.190E-6</d_gn>
<sigma_g type="real">1.410E+0</sigma_g>
<rho type="real">2.600E+3</rho>
</aerosol>
</modes>
</syntaxhighlight>

Here is an example for a more complex .xml using [[AERODYN]]:
This is an example for the number and mass concentration of dust in a tracer .xml.
<syntaxhighlight lang=xml line>
<tracers>
<aerosol id="nmb">
<moment type="int">0</moment>
<mode type="char">insol_acc,insol_coa</mode>
<unit type="char">kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="dust">
<moment type="int">3</moment>
<mode type="char">insol_acc,insol_coa</mode>
<sol type="real">1.0</sol>
<mol_weight type="real">50.00E-3</mol_weight>
<rho type="real">2.650E3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
</tracers>
</syntaxhighlight>

Further .xml examples can be fount in <code>/your_ART_Directory/runctrl_examples/xml_ctrl. </code>

== Input Data ==

Depending on the type of simulation there might be additional input files required. These are essential files that are not contained in classical ICON initialisation data. For example for the emission of mineral dust there has to be information about the soil types supplied to the model. The additional input files should be renamed to a netcdf file and follow the naming convention shown in fig [[#input-names|1.1]]. Please note that the XXX has to be replaced by one of the indices mentionend in Table [[#tab:input-init|input-init]] and [[#tab:input-emissions|input-emissions]].

[[File:input-names.png|800px|none|alt=ICON-ART input file naming convention.|ICON-ART input file naming convention.]]

<div id="tab:input-init">

{| class="wikitable" style="text-align:left;"
|+ Additional input files for the initialisation
! Species
! Namelist switch
! Options
! XXX
|-
| Gas
| <code>iart_init_gas</code>
| 0 (cold start), 5 (from file)
| IGX
|-
| Aerosol
| <code>iart_init_aero</code>
| 0 (cold start), 5 (from file)
| IAE
|}

</div>
<div id="tab:input-emissions">

=== Obtaining Input Data ===

The 2 ways of obtaining input data are to generate it yourself or download it.

[[How to generate input data]]

[[How to download input data]]

=== Emission Data ===
In every ICON-ART, there is the possibility to add additional input data like emission data that correspond with different sources. A quick overview about them can be found below.
{| class="wikitable" style="text-align:left;"
|+ Additional input files for the emissions
! Type
! Data
! XXX
|-
| Point souces
| XML-file
| -
|-
| Sea salt
| no extra data necessary
| -
|-
| Mineral Dust
| Soil Type Data
| ART_STY
|-
| Biogenic VOCs
| Emissions or Vegetatiom
| ART_STY
|-
| Athropogenic emissions
| Emission data sets
| ART_BIO ART_ANT
|-
| Biomass burning
| Satellite data
| ART_BCF
|}

==== Chemical Tracer ====
Emission data can be obtained from several sources. The following table should give an overview about which emission data are available for a corresponding tracer. To find out when to use which emission data type we recommend respective further reading.
{| class="mw-collapsible mw-collapsed wikitable" style="text-align:left;"
|+ style=white-space:nowrap | <span id="Emission Data Sources">Emission Data Sources</span>
|-
! Tracer !! Emission Type !! Emission !! Resolution
|-
| C2H6 || ANT || EDGA || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || BBE || GFED || R2B04_ECHAM
|-
| || || GFED3 || R2B06_0024
|-
| || || || R2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || BIO || MEGA || R2B04_ECHAM
|-
| || || MEGAN-MACC || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0024
|-
| || || || R2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| C3H8 || ANT || EDGA || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR432-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BBE || GFED || R2B04_ECHAM
|-
| || || GFED3 || R2B06_0024
|-
| || || || R2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || BIO || MEGA || R2B04_ECHAM
|-
| || || MEGAN-MACC || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0024
|-
| || || || R2B07_nest
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| C5H8 || BBE || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BIO || MEGAN-MACC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CFCl3 || ANT || GEIA || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH3CN || BBE || GFED.1s_Akagi_daily || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || GFED.1s_Akagi_monthly || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || GFED.1s_Stockwell_daily || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || GFED.1s_Stockwell_monthly || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH3COCH3 || ANT || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || || R3B07_0026
|-
| || || POET || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GICC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GUESS-ES || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || POET || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || RETRO || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BIO || MEGAN-MACC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MEGANv2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH3I || BIO || Bell || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || Ziska || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CH4 || ANT || EDGA || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR432-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || || R3B08_0049
|-
| || || EDGARv4.2 || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || || EDGARv4.3.1 || R2B04_0012
|-
| || || || R2B05_0014
|-
| CHBr3 || BIO || Liang || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || Ordonez || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || Ziska || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CO || ANT || EDGAR || R2B04_ECHAM
|-
| || || EDGAR+ || R2B05_0014
|-
| || || EDGAR432-monthly || R2B05_0014
|-
| || || || R2B06_0024
|-
| || || || B2B06_EU_nest
|-
| || || || R2B07_nest
|-
| || || || R2B07_NSR_D_nest
|-
| || || || R3B07_0026
|-
| || || EDGARv4.2 || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || POET || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GICC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || GUESS-ES || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || RETRO || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || BIO || MEGAN || R2B04_ECHAM
|-
| || || MEGAN-MACC || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MEGANv2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| CO2 || ANT || EDGARv4.2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0024
|-
| || || || R2B06_0031_nest
|-
| || || || R2B07_0018
|-
| || || || R2B07_nest
|-
| || || || R3B07_0022
|-
| || || || B3B07_0026
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| DMS || BBE || GFED3 || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MACCity || R2B04_0012
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| N2O || ANT || EDGARv4.2 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| || BBE || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_polstrNest
|-
| || || || R3B07_0022
|-
| NH3 || ANT || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| || || MACCity_agriculture || R2B06_0016
|-
| || || MACCity_agric_waste || R2B06_0016
|-
| || || MACCity_energy || R2B06_0016
|-
| || || MACCity_industrial || R2B06_0016
|-
| || || MACCity_residential || R2B06_0016
|-
| || || MACCity_transport || R2B06_0016
|-
| SF6 || ANT || EDGARLevin || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|-
| SO2 || ANT || MACCity || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R2B07_nest
|-
| || || || R3B07_0022
|-
| || || GFED3 || R2B04_0012
|-
| || || || R2B04_ECHAM
|-
| || || || R2B05_0014
|-
| || || || R2B06_0016
|-
| || || || R2B07_0018
|-
| || || || R3B07_0022
|}

==== Remapping Emissions ====
A comprehensible manual can be found [https://gitlab.dkrz.de/art/kit-wiki/-/wikis/uploads/44c05db0bff5da516d7812292de3dff8/MECCA_Emissions01.pdf here]. The document was provided by M. Weimer (June 2019).
This document provides an overview of the workflow to be done in order to remap a set of emission data onto your own ICON grid. The raw emission data can be taken from emission inventories such as Edgar, MACCity, etc. (see above).
The desired files can be copied to an own directory where they serve as input for the remap procedure described in the manual.
Additional remarks:
* The mentioned workflow was initially designed for FH2. Should be tested on other machines as well
* The automatic addition of emission tags to the tracer.xml (add_emissions_to_tracer_xml.py) is very sensitive to tracer names and emission species. In doubt, add emission tags manually.
* In any case, double check if all emission tags have been assigned to the correct tracer

== Running a Limited Area Meteorology (LAM) Simulation ==

=== General ===

Here are some notes on setting up an ICON-ART LAM simulation. Theses settings are important if you use initial data and boundary data from different sources. It is preferable to use data from the same source to be consistent. However, in certain situations this is not possible due to limitations of the model (e.g. initialization routines).
=== Required data for LAM domain ===

Grid of LAM domain external parameters of LAM domain external parameters containing soil parameters (only necessary for dust simulations) initial data (ICON-ART or IFS)

=== Required data for LAM boundaries ===

Auxiliary grid (grid containing boundary area of the LAM domain, generated during remapping process with ICONtools) forcing data for the boundaries

===Initialization ===

There are two different possible methods to read in the dust during initialization. You can either pass a file containing meteorological variables and a second file containing dust data. The vertical levels may differ between these two files and the dust must be delivered as ART_IAE file. The corresponding namelist setting in<code> &art_nml </code> is <code> iart_init_aero=5 </code>

The other possibility is to pass all variables required for the initialization in a single file. The vertical levels must all be consistent and the corresponding namelist setting in <code> &art_nml</code> is <code> iart_init_aero=0 </code>. Furthermore you have to add file in the tracer xml file.

=== Boundary Data ===

The boundary data can only be passed to the model as one single file per time step. The vertical levels for all time steps must be the same. Otherwise an error occurs. If you use data from a different source than the one used for initialization, it is crucial to decouple the reading of the boundary data from the reading of initial data. During the start of the simulation it is possible to read the first boundary data from the initial data when using ICON-ART data. To prevent this and to read the boundary data from a separate file during initialization, set <code>init_latbc_from_fg = .FALSE. </code>in <code> &limarea_nml </code>. Additionally you have to add file in the tracer xml file.

------
[[#Aerosol Tracers]]

== Creating A Nested ICON-Grid ==
There are four steps to create a grid. The steps have to be run separately as they are dependent on each other.
=== Graph Generation ===
The first step is creating the graph. Ensure that you specify R and B for the finest nest! I.e. if you plan a global R2B6 grid with a R2B7 Nest, you have to set R=2 and B=7.
An example namelist in a runscript looks like this:
<pre>cat > NAMELIST_GRAPH << EOF
&graph_ini
nroot = ${R}
grid_levels = ${B}
/
EOF

echo global_graph_generator null > $commandFile
job_submit -c p -p 1 -t 60 -m 32000 ${run_commmand}</pre>
You might have to increase the allocated memory (-m 64000) if the process crashes (-p 1 is maximum because of lacking parallelization)
=== Grid generation ===
The second step is the generation of all (global) grids. I.e., if you choose R=2 and B=7, you get global grids for R2B1, R2B2, R2B3, R2B4, R2B5, R2B6 and R2B7. Even if you want R2B7 to be your nest, you have to do this step down to R2B7! This means, that the values of R and B for step 1 and 2 must not differ.
As far as I could figure out, the spring dynamics optimization is the one to choose. Therefore, you should choose @itype_optimize = 4@.
An example namelist looks like this:
<pre>cat > NAMELIST_GRID << EOF
&grid_ini
nroot = ${R}
grid_levels = ${B}
/
&grid_options
itype_optimize = 4
/
EOF

echo global_grid_generator null > $commandFile
job_submit -c p -p 1 -t 600 -m 32000 ${run_commmand}</pre>
If you want to generate a finer grid (e.g. R2B10) you might have to increase the allocated memory (-m 256000).
=== Modify the filenames ===
The spring-dynamics-optimized files carry this information within their filename. In order to continue, the names have to be changed to the standard names of grids. This can be done within a script as shown in the following. (@maxlev_optim@ is a parameter, that specifies the maximum level to which optimizations are applied. This is set in the previous step within the grid_options namelist. As the default is 100, there is usually no need to change this. You just have to set the variable @maxlev_optim@ within the script for the copying):
<pre>level=1;
while [[ $level -le $maxlev_optim ]] ; do
cp iconR${R}B0${level}-grid_spr0.90.nc iconR${R}B0${level}-grid.nc
((level=$level+1))
done</pre>
=== Nested grid creation ===
As a last step, you have to specify the nests. In the following example, three nests are added to a global R2B6 grid. Therefore, @start_lev = 6@ and @n_dom = 4@. As these nests are subsequent, the @parent_id@ of each nest is the one of the domain with one rank higher. I.e. R2B6 has the ID 1, therefore the @parent_id@ of R2B7 is 1. R2B8 has the ID of the R2B7 as @parent_id@ and therefore 2. The different domains are seperated by commas in the namelist. The global domain does of course not show up (you produced the global grid files in step 2).
In this example, the further namelist variables mean the following:
@l_circ@ gives the nests a circular instead of an rectangular shape.
@l_plot@ provides output which can be used to plot the grids with GMT scripts.
@radius, center_lon, center_lat@ define the location of the nests.
With @lsep_gridref_info = .true.@ the grid information is stored within an additional grid description file. This needs then to be specified within ICON!
<pre>cat > NAMELIST_GRIDREF << EOF
&gridref_ini
grid_root = 2
start_lev = 6
n_dom = 4
parent_id = 1,2,3
l_circ = .true.
l_plot = .true.
radius = 20.,12.,12.
center_lon = 10.,5.,5.
center_lat = 40.,47.5,47.5
bdy_indexing_depth = 14
lsep_gridref_info = .false.
/
EOF

echo global_grid_refine null > $commandFile
job_submit -c p -p 64 -t 60 -m 16000 ${run_commmand}</pre>
If you want to generate a finer grid (e.g. R2B10) you might have to increase the allocated memory (-m 64000).

Tutorial Examples

2023-08-31T08:35:32Z

Luca R: removed old Training slides pngs,uploaded new Training PDF

File:ICON-ART-TC2023.pdf

2023-08-31T08:33:24Z

Luca R: Summary Slides of The ICON-ART Training course 2023

== Summary ==
Summary Slides of The ICON-ART Training course 2023

AERODYN

2023-08-09T08:41:54Z

Luca R: syntaxhighlight for xmls

AERODYN (AEROsol DYNamics) is a new aerosol dynamics module used in ICON-ART. With AERODYN, secondary aerosol formations (e.g. the accumulation of SO2 on a volcanic ash particle) can be accounted for.
Aerosol dynamic processes are the processes which change an aerosols properties after it has been released into the atmosphere. This change of properties (omitting chemical reactions) occurs either via condensation or coagulation.
Since it is not feasible to account for every individual aerosol to calculate their movement, coagulation and condensation, AERODYN uses a total of 12 categories for aerosol modes.

== Aerosol Modes in AERODYN ==

Traditionally, There are three categories to account for the aerosol sizes [https://www.sciencedirect.com/science/article/abs/pii/0004698178901968 Whitby et. Al. (1978)] : Aitken mode, accumulation mode and coarse mode. In AERODYN there is also a fourth mode, called giant mode, which accounts for especially large Aerosols, e.g. some species of pollen.

Since condensation can occur on a solid aerosol particle, as well as for a number of other reasons, it is possible to obtain aerosol modes which are a mixture of solid and liquid modes. To account for this, AERODYN combines 4 modes of aerosol size with the 3 modes soluble, insoluble and mixed.
A Visual Representation is given here:

[[File:AERODYNModes.png]]

when also adding the giant mode, the following configurations are possible:

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ AERODYN mode configurations
! || aitken || accumulation || coarse || giant
|-
! soluble
| sol_aitken || sol_acc || sol_coarse || sol_giant
|-
! insoluble
|insol_aitken || insol_acc || insol_coarse || insol_giant
|-
! mixed
| mixed_aitken ||mixed_acc || mixed_coarse || mixed_giant
|-
|}

== Condensation, Nucleation and Accumulation Processes ==
Here is a quick reminder of the different aerosol processes:

* Condensation is the process of a gas turning into a liquid. In practice, some sort of surface, or condensation nucleus is needed to enable the phase change.
* Nucleation is the process of an aerosol particle forming from a gaseous phase by random collisions of the gas molecules. The produced aerosol sizes are in the Aitken mode.
* Accumulation in the Process of the initial nucleus growing in size through collision and deposition Processes.

A schematic overview of the different microphysical aerosol processes can be found in the following Figure from [https://doi.org/10.1016/S1352-2310(00)00239-9 Raes et. Al(2000)]:

[[File:Raes2000.jpg| Overview of Aerosol Processes.]]

In AERODYN, some modes are treated differently with respect to the dynamic processes: Nucleation is only considered for the Aitken mode, and the Giant mode is not affected by condensation. This is done to save computing power while still retaining realism, and since Nucleation of larger Particles without first having Aitken mode Aerosols is highly unlikely, this is a good approximation. Also, as condensation highly prefers smaller modes to the giant mode, this also is a reasonable way to save computing power. For details how the processes are calculated, see [https://doi.org/10.1029/2018RG000615 Riemer et Al. (2003)] and [https://doi.org/10.5194/acp-9-8661-2009 Vogel et. Al. (2009)].

== Defining Aerosols with AERODYN ==

AERODYN is meant to be as flexible as possible, allowing the user to quickly define their own aerosols and then use ICON-ART to simulate their dynamics.
Aerosols are defined in .xml files, which can then be included in the run by specifying their path in the runfile (see [[Input]] for more details)

There are several Types of .xml files for different applications. The most important/most used .xml files are compiled in the following paragraphs.

=== Defining Tracers ===

Below is the tracers.xml, where the aerosols as well as their properties are defined. The properties to be defined are listed in the following table.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ tracers.xml Parameters
! Parameter || type || description || example
|-
! aerosol id
| char ||specifies name of the aerosol || "my_first_aerosol"
|-
! moment
| int || this is the moment of the log-normal distribution (0 for number concentration and 3 for mass mixing ratio)|| 3
|-
! mode
| char || mode name from AERODYN mode configurations (multiple modes can be given, seperated by a comma) || insol_giant,sol_aitken
|-

! mol_weight
| real || specify the molar mass || 30E-3
|-
! rho
| real || specify the density in terms of <code> unit </code> || 2e3

|-
! unit
| char || specify what unit the aerosol will have || kg-1
|-
! transport
| char || specify transport scheme to be used || hadv52aero
|-
|}

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of tracers.xml: ====
<syntaxhighlight class="mw-collapsible-content" lang=xml line>
<tracers>
<aerosol id="nmb">
<moment type="int">0</moment>
<mode type="char">insol_acc,insol_coa,sol_acc,sol_coa,mixed_acc,mixed_coa</mode>
<unit type="char">kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="dust">
<moment type="int">3</moment>
<mode type="char">insol_acc,insol_coa,mixed_acc,mixed_coa</mode>
<sol type="real">0.0</sol>
<mol_weight type="real">50.00E-3</mol_weight>
<rho type="real">2.650E3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="na">
<moment type="int">3</moment>
<mode type="char">sol_acc,sol_coa,mixed_acc,mixed_coa</mode>
<sol type="real">1.</sol>
<mol_weight type="real">22.9898E-3</mol_weight>
<rho type="real">2.2E+3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="cl">
<moment type="int">3</moment>
<mode type="char">sol_acc,sol_coa,mixed_acc,mixed_coa</mode>
<sol type="real">1.</sol>
<mol_weight type="real">35.453E-3</mol_weight>
<rho type="real">2.2E+3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="h2o">
<moment type="int">3</moment>
<mode type="char">sol_acc,sol_coa</mode>
<sol type="real">1.</sol>
<mol_weight type="real">18.01528E-3</mol_weight>
<rho type="real">1.E3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">off</transport>
</aerosol>
</tracers>
</syntaxhighlight>
</div>

=== Defining Modes ===

In the modes.xml the modes which are going to be used in the simulation can be defined. The following .xml contains the same information as Figure 3.1, but in a way the computer can interpret it.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ modes.xml Parameters
! Parameter || type || description || example
|-
! aerosol id
| char ||specifies name of the aerosol || "my_first_aerosol"
|-
! kind
| char || specifies which properties are to be calculated in the model (eg 2mom -> mass mixing ratio and number concentration) || 2mom
|-
! d_gn
| real || median diameter of the lognormal aerosol distribution || 2.0E-6
|-
! sigma_g
| real || standard deviation of the lognormal aerosol distribution || 2.2
|-
! condensation
| int || switch to control if condensation can occur || 0
|-
! icoag
| int || switch to control if coagulation can occur || 0
|-
! shift2larger
| char || prevents mode from getting to large, then shifts size to given mode (mostly used for Aitken mode) || insol_acc
|-
! shift_diam
| real || defines the diameter where the aerosol is shifted to the given mode || 0.08E-6
|-
! shift2mixed
| char || shift to given mode when 5% of mass is liquid || mixed_coa
|-
|}

Sometimes it is necessary for the type of an aerosol to change its given mode if its diameter gets too large or a large amount of condensation occurs. The mass can then be transferred to defined modes, as seen in the table above.

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of modes.xml: ====
<syntaxhighlight class="mw-collapsible-content" lang=xml line>
<modes>
<aerosol id="sol_acc">
<kind type="char">2mom</kind>
<d_gn type="real">0.2E-6</d_gn>
<sigma_g type="real">2.0</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="insol_acc">
<kind type="char">2mom</kind>
<d_gn type="real">6.445E-7</d_gn>
<sigma_g type="real">1.700E+0</sigma_g>
<condensation type="int">0</condensation>
<shift2mixed type="char">mixed_acc</shift2mixed>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="mixed_acc">
<kind type="char">2mom</kind>
<d_gn type="real">6.445E-7</d_gn>
<sigma_g type="real">1.700E+0</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="sol_coa">
<kind type="char">2mom</kind>
<d_gn type="real">2.0E-6</d_gn>
<sigma_g type="real">2.2</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="insol_coa">
<kind type="char">2mom</kind>
<d_gn type="real">2.0E-6</d_gn>
<sigma_g type="real">2.2</sigma_g>
<condensation type="int">0</condensation>
<shift2mixed type="char">mixed_coa</shift2mixed>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="mixed_coa">
<kind type="char">2mom</kind>
<d_gn type="real">2.0E-6</d_gn>
<sigma_g type="real">2.2</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
</modes>
</syntaxhighlight>
</div>

=== Defining Coagulation ===

Here the behaviour of collision and subsequent merging (coagulation) between aerosols can be defined. If for example a solid (insol) aerosol merges with a liquid (sol) aerosol, the result should logically be a mixed aerosol, and so on. This can further be defined for all sizes of aerosols as seen in the example .xml for accumulation and coarse modes.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ coagulate.xml Parameters
! Parameter || type || description || example
|-
!smallmode id
| char || The aerosol_mode which is gonig to coagulate || insol_acc
|-
!nmodes
| int || Defines the number of modes for which the coagulation preocesses are defined || 5
|-
!bigmode id
| char || Inside angle brackets : coagulation partner of the given aerosol mode. Outside angle brackets: the new mode of the cagulated aerosols. || <bigmode id="sol_acc" type="char">mixed_acc</bigmode>
|-
|}

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of coag.xml: ====
<syntaxhighlight class="mw-collapsible-content" lang=xml line>
<coagulate>
<smallmode id="sol_acc">
<nmodes type="int">6</nmodes>
<bigmode id="sol_acc" type="char">sol_acc</bigmode>
<bigmode id="insol_acc" type="char">mixed_acc</bigmode>
<bigmode id="mixed_acc" type="char">mixed_acc</bigmode>
<bigmode id="sol_coa" type="char">sol_coa</bigmode>
<bigmode id="insol_coa" type="char">mixed_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="insol_acc">
<nmodes type="int">5</nmodes>
<bigmode id="insol_acc" type="char">insol_acc</bigmode>
<bigmode id="mixed_acc" type="char">mixed_acc</bigmode>
<bigmode id="sol_coa" type="char">mixed_coa</bigmode>
<bigmode id="insol_coa" type="char">insol_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="mixed_acc">
<nmodes type="int">4</nmodes>
<bigmode id="mixed_acc" type="char">mixed_acc</bigmode>
<bigmode id="sol_coa" type="char">mixed_coa</bigmode>
<bigmode id="insol_coa" type="char">mixed_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="sol_coa">
<nmodes type="int">3</nmodes>
<bigmode id="sol_coa" type="char">sol_coa</bigmode>
<bigmode id="insol_coa" type="char">mixed_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="insol_coa">
<nmodes type="int">2</nmodes>
<bigmode id="insol_coa" type="char">insol_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="mixed_coa">
<nmodes type="int">1</nmodes>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
</coagulate>
</syntaxhighlight>
</div>
=== Defining the Model Output ===
In the diagnostics.xml the diagnostic variables that can be written out are defined. Below are some examples, however this is a heavily shortened version of the diagnostics.xml used in ICON-ART.

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of diagnostics.xml: ====
<syntaxhighlight class="mw-collapsible-content" lang=xml line>
<diagnostics>
<aerosol id="diam_insol_acc">
<productDefinitionTemplate type="int">40</productDefinitionTemplate>
<discipline type="int">0</discipline>
<parameterCategory type="int">254</parameterCategory>
<parameterNumber type="int">201</parameterNumber>
<constituentType type="int">62001</constituentType>
</aerosol>
<aerosol id="acc_drydepo">
<tracername_list type="char">dust_insol_acc,dust_insol_coa,na_sol_acc,na_sol_coa,cl_sol_acc,cl_sol_coa,nmb_sol_acc,nmb_sol_coa,nmb_insol_acc,nmb_insol_coa</tracername_list>
<parameterNumber_list type="char">6,6,6,6,6,6,192,192,192,192</parameterNumber_list>
<productDefinitionTemplate_list type="char">67,67,67,67,67,67,67,67,67,67</productDefinitionTemplate_list>
<typeOfStatisticalProcessing type="int">1</typeOfStatisticalProcessing>
<bitsPerValue type="int">32</bitsPerValue>
<typeOfDistributionFunction type="int">8</typeOfDistributionFunction>
<numberOfModeOfDistribution type="int">3</numberOfModeOfDistribution>
<numberOfDistributionFunctionParameters type="int">2</numberOfDistributionFunctionParameters>
</aerosol>
<aerosol id="acc_emiss">
<tracername_list type="char">dust_insol_acc,dust_insol_coa,na_sol_acc,na_sol_coa,cl_sol_acc,cl_sol_coa,nmb_sol_acc,nmb_sol_coa,nmb_insol_acc,nmb_insol_coa</tracername_list>
<parameterNumber_list type="char">3,3,3,3,3,3,3,3,3,3</parameterNumber_list>
<productDefinitionTemplate_list type="char">67,67,67,67,67,67,67,67,67,67</productDefinitionTemplate_list>
<typeOfStatisticalProcessing type="int">1</typeOfStatisticalProcessing>
<bitsPerValue type="int">32</bitsPerValue>
<typeOfDistributionFunction type="int">8</typeOfDistributionFunction>
<numberOfModeOfDistribution type="int">3</numberOfModeOfDistribution>
<numberOfDistributionFunctionParameters type="int">2</numberOfDistributionFunctionParameters>
</aerosol>
</diagnostics>
</syntaxhighlight>
</div>
=== Defining the Aerosol Emissions ===

Aerosol emissions can be defined via the aero_emiss.xml. Here you do not define the modes themselves but the distribution of the aerosols. ICON-ART then maps the given distribution on to the best fitting modes.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ aero_emiss.xml Parameters
! Parameter || type || description || example
|-
! routine id
| char || define the name of the emission routine || dust
|-
! nmodes
| int || number of aerosol modes to be emitted. The "X" in the following is to replaced by the numbers 1 to nmodes. || 3
|-
! d_g0_X
| real|| median diameter weighted by number concentration. Has to be defined for every mode given by nmodes. ||6.445E-7
|-
! d_g3_X
| real || median diameter weighted by mass concentration. Has to be defined for every mode given by nmodes. || 1.500E-6
|-
! rho
| real || Density of the emitted aerosol || 2.2E3
|-
! substances
| char || name of the tracers about to be emitted|| na,cl
|-
|}

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of aero_emiss.xml: ====
<syntaxhighlight class="mw-collapsible-content" lang=xml line>
<emiss>
<routine id="dust">
<nmodes type="int">3</nmodes>
<d_g0_1 type="real">6.445E-7</d_g0_1>
<d_g3_1 type="real">1.500E-6</d_g3_1>
<sigma_g_1 type="real">1.700E+0</sigma_g_1>
<d_g0_2 type="real">3.454E-6</d_g0_2>
<d_g3_2 type="real">6.700E-6</d_g3_2>
<sigma_g_2 type="real">1.600E+0</sigma_g_2>
<d_g0_3 type="real">8.672E-6</d_g0_3>
<d_g3_3 type="real">1.420E-5</d_g3_3>
<sigma_g_3 type="real">1.500E+0</sigma_g_3>
<rho type="real">2.650E3</rho>
<substances type="char">dust</substances>
</routine>
<routine id="seas_mode1">
<nmodes type="int">1</nmodes>
<d_g0_1 type="real">0.100E-6</d_g0_1>
<d_g3_1 type="real">0.433E-6</d_g3_1>
<sigma_g_1 type="real">1.900E+0</sigma_g_1>
<rho type="real">2.2E3</rho>
<substances type="char">na,cl</substances>
</routine>
<routine id="seas_mode2">
<nmodes type="int">1</nmodes>
<d_g0_1 type="real">3.000E-6</d_g0_1>
<d_g3_1 type="real">1.268E-5</d_g3_1>
<sigma_g_1 type="real">2.000E+0</sigma_g_1>
<rho type="real">2.2E3</rho>
<substances type="char">na,cl</substances>
</routine>
</emiss>
</syntaxhighlight>
</div>

Input

2023-08-09T08:33:15Z

Luca R: syntaxhighlight for xmls

File:Nature14956.pdf

2023-08-03T10:35:33Z

Luca R: test pdf upload

== Summary ==
test pdf upload

Input

2023-07-20T14:36:09Z

Luca R: added column ""Further explainations" in xml table

Input

2023-07-20T14:20:19Z

Luca R: moved xml tables around and testing for working in-site links

Input

2023-07-20T08:58:25Z

Luca R: Added Julia's XML tables at the end of the file

Tutorial Examples

2023-07-13T12:11:09Z

Luca R: added Pointsource , hid tutorial slides

- under construction! -

In this article several examples from different application areas in ICON-ART are presented. Just click on the highlighted text to get to the page of the respecting example. For every example a configuration is explained to try out typical cases in ICON-ART yourselves. Thereby, these pages can be viewed as a Tutorial in ICON-ART.

== Overview of topics ==

=== Atmospheric Chemistry Simulations ===
==== Simplified Chemistry ====
* here a simulation of the Hydroxylradical (OH) is performed with simplified chemistry
* this example teaches you...
** the basics of setting up an ICON runscript with ICON-ART settings
** the use of the most simple ICON-ART namelist parameter
** the implementation of the desired chemical species in a simulation by setting up a chemtracer xml-data for simplified chemistry simulations
** the implementation of emission data in a simulation
* click here to see more details

==== MECCA-based (full) Chemistry ====
* here a simulation with (full) MECCA-based chemistry is performed, for that a complete reaction mechanism is created and transferred to ICON-ART
* this example teaches you...
** the implementation of (full) MECCA-based chemistry in ICON-ART
** the creation of a chemical mechanism and the selection of the respective desired chemical species and their reactions in MECCA
** the implementation of not yet in MECCA implemented reactions in your mechanism
** the creation of the to your mechanism belonging Mecca-xml data to link MECCA calculations with ICON-ART
* click here to see more details

==== Stratospheric Ozone chemistry (with linearized Ozone, LINOZ) ====
* Stratospheric Ozone chemistry (with linearized Ozone, LINOZ)

==== Lifetime tracer Simulation ====
* Lifetime tracer Simulation

---------
=== Aerosol Simulations ===
==== [[Simulating a Point Source]] ====
* In this example a Point Source is installed at the location of the volcano Raikoke and a constant emission of SO2 and Ash is emitted into the atmosphere.
* this example teaches you...
** How to setup and run your first Simulation
** The setup of the directory structure when doing an ICON-ART Simulation
** How to modify a runscript to suit your needs
** How to modify .xml data to set up a Simulation
** A possible approach to visualise the output data

------
<div class="toccolours mw-collapsible mw-collapsed">
== Slides ICON-Training ==
<div class="mw-collapsible-content">

[[File:TD_slides_A.png|left]]
[[File:TD_slides_B.png|left]]

</div>
</div>

Talk:AERODYN

2023-07-13T07:57:39Z

Luca R: Restructure announcement

== Planned Restructuring Announcement==
As more Articles containing example .xml's are written, the .xml's in the dropdowns here will be obsolete and are then going to be deleted.

== table 'defining tracers' ==
moment: this is the moment of the log-normal distribution (0 for number concentration and 3 for mass mixing ratio) --[[User:Julia B|Julia B]] ([[User talk:Julia B|talk]]) 17:47, 20 June 2023 (CEST)

AERODYN

2023-07-13T07:54:39Z

Luca R: /* Defining Tracers */ added Julias comment about moment

AERODYN (AEROsol DYNamics) is a new aerosol dynamics module used in ICON-ART. With AERODYN, secondary aerosol formations (e.g. the accumulation of SO2 on a volcanic ash particle) can be accounted for.
Aerosol dynamic processes are the processes which change an aerosols properties after it has been released into the atmosphere. This change of properties (omitting chemical reactions) occurs either via condensation or coagulation.
Since it is not feasible to account for every individual aerosol to calculate their movement, coagulation and condensation, AERODYN uses a total of 12 categories for aerosol modes.

== Aerosol Modes in AERODYN ==

Traditionally, There are three categories to account for the aerosol sizes [https://www.sciencedirect.com/science/article/abs/pii/0004698178901968 Whitby et. Al. (1978)] : Aitken mode, accumulation mode and coarse mode. In AERODYN there is also a fourth mode, called giant mode, which accounts for especially large Aerosols, e.g. some species of pollen.

Since condensation can occur on a solid aerosol particle, as well as for a number of other reasons, it is possible to obtain aerosol modes which are a mixture of solid and liquid modes. To account for this, AERODYN combines 4 modes of aerosol size with the 3 modes soluble, insoluble and mixed.
A Visual Representation is given here:

[[File:AERODYNModes.png]]

when also adding the giant mode, the following configurations are possible:

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ AERODYN mode configurations
! || aitken || accumulation || coarse || giant
|-
! soluble
| sol_aitken || sol_acc || sol_coarse || sol_giant
|-
! insoluble
|insol_aitken || insol_acc || insol_coarse || insol_giant
|-
! mixed
| mixed_aitken ||mixed_acc || mixed_coarse || mixed_giant
|-
|}

== Condensation, Nucleation and Accumulation Processes ==
Here is a quick reminder of the different aerosol processes:

* Condensation is the process of a gas turning into a liquid. In practice, some sort of surface, or condensation nucleus is needed to enable the phase change.
* Nucleation is the process of an aerosol particle forming from a gaseous phase by random collisions of the gas molecules. The produced aerosol sizes are in the Aitken mode.
* Accumulation in the Process of the initial nucleus growing in size through collision and deposition Processes.

A schematic overview of the different microphysical aerosol processes can be found in the following Figure from [https://doi.org/10.1016/S1352-2310(00)00239-9 Raes et. Al(2000)]:

[[File:Raes2000.jpg| Overview of Aerosol Processes.]]

In AERODYN, some modes are treated differently with respect to the dynamic processes: Nucleation is only considered for the Aitken mode, and the Giant mode is not affected by condensation. This is done to save computing power while still retaining realism, and since Nucleation of larger Particles without first having Aitken mode Aerosols is highly unlikely, this is a good approximation. Also, as condensation highly prefers smaller modes to the giant mode, this also is a reasonable way to save computing power. For details how the processes are calculated, see [https://doi.org/10.1029/2018RG000615 Riemer et Al. (2003)] and [https://doi.org/10.5194/acp-9-8661-2009 Vogel et. Al. (2009)].

== Defining Aerosols with AERODYN ==

AERODYN is meant to be as flexible as possible, allowing the user to quickly define their own aerosols and then use ICON-ART to simulate their dynamics.
Aerosols are defined in .xml files, which can then be included in the run by specifying their path in the runfile (see [[Input]] for more details)

There are several Types of .xml files for different applications. The most important/most used .xml files are compiled in the following paragraphs.

=== Defining Tracers ===

Below is the tracers.xml, where the aerosols as well as their properties are defined. The properties to be defined are listed in the following table.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ tracers.xml Parameters
! Parameter || type || description || example
|-
! aerosol id
| char ||specifies name of the aerosol || "my_first_aerosol"
|-
! moment
| int || this is the moment of the log-normal distribution (0 for number concentration and 3 for mass mixing ratio)|| 3
|-
! mode
| char || mode name from AERODYN mode configurations (multiple modes can be given, seperated by a comma) || insol_giant,sol_aitken
|-

! mol_weight
| real || specify the molar mass || 30E-3
|-
! rho
| real || specify the density in terms of <code> unit </code> || 2e3

|-
! unit
| char || specify what unit the aerosol will have || kg-1
|-
! transport
| char || specify transport scheme to be used || hadv52aero
|-
|}

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of tracers.xml: ====
<pre class="mw-collapsible-content">
<tracers>
<aerosol id="nmb">
<moment type="int">0</moment>
<mode type="char">insol_acc,insol_coa,sol_acc,sol_coa,mixed_acc,mixed_coa</mode>
<unit type="char">kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="dust">
<moment type="int">3</moment>
<mode type="char">insol_acc,insol_coa,mixed_acc,mixed_coa</mode>
<sol type="real">0.0</sol>
<mol_weight type="real">50.00E-3</mol_weight>
<rho type="real">2.650E3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="na">
<moment type="int">3</moment>
<mode type="char">sol_acc,sol_coa,mixed_acc,mixed_coa</mode>
<sol type="real">1.</sol>
<mol_weight type="real">22.9898E-3</mol_weight>
<rho type="real">2.2E+3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="cl">
<moment type="int">3</moment>
<mode type="char">sol_acc,sol_coa,mixed_acc,mixed_coa</mode>
<sol type="real">1.</sol>
<mol_weight type="real">35.453E-3</mol_weight>
<rho type="real">2.2E+3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">hadv52aero</transport>
</aerosol>
<aerosol id="h2o">
<moment type="int">3</moment>
<mode type="char">sol_acc,sol_coa</mode>
<sol type="real">1.</sol>
<mol_weight type="real">18.01528E-3</mol_weight>
<rho type="real">1.E3</rho>
<unit type="char">mug kg-1</unit>
<transport type="char">off</transport>
</aerosol>
</tracers>
</pre>
</div>

=== Defining Modes ===

In the modes.xml the modes which are going to be used in the simulation can be defined. The following .xml contains the same information as Figure 3.1, but in a way the computer can interpret it.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ modes.xml Parameters
! Parameter || type || description || example
|-
! aerosol id
| char ||specifies name of the aerosol || "my_first_aerosol"
|-
! kind
| char || specifies which properties are to be calculated in the model (eg 2mom -> mass mixing ratio and number concentration) || 2mom
|-
! d_gn
| real || median diameter of the lognormal aerosol distribution || 2.0E-6
|-
! sigma_g
| real || standard deviation of the lognormal aerosol distribution || 2.2
|-
! condensation
| int || switch to control if condensation can occur || 0
|-
! icoag
| int || switch to control if coagulation can occur || 0
|-
! shift2larger
| char || prevents mode from getting to large, then shifts size to given mode (mostly used for Aitken mode) || insol_acc
|-
! shift_diam
| real || defines the diameter where the aerosol is shifted to the given mode || 0.08E-6
|-
! shift2mixed
| char || shift to given mode when 5% of mass is liquid || mixed_coa
|-
|}

Sometimes it is necessary for the type of an aerosol to change its given mode if its diameter gets too large or a large amount of condensation occurs. The mass can then be transferred to defined modes, as seen in the table above.

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of modes.xml: ====
<pre class="mw-collapsible-content">
<modes>
<aerosol id="sol_acc">
<kind type="char">2mom</kind>
<d_gn type="real">0.2E-6</d_gn>
<sigma_g type="real">2.0</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="insol_acc">
<kind type="char">2mom</kind>
<d_gn type="real">6.445E-7</d_gn>
<sigma_g type="real">1.700E+0</sigma_g>
<condensation type="int">0</condensation>
<shift2mixed type="char">mixed_acc</shift2mixed>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="mixed_acc">
<kind type="char">2mom</kind>
<d_gn type="real">6.445E-7</d_gn>
<sigma_g type="real">1.700E+0</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="sol_coa">
<kind type="char">2mom</kind>
<d_gn type="real">2.0E-6</d_gn>
<sigma_g type="real">2.2</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="insol_coa">
<kind type="char">2mom</kind>
<d_gn type="real">2.0E-6</d_gn>
<sigma_g type="real">2.2</sigma_g>
<condensation type="int">0</condensation>
<shift2mixed type="char">mixed_coa</shift2mixed>
<icoag type="int">1</icoag>
</aerosol>
<aerosol id="mixed_coa">
<kind type="char">2mom</kind>
<d_gn type="real">2.0E-6</d_gn>
<sigma_g type="real">2.2</sigma_g>
<condensation type="int">0</condensation>
<icoag type="int">1</icoag>
</aerosol>
</modes>
</pre>
</div>

=== Defining Coagulation ===

Here the behaviour of collision and subsequent merging (coagulation) between aerosols can be defined. If for example a solid (insol) aerosol merges with a liquid (sol) aerosol, the result should logically be a mixed aerosol, and so on. This can further be defined for all sizes of aerosols as seen in the example .xml for accumulation and coarse modes.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ coagulate.xml Parameters
! Parameter || type || description || example
|-
!smallmode id
| char || The aerosol_mode which is gonig to coagulate || insol_acc
|-
!nmodes
| int || Defines the number of modes for which the coagulation preocesses are defined || 5
|-
!bigmode id
| char || Inside angle brackets : coagulation partner of the given aerosol mode. Outside angle brackets: the new mode of the cagulated aerosols. || <bigmode id="sol_acc" type="char">mixed_acc</bigmode>
|-
|}

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of coag.xml: ====
<pre class="mw-collapsible-content">
<coagulate>
<smallmode id="sol_acc">
<nmodes type="int">6</nmodes>
<bigmode id="sol_acc" type="char">sol_acc</bigmode>
<bigmode id="insol_acc" type="char">mixed_acc</bigmode>
<bigmode id="mixed_acc" type="char">mixed_acc</bigmode>
<bigmode id="sol_coa" type="char">sol_coa</bigmode>
<bigmode id="insol_coa" type="char">mixed_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="insol_acc">
<nmodes type="int">5</nmodes>
<bigmode id="insol_acc" type="char">insol_acc</bigmode>
<bigmode id="mixed_acc" type="char">mixed_acc</bigmode>
<bigmode id="sol_coa" type="char">mixed_coa</bigmode>
<bigmode id="insol_coa" type="char">insol_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="mixed_acc">
<nmodes type="int">4</nmodes>
<bigmode id="mixed_acc" type="char">mixed_acc</bigmode>
<bigmode id="sol_coa" type="char">mixed_coa</bigmode>
<bigmode id="insol_coa" type="char">mixed_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="sol_coa">
<nmodes type="int">3</nmodes>
<bigmode id="sol_coa" type="char">sol_coa</bigmode>
<bigmode id="insol_coa" type="char">mixed_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="insol_coa">
<nmodes type="int">2</nmodes>
<bigmode id="insol_coa" type="char">insol_coa</bigmode>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
<smallmode id="mixed_coa">
<nmodes type="int">1</nmodes>
<bigmode id="mixed_coa" type="char">mixed_coa</bigmode>
</smallmode>
</coagulate>
</pre>
</div>
=== Defining the Model Output ===
In the diagnostics.xml the diagnostic variables that can be written out are defined. Below are some examples, however this is a heavily shortened version of the diagnostics.xml used in ICON-ART.

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of diagnostics.xml: ====
<pre class="mw-collapsible-content">
<diagnostics>
<aerosol id="diam_insol_acc">
<productDefinitionTemplate type="int">40</productDefinitionTemplate>
<discipline type="int">0</discipline>
<parameterCategory type="int">254</parameterCategory>
<parameterNumber type="int">201</parameterNumber>
<constituentType type="int">62001</constituentType>
</aerosol>
<aerosol id="acc_drydepo">
<tracername_list type="char">dust_insol_acc,dust_insol_coa,na_sol_acc,na_sol_coa,cl_sol_acc,cl_sol_coa,nmb_sol_acc,nmb_sol_coa,nmb_insol_acc,nmb_insol_coa</tracername_list>
<parameterNumber_list type="char">6,6,6,6,6,6,192,192,192,192</parameterNumber_list>
<productDefinitionTemplate_list type="char">67,67,67,67,67,67,67,67,67,67</productDefinitionTemplate_list>
<typeOfStatisticalProcessing type="int">1</typeOfStatisticalProcessing>
<bitsPerValue type="int">32</bitsPerValue>
<typeOfDistributionFunction type="int">8</typeOfDistributionFunction>
<numberOfModeOfDistribution type="int">3</numberOfModeOfDistribution>
<numberOfDistributionFunctionParameters type="int">2</numberOfDistributionFunctionParameters>
</aerosol>
<aerosol id="acc_emiss">
<tracername_list type="char">dust_insol_acc,dust_insol_coa,na_sol_acc,na_sol_coa,cl_sol_acc,cl_sol_coa,nmb_sol_acc,nmb_sol_coa,nmb_insol_acc,nmb_insol_coa</tracername_list>
<parameterNumber_list type="char">3,3,3,3,3,3,3,3,3,3</parameterNumber_list>
<productDefinitionTemplate_list type="char">67,67,67,67,67,67,67,67,67,67</productDefinitionTemplate_list>
<typeOfStatisticalProcessing type="int">1</typeOfStatisticalProcessing>
<bitsPerValue type="int">32</bitsPerValue>
<typeOfDistributionFunction type="int">8</typeOfDistributionFunction>
<numberOfModeOfDistribution type="int">3</numberOfModeOfDistribution>
<numberOfDistributionFunctionParameters type="int">2</numberOfDistributionFunctionParameters>
</aerosol>
</diagnostics>
</pre>
</div>
=== Defining the Aerosol Emissions ===

Aerosol emissions can be defined via the aero_emiss.xml. Here you do not define the modes themselves but the distribution of the aerosols. ICON-ART then maps the given distribution on to the best fitting modes.

<div id="OutputTable"></div>
::{| class="wikitable" style="text-align:left;"
|+ aero_emiss.xml Parameters
! Parameter || type || description || example
|-
! routine id
| char || define the name of the emission routine || dust
|-
! nmodes
| int || number of aerosol modes to be emitted. The "X" in the following is to replaced by the numbers 1 to nmodes. || 3
|-
! d_g0_X
| real|| median diameter weighted by number concentration. Has to be defined for every mode given by nmodes. ||6.445E-7
|-
! d_g3_X
| real || median diameter weighted by mass concentration. Has to be defined for every mode given by nmodes. || 1.500E-6
|-
! rho
| real || Density of the emitted aerosol || 2.2E3
|-
! substances
| char || name of the tracers about to be emitted|| na,cl
|-
|}

<div class="toccolours mw-collapsible mw-collapsed">
==== Example File of aero_emiss.xml: ====
<pre class="mw-collapsible-content">
<emiss>
<routine id="dust">
<nmodes type="int">3</nmodes>
<d_g0_1 type="real">6.445E-7</d_g0_1>
<d_g3_1 type="real">1.500E-6</d_g3_1>
<sigma_g_1 type="real">1.700E+0</sigma_g_1>
<d_g0_2 type="real">3.454E-6</d_g0_2>
<d_g3_2 type="real">6.700E-6</d_g3_2>
<sigma_g_2 type="real">1.600E+0</sigma_g_2>
<d_g0_3 type="real">8.672E-6</d_g0_3>
<d_g3_3 type="real">1.420E-5</d_g3_3>
<sigma_g_3 type="real">1.500E+0</sigma_g_3>
<rho type="real">2.650E3</rho>
<substances type="char">dust</substances>
</routine>
<routine id="seas_mode1">
<nmodes type="int">1</nmodes>
<d_g0_1 type="real">0.100E-6</d_g0_1>
<d_g3_1 type="real">0.433E-6</d_g3_1>
<sigma_g_1 type="real">1.900E+0</sigma_g_1>
<rho type="real">2.2E3</rho>
<substances type="char">na,cl</substances>
</routine>
<routine id="seas_mode2">
<nmodes type="int">1</nmodes>
<d_g0_1 type="real">3.000E-6</d_g0_1>
<d_g3_1 type="real">1.268E-5</d_g3_1>
<sigma_g_1 type="real">2.000E+0</sigma_g_1>
<rho type="real">2.2E3</rho>
<substances type="char">na,cl</substances>
</routine>
</emiss>
</pre>
</div>