Solution to UKCA Chemistry and Aerosol Tutorial 8 Task 8.1

From UKCA
Revision as of 14:19, 15 December 2015 by Nla27 (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

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_aer 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