Difference between revisions of "UKCA Chemistry and Aerosol UMvn13.0 Tutorial 1"

From UKCA
 
(36 intermediate revisions by the same user not shown)
Line 7: Line 7:
 
| Time to Complete || '''Under 1 hour'''
 
| Time to Complete || '''Under 1 hour'''
 
|-
 
|-
| Video instructions ||
+
| Video instructions || [https://www.youtube.com/watch?v=n3yE-oYy-gU&t=89s Walkthrough (YouTube)]
 
|}
 
|}
   
'''Remember to run <tt>mosrs-cache-password</tt> in your Rose suite and UM branch terminals.'''
+
'''Remember to run <tt>mosrs-cache-password</tt>.'''
  +
  +
==What you will learn in this tutorial==
  +
  +
In this tutorial you will learn how to copy a UKCA Box Model suite, run it, and plot the output.
   
 
==Copying and Running an Existing Rose Suite==
 
==Copying and Running an Existing Rose Suite==
Line 23: Line 27:
 
|}
 
|}
   
If you are using your own '''Met Office Virtual Machine''', you will need to login to the VM, e.g.
+
You will need to login to the VM using SSH via a Terminal
   
  +
ssh -X ukca
vagrant ssh
 
   
  +
When you first login to your virtual machine you should then be asked for your Met Office Science Repository Service (MOSRS) password. Once you enter this, if this is the first time you are logging in, you will also be asked for your MOSRS username. This will likely be your name in lower case, and is '''not''' your username on the VM.
If you are using a VM provided for you, you should connect via X2Go or via a Terminal.
 
 
You should then be asked for your Met Office Science Repository Service (MOSRS) password. Once you enter this, if this is the first time you are logging in, you will also be asked for your MOSRS username. This will likely be your name in lower case, and is '''not''' your username on the VM.
 
   
 
Then copy the rose suite by:
 
Then copy the rose suite by:
Line 62: Line 64:
 
rose suite-run
 
rose suite-run
   
in the terminal.
+
in the terminal. For this suite to run you need the UKCA branch
   
  +
fcm:ukca.x_br/dev/lukeabraham/um13.0_ukca_box
If you are using your own VM, times will vary depending on the specifications of the host. It could take about 10 minutes to compile the code (''fcm_make'') for the first time (and about 1-2 minutes when recompiling), followed by about 1 minute to run the reconfiguration step (''recon''), and then around 10 minutes to run the UM itself (''atmos''). If you have been provided a VM the compile should take 1-2 minutes.
 
  +
  +
checked-out into your <tt>home/</tt> directory. This has been done for you.
  +
  +
If you are using your own VM, times will vary depending on the specifications of the host. It may take 2 or 3 minutes to compile the code (''fcm_make'') for the first time (and a few seconds when recompiling), followed by another few seconds to run UKCA itself (''ukca''). Any unsuccessful jobs will be highlighted in red and have a status of "failed".
   
 
When the suite has finished successfully it will then become blank with the message '''stopped with 'succeeded'''' in the bottom-left corner.
 
When the suite has finished successfully it will then become blank with the message '''stopped with 'succeeded'''' in the bottom-left corner.
Line 70: Line 76:
 
==Version Control==
 
==Version Control==
   
