Solution to UKCA & UMUI Tutorial 7 Task 7.1

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