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:

$\ k(298)\$	$\ -\left({\Delta H}/R\right)\$	$\ k(298)$ for the 1st dissociation	$\ -\left({\Delta H}/R\right)$ for the 1st dissociation	$\ k(298)$ for the 2nd dissociation	$\ -\left({\Delta H}/R\right)$ for the 2nd dissociation
$\ 0.21\times 10^{+06}\$	$\ 0.87\times 10^{+04}\$	$\ 0.2\times 10^{+02}\$	$\ 0.0\$	$\ 0.0\$	$\ 0.0\$

Solution

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).

Hand-edit

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

/home/ukca/hand_edits/VN8.2/Tutorial/Task8.1_incr_JPvals.ed

Output

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

/work/n02/n02/ukca/Tutorial/sample_output/Task8.1/

on HECToR, and in

/projects/ukca/Tutorial/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 -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 2014