Difference between revisions of "Solution to UKCA Chemistry and Aerosol Tutorial 10 Task 10.1"
From UKCA
(Created page with " Back to UKCA Chemistry and Aerosol Tutorials [[UKCA Chemistry and Aerosol Tutorial 10#Task 10.1: Add wet deposition of a species | Bac…") |
(Blanked the page) |
||
Line 1: | Line 1: | ||
− | [[UKCA Chemistry and Aerosol Tutorials | Back to UKCA Chemistry and Aerosol Tutorials]] |
||
− | |||
− | [[UKCA Chemistry and Aerosol Tutorial 10#Task 10.1: Add wet deposition of a species | Back to the adding wet deposition tutorial]] |
||
− | |||
− | ==Task== |
||
− | |||
− | You were asked to |
||
− | |||
− | <blockquote> |
||
− | Add in wet deposition for '''BOB''', using the following values: |
||
− | |||
− | {| border="1" |
||
− | ! <math>\ k(298)\ </math> || <math>\ -\left({\Delta H}/R\right)\ </math> || <math>\ k(298)</math> for the 1st dissociation || <math>\ -\left({\Delta H}/R\right)</math> for the 1st dissociation || <math>\ k(298)</math> for the 2nd dissociation || <math>\ -\left({\Delta H}/R\right)</math> for the 2nd dissociation |
||
− | |- |
||
− | | <math>\ 0.21 \times 10^{+06}\ </math> || <math>\ 0.87 \times 10^{+04}\ </math> || <math>\ 0.2 \times 10^{+02}\ </math> || <math>\ 0.0\ </math> || <math>\ 0.0\ </math> || <math>\ 0.0\ </math> |
||
− | |} |
||
− | </blockquote> |
||
− | |||
− | ==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 '''[[UKCA & UMUI Tutorials: Things to know before you start#Archiving|archive]]''' directory) after these changes, and compare it to the equivalent field after the [[Solution to UKCA & UMUI Tutorial 7 Task 7.1 | previous Task (7.1)]] you will see that there are some differences. However, we will not be able to quantify these differences correctly until [[UKCA & UMUI Tutorial 9 | 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 [[User:Nla27 | Luke Abraham]] 2014'' |