Solution to UKCA & UMUI Tutorial 8 Task 8.1
You were asked to
Add in wet deposition for BOB, using the following values:
for the 1st dissociation for the 1st dissociation for the 2nd dissociation for the 2nd dissociation
Changes to ukca_chem_strattrop.F90
First, you will need to increase the value of nwet_strattrop from 29 to 30.
Then, you will need to change the 0 in the 7th column of the BOB entry in the chch_defs_strattrop_chem array to 1, e.g.
chch_t( 77,'BOB ', 1,'TR ',' ', 0, 1, 0) & ! 77
Finally, you will need to add the following to the end of the henry_defs_strattrop_chem array
0.2100E+06, 0.8700E+04, 0.2000E+02, 0.0000E+00, 0.0000E+00, 0.0000E+00 & ! 30 BOB
(remembering to add a comma to the end of the preceding line).
You will need to make a hand-edit containing
ed CNTLATM<<\EOF /JPDW/ d i JPDW = 30, . w q EOF
and add this to the UMUI in the Model Selection → Input/Output Control and Resource → User hand edit files by placing it in the table and putting a Y in the Include Y/N column.
An example hand-edit can be found at
If you view the BOB (34065) field in the pb file (located in your archive directory) after these changes, and compare it to the equivalent field after the previous Task (7.1) you will see that there are some differences. However, we will not be able to quantify these differences correctly until Tutorial 9: Adding new UKCA diagnostics.
Sample output from this job can be found in
on HECToR, and in
There is a worked solution to this problem in the UMUI Tutorial experiment. This is job i: Tutorial: solution to Task 8.1 - add new wet deposition .
The code changes can be viewed by using the following FCM command
fcm diff -g fcm:um_br/dev/luke/vn8.2_UKCA_Tutorial_Solns@12214 fcm:um_br/dev/luke/vn8.2_UKCA_Tutorial_Solns@12227
This gives the following (non-graphical) output:
Index: src/atmosphere/UKCA/ukca_chem_strattrop.F90 =================================================================== --- src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 12214) +++ src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 12227) @@ -60,7 +60,7 @@ INTEGER, PARAMETER, PUBLIC :: ndry_st_aer = 5 ! Aerosol chemistry ! No of wet deposited species -INTEGER, PARAMETER, PUBLIC :: nwet_strattrop = 29 ! Stratospheric chemistry +INTEGER, PARAMETER, PUBLIC :: nwet_strattrop = 30 ! Stratospheric chemistry INTEGER, PARAMETER, PUBLIC :: nwet_st_aer = 34 ! Aerosol chemistry @@ -142,7 +142,7 @@ chch_t( 74,'O2 ', 1,'CT ',' ', 0, 0, 0), & ! 74 chch_t( 75,'N2 ', 1,'CT ',' ', 0, 0, 0), & ! 75 chch_t( 76,'ALICE ', 1,'TR ',' ', 1, 0, 0), & ! 76 -chch_t( 77,'BOB ', 1,'TR ',' ', 0, 0, 0) & ! 77 +chch_t( 77,'BOB ', 1,'TR ',' ', 0, 1, 0) & ! 77 /) TYPE(CHCH_T), DIMENSION( 87), PUBLIC :: chch_defs_strattrop_aer=(/ & @@ -689,7 +689,8 @@ 0.6900E+04, 0.5600E+04, 0.1800E-03,-0.1510E+04, 0.0000E+00, 0.0000E+00,& ! 26 HCOOH 0.7500E+03, 0.5300E+04, 0.6300E-08, 0.0000E+00, 0.0000E+00, 0.0000E+00,& ! 27 MeCO3H 0.4700E+04, 0.6000E+04, 0.1800E-04, 0.0000E+00, 0.0000E+00, 0.0000E+00,& ! 28 MeCO2H - 0.2300E+03, 0.4900E+04, 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00 & ! 29 MeOH + 0.2300E+03, 0.4900E+04, 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,& ! 29 MeOH + 0.2100E+06, 0.8700E+04, 0.2000E+02, 0.0000E+00, 0.0000E+00, 0.0000E+00 & ! 30 BOB /),(/ 6, nwet_strattrop/) )
Written by Luke Abraham 2013