Difference between revisions of "Solution to UKCA & UMUI Tutorial 7 Task 7.1"
Line 33: | Line 33: | ||
You will need to make changes to 3 UKCA routines and make 1 hand-edit. |
You will need to make changes to 3 UKCA routines and make 1 hand-edit. |
||
− | ===Changes to |
+ | ===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 preceeding 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=== |
||
+ | |||
+ | You will need to add a hand-edit to UMUI in the '''Model Selectiom → 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 |
||
+ | |||
+ | ==Worked Solution== |
||
+ | |||
+ | |||
---- |
---- |
Revision as of 14:51, 24 June 2013
Back to the adding dry deposition tutorial
Task
You were asked to
You should now add in the dry deposition of ALICE. This species deposits in a similar way to CO. You will need to adjust depvel_defs_scheme as well ukca_aerod.F90 and ukca_surfddr.F90. The values for depvel_defs_scheme are:
0.00 0.00 0.00 0.00 0.00 0.00 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.03 0.00 0.00 0.00 0.00 0.00 0.00 i.e. the same as for CO.
and were given the hint
The changes to ukca_surfddr.F90 can be made very easily by adding ALICE to the CO block in the CASE statement.
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 chch_defs_strattrop_chem 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 preceeding 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 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 'ALICE ' 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
You will need to add a hand-edit to UMUI in the Model Selectiom → 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
Worked Solution
Written by Luke Abraham 2013