UKCA Chemistry and Aerosol vn10.9 Tutorial 10
UKCA Chemistry and Aerosol Tutorials at vn10.9
What you will learn in this Tutorial
Task 10.1: Output aerosol diagnostics
TASK 10.1: Output the following aerosol and radiation diagnostics to the UPA output stream. You should make a new time profile (called TRAD) to only output these on radiation timesteps. You will also need to make a new domain profile (called D3DAR) for s02i530 and s02i540.
Hint |
---|
You can find the number of radiation timesteps in a 24-hour period at: um namelist UM Science Settings Section 01 - 02 Radiation Shortwave. |
STASH Section | STASH Item | STASH Name |
---|---|---|
1 | 207 | INCOMING SW RAD FLUX (TOA): ALL TSS |
1 | 208 | OUTGOING SW RAD FLUX (TOA) |
2 | 205 | OUTGOING LW RAD FLUX (TOA) |
2 | 285 | MINERAL DUST OPTICAL DEPTH IN RADN. |
2 | 300 | AITKEN MODE (SOLUBLE) OPTICAL DEPTH |
2 | 301 | ACCUM MODE (SOLUBLE) OPTICAL DEPTH |
2 | 302 | COARSE MODE (SOLUBLE) OPTICAL DEPTH |
2 | 303 | AITKEN MODE (INSOL) OPTICAL DEPTH |
2 | 304 | ACCUM MODE (INSOL) OPTICAL DEPTH |
2 | 305 | COARSE MODE (INSOL) OPTICAL DEPTH |
2 | 585 | MINERAL DUST ABS. OPICAL DEPTH |
2 | 240 | AITKEN (SOLUBLE) ABS OPTICAL DEPTH |
2 | 241 | ACCUM (SOLUBLE) ABS OPTICAL DEPTH |
2 | 242 | COARSE (SOLUBLE) ABS OPTICAL DEPTH |
2 | 243 | AITKEN (INSOL) ABS OPTICAL DEPTH |
2 | 244 | ACCUM (INSOL) ABS OPTICAL DEPTH |
2 | 245 | COARSE (INSOL) ABS OPTICAL DEPTH |
2 | 530 | UKCA 3D AEROSOL EXTINCTION |
2 | 540 | CLASSIC 3D AEROSOL EXTINCTION |
Solution to Task 10.1
You were given the task
- Output the following aerosol and radiation diagnostics to the UPA output stream. You should make a new time profile (called TRAD) to only output these on radiation timesteps. You will also need to make a new domain profile (called D3DAR) for s02i530 and s02i540.
STASH Section | STASH Item | STASH Name |
---|---|---|
1 | 207 | INCOMING SW RAD FLUX (TOA): ALL TSS |
1 | 208 | OUTGOING SW RAD FLUX (TOA) |
2 | 205 | OUTGOING LW RAD FLUX (TOA) |
2 | 285 | MINERAL DUST OPTICAL DEPTH IN RADN. |
2 | 300 | AITKEN MODE (SOLUBLE) OPTICAL DEPTH |
2 | 301 | ACCUM MODE (SOLUBLE) OPTICAL DEPTH |
2 | 302 | COARSE MODE (SOLUBLE) OPTICAL DEPTH |
2 | 303 | AITKEN MODE (INSOL) OPTICAL DEPTH |
2 | 304 | ACCUM MODE (INSOL) OPTICAL DEPTH |
2 | 305 | COARSE MODE (INSOL) OPTICAL DEPTH |
2 | 585 | MINERAL DUST ABS. OPICAL DEPTH |
2 | 240 | AITKEN (SOLUBLE) ABS OPTICAL DEPTH |
2 | 241 | ACCUM (SOLUBLE) ABS OPTICAL DEPTH |
2 | 242 | COARSE (SOLUBLE) ABS OPTICAL DEPTH |
2 | 243 | AITKEN (INSOL) ABS OPTICAL DEPTH |
2 | 244 | ACCUM (INSOL) ABS OPTICAL DEPTH |
2 | 245 | COARSE (INSOL) ABS OPTICAL DEPTH |
2 | 530 | UKCA 3D AEROSOL EXTINCTION |
2 | 540 | CLASSIC 3D AEROSOL EXTINCTION |
and were given the hint
- You can find the number of radiation timesteps in a 24-hour period at: um namelist UM Science Settings Section 01 - 02 Radiation Shortwave.
For a working Rose suite that has completed this task, please see
- ARCHER:
u-as292@62651
- vm:
u-as297@62631
The specific Rose changes made are:
- ARCHER: https://code.metoffice.gov.uk/trac/roses-u/changeset/62651/a/s/2/9/2/trunk
- vm: https://code.metoffice.gov.uk/trac/roses-u/changeset/62631/a/s/2/9/7/trunk
The specific Rose changes made are:
ARCHER:
Index: app/um/rose-app.conf =================================================================== --- app/um/rose-app.conf (revision 60289) +++ app/um/rose-app.conf (revision 62651) @@ -2744,6 +2744,40 @@ precip_segment_size=32 ukca_mode_seg_size=4 +[namelist:umstash_domain(d3dar_72578706)] +dom_name='D3DAR' +!!iest=0 +ilevs=1 +imn=0 +imsk=1 +!!inth=0 +iopa=1 +iopl=2 +!!isth=0 +!!iwst=0 +iwt=0 +!!l_spml_ts=.false. +levb=01 +!!levlst=0 +levt=85 +plt=4 +pslist=1,2,3,4,5,6 +!!rlevlst=0 +!!spml_bot=0 +!!spml_ew=0 +!!spml_ns=0 +!!spml_top=0 +!!tblim=0 +!!tblimr=0 +!!telim=0 +!!tnlim=0 +ts=.false. +!!tslim=0 +!!tsnum=0 +!!ttlim=0 +!!ttlimr=0 +!!twlim=0 + [namelist:umstash_domain(dallrh_0496a967)] dom_name='DALLRH' !!iest=0 @@ -3025,6 +3059,22 @@ tim_name='TALLTS' use_name='UPUKCA' +[namelist:umstash_streq(01207_ed72c304)] +dom_name='DIAG' +isec=1 +item=207 +package= +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(01208_83140cd8)] +dom_name='DIAG' +isec=1 +item=208 +package= +tim_name='TRAD' +use_name='UPA' + [namelist:umstash_streq(01235_3511dd9f)] dom_name='DIAG' isec=1 @@ -3033,14 +3083,142 @@ tim_name='TALLTS' use_name='UPUKCA' -[namelist:umstash_streq(02301_0f7c5f4a)] +[namelist:umstash_streq(02205_357bf644)] +dom_name='DIAG' +isec=2 +item=205 +package= +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02240_d97aaab7)] dom_name='DIAGAOT' isec=2 +item=240 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02241_8cda3169)] +dom_name='DIAGAOT' +isec=2 +item=241 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02242_91e371db)] +dom_name='DIAGAOT' +isec=2 +item=242 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02243_21bc5f11)] +dom_name='DIAGAOT' +isec=2 +item=243 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02244_32fe0790)] +dom_name='DIAGAOT' +isec=2 +item=244 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02245_da00b6ef)] +dom_name='DIAGAOT' +isec=2 +item=245 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02285_1d9800f0)] +dom_name='DIAGAOT' +isec=2 +item=285 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02300_8b9907b5)] +dom_name='DIAGAOT' +isec=2 +item=300 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02301_b8711d23)] +dom_name='DIAGAOT' +isec=2 item=301 -package='UKCA Testing' -tim_name='T3HMN' +package=' ' +tim_name='TRAD' use_name='UPA' +[namelist:umstash_streq(02302_420d0ec7)] +dom_name='DIAGAOT' +isec=2 +item=302 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02303_ad5c3af4)] +dom_name='DIAGAOT' +isec=2 +item=303 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02304_8c1869b6)] +dom_name='DIAGAOT' +isec=2 +item=304 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02305_9ecd020a)] +dom_name='DIAGAOT' +isec=2 +item=305 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02530_7a218781)] +dom_name='D3DAR' +isec=2 +item=530 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02540_56484a28)] +dom_name='D3DAR' +isec=2 +item=540 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02585_38287367)] +dom_name='DIAGAOT' +isec=2 +item=585 +package=' ' +tim_name='TRAD' +use_name='UPA' + [namelist:umstash_streq(03025_c8768f77)] dom_name='DIAG' isec=3 @@ -3551,6 +3729,25 @@ !!unt2=2 unt3=1 +[namelist:umstash_time(trad_4c3a45a7)] +!!iedt=0 +iend=-1 +ifre=3 +!!intv=0 +!!ioff=0 +iopt=1 +!!isam=0 +!!isdt=0 +!!iser=0 +istr=1 +!!itimes=0 +ityp=1 +!!lts0=.false. +tim_name='TRAD' +!!unt1=2 +!!unt2=2 +unt3=1 + [namelist:umstash_use(upa_ffb3f00b)] file_id='pp0' locn=3
These differences can be found in the file /home/ukca/Tutorial/vn10.9/worked_solutions/Task10.1/Task10.1_rose.patch
on PUMA.
vm:
Index: app/um/rose-app.conf =================================================================== --- app/um/rose-app.conf (revision 60286) +++ app/um/rose-app.conf (revision 62631) @@ -2849,6 +2849,40 @@ !!ttlimr=0.0 !!twlim=0 +[namelist:umstash_domain(d3dar_72578706)] +dom_name='D3DAR' +!!iest=0 +ilevs=1 +imn=0 +imsk=1 +!!inth=0 +iopa=1 +iopl=2 +!!isth=0 +!!iwst=0 +iwt=0 +!!l_spml_ts=.false. +levb=01 +!!levlst=0 +levt=85 +plt=4 +pslist=1,2,3,4,5,6 +!!rlevlst=0 +!!spml_bot=0 +!!spml_ew=0 +!!spml_ns=0 +!!spml_top=0 +!!tblim=0 +!!tblimr=0 +!!telim=0 +!!tnlim=0 +ts=.false. +!!tslim=0 +!!tsnum=0 +!!ttlim=0 +!!ttlimr=0 +!!twlim=0 + [namelist:umstash_domain(dallrh_0496a967)] dom_name='DALLRH' !!iest=0 @@ -3617,6 +3651,22 @@ tim_name='TALLTS' use_name='UPUKCA' +[namelist:umstash_streq(01207_ed72c304)] +dom_name='DIAG' +isec=1 +item=207 +package= +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(01208_83140cd8)] +dom_name='DIAG' +isec=1 +item=208 +package= +tim_name='TRAD' +use_name='UPA' + [namelist:umstash_streq(01235_3511dd9f)] dom_name='DIAG' isec=1 @@ -3625,14 +3675,142 @@ tim_name='TALLTS' use_name='UPUKCA' -[namelist:umstash_streq(02301_0f7c5f4a)] +[namelist:umstash_streq(02205_357bf644)] +dom_name='DIAG' +isec=2 +item=205 +package= +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02240_d97aaab7)] dom_name='DIAGAOT' isec=2 +item=240 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02241_8cda3169)] +dom_name='DIAGAOT' +isec=2 +item=241 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02242_91e371db)] +dom_name='DIAGAOT' +isec=2 +item=242 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02243_21bc5f11)] +dom_name='DIAGAOT' +isec=2 +item=243 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02244_32fe0790)] +dom_name='DIAGAOT' +isec=2 +item=244 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02245_da00b6ef)] +dom_name='DIAGAOT' +isec=2 +item=245 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02285_1d9800f0)] +dom_name='DIAGAOT' +isec=2 +item=285 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02300_8b9907b5)] +dom_name='DIAGAOT' +isec=2 +item=300 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02301_b8711d23)] +dom_name='DIAGAOT' +isec=2 item=301 -package='UKCA Testing' -tim_name='T3HMN' +package=' ' +tim_name='TRAD' use_name='UPA' +[namelist:umstash_streq(02302_420d0ec7)] +dom_name='DIAGAOT' +isec=2 +item=302 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02303_ad5c3af4)] +dom_name='DIAGAOT' +isec=2 +item=303 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02304_8c1869b6)] +dom_name='DIAGAOT' +isec=2 +item=304 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02305_9ecd020a)] +dom_name='DIAGAOT' +isec=2 +item=305 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02530_7a218781)] +dom_name='D3DAR' +isec=2 +item=530 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02540_56484a28)] +dom_name='D3DAR' +isec=2 +item=540 +package=' ' +tim_name='TRAD' +use_name='UPA' + +[namelist:umstash_streq(02585_38287367)] +dom_name='DIAGAOT' +isec=2 +item=585 +package=' ' +tim_name='TRAD' +use_name='UPA' + [namelist:umstash_streq(03025_c8768f77)] dom_name='DIAG' isec=3 @@ -4504,6 +4682,25 @@ unt2=2 unt3=4 +[namelist:umstash_time(trad_4c3a45a7)] +!!iedt=0 +iend=-1 +ifre=3 +!!intv=0 +!!ioff=0 +iopt=1 +!!isam=0 +!!isdt=0 +!!iser=0 +istr=1 +!!itimes=0 +ityp=1 +!!lts0=.false. +tim_name='TRAD' +!!unt1=2 +!!unt2=2 +unt3=1 + [namelist:umstash_time(traddm_fa7c24ce)] !!iedt=0 iend=-1
If you open the .pa file in Xconv, you should see the following additional fields:
0 : 96 72 1 2 field200: INCOMING SW RAD FLUX (TOA): ALL TSS 1 : 96 72 1 2 field201: OUTGOING SW RAD FLUX (TOA) 2 : 96 72 1 2 olr: OUTGOING LW RAD FLUX (TOA) 3 : 96 72 6 2 unspecified: Stash code = 2240 4 : 96 72 6 2 unspecified: Stash code = 2241 5 : 96 72 6 2 unspecified: Stash code = 2242 6 : 96 72 6 2 unspecified: Stash code = 2243 7 : 96 72 6 2 unspecified: Stash code = 2244 8 : 96 72 6 2 unspecified: Stash code = 2245 9 : 96 72 6 2 unspecified: MINERAL DUST OPTICAL DEPTH IN RADN. 10 : 96 72 6 2 unspecified: AITKEN MODE (SOLUBLE) OPTICAL DEPTH 11 : 96 72 6 2 unspecified: ACCUM MODE (SOLUBLE) OPTICAL DEPTH 12 : 96 72 6 2 unspecified: COARSE MODE (SOLUBLE) OPTICAL DEPTH 13 : 96 72 6 2 unspecified: AITKEN MODE (INSOL) OPTICAL DEPTH 14 : 96 72 6 2 unspecified: ACCUM MODE (INSOL) OPTICAL DEPTH 15 : 96 72 6 2 unspecified: COARSE MODE (INSOL) OPTICAL DEPTH 16 : 96 72 38 2 unspecified: Stash code = 2530 17 : 96 72 38 2 unspecified: Stash code = 2530 18 : 96 72 38 2 unspecified: Stash code = 2530 19 : 96 72 38 2 unspecified: Stash code = 2530 20 : 96 72 38 2 unspecified: Stash code = 2530 21 : 96 72 38 2 unspecified: Stash code = 2530 22 : 96 72 38 2 unspecified: Stash code = 2540 23 : 96 72 38 2 unspecified: Stash code = 2540 24 : 96 72 38 2 unspecified: Stash code = 2540 25 : 96 72 38 2 unspecified: Stash code = 2540 26 : 96 72 38 2 unspecified: Stash code = 2540 27 : 96 72 38 2 unspecified: Stash code = 2540 28 : 96 72 6 2 unspecified: Stash code = 2585
Sample output from this task can be found at /work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.1/atmosa.pa19810901_00
on ARCHER.
Task 10.2: Calculate aerosol optical depth
TASK 10.2:
Solution to Task 10.2
Sample output from this task can be found in the /work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.2/
directory on ARCHER, containing the following:
Task102_AOD.nc write_AOD.py
Task 10.3: Calculate the single-scattering albedo
TASK 10.3:
Solution to Task 10.3
Sample output from this task can be found in the /work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.3/
directory on ARCHER, containing the following:
Task103_SSA.nc write_SSA.py
Task 10.4: Calculate the top of the atmosphere net downward radiative flux
TASK 10.4:
Solution to Task 10.4
Sample output from this task can be found in the /work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.4/
directory on ARCHER, containing the following:
Task104_TOA.nc write_TOA.py
Task 10.5: Calculate aerosol optical depth from the 3D aerosol extinction
TASK 10.5:
Solution to Task 10.5
Sample output from this task can be found in the /work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.5/
directory on ARCHER, containing the following:
Task105_AOD.nc Task105_AOD_diff.nc calc_AOD.py diff_AOD_methods.py
Task 10.6: Calculate the difference in aerosol impacts when Sec_Org
is no longer formed from ALICE
TASK 10.6:
Solution to Task 10.5
For a working Rose suite that has completed this task, please see
- ARCHER:
u-as292@62669
- vm:
u-as297@62632
The specific Rose changes made are:
- ARCHER: https://code.metoffice.gov.uk/trac/roses-u/changeset/62669/a/s/2/9/2/trunk
- vm: https://code.metoffice.gov.uk/trac/roses-u/changeset/62632/a/s/2/9/7/trunk
The specific Rose changes made are:
ARCHER:
Index: app/fcm_make/rose-app.conf =================================================================== --- app/fcm_make/rose-app.conf (revision 62651) +++ app/fcm_make/rose-app.conf (revision 62669) @@ -42,4 +42,4 @@ stash_version=1A timer_version=3A um_rev=vn10.9 -um_sources=branches/dev/lukeabraham/vn10.9_UKCA_Tutorial_Solns@46718 +um_sources=branches/dev/lukeabraham/vn10.9_UKCA_Tutorial_Solns@47380
These differences can be found in the file /home/ukca/Tutorial/vn10.9/worked_solutions/Task10.6/Task10.6_rose.patch
on PUMA.
vm:
Index: app/fcm_make/rose-app.conf =================================================================== --- app/fcm_make/rose-app.conf (revision 62631) +++ app/fcm_make/rose-app.conf (revision 62632) @@ -42,4 +42,4 @@ stash_version=1A timer_version=3A um_rev=vn10.9 -um_sources=branches/dev/lukeabraham/vn10.9_UKCA_Tutorial_Solns@46718 +um_sources=branches/dev/lukeabraham/vn10.9_UKCA_Tutorial_Solns@47380
The specific UM changes made are:
Index: src/atmosphere/UKCA/asad_flux_dat.F90 =================================================================== --- src/atmosphere/UKCA/asad_flux_dat.F90 (revision 46718) +++ src/atmosphere/UKCA/asad_flux_dat.F90 (revision 47380) @@ -1287,9 +1287,9 @@ TYPE(asad_flux_defn), PARAMETER, PUBLIC :: & ukca_tutorial_fluxes(3) = (/ & -asad_flux_defn('RXN',50134,'B',.FALSE.,0,4, & +asad_flux_defn('RXN',50134,'B',.FALSE.,0,3, & (/'ALICE ','OH '/), & -(/'BOB ','Sec_Org ',' ',' '/)), & +(/'BOB ',' ',' ',' '/)), & asad_flux_defn('DEP',50135,'D',.FALSE.,0,1, & (/'ALICE ',' '/), & (/' ',' ',' ',' '/)), & Index: src/atmosphere/UKCA/ukca_chem_master.F90 =================================================================== --- src/atmosphere/UKCA/ukca_chem_master.F90 (revision 46718) +++ src/atmosphere/UKCA/ukca_chem_master.F90 (revision 47380) @@ -2156,7 +2156,7 @@ 'HCHO ',1.00e-12, 0.00, 0.00, 1.00, 0.75, 0.25, 2.75, TI,0,0,107),& ratb_t1(277,'MACRO2 ','MeOO ','HO2 ','CO ',' ',& ' ',1.00e-12, 0.00, 0.00, 1.17, 0.25, 0.00, 0.00, TI,0,0,107),& -ratb_t1(278,'ALICE ','OH ','BOB ','Sec_Org ',' ',& +ratb_t1(278,'ALICE ','OH ','BOB ',' ',' ',& ' ',2.70E-11, 0.00, -390.00, 0.00, 0.00, 0.00, 0.00, ST,0,0,107) /) !----------------------------------------------------------------------
These differences can be found in the file /home/ukca/Tutorial/vn10.9/worked_solutions/Task10.6/Task10.6_code.patch
on PUMA.
Sample output from this task can be found in the /work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.6/
directory on ARCHER, containing the following:
Task106_AOD.nc Task106_AOD_diff.nc Task106_SSA.nc Task106_SSA_diff.nc Task106_TOA.nc Task106_TOA_diff.nc atmosa.pa19810901_00 diff_AOD_rxn.py diff_SSA_rxn.py diff_TOA_rxn.py
Checklist
- ☐
Written by Luke Abraham, Nicolas Bellouin, & Anja Schmidt 2017