Difference between revisions of "Solution to UKCA Chemistry and Aerosol Tutorial 9 Task 9.1"
From UKCA
m (Protected "Solution to UKCA Chemistry and Aerosol Tutorial 9 Task 9.1" ([edit=sysop] (indefinite) [move=sysop] (indefinite)) [cascading]) |
(Blanked the page) |
||
| Line 1: | Line 1: | ||
| − | [[UKCA Chemistry and Aerosol Tutorials | Back to UKCA Chemistry and Aerosol Tutorials]] |
||
| − | |||
| − | [[UKCA Chemistry and Aerosol Tutorial 9#Task 9.1: adding new dry deposition values | Back to the adding dry deposition tutorial]] |
||
| − | |||
| − | ==Task== |
||
| − | |||
| − | You were asked to |
||
| − | |||
| − | <blockquote> |
||
| − | You should now add in the dry deposition of '''ALICE'''. This species deposits in a similar way to '''CO'''. The values for '''depvel_defs_strattrop''' are: |
||
| − | |||
| − | {| border="1" |
||
| − | ! Surface Type || Summer (day) || Summer (night) || Summer (24h ave) || Winter (day) || Winter (night) || Winter (24h ave) |
||
| − | |- |
||
| − | | '''Water''' || 0.00 || 0.00 || 0.00 || 0.00 || 0.00 || 0.00 |
||
| − | |- |
||
| − | | '''Forest''' || 0.03 || 0.03 || 0.03 || 0.03 || 0.03 || 0.03 |
||
| − | |- |
||
| − | | '''Grass''' || 0.03 || 0.03 || 0.03 || 0.03 || 0.03 || 0.03 |
||
| − | |- |
||
| − | | '''Desert''' || 0.03 || 0.03 || 0.03 || 0.03 || 0.03 || 0.03 |
||
| − | |- |
||
| − | | '''Ice''' || 0.00 || 0.00 || 0.00 || 0.00 || 0.00 || 0.00 |
||
| − | |} |
||
| − | |||
| − | i.e. the same as for CO. |
||
| − | </blockquote> |
||
| − | |||
| − | and were given the hint |
||
| − | |||
| − | <blockquote> |
||
| − | You will need to make changes for both dry deposition schemes. The changes to '''ukca_surfddr.F90''' can be made very easily by adding '''ALICE''' to the '''CO''' block in the <tt>CASE</tt> statement. |
||
| − | </blockquote> |
||
| − | |||
| − | ==Solution== |
||
| − | |||
| − | You will need to make changes to 3 UKCA routines and make 1 hand-edit. |
||
| − | |||
| − | ===Changes to ukca_chem_strattrop.F90=== |
||
| − | |||
| − | ====chch_defs array==== |
||
| − | |||
| − | First, you will need to change the '''0''' in the 6th column of the <tt>chch_defs_strattrop_chem</tt> to '''1''': |
||
| − | |||
| − | chch_t( 76,'ALICE ', 1,'TR ',' ', 1, 0, 0), & ! 76 |
||
| − | |||
| − | ====Additions for depvel_defs==== |
||
| − | |||
| − | To add in the values for the deposition velocity into the '''depvel_defs_strattrop''' array you will need to |
||
| − | |||
| − | * Increase the value of '''ndry_strattrop''' from 36 to '''37''' |
||
| − | * Increase the size of '''depvel_defs_strattrop03''' from 360 to '''390''' |
||
| − | * Add this code block to end of the '''depvel_defs_strattrop03''' array (remembering to add a comma to end of the line preceding this) |
||
| − | |||
| − | 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, & ! 37.1 |
||
| − | 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 37.2 |
||
| − | 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 37.3 |
||
| − | 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 37.4 |
||
| − | 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 & ! 37.5 |
||
| − | |||
| − | ===Changes to ukca_aerod.F90=== |
||
| − | |||
| − | You will need to add the following block to the '''CASE''' statement within this routine |
||
| − | |||
| − | CASE ('ALICE ') |
||
| − | d0(j) = d_h2o * SQRT(m_h2o / m_ALICE) |
||
| − | |||
| − | the value of m_ALICE having been set in [[Solution to UKCA & UMUI Tutorial 5 Task 5.2#Code changes to ukca_constants.F90 | Task 5.2 - adding new emissions to UKCA]]. |
||
| − | |||
| − | ===Changes to ukca_surfddr.F90=== |
||
| − | |||
| − | As has been noted in the ''hint'' for this task, the only change required to '''ukca_surfddr.F90''' is to add <tt>'ALICE '</tt> to the '''CO''' section of the '''CASE''' statement: |
||
| − | |||
| − | CASE ('CO ','ALICE ') |
||
| − | rsurf(:,n)=(/3700.,7300.,4550.,1960.,4550.0,r_null,r_null, & |
||
| − | 4550.0,r_null /) ! Shrub+bare soil set to C3 grass (guess) |
||
| − | |||
| − | ===Hand-edit to increase the value of JPDD=== |
||
| − | |||
| − | You will need to add a hand-edit to 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. |
||
| − | |||
| − | It should contain the following: |
||
| − | |||
| − | ed CNTLATM<<\EOF |
||
| − | /JPDD/ |
||
| − | d |
||
| − | i |
||
| − | JPDD = 37, |
||
| − | . |
||
| − | w |
||
| − | q |
||
| − | EOF |
||
| − | |||
| − | An example hand-edit can be found at |
||
| − | |||
| − | /home/ukca/hand_edits/VN8.2/Tutorial/Task7.1_incr_JPvals.ed |
||
| − | |||
| − | ==Output== |
||
| − | |||
| − | If you view the ALICE (34064) 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 6 Task 6.1 | previous Task (6.1)]] you will see that there are some difference. 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/Task7.1/ |
||
| − | |||
| − | on HECToR, and in |
||
| − | |||
| − | /projects/ukca/Tutorial/sample_ouput/Task7.1/ |
||
| − | |||
| − | on MONSooN. |
||
| − | |||
| − | ==Worked Solution== |
||
| − | |||
| − | There is a worked solution to this problem in the UMUI Tutorial experiment. This is job '''h''': ''Tutorial: solution to Task 7.1 - add new dry 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@12170 fcm:um_br/dev/luke/vn8.2_UKCA_Tutorial_Solns@12214 |
||
| − | |||
| − | This gives the following (non-graphical) output: |
||
| − | |||
| − | Index: src/atmosphere/UKCA/ukca_surfddr.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_surfddr.F90 (revision 12170) |
||
| − | +++ src/atmosphere/UKCA/ukca_surfddr.F90 (revision 12214) |
||
| − | @@ -277,7 +277,7 @@ |
||
| − | 500.0, 12500. /) |
||
| − | CASE ('NH3 ') |
||
| − | rsurf(:,n)=hno3 |
||
| − | - CASE ('CO ') |
||
| − | + CASE ('CO ','ALICE ') |
||
| − | rsurf(:,n)=(/3700.,7300.,4550.,1960.,4550.0,r_null,r_null, & |
||
| − | 4550.0,r_null /) ! Shrub+bare soil set to C3 grass (guess) |
||
| − | CASE ('CH4 ') |
||
| − | Index: src/atmosphere/UKCA/ukca_chem_strattrop.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 12170) |
||
| − | +++ src/atmosphere/UKCA/ukca_chem_strattrop.F90 (revision 12214) |
||
| − | @@ -56,7 +56,7 @@ |
||
| − | INTEGER, PARAMETER, PUBLIC :: nhet_st_tpht = 2 ! trophet rxns |
||
| − | |||
| − | ! No of dry deposited species |
||
| − | -INTEGER, PARAMETER, PUBLIC :: ndry_strattrop = 36 ! Stratospheric chemistry |
||
| − | +INTEGER, PARAMETER, PUBLIC :: ndry_strattrop = 37 ! Stratospheric chemistry |
||
| − | INTEGER, PARAMETER, PUBLIC :: ndry_st_aer = 5 ! Aerosol chemistry |
||
| − | |||
| − | ! No of wet deposited species |
||
| − | @@ -141,7 +141,7 @@ |
||
| − | chch_t( 73,'CO2 ', 1,'CT ',' ', 0, 0, 0), & ! 73 |
||
| − | 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 ',' ', 0, 0, 0), & ! 76 |
||
| − | +chch_t( 76,'ALICE ', 1,'TR ',' ', 1, 0, 0), & ! 76 |
||
| − | chch_t( 77,'BOB ', 1,'TR ',' ', 0, 0, 0) & ! 77 |
||
| − | /) |
||
| − | |||
| − | @@ -759,7 +759,7 @@ |
||
| − | |||
| − | REAL, DIMENSION( 360) :: depvel_defs_strattrop01 |
||
| − | REAL, DIMENSION( 360) :: depvel_defs_strattrop02 |
||
| − | -REAL, DIMENSION( 360) :: depvel_defs_strattrop03 |
||
| − | +REAL, DIMENSION( 390) :: depvel_defs_strattrop03 |
||
| − | |||
| − | INTEGER(KIND=jpim), PARAMETER :: zhook_in = 0 |
||
| − | INTEGER(KIND=jpim), PARAMETER :: zhook_out = 1 |
||
| − | @@ -2120,7 +2120,13 @@ |
||
| − | 0.83, 0.04, 0.44, 0.06, 0.05, 0.06, & ! 36.2 |
||
| − | 0.63, 0.06, 0.35, 0.08, 0.06, 0.07, & ! 36.3 |
||
| − | 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 36.4 |
||
| − | - 0.01, 0.01, 0.01, 0.01, 0.01, 0.01 & ! 36.5 |
||
| − | + 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, & ! 36.5 |
||
| − | +! 37 ALICE - same as for CO |
||
| − | + 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, & ! 37.1 |
||
| − | + 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 37.2 |
||
| − | + 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 37.3 |
||
| − | + 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, & ! 37.4 |
||
| − | + 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 & ! 37.5 |
||
| − | /) |
||
| − | |||
| − | |||
| − | Index: src/atmosphere/UKCA/ukca_aerod.F90 |
||
| − | =================================================================== |
||
| − | --- src/atmosphere/UKCA/ukca_aerod.F90 (revision 12170) |
||
| − | +++ src/atmosphere/UKCA/ukca_aerod.F90 (revision 12214) |
||
| − | @@ -222,6 +222,8 @@ |
||
| − | d0(j) = d_h2o * SQRT(m_h2o / m_meoh) |
||
| − | CASE ('Monoterp') |
||
| − | d0(j) = d_h2o * SQRT(m_h2o / m_monoterp) |
||
| − | + CASE ('ALICE ') |
||
| − | + d0(j) = d_h2o * SQRT(m_h2o / m_ALICE) |
||
| − | END SELECT |
||
| − | END DO |
||
| − | ! |
||
| − | |||
| − | ---- |
||
| − | ''Written by [[User:Nla27 | Luke Abraham]] 2014'' |
||
