Solution to UKCA Chemistry and Aerosol Tutorial 8 Task 8.1
Back to UKCA Chemistry and Aerosol Tutorials
Back to the adding wet deposition tutorial
Task
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
Solution
Changes to ukca_chem_strattrop.F90
First, you will need to increase the value of nwet_st_aer from 34 to 35.
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.
! 77 chch_t( 77,'BOB ', 1,'TR ',' ', 0, 1, 0), &
Finally, you will need to add the following to the end of the henry_defs_strattrop_aer array
! 31 BOB 0.2100E+06, 0.8700E+04, 0.2000E+02, 0.0000E+00, 0.0000E+00, 0.0000E+00,&
(there is a separate and complete array for when using just the CheST/StratTrop chemistry, henry_defs_strattrop_chem)
Hand-edit
You will need to make a hand-edit containing
ed CNTLATM<<\EOF /JPDW/ d i JPDW = 35, . 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
/home/ukca/hand_edits/VN8.4/Tutorial/Task8.1_incr_JPvals.ed
Output
If you view the BOB (34065) field in the pa 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
/work/n02/n02/ukca/Tutorial/vn8.4/sample_output/Task8.1/
on ARCHER, and in
/projects/ukca/Tutorial/vn8.4/sample_ouput/Task8.1/
on MONSooN.
Worked Solution
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 fcm:um_br/dev/luke/vn8.4_UKCA_Tutorial_Solns@14711 fcm:um_br/dev/luke/vn8.4_UKCA_Tutorial_Solns@14713
This gives the following output:
Index: src/atmosphere/UKCA/ukca_chem_strattrop.F90 =================================================================== --- src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 14711) +++ src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 14713) @@ -61,7 +61,7 @@ ! No of wet deposited species INTEGER, PARAMETER, PUBLIC :: nwet_strattrop = 29 ! Stratos chemistry -INTEGER, PARAMETER, PUBLIC :: nwet_st_aer = 34 ! Aerosol chemistry +INTEGER, PARAMETER, PUBLIC :: nwet_st_aer = 35 ! Aerosol chemistry ! ATA NLA CheST Chemistry v1.2 @@ -371,19 +371,19 @@ chch_t( 75,'N2 ', 1,'CT ',' ', 0, 0, 0), & ! 76 DD:37 chch_t( 76,'ALICE ', 1,'TR ',' ', 1, 0, 0), & -! 77 -chch_t( 77,'BOB ', 1,'TR ',' ', 0, 0, 0), & +! 77 WD:31 +chch_t( 77,'BOB ', 1,'TR ',' ', 0, 1, 0), & ! 78 chch_t( 78,'DMS ', 1,'TR ',' ', 0, 0, 0), & -! 79 DD:38,WD:31,EM:10 +! 79 DD:38,WD:32,EM:10 chch_t( 79,'SO2 ', 1,'TR ',' ', 1, 1, 1), & ! 80 chch_t( 80,'H2SO4 ', 1,'TR ',' ', 0, 0, 0), & ! 81 chch_t( 81,'MSA ', 1,'TR ',' ', 0, 0, 0), & -! 82 DD:39,WD:32 +! 82 DD:39,WD:33 chch_t( 82,'DMSO ', 1,'TR ',' ', 1, 1, 0), & -! 83 DD:40,WD:33,EM:11 +! 83 DD:40,WD:34,EM:11 chch_t( 83,'NH3 ', 1,'TR ',' ', 1, 1, 1), & ! 84 chch_t( 84,'CS2 ', 1,'TR ',' ', 0, 0, 0), & @@ -395,7 +395,7 @@ chch_t( 87,'SO3 ', 1,'TR ',' ', 0, 0, 0), & ! 88 DD:41, EM: 12 chch_t( 88,'Monoterp ', 1,'TR ',' ', 1, 0, 1), & -! 89 DD:42,WD:34 +! 89 DD:42,WD:35 chch_t( 89,'Sec_Org ', 1,'TR ',' ', 1, 1, 0) & /) @@ -975,13 +975,15 @@ 0.4700E+04, 0.6000E+04, 0.1800E-04, 0.0000E+00, 0.0000E+00, 0.0000E+00,& ! 30 MeOH 0.2300E+03, 0.4900E+04, 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,& -! 31 SO2 +! 31 BOB + 0.2100E+06, 0.8700E+04, 0.2000E+02, 0.0000E+00, 0.0000E+00, 0.0000E+00,& +! 32 SO2 0.1230E+01, 0.3020E+04, 0.1230E-01, 0.2010E+04, 0.6000E-07, 0.1120E+04,& -! 32 DMSO +! 33 DMSO 0.5000E+05, 0.6425E+04, 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,& -! 33 NH3 (H*) +! 34 NH3 (H*) 0.1000E+07, 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00, 0.0000E+00,& -! 34 SEC_ORG +! 35 SEC_ORG 0.1000E+06, 0.1200E+02, 0.0000E-00, 0.0000E+00, 0.0000E+00, 0.0000E+00 & /),(/ 6, nwet_st_aer/) )
Written by Luke Abraham 2014