Difference between revisions of "UKCA Chemistry and Aerosol vn10.9 Tutorial 10"
Line 59: | Line 59: | ||
|} |
|} |
||
− | ==Solution to Task 10.1 |
+ | ==Solution to Task 10.1== |
− | |||
You were given the task |
You were given the task |
||
Line 644: | Line 643: | ||
Sample output from this task can be found at '''<code>/work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.1/atmosa.pa19810901_00</code>''' on ARCHER. |
Sample output from this task can be found at '''<code>/work/n02/n02/ukca/Tutorial/vn10.9/sample_output/Task10.1/atmosa.pa19810901_00</code>''' on ARCHER. |
||
+ | |||
+ | ==Task 10.2: Calculate Aerosol Optical Depth== |
||
+ | |||
+ | <span style="color:green">'''TASK 10.1:''' </span> |
||
+ | |||
+ | ==Solution to Task 10.2 == |
||
==Checklist== |
==Checklist== |
Revision as of 15:48, 12 December 2017
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.1:
Solution to Task 10.2
Checklist
- ☐
Written by Luke Abraham, Nicolas Bellouin, & Anja Schmidt 2017