Rose suites are all held under version control, using [http://cms.ncas.ac.uk/wiki/Fcm fcm]. When making changes to a suite, you will need to '''save''' it before you can '''run''' the suite. Once you are happy with the settings, you can also '''commit''' these changes back to the repository - to do this change directory to the
+
Rose suites are all held under version control, using [https://www.metoffice.gov.uk/research/weather/weather-science-it/fcm fcm]. When making changes to a suite, you will need to '''save''' it before you can '''run''' the suite. Once you are happy with the settings, you can also '''commit''' these changes back to the repository - to do this change directory to the
   
 
/home/'''$USER'''/roses/'''[SUITE-ID]'''
 
/home/'''$USER'''/roses/'''[SUITE-ID]'''
Line 90: Line 96:
 
</gallery>
 
</gallery>
   
The output directory structure of rose suites are rather complex. The schematic gives a broad overview of the general structure. Everything can be found within the '''<tt>[SUITE-ID]</tt>''' (e.g. ''u-cq782'' etc.) directory, which can be found within your '''<tt>$HOME/cylc-run</tt>''' directory. On the VM, everything in held in this single place, but things are more complicated when using other machines, such as PUMA & ARCHER2, as some files are held on the ARCHER2 ''<tt>/home</tt>'' filesystem, some on the ''<tt>/work</tt>'' filesystem, and some files are synced back to PUMA as well. On Monsoon2, everything is held on the same machine, but with some files in <tt>/projects</tt> and some in <tt>/home</tt>.
+
The output directory structure of rose suites are rather complex. The schematic gives a broad overview of the general structure. Everything can be found within the '''<tt>[SUITE-ID]</tt>''' (e.g. ''u-cq782'' etc.) directory, which can be found within your '''<tt>$HOME/cylc-run</tt>''' directory.
   
 
Within this directory there are several directories, including:
 
Within this directory there are several directories, including:
   
 
* '''<tt>log</tt>''': a symbolic link to a directory called something like <tt>log.20220914T092709Z</tt>, which contains all the job.out (containing output from UM <tt>WRITE</tt> statements from the <tt>umPrint</tt> subroutine) and job.err files etc., as well as the script used to run the job.
 
* '''<tt>log</tt>''': a symbolic link to a directory called something like <tt>log.20220914T092709Z</tt>, which contains all the job.out (containing output from UM <tt>WRITE</tt> statements from the <tt>umPrint</tt> subroutine) and job.err files etc., as well as the script used to run the job.
* '''<tt>work</tt>''': this contains the directories used when the job actually runs. Real-time output will be sent to files here (held in a <tt>pe_output</tt> directory. Sometimes model output will also be here. On ARCHER this directory exists on the <tt>/work</tt> filesystem.
+
* '''<tt>work</tt>''': this contains the directories used when the job actually runs. Real-time output will be sent to files here (held in a <tt>pe_output</tt> directory. Sometimes model output will also be here.
* '''<tt>share</tt>''': This isn't used by output from the UKCA Training Suites, but for UM climate jobs model output will usually be held in this directory. On ARCHER2 this directory exists on the <tt>/work</tt> filesystem.
 
   
  +
Within these directories there will be many sub-directories. Some of these will be named from the '''cycle-point''' (labelled '''<tt>[CYCLE]</tt>''' in the graphic. For the UKCA Training Suite this will be '''1''', but under normal climate operation this will be a date string, e.g. '''1''' or '''19880901T0000Z'''. The sub-directories will (eventually) be named after the '''app''' (labelled by '''<tt>[JOB NAME]</tt>''' in the table below) that the output is from, e.g. <tt>fcm_make</tt>, <tt>ukca</tt>, <tt>recon</tt>, <tt>atmos</tt> etc.
 
  +
Within these directories there will be many sub-directories. Some of these will be named from the '''cycle-point''' (labelled '''<tt>[CYCLE]</tt>''' in the graphic. For the UKCA Training Suite this will be '''1'''. The sub-directories will (eventually) be named after the '''app''' (labelled by '''<tt>[JOB NAME]</tt>''' in the table below) that the output is from, e.g. <tt>fcm_make</tt>, <tt>ukca</tt> etc.
   
 
{| class="wikitable" border = "1px"
 
{| class="wikitable" border = "1px"
Line 105: Line 111:
 
!style="text-align:top;|Example
 
!style="text-align:top;|Example
 
|-
 
|-
| Most recent ''job.out'' files || '''vm:''' <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/log/job/1/'''[JOB NAME]'''/NN</tt> || <tt>/home/vagrant/cylc-run/u-cq782/log/job/1/ukca/NN</tt>
+
| Most recent ''job.out'' files || <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/log/job/1/'''[JOB NAME]'''/NN</tt> || <tt>/home/vagrant/cylc-run/u-cq782/log/job/1/ukca/NN</tt>
 
|-
 
|-
| Processor output (while running) || '''vm:''' <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/work/1/'''[JOB NAME]'''/pe_output</tt> || <tt>/home/vagrant/cylc-run/u-cq782/work/1/ukca/pe_output</tt>
+
| Processor output (while running) || <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/work/1/'''[JOB NAME]'''/pe_output</tt> || <tt>/home/vagrant/cylc-run/u-cq782/work/1/ukca/pe_output</tt>
 
|-
 
|-
| UKCA Box Model output text files in ASCII format<br/>e.g. <tt>tracer_out.csv</tt> and <tt>flux_out.csv</tt> || '''vm:''' <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/work/1/ukca</tt> || <tt>/home/vagrant/cylc-run/u-cq782/work/1/ukca</tt>
+
| UKCA Box Model output text files in ASCII format<br/><tt>tracer_out.csv</tt><br/><tt>flux_out.csv</tt> || <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/work/1/ukca</tt> || <tt>/home/vagrant/cylc-run/u-cq782/work/1/ukca</tt>
|-
 
| Output files in 64-bit fieldsfile format || '''vm:''' <tt>/home/vagrant/cylc-run/'''[SUITE-ID]'''/work/1/atmos</tt> || <tt>/home/vagrant/cylc-run/u-cq148/work/1/atmos</tt>
 
 
|}
 
|}
   
Line 117: Line 121:
   
 
* The <tt>job.out</tt> files can also be viewed through the Gcylc GUI right-click menu from each job.
 
* The <tt>job.out</tt> files can also be viewed through the Gcylc GUI right-click menu from each job.
* UKCA Box Model output file will be <tt>tracer_out.csv</tt> for tracers and <tt>flux_out.csv</tt> for chemical reaction fluxes. These can be renamed in Rose in the '''ukca <math>\rightarrow</math> Model Input and Output''' panel.
+
* UKCA Box Model output file will be <tt>tracer_out.csv</tt> for tracers and <tt>flux_out.csv</tt> for chemical reaction fluxes. These can be renamed in Rose in the '''ukca <math>\rightarrow</math> namelist <math>\rightarrow</math> Model Input and Output''' panel.
  +
* A handy command to check the progress of the <tt>atmos</tt> job is
 
  +
==Exploring Rose==
: '''vm:''' <tt>tail -1000f atmos.fort6.pe0 | grep Atm_Step</tt>
 
  +
* Other suites may put the output fieldsfiles into a different location, e.g.
 
  +
Take a look at your suite, and click through the options. There are 2 key sections that you will need to consider:
:: <tt>/work/n02/n02/'''$USER'''/cylc-run/'''[SUITE-ID]'''/share/data/History_Data</tt>
 
  +
: and there are also post-processing settings that can be used to copy data to the RDF or JASMIN as the suite runs.
 
  +
===<tt>fcm_make</tt>===
* Output fieldsfiles can have various naming conventions. For these suites, they will be '''atmosa.pa19810901_00''', but more generally they are likely to be of the form
 
  +
: <tt>''[SHORT SUITE-ID]''a.p'''[abcdefghijkmsyx]'''YYYYMMDD</tt>
 
  +
This section controls which fcm branches are used by the model, and also what compiler settings are used etc. You should go to '''fcm_make <math>\rightarrow</math> env <math>\rightarrow</math> Sources''' to see what branches are being used. During the course of these tutorials you will make a branch and add it here to complete the tutorial tasks. You can add both relative paths to the repository (starting with <tt>branches/</tt>) and absolute paths to a working copy (starting with <tt>/home/</tt>).
:: e.g. ag308a.pk19880901
 
  +
: <tt>''[SHORT SUITE-ID]''a.p'''[abcdefghijkm]'''YYYYmon</tt>
 
  +
===<tt>ukca</tt>===
:: e.g. ak468a.pe1989oct
 
* Climate-mean files ('''<tt>.pm</tt>''') will often also have seasonal ('''<tt>.ps</tt>'''), annual ('''<tt>.py</tt>'''), and decadal ('''<tt>.px</tt>''') equivalents.
 
   
  +
Click the arrow by '''ukca''' to view the list of sections. We'll be most interested in the '''namelist''' section, which contains all the settings that configure the science that can be done by UKCA. Take a look through the panels and sub-panels in this section to get a feel for how it is organised.
'''Whilst this suite is running, take a look at [[UKCA Chemistry and Aerosol vn11.8 Tutorial 2|Tutorial 2: Exploring Rose]].'''
 
   
 
==Viewing Output==
 
==Viewing Output==
   
===UKCA Box Model===
+
===gnuplot===
   
  +
<gallery widths=400px heights=350px>
  +
Image:T01_gnuplot.png|Plotting O3 concentrations using Gnuplot
  +
Image:T01_o3_002.png|O3 mass-mixing ratio plotted using gnuplot
  +
</gallery>
   
  +
The quickest and easiest way to plot UKCA box model output is using '''[http://www.gnuplot.info/ Gnuplot]''', which is a relatively simple command-line plotting tool. The box model outputs tracer concentrations and chemical reaction fluxes every timestep to the files '''<tt>tracer_out.csv</tt>''' and '''<tt>flux_out.csv</tt>''' respectively.
===Met Office Unified Model===
 
   
  +
To plot these files, change into your suite's cylc-run directory and navigate to the ukca directory in work, e.g.
[[Image:Vn10.9_base_output.png|thumb|300px|right|Example output from the UKCA training suite.]]
 
   
  +
cd /home/vagrant/cylc-run/u-cq782/work/1/ukca
For more detailed plotting, tools such as the [http://scitools.org.uk/iris/ Iris] and [http://cfpython.bitbucket.org/ cf-python] libraries can be used to view UM file formats directly. However, for quick viewing, [http://cms.ncas.ac.uk/documents/xconv/ Xconv] is a very useful tool. Information on how to install Xconv on the VM can be found [[UKCA_Chemistry_and_Aerosol_Tutorials:_Things_to_know_before_you_start_the_vn11.8_practicals#Xconv|here]].
 
   
  +
then load gnuplot by typing
To view these files, do
 
   
  +
gnuplot
'''xconv atmosa.pa19810901_00'''
 
   
  +
To plot the tracer file, type
As well as viewing files, you can use Xconv to convert these files to netCDF, by filling in the '''''Output file name:''''' box (e.g. ''foo.nc''), and then clicking '''convert'''. If no path is defined, this will save the file in the same directory that you opened Xconv from.
 
   
  +
plot "tracer_out.csv" using 1:3 with linespoints
Example output from the UKCA training suite can be found at
 
  +
  +
which can be shortened to
  +
  +
p "tracer_out.csv" u 1:3 w lp
  +
  +
The file can be saved from the export menu on the plot itself. Additional options for the [http://gnuplot.info/docs_5.5/loc9418.html <tt>set</tt> command] can be used to define labels, titles, change how the axis labels look, e.g.
  +
  +
gnuplot> set title "O3 MMR"
  +
gnuplot> set xlabel "timestep number"
  +
gnuplot> set ylabel "mass mixing ratio /kg kg^{-1}"
  +
gnuplot> set format y "%.2t*10^{%+03T}"
  +
gnuplot> unset key
  +
gnuplot> p 'tracer_out.csv' u 1:3 w lp
  +
  +
The ordering of the species in the <tt>tracer_out.csv</tt> can be seen by opening it in e.g. a simple spreadsheet program like '''[http://www.gnumeric.org/ Gnumeric]'''. You may need to install this first using
  +
  +
sudo apt install gnumeric
  +
  +
before opening the file
  +
  +
gnumeric tracer_out.csv
  +
  +
The column number is given in the first line with the tracer name in the line below. You could also open the file in <tt>emacs</tt> and then <tt>M-x toggle-truncate-lines</tt>. An example of this header for the StratTrop chemistry with GLOMAP-mode aerosols is
  +
  +
<pre style="white-space: pre;">
  +
# COL 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96
  +
# Timestep, O(3P), O3, N, NO, NO3, NO2, N2O5, HO2NO2, HONO2, H2O2, CH4, CO, HCHO, MeOOH, H, H2O, OH, HO2, Cl, Cl2O2, ClO, OClO, Br, BrO, BrCl, BrONO2, N2O, HCl, HOCl, HBr, HOBr, ClONO2, CFCl3, CF2Cl2, MeBr, HONO, C2H6, EtOOH, MeCHO, PAN, C3H8, n-PrOOH, i-PrOOH, EtCHO, Me2CO, MeCOCH2OOH, PPAN, MeONO2, C5H8, ISOOH, ISON, MACR, MACROOH, MPAN, HACET, MGLY, NALD, HCOOH, MeCO3H, MeCO2H, H2, MeOH, DMS, SO2, H2SO4, MSA, DMSO, NH3, CS2, COS, H2S, SO3, Monoterp, Sec_Org, Nuc_SOL_N, Nuc_SOL_SU, Ait_SOL_N, Ait_SOL_SU, Ait_SOL_BC, Ait_SOL_OM, Acc_SOL_N, Acc_SOL_SU, Acc_SOL_BC, Acc_SOL_OM, Acc_SOL_SS, Cor_SOL_N, Cor_SOL_SU, Cor_SOL_BC, Cor_SOL_OM, Cor_SOL_SS, Ait_INS_N, Ait_INS_BC, Ait_INS_OM, Nuc_SOL_OM, PASSIVE O3
  +
</pre>
  +
  +
The <tt>flux_out.csv</tt> file has a more complicated header as it also gives all the chemical reactants and products, e.g.
  +
  +
<pre style="white-space: pre;">
  +
# Timestep, UKCA reaction fluxes (units = molecules.cm^-3.s^-1)
  +
# COL 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299, 300, 301, 302, 303, 304, 305, 306, 307, 308, 309, 310
  +
# PROD 1, Cl2O2, HO2NO2, PAN, PPAN, MPAN, N2O5, EtOOH, H2O2, HCHO, HCHO, HO2NO2, HONO2, MeCHO, MeCHO, MeOOH, N2O5, NO2, NO3, NO3, O2, O3, O3, PAN, HONO, EtCHO, Me2CO, n-PrOOH, i-PrOOH, MeCOCH2OOH, PPAN, MeONO2, ISOOH, ISON, MACR, MPAN, MACROOH, MACROOH, HACET, MGLY, NALD, MeCO3H, BrCl, BrO, BrONO2, BrONO2, O2, OClO, NO, HOBr, N2O, H2O, ClONO2, ClONO2, HCl, HOCl, Cl2O2, CFCl3, CF2Cl2, MeBr, CH4, CO2, HO2NO2, CS2, COS, H2SO4, SO3, Br, Br, Br, Br, Br, BrO, BrO, BrO, BrO, BrO, BrO, BrO, CF2Cl2, CFCl3, Cl, Cl, Cl, Cl, Cl, Cl, Cl, Cl, Cl, Cl, Cl, Cl, Cl, ClO, ClO, ClO, ClO, ClO, ClO, ClO, EtCO3, EtCO3, EtOO, EtOO, H, H, H, H, H, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, HO2, MeBr, MeBr, MeBr, MeCO3, MeCO3, MeCOCH2OO, MeCOCH2OO, MeOO, MeOO, MeOO, N, N, N, N2O5, NO, NO, NO, NO, NO, NO, NO2, NO2, NO3, NO3, NO3, NO3, NO3, NO3, NO3, O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(1D), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O(3P), O3, O3, O3, O3, O3, O3, O3, OClO, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, OH, i-PrOO, i-PrOO, n-PrOO, n-PrOO, CS2, CS2, DMS, DMS, DMSO, H2S, H2S, COS, COS, SO2, SO3, Monoterp, Monoterp, Monoterp, MeOO, MeOO, EtOO, EtOO, n-PrOO, n-PrOO, i-PrOO, i-PrOO, MeCO3, EtCO3, MeCOCH2OO, MeCOCH2OO, ISO2, MACRO2, MACRO2, HO2, O(3P), O(3P), O(3P), O(1D), BrO, ClO, ClO, H, HO2, HO2, OH, OH, OH, MeCO3, EtCO3, MACRO2, NO2, NO, SO2, DMS, DMS, ClONO2, ClONO2, HOCl, N2O5, N2O5, SO2, SO2, SO2
  +
# PROD 2, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Cl2O2, HCHO, HO2, O3, OClO, BrO, ClO, ClO, ClO, HO2, NO, OH, O(1D), O(1D), CH4, Cl2O2, ClONO2, H2, H2O2, HCHO, HO2, HO2, HOCl, MeOOH, NO3, O3, OClO, ClO, ClO, ClO, HO2, MeOO, NO, NO3, NO, NO3, NO, NO3, HO2, HO2, HO2, NO2, O3, EtCO3, EtCO3, EtOO, HO2, ISO2, MACRO2, MeCO3, MeCO3, MeCO3, MeCOCH2OO, MeOO, MeOO, NO, NO3, O3, i-PrOO, n-PrOO, Cl, O(1D), OH, NO, NO3, NO, NO3, NO, NO, NO3, NO, NO2, O2, H2O, ISO2, ISO2, MACRO2, MACRO2, NO3, O3, NO3, O3, Br, C5H8, EtCHO, HCHO, MGLY, Me2CO, MeCHO, CH4, CH4, CH4, CO2, H2, H2O, HBr, HBr, HCl, HCl, HCl, N2, N2O, N2O, O2, O3, O3, BrO, ClO, ClONO2, H2O2, HBr, HCl, HO2, HOCl, NO2, NO3, O3, OClO, OH, C5H8, C5H8, C5H8, MACR, MACR, MACR, MACR, NO, C2H6, C3H8, C3H8, C5H8, CH4, CO, ClO, ClO, ClONO2, EtCHO, EtOOH, EtOOH, H2, H2O2, HACET, HBr, HCHO, HCOOH, HCl, HO2, HO2NO2, HOCl, HONO, HONO2, ISON, ISOOH, MACR, MACR, MACROOH, MGLY, MPAN, Me2CO, Me2CO, MeCHO, MeCO2H, MeCO3H, MeCOCH2OOH, MeCOCH2OOH, MeOH, MeONO2, MeOOH, MeOOH, NALD, NO3, O3, OClO, OH, PAN, PPAN, i-PrOOH, i-PrOOH, n-PrOOH, n-PrOOH, NO, NO3, NO, NO3, O(3P), OH, OH, NO3, OH, O(3P), OH, O(3P), OH, O3, H2O, OH, O3, NO3, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, RO2, EtCO3, O2, NO, NO2, N2, NO2, ClO, NO2, O2, HO2, NO2, NO, NO2, OH, NO2, NO2, NO2, NO3, NO, OH, OH, OH, H2O, HCl, HCl, H2O, HCl, H2O2, O3, O3
  +
# REACT 1, ClO, NO2, MeCO3, EtCO3, MACRO2, NO2, MeCHO, OH, HO2, H2, HO2, OH, MeOO, CH4, HO2, NO3, NO, NO, NO2, O(3P), O2, O2, MeCO3, OH, EtOO, MeCO3, EtCHO, Me2CO, MeCO3, EtCO3, HO2, OH, NO2, MeCO3, MACRO2, OH, HACET, MeCO3, MeCO3, HCHO, MeOO, Br, Br, Br, BrO, O(3P), O(3P), N, OH, N2, OH, Cl, ClO, H, OH, Cl, Cl, Cl, Br, MeOO, CO, OH, COS, CO, SO3, SO2, BrCl, HBr, HBr, BrO, BrO, Br, Br, Br, BrCl, HOBr, Br, Br, Cl, Cl, HCl, Cl, Cl, HCl, HCl, HCl, ClO, HCl, Cl, HCl, ClO, ClO, ClO, Cl, Cl, Cl, HOCl, Cl, Cl, Cl, EtOO, EtOO, MeCHO, MeCHO, H2, O(3P), OH, OH, OH, O2, O3, EtOOH, H2O2, ISOOH, MACROOH, MeCO2H, MeCO3H, OH, MeCOCH2OOH, HCHO, MeOOH, OH, OH, OH, i-PrOOH, n-PrOOH, Br, Br, Br, MeOO, MeOO, MeCO3, MeCO3, HO2, MeONO2, HO2, N2, N2O, NO, HONO2, ISON, NO2, MGLY, NO2, NO2, NO2, NO, NO3, BrO, ISON, HONO2, HONO2, MeCO3, HONO2, HONO2, HCHO, HCHO, OH, O(3P), OH, OH, HBr, OH, H, O(3P), OH, O(3P), N2, NO, O(3P), O2, O2, O2, Cl, ClO, OH, OH, OH, OH, OH, NO, O2, O2, O2, O2, HO2, MACR, MeOO, MGLY, MGLY, OH, OH, NO2, H2O, i-PrOO, n-PrOO, ISO2, H2O, H, HCl, HO2, HOCl, H2O, H2O, H2O, H2O, H2O, MGLY, H2O, H2O, CO2, H2O, H2O, H2O, ClO, H2O, H2O, HACET, MACR, MACRO2, MACRO2, MACRO2, MeCO3, HACET, H2O, H2O, H2O, MeOO, MeCO3, H2O, OH, HO2, HCHO, H2O, H2O, HCHO, HO2, HO2, HOCl, H2O, HCHO, MeCHO, Me2CO, i-PrOO, EtCHO, n-PrOO, Me2CO, Me2CO, EtCHO, EtCHO, CO, COS, SO2, SO2, SO2, OH, SO2, CO, CO2, SO3, H2SO4, Sec_Org, Sec_Org, Sec_Org, HO2, MeOH, HO2, MeCHO, HO2, EtCHO, HO2, Me2CO, MeCO2H, , HCHO, HACET, MACR, HACET, HO2, OH, O3, NO2, NO3, N2O, BrONO2, Cl2O2, ClONO2, HO2, H2O2, HO2NO2, HONO, HONO2, H2O2, PAN, PPAN, MPAN, N2O5, NO2, SO3, DMSO, SO2, HOCl, Cl, Cl, HONO2, Cl, , ,
  +
# REACT 2, ClO, HO2, NO2, NO2, NO2, NO3, HO2, OH, HO2, CO, NO2, NO2, HO2, CO, HCHO, NO2, O(3P), O2, O(3P), O(3P), O(1D), O(3P), NO2, NO, HO2, MeOO, HO2, HO2, HCHO, NO2, HCHO, MACR, MACR, HCHO, NO2, HO2, CO, HCHO, CO, CO, OH, Cl, O(3P), NO3, NO2, O(1D), ClO, O(3P), Br, O(1D), H, NO3, NO2, Cl, Cl, Cl, Cl, Cl, H, H, O(3P), NO3, SO2, SO2, OH, O(3P), Cl, CO, O2, O2, ClO, Br, Cl, OClO, O2, O2, NO2, HO2, ClO, Cl, MeOO, Cl, Cl, H, HO2, CO, OH, O2, Cl, MeOO, NO2, O2, ClO, Cl, Cl, OClO, O2, HCHO, NO2, O2, CO2, CO2, HO2, HO2, O2, H2O, OH, NO, O2, , , , , , , O3, , MeOO, , , , NO2, NO2, O2, , , HCl, OH, H2O, CO2, CO2, HCHO, HCHO, HCHO, , HCHO, O(3P), O(3P), O(3P), HONO2, , MACR, HCHO, MeCO3, NO2, , NO2, , NO2, , EtCO3, HO2, CO, MeCOCH2OO, MeCO3, H2, HO2, MeOO, CO2, H, OH, O(3P), Br, ClO, HCl, Cl, N2, O2, NO, O2, O(3P), O2, Br, O2, NO3, HO2, Br, Cl, O2, ClO, O2, NO2, O2, ClO, H, OH, HCHO, HCOOH, HCOOH, HCOOH, MeCO3, MeCO3, ClO, EtOO, H2O, H2O, , MeOO, CO2, O2, Cl, NO3, EtCO3, EtOO, MeCHO, H, HO2, HO2, Br, HO2, H2O, Cl, , NO2, H2O, NO2, NO3, NALD, OH, , , , CO, NO2, MeCOCH2OO, MeCOCH2OO, MeCO3, , , MeCOCH2OO, MGLY, HCHO, NO2, HCHO, MeOO, CO, NO2, O2, O2, O(3P), NO2, NO2, OH, H2O, H2O, H2O, HO2, HO2, HO2, HO2, SO2, SO2, , , MSA, SO2, H2O, SO2, SO2, , H2O, , , , HCHO, HCHO, MeCHO, , EtCHO, , Me2CO, , MeOO, EtOO, MeCO3, MGLY, HCHO, MGLY, , EtOO, , , , , , , , , O2, , , , , , , , , NO2, HO2, HO2, MeOO, HONO2, Cl, Cl, HONO2, NO2, , ,
  +
# REACT 3, , , , , , , OH, , CO, , , , CO, , OH, , , , , , , , , , CO, , OH, OH, OH, , NO2, HCHO, HCHO, CO, , OH, MGLY, HO2, HO2, NO2, , , , , , , , , , , , , , , , O2, Cl, , , , , , , , H, , O2, HO2, , , , O2, O2, , , , , , , ClO, , Cl, NO3, , , HO2, , , OH, , , , , O2, O2, , , HO2, , NO2, NO2, NO2, NO2, NO2, , , , , , , , , , , , , , , , , , , , , , , , , , NO2, NO2, NO2, NO2, NO2, , NO2, , , , , , HCHO, HO2, HACET, , , O2, , , , , CO, HONO2, , , , HO2, , , , , , , , , , , , , , O(3P), , , , , , , , , , , , , , , , MACRO2, CO, HO2, HO2, , , , , , , , , , , , , , , OH, , , , , CO, H, , , O2, , , , , , , , , , , , , , , , , , , H2O, OH, , NO2, , , , , H2O, H2O, , , OH, , NO2, NO2, NO2, NO2, SO2, , , , , OH, , , , , , , , , , , , , , , , , , , , , HO2, HCHO, , CO2, , , , , , , , , , , , , , , , , , , , , , , HONO2, H2O, , HONO2, , ,
  +
# REACT 4, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , HO2, HO2, HO2, , HO2, HCHO, , , HO2, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , HO2, , CO, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , MeCO3, H2O2, CO, CO, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , O(3P), , , , , , , , , , , , , , , , , , , , , , CO, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ,
  +
</pre>
  +
  +
===Jupyter Notebook===
  +
  +
<gallery widths=400px heights=350px>
  +
Image:T01_o3_001.png|Ozone mass-mixing ration
  +
Image:T01_f200_001.png|Flux through the OH+ClO->HO2+Cl reaction
  +
</gallery>
  +
  +
An example Jupyter Notebook which plots both tracers and chemical reaction fluxes can be found in
  +
  +
Tutorials/UMvn13.0/notebooks/BoxModel_T01-03.ipynb
  +
  +
and can be viewed online here
  +
  +
* [https://github.com/theabro/ukca/blob/master/Tutorials/UMvn13.0/notebooks/BoxModel_T01-03.ipynb BoxModel_T01-03.ipynb]
  +
  +
  +
You can try changing the chemical species and reaction flux selected.
  +
  +
Connect to the notebook using
  +
  +
* <div style="font-size: 140%;">'''http://localhost:4801/?token=UKCATraining'''</div>
   
 
==Checklist==
 
==Checklist==
Line 156: Line 226:
 
: <span style="font-size:20px">☐</span> Copy suites using the right-click menu
 
: <span style="font-size:20px">☐</span> Copy suites using the right-click menu
 
: <span style="font-size:20px">☐</span> Run suites using the play button
 
: <span style="font-size:20px">☐</span> Run suites using the play button
  +
: <span style="font-size:20px">☐</span> Plot output using gnuplot or python
  +
   
  +
[[UKCA Chemistry and Aerosol UMvn13.0 Tutorial 2|Tutorial 2]]
   
[[UKCA Chemistry and Aerosol vn11.8 Tutorial 2|Tutorial 2]]
+
[[UKCA Chemistry and Aerosol Tutorials at UMvn13.0]]
 
----
 
----
 
''Written by [[User:Nla27 | Luke Abraham]] 2022''
 
''Written by [[User:Nla27 | Luke Abraham]] 2022''

Latest revision as of 15:17, 29 January 2024

UKCA Chemistry and Aerosol Tutorials at UMvn13.0

Difficulty EASY
Time to Complete Under 1 hour
Video instructions Walkthrough (YouTube)

Remember to run mosrs-cache-password.

What you will learn in this tutorial

In this tutorial you will learn how to copy a UKCA Box Model suite, run it, and plot the output.

Copying and Running an Existing Rose Suite

Machine UM Version/Configuration Suite ID
vm UMvn13.0 UKCA Box Model u-cq774

You will need to login to the VM using SSH via a Terminal

ssh -X ukca

When you first login to your virtual machine you should then be asked for your Met Office Science Repository Service (MOSRS) password. Once you enter this, if this is the first time you are logging in, you will also be asked for your MOSRS username. This will likely be your name in lower case, and is not your username on the VM.

Then copy the rose suite by:

rosie copy u-cq774
UKCA Box Model opened with Rose and GCylc running on the VM via X2Go.

This will open an editor containing text similar to this:

description=Copy of u-cq774/trunk@236265
owner=lukeabraham
project=ukca
title=UMvn13.0 UKCA Box Model for UKCA Tutorials

# Make changes ABOVE these lines.
# The "owner", "project" and "title" fields are compulsory.
# Any KEY=VALUE pairs can be added. Known fields include:
# "access-list", "description" and "sub-project".

You can edit the title to make it more descriptive if you like. Then save (e.g. Ctrl-s) and close the file (e.g. Ctrl-q). The press y on the keyboard when prompted and press return.

The suite will now copy and checkout to your /home/$USER/roses directory, into a sub-directory that is the same as the SUITE-ID. You should change into this directory now using cd, e.g.

cd roses/u-ab123

then open the suite by typing

rose edit &

and pressing return. You can either run the suite by pressing the play button that looks like a triangle to the top right of the options bar, or by typing

rose suite-run

in the terminal. For this suite to run you need the UKCA branch

fcm:ukca.x_br/dev/lukeabraham/um13.0_ukca_box

checked-out into your home/ directory. This has been done for you.

If you are using your own VM, times will vary depending on the specifications of the host. It may take 2 or 3 minutes to compile the code (fcm_make) for the first time (and a few seconds when recompiling), followed by another few seconds to run UKCA itself (ukca). Any unsuccessful jobs will be highlighted in red and have a status of "failed".

When the suite has finished successfully it will then become blank with the message stopped with 'succeeded' in the bottom-left corner.

Version Control

Rose suites are all held under version control, using fcm. When making changes to a suite, you will need to save it before you can run the suite. Once you are happy with the settings, you can also commit these changes back to the repository - to do this change directory to the

/home/$USER/roses/[SUITE-ID]

and then type

fcm commit

a text editor will then open, and you should type a short message describing what the changes you have made do. You should then close the editor and type y in the terminal. It is recommended that you commit frequently (even on configurations that aren't working) as this protects you against mistakes and accidental deletions etc.

These suites can be viewed on the SRS here: https://code.metoffice.gov.uk/trac/roses-u (password required)

It is recommended that you commit your suites regularly.

Output Directory Structure

The output directory structure of rose suites are rather complex. The schematic gives a broad overview of the general structure. Everything can be found within the [SUITE-ID] (e.g. u-cq782 etc.) directory, which can be found within your $HOME/cylc-run directory.

Within this directory there are several directories, including:

  • log: a symbolic link to a directory called something like log.20220914T092709Z, which contains all the job.out (containing output from UM WRITE statements from the umPrint subroutine) and job.err files etc., as well as the script used to run the job.
  • work: this contains the directories used when the job actually runs. Real-time output will be sent to files here (held in a pe_output directory. Sometimes model output will also be here.


Within these directories there will be many sub-directories. Some of these will be named from the cycle-point (labelled [CYCLE] in the graphic. For the UKCA Training Suite this will be 1. The sub-directories will (eventually) be named after the app (labelled by [JOB NAME] in the table below) that the output is from, e.g. fcm_make, ukca etc.

Files General Path Example
Most recent job.out files /home/vagrant/cylc-run/[SUITE-ID]/log/job/1/[JOB NAME]/NN /home/vagrant/cylc-run/u-cq782/log/job/1/ukca/NN
Processor output (while running) /home/vagrant/cylc-run/[SUITE-ID]/work/1/[JOB NAME]/pe_output /home/vagrant/cylc-run/u-cq782/work/1/ukca/pe_output
UKCA Box Model output text files in ASCII format
tracer_out.csv
flux_out.csv
/home/vagrant/cylc-run/[SUITE-ID]/work/1/ukca /home/vagrant/cylc-run/u-cq782/work/1/ukca

Notes

  • The job.out files can also be viewed through the Gcylc GUI right-click menu from each job.
  • UKCA Box Model output file will be tracer_out.csv for tracers and flux_out.csv for chemical reaction fluxes. These can be renamed in Rose in the ukca namelist Model Input and Output panel.

Exploring Rose

Take a look at your suite, and click through the options. There are 2 key sections that you will need to consider:

fcm_make

This section controls which fcm branches are used by the model, and also what compiler settings are used etc. You should go to fcm_make env Sources to see what branches are being used. During the course of these tutorials you will make a branch and add it here to complete the tutorial tasks. You can add both relative paths to the repository (starting with branches/) and absolute paths to a working copy (starting with /home/).

ukca

Click the arrow by ukca to view the list of sections. We'll be most interested in the namelist section, which contains all the settings that configure the science that can be done by UKCA. Take a look through the panels and sub-panels in this section to get a feel for how it is organised.

Viewing Output

gnuplot

The quickest and easiest way to plot UKCA box model output is using Gnuplot, which is a relatively simple command-line plotting tool. The box model outputs tracer concentrations and chemical reaction fluxes every timestep to the files tracer_out.csv and flux_out.csv respectively.

To plot these files, change into your suite's cylc-run directory and navigate to the ukca directory in work, e.g.

cd /home/vagrant/cylc-run/u-cq782/work/1/ukca

then load gnuplot by typing

gnuplot

To plot the tracer file, type

plot "tracer_out.csv" using 1:3 with linespoints

which can be shortened to

p "tracer_out.csv" u 1:3 w lp

The file can be saved from the export menu on the plot itself. Additional options for the set command can be used to define labels, titles, change how the axis labels look, e.g.

gnuplot> set title "O3 MMR"
gnuplot> set xlabel "timestep number"
gnuplot> set ylabel "mass mixing ratio /kg kg^{-1}"
gnuplot> set format y "%.2t*10^{%+03T}"
gnuplot> unset key
gnuplot> p 'tracer_out.csv' u 1:3 w lp

The ordering of the species in the tracer_out.csv can be seen by opening it in e.g. a simple spreadsheet program like Gnumeric. You may need to install this first using

sudo apt install gnumeric

before opening the file

gnumeric tracer_out.csv

The column number is given in the first line with the tracer name in the line below. You could also open the file in emacs and then M-x toggle-truncate-lines. An example of this header for the StratTrop chemistry with GLOMAP-mode aerosols is

 # COL    1,               2,               3,               4,               5,               6,               7,               8,               9,              10,              11,              12,              13,              14,              15,              16,              17,              18,              19,              20,              21,              22,              23,              24,              25,              26,              27,              28,              29,              30,              31,              32,              33,              34,              35,              36,              37,              38,              39,              40,              41,              42,              43,              44,              45,              46,              47,              48,              49,              50,              51,              52,              53,              54,              55,              56,              57,              58,              59,              60,              61,              62,              63,              64,              65,              66,              67,              68,              69,              70,              71,              72,              73,              74,              75,              76,              77,              78,              79,              80,              81,              82,              83,              84,              85,              86,              87,              88,              89,              90,              91,              92,              93,              94,              95,              96
 # Timestep,           O(3P),              O3,               N,              NO,             NO3,             NO2,            N2O5,          HO2NO2,           HONO2,            H2O2,             CH4,              CO,            HCHO,           MeOOH,               H,             H2O,              OH,             HO2,              Cl,           Cl2O2,             ClO,            OClO,              Br,             BrO,            BrCl,          BrONO2,             N2O,             HCl,            HOCl,             HBr,            HOBr,          ClONO2,           CFCl3,          CF2Cl2,            MeBr,            HONO,            C2H6,           EtOOH,           MeCHO,             PAN,            C3H8,         n-PrOOH,         i-PrOOH,           EtCHO,           Me2CO,      MeCOCH2OOH,            PPAN,          MeONO2,            C5H8,           ISOOH,            ISON,            MACR,         MACROOH,            MPAN,           HACET,            MGLY,            NALD,           HCOOH,          MeCO3H,          MeCO2H,              H2,            MeOH,             DMS,             SO2,           H2SO4,             MSA,            DMSO,             NH3,             CS2,             COS,             H2S,             SO3,        Monoterp,         Sec_Org,       Nuc_SOL_N,      Nuc_SOL_SU,       Ait_SOL_N,      Ait_SOL_SU,      Ait_SOL_BC,      Ait_SOL_OM,       Acc_SOL_N,      Acc_SOL_SU,      Acc_SOL_BC,      Acc_SOL_OM,      Acc_SOL_SS,       Cor_SOL_N,      Cor_SOL_SU,      Cor_SOL_BC,      Cor_SOL_OM,      Cor_SOL_SS,       Ait_INS_N,      Ait_INS_BC,      Ait_INS_OM,      Nuc_SOL_OM,      PASSIVE O3

The flux_out.csv file has a more complicated header as it also gives all the chemical reactants and products, e.g.

# Timestep, UKCA reaction fluxes (units = molecules.cm^-3.s^-1)
# COL    1,               2,               3,               4,               5,               6,               7,               8,               9,              10,              11,              12,              13,              14,              15,              16,              17,              18,              19,              20,              21,              22,              23,              24,              25,              26,              27,              28,              29,              30,              31,              32,              33,              34,              35,              36,              37,              38,              39,              40,              41,              42,              43,              44,              45,              46,              47,              48,              49,              50,              51,              52,              53,              54,              55,              56,              57,              58,              59,              60,              61,              62,              63,              64,              65,              66,              67,              68,              69,              70,              71,              72,              73,              74,              75,              76,              77,              78,              79,              80,              81,              82,              83,              84,              85,              86,              87,              88,              89,              90,              91,              92,              93,              94,              95,              96,              97,              98,              99,             100,             101,             102,             103,             104,             105,             106,             107,             108,             109,             110,             111,             112,             113,             114,             115,             116,             117,             118,             119,             120,             121,             122,             123,             124,             125,             126,             127,             128,             129,             130,             131,             132,             133,             134,             135,             136,             137,             138,             139,             140,             141,             142,             143,             144,             145,             146,             147,             148,             149,             150,             151,             152,             153,             154,             155,             156,             157,             158,             159,             160,             161,             162,             163,             164,             165,             166,             167,             168,             169,             170,             171,             172,             173,             174,             175,             176,             177,             178,             179,             180,             181,             182,             183,             184,             185,             186,             187,             188,             189,             190,             191,             192,             193,             194,             195,             196,             197,             198,             199,             200,             201,             202,             203,             204,             205,             206,             207,             208,             209,             210,             211,             212,             213,             214,             215,             216,             217,             218,             219,             220,             221,             222,             223,             224,             225,             226,             227,             228,             229,             230,             231,             232,             233,             234,             235,             236,             237,             238,             239,             240,             241,             242,             243,             244,             245,             246,             247,             248,             249,             250,             251,             252,             253,             254,             255,             256,             257,             258,             259,             260,             261,             262,             263,             264,             265,             266,             267,             268,             269,             270,             271,             272,             273,             274,             275,             276,             277,             278,             279,             280,             281,             282,             283,             284,             285,             286,             287,             288,             289,             290,             291,             292,             293,             294,             295,             296,             297,             298,             299,             300,             301,             302,             303,             304,             305,             306,             307,             308,             309,             310
#   PROD 1,           Cl2O2,          HO2NO2,             PAN,            PPAN,            MPAN,            N2O5,           EtOOH,            H2O2,            HCHO,            HCHO,          HO2NO2,           HONO2,           MeCHO,           MeCHO,           MeOOH,            N2O5,             NO2,             NO3,             NO3,              O2,              O3,              O3,             PAN,            HONO,           EtCHO,           Me2CO,         n-PrOOH,         i-PrOOH,      MeCOCH2OOH,            PPAN,          MeONO2,           ISOOH,            ISON,            MACR,            MPAN,         MACROOH,         MACROOH,           HACET,            MGLY,            NALD,          MeCO3H,            BrCl,             BrO,          BrONO2,          BrONO2,              O2,            OClO,              NO,            HOBr,             N2O,             H2O,          ClONO2,          ClONO2,             HCl,            HOCl,           Cl2O2,           CFCl3,          CF2Cl2,            MeBr,             CH4,             CO2,          HO2NO2,             CS2,             COS,           H2SO4,             SO3,              Br,              Br,              Br,              Br,              Br,             BrO,             BrO,             BrO,             BrO,             BrO,             BrO,             BrO,          CF2Cl2,           CFCl3,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,              Cl,             ClO,             ClO,             ClO,             ClO,             ClO,             ClO,             ClO,           EtCO3,           EtCO3,            EtOO,            EtOO,               H,               H,               H,               H,               H,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,             HO2,            MeBr,            MeBr,            MeBr,           MeCO3,           MeCO3,       MeCOCH2OO,       MeCOCH2OO,            MeOO,            MeOO,            MeOO,               N,               N,               N,            N2O5,              NO,              NO,              NO,              NO,              NO,              NO,             NO2,             NO2,             NO3,             NO3,             NO3,             NO3,             NO3,             NO3,             NO3,           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(1D),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),           O(3P),              O3,              O3,              O3,              O3,              O3,              O3,              O3,            OClO,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,              OH,          i-PrOO,          i-PrOO,          n-PrOO,          n-PrOO,             CS2,             CS2,             DMS,             DMS,            DMSO,             H2S,             H2S,             COS,             COS,             SO2,             SO3,        Monoterp,        Monoterp,        Monoterp,            MeOO,            MeOO,            EtOO,            EtOO,          n-PrOO,          n-PrOO,          i-PrOO,          i-PrOO,           MeCO3,           EtCO3,       MeCOCH2OO,       MeCOCH2OO,            ISO2,          MACRO2,          MACRO2,             HO2,           O(3P),           O(3P),           O(3P),           O(1D),             BrO,             ClO,             ClO,               H,             HO2,             HO2,              OH,              OH,              OH,           MeCO3,           EtCO3,          MACRO2,             NO2,              NO,             SO2,             DMS,             DMS,          ClONO2,          ClONO2,            HOCl,            N2O5,            N2O5,             SO2,             SO2,             SO2
#   PROD 2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,           Cl2O2,            HCHO,             HO2,              O3,            OClO,             BrO,             ClO,             ClO,             ClO,             HO2,              NO,              OH,           O(1D),           O(1D),             CH4,           Cl2O2,          ClONO2,              H2,            H2O2,            HCHO,             HO2,             HO2,            HOCl,           MeOOH,             NO3,              O3,            OClO,             ClO,             ClO,             ClO,             HO2,            MeOO,              NO,             NO3,              NO,             NO3,              NO,             NO3,             HO2,             HO2,             HO2,             NO2,              O3,           EtCO3,           EtCO3,            EtOO,             HO2,            ISO2,          MACRO2,           MeCO3,           MeCO3,           MeCO3,       MeCOCH2OO,            MeOO,            MeOO,              NO,             NO3,              O3,          i-PrOO,          n-PrOO,              Cl,           O(1D),              OH,              NO,             NO3,              NO,             NO3,              NO,              NO,             NO3,              NO,             NO2,              O2,             H2O,            ISO2,            ISO2,          MACRO2,          MACRO2,             NO3,              O3,             NO3,              O3,              Br,            C5H8,           EtCHO,            HCHO,            MGLY,           Me2CO,           MeCHO,             CH4,             CH4,             CH4,             CO2,              H2,             H2O,             HBr,             HBr,             HCl,             HCl,             HCl,              N2,             N2O,             N2O,              O2,              O3,              O3,             BrO,             ClO,          ClONO2,            H2O2,             HBr,             HCl,             HO2,            HOCl,             NO2,             NO3,              O3,            OClO,              OH,            C5H8,            C5H8,            C5H8,            MACR,            MACR,            MACR,            MACR,              NO,            C2H6,            C3H8,            C3H8,            C5H8,             CH4,              CO,             ClO,             ClO,          ClONO2,           EtCHO,           EtOOH,           EtOOH,              H2,            H2O2,           HACET,             HBr,            HCHO,           HCOOH,             HCl,             HO2,          HO2NO2,            HOCl,            HONO,           HONO2,            ISON,           ISOOH,            MACR,            MACR,         MACROOH,            MGLY,            MPAN,           Me2CO,           Me2CO,           MeCHO,          MeCO2H,          MeCO3H,      MeCOCH2OOH,      MeCOCH2OOH,            MeOH,          MeONO2,           MeOOH,           MeOOH,            NALD,             NO3,              O3,            OClO,              OH,             PAN,            PPAN,         i-PrOOH,         i-PrOOH,         n-PrOOH,         n-PrOOH,              NO,             NO3,              NO,             NO3,           O(3P),              OH,              OH,             NO3,              OH,           O(3P),              OH,           O(3P),              OH,              O3,             H2O,              OH,              O3,             NO3,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,             RO2,           EtCO3,              O2,              NO,             NO2,              N2,             NO2,             ClO,             NO2,              O2,             HO2,             NO2,              NO,             NO2,              OH,             NO2,             NO2,             NO2,             NO3,              NO,              OH,              OH,              OH,             H2O,             HCl,             HCl,             H2O,             HCl,            H2O2,              O3,              O3
#  REACT 1,             ClO,             NO2,           MeCO3,           EtCO3,          MACRO2,             NO2,           MeCHO,              OH,             HO2,              H2,             HO2,              OH,            MeOO,             CH4,             HO2,             NO3,              NO,              NO,             NO2,           O(3P),              O2,              O2,           MeCO3,              OH,            EtOO,           MeCO3,           EtCHO,           Me2CO,           MeCO3,           EtCO3,             HO2,              OH,             NO2,           MeCO3,          MACRO2,              OH,           HACET,           MeCO3,           MeCO3,            HCHO,            MeOO,              Br,              Br,              Br,             BrO,           O(3P),           O(3P),               N,              OH,              N2,              OH,              Cl,             ClO,               H,              OH,              Cl,              Cl,              Cl,              Br,            MeOO,              CO,              OH,             COS,              CO,             SO3,             SO2,            BrCl,             HBr,             HBr,             BrO,             BrO,              Br,              Br,              Br,            BrCl,            HOBr,              Br,              Br,              Cl,              Cl,             HCl,              Cl,              Cl,             HCl,             HCl,             HCl,             ClO,             HCl,              Cl,             HCl,             ClO,             ClO,             ClO,              Cl,              Cl,              Cl,            HOCl,              Cl,              Cl,              Cl,            EtOO,            EtOO,           MeCHO,           MeCHO,              H2,           O(3P),              OH,              OH,              OH,              O2,              O3,           EtOOH,            H2O2,           ISOOH,         MACROOH,          MeCO2H,          MeCO3H,              OH,      MeCOCH2OOH,            HCHO,           MeOOH,              OH,              OH,              OH,         i-PrOOH,         n-PrOOH,              Br,              Br,              Br,            MeOO,            MeOO,           MeCO3,           MeCO3,             HO2,          MeONO2,             HO2,              N2,             N2O,              NO,           HONO2,            ISON,             NO2,            MGLY,             NO2,             NO2,             NO2,              NO,             NO3,             BrO,            ISON,           HONO2,           HONO2,           MeCO3,           HONO2,           HONO2,            HCHO,            HCHO,              OH,           O(3P),              OH,              OH,             HBr,              OH,               H,           O(3P),              OH,           O(3P),              N2,              NO,           O(3P),              O2,              O2,              O2,              Cl,             ClO,              OH,              OH,              OH,              OH,              OH,              NO,              O2,              O2,              O2,              O2,             HO2,            MACR,            MeOO,            MGLY,            MGLY,              OH,              OH,             NO2,             H2O,          i-PrOO,          n-PrOO,            ISO2,             H2O,               H,             HCl,             HO2,            HOCl,             H2O,             H2O,             H2O,             H2O,             H2O,            MGLY,             H2O,             H2O,             CO2,             H2O,             H2O,             H2O,             ClO,             H2O,             H2O,           HACET,            MACR,          MACRO2,          MACRO2,          MACRO2,           MeCO3,           HACET,             H2O,             H2O,             H2O,            MeOO,           MeCO3,             H2O,              OH,             HO2,            HCHO,             H2O,             H2O,            HCHO,             HO2,             HO2,            HOCl,             H2O,            HCHO,           MeCHO,           Me2CO,          i-PrOO,           EtCHO,          n-PrOO,           Me2CO,           Me2CO,           EtCHO,           EtCHO,              CO,             COS,             SO2,             SO2,             SO2,              OH,             SO2,              CO,             CO2,             SO3,           H2SO4,         Sec_Org,         Sec_Org,         Sec_Org,             HO2,            MeOH,             HO2,           MeCHO,             HO2,           EtCHO,             HO2,           Me2CO,          MeCO2H,                ,            HCHO,           HACET,            MACR,           HACET,             HO2,              OH,              O3,             NO2,             NO3,             N2O,          BrONO2,           Cl2O2,          ClONO2,             HO2,            H2O2,          HO2NO2,            HONO,           HONO2,            H2O2,             PAN,            PPAN,            MPAN,            N2O5,             NO2,             SO3,            DMSO,             SO2,            HOCl,              Cl,              Cl,           HONO2,              Cl,                ,                ,
#  REACT 2,             ClO,             HO2,             NO2,             NO2,             NO2,             NO3,             HO2,              OH,             HO2,              CO,             NO2,             NO2,             HO2,              CO,            HCHO,             NO2,           O(3P),              O2,           O(3P),           O(3P),           O(1D),           O(3P),             NO2,              NO,             HO2,            MeOO,             HO2,             HO2,            HCHO,             NO2,            HCHO,            MACR,            MACR,            HCHO,             NO2,             HO2,              CO,            HCHO,              CO,              CO,              OH,              Cl,           O(3P),             NO3,             NO2,           O(1D),             ClO,           O(3P),              Br,           O(1D),               H,             NO3,             NO2,              Cl,              Cl,              Cl,              Cl,              Cl,               H,               H,           O(3P),             NO3,             SO2,             SO2,              OH,           O(3P),              Cl,              CO,              O2,              O2,             ClO,              Br,              Cl,            OClO,              O2,              O2,             NO2,             HO2,             ClO,              Cl,            MeOO,              Cl,              Cl,               H,             HO2,              CO,              OH,              O2,              Cl,            MeOO,             NO2,              O2,             ClO,              Cl,              Cl,            OClO,              O2,            HCHO,             NO2,              O2,             CO2,             CO2,             HO2,             HO2,              O2,             H2O,              OH,              NO,              O2,                ,                ,                ,                ,                ,                ,              O3,                ,            MeOO,                ,                ,                ,             NO2,             NO2,              O2,                ,                ,             HCl,              OH,             H2O,             CO2,             CO2,            HCHO,            HCHO,            HCHO,                ,            HCHO,           O(3P),           O(3P),           O(3P),           HONO2,                ,            MACR,            HCHO,           MeCO3,             NO2,                ,             NO2,                ,             NO2,                ,           EtCO3,             HO2,              CO,       MeCOCH2OO,           MeCO3,              H2,             HO2,            MeOO,             CO2,               H,              OH,           O(3P),              Br,             ClO,             HCl,              Cl,              N2,              O2,              NO,              O2,           O(3P),              O2,              Br,              O2,             NO3,             HO2,              Br,              Cl,              O2,             ClO,              O2,             NO2,              O2,             ClO,               H,              OH,            HCHO,           HCOOH,           HCOOH,           HCOOH,           MeCO3,           MeCO3,             ClO,            EtOO,             H2O,             H2O,                ,            MeOO,             CO2,              O2,              Cl,             NO3,           EtCO3,            EtOO,           MeCHO,               H,             HO2,             HO2,              Br,             HO2,             H2O,              Cl,                ,             NO2,             H2O,             NO2,             NO3,            NALD,              OH,                ,                ,                ,              CO,             NO2,       MeCOCH2OO,       MeCOCH2OO,           MeCO3,                ,                ,       MeCOCH2OO,            MGLY,            HCHO,             NO2,            HCHO,            MeOO,              CO,             NO2,              O2,              O2,           O(3P),             NO2,             NO2,              OH,             H2O,             H2O,             H2O,             HO2,             HO2,             HO2,             HO2,             SO2,             SO2,                ,                ,             MSA,             SO2,             H2O,             SO2,             SO2,                ,             H2O,                ,                ,                ,            HCHO,            HCHO,           MeCHO,                ,           EtCHO,                ,           Me2CO,                ,            MeOO,            EtOO,           MeCO3,            MGLY,            HCHO,            MGLY,                ,            EtOO,                ,                ,                ,                ,                ,                ,                ,                ,              O2,                ,                ,                ,                ,                ,                ,                ,                ,             NO2,             HO2,             HO2,            MeOO,           HONO2,              Cl,              Cl,           HONO2,             NO2,                ,                ,
#  REACT 3,                ,                ,                ,                ,                ,                ,              OH,                ,              CO,                ,                ,                ,              CO,                ,              OH,                ,                ,                ,                ,                ,                ,                ,                ,                ,              CO,                ,              OH,              OH,              OH,                ,             NO2,            HCHO,            HCHO,              CO,                ,              OH,            MGLY,             HO2,             HO2,             NO2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,              O2,              Cl,                ,                ,                ,                ,                ,                ,                ,               H,                ,              O2,             HO2,                ,                ,                ,              O2,              O2,                ,                ,                ,                ,                ,                ,             ClO,                ,              Cl,             NO3,                ,                ,             HO2,                ,                ,              OH,                ,                ,                ,                ,              O2,              O2,                ,                ,             HO2,                ,             NO2,             NO2,             NO2,             NO2,             NO2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,             NO2,             NO2,             NO2,             NO2,             NO2,                ,             NO2,                ,                ,                ,                ,                ,            HCHO,             HO2,           HACET,                ,                ,              O2,                ,                ,                ,                ,              CO,           HONO2,                ,                ,                ,             HO2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,           O(3P),                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,          MACRO2,              CO,             HO2,             HO2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,              OH,                ,                ,                ,                ,              CO,               H,                ,                ,              O2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,             H2O,              OH,                ,             NO2,                ,                ,                ,                ,             H2O,             H2O,                ,                ,              OH,                ,             NO2,             NO2,             NO2,             NO2,             SO2,                ,                ,                ,                ,              OH,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,             HO2,            HCHO,                ,             CO2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,           HONO2,             H2O,                ,           HONO2,                ,                ,
#  REACT 4,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,             HO2,             HO2,             HO2,                ,             HO2,            HCHO,                ,                ,             HO2,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,             HO2,                ,              CO,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,           MeCO3,            H2O2,              CO,              CO,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,           O(3P),                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,              CO,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,                ,

Jupyter Notebook

An example Jupyter Notebook which plots both tracers and chemical reaction fluxes can be found in

Tutorials/UMvn13.0/notebooks/BoxModel_T01-03.ipynb

and can be viewed online here


You can try changing the chemical species and reaction flux selected.

Connect to the notebook using

Checklist

List suites using rosie go
Copy suites using the right-click menu
Run suites using the play button
Plot output using gnuplot or python


Tutorial 2

UKCA Chemistry and Aerosol Tutorials at UMvn13.0


Written by Luke Abraham 2022