Minutes of the code management group meeting 2021-02-26
Teleconference Details
Friday 26th February 2021, 10.00am - 11.30am GMT
Meeting via Zoom - details have been emailed
Met Office: AJH, NC, CH, FM, JH, MG, AH, MD
Cambridge: NLA, SAN
Leeds: MR, DG, GM
Reading: NB
NIWA: OM
CSIRO: AL
Apologies: FOC, MS
Actions from this meeting
- (ongoing) LA to put
C_SPECIES
definitions in the Redundant Code table (short-term), and switch to useM_SPECIES
instead. Talk with JH about AQUM tests and requirements. Note that there will be problems with exoplanets and >100km top configs where M_air not constant. LA to make a branch and progress. - (ongoing) Monsoon users MASS data Management. Points-of-contact identified and list of UMUI jobs created up to Dec-2014 circulated to most centres for review.
- (ongoing) GM (supported by KP/MY/AH) has made a table on UKCA wiki documenting GLOMAP versions. ALL to read, GM & others to update table.
- (ongoing) JH to update API document to reflect the UMvn11.9 release
- (new) AJH to investigate which UKCA (and GLOMAP_CLIM) options not protected under rose-stem.
- (new) NLA to chat with OM offline regarding Solar Cycle routine
- (new) NLA to contact JM regarding bug #128 (Bug in UKCA stratospheric chemistry diagnostics)
- (new) AJH to take on bug #129 (warning messages in ukca_radaer_lut_in) as has to change routine anyway for LFRic port
- (new) NLA to make new table for code not tested listing tickets for rose-stem tests that need to be added.
- (new) MR to contact Leighton Regayre regarding flight-track work. Duncan Watson-Paris might have something too.
- (new) NLA to send round Doodle Poll for end of April.
Minutes
Actions from the Last Meeting (LA)
- (ongoing) LA to put
C_SPECIES
definitions in the Redundant Code table (short-term), and switch to useM_SPECIES
instead. Talk with JH about AQUM tests and requirements. Note that there will be problems with exoplanets and >100km top configs where M_air not constant. LA to make a branch and progress for vn11.6.- Ongoing
- (ongoing) Monsoon users MASS data Management
- Status: Points-of-contact identified and list of UMUI jobs created up to Dec-2014 circulated to most centres for review.
- Leave open until first round completed
- (ongoing) GM (supported by KP/MY/AH) has made a table on UKCA wiki documenting GLOMAP versions. ALL to read, GM & others to update table.
- GM has a ticket to take some code back.
- There is also a GLOMAP Dust ticket, necessary for UKCA consolidation away from CLASSIC, as only UKCA will be ported to LFRic.
- Leave open
- (ongoing) NLA to make a ticket to add N48 Linux and update EXP test
- Closed
- (ongoing) NLA to contact AJH/GM offline regarding coupling StratTrop to GLOMAP-CLIM
- Closed. Might align with work of GM/SD
- (ongoing) NLA to contact GF about bug #109 (Bug in dry deposition of CH4 for 13/17/27 tiles)
- Closed
- (new) NLA to remove Old (non-interactive) Dry deposition scheme from the Redundant Code table
- Closed
- (ongoing) JH, NC, LA, FOC, OM (and possibly others) to discuss the need for data from 2-D Cambridge model as well as other top-boundary options for tropospheric configurations. NLA to check with FOC & OM
- FastJX rates have been updated as part of StratO3 work. FJX formatting differences to get in the UM. Cloud-J is also available to be incorporated if desired.
- No-one can run the 2D model any more. 2D model not working “off the shelf”. API work can still progress. NLA to open ticket for rose-stem test, NC to progress API ticket.
- Closed
- (new) AJH to take ticket UM:#5918 (bug #129, "Warning messages in ukca_radaer_lut_in")
- Closed
- (new) JH to update API document to reflect the UMvn11.8 release
- Ongoing - do for UMvn11.9 release
- (new) NLA to email Adrian Lock or Ian Boutle regarding Bug #131 ("Invalid surface roughness input to aerosol dry deposition over land")
- MD has opened a ticket and will progress with suggested replacement STASH code. Should cc GM
- Closed
- (new) AJH to re-write descriptions of Priorities #1 & #7.
- Ongoing
- (new) NLA to put Dry Deposition down as a high(ish) priority.
- Closed
- (new) NLA to send round Doodle poll for February 2021.
- Closed
Standing actions
- (ongoing - standing item) CMG to monitor list of bugs and circulate between group for owners.
- (ongoing - standing item) At CMG future meetings there should be included an appendix listing the open bugs. If a ticket is included it is not an open bug (as the owner of ticket is person doing lodging).
- (ongoing - standing item) ALL with access to keep monitoring ukca_bugs email list.
- (ongoing - standing item) ALL to think about if code changes will change KGO and let systems team know in good time, and remember to include the UKCA keywords to all UKCA tickets. This should be a standard part of our working practice. Note also the required keywords.
- (ongoing - standing item) ALL to note, review, and remove obsolete/redundant code
- (ongoing - standing item) ALL closed ticket owners to produce a summary of their change for the UKCA newsletter for special occasions, e.g. to coincide with Composition-Climate Meeting.
- (ongoing - standing item) Code Owner to produce draft newsletter for CMG to approve.
- (ongoing - standing item) ALL to consider other rose-stem tests that could be included
- (ongoing - standing item) ALL to email LA with suggestions for next set of tutorials
Bugs/Bugfixes (LA)
Open Bugs
# | Description | Ticket | Owner | Comments | Comments from this meeting |
---|---|---|---|---|---|
35 | Error in labelling: Surface resistance <--> stomatal resistance ? | UM:#1409 | Fiona O'Connor | Leave open. Garry Hayman started work in offline JULES. Still ongoing. JM found bug in SO2 deposition, where there is an inconsistency when changing tiles. Making progress. Leave open. GH gave update, has done a lot of work. Concerned it might not be bit-reproducable, but could be solvable. Encouraging to lodge existing. Resistances done as metadata. Not clear as to timescales. FOC to encourage Garry Hayman to start lodging JULES code piecemeal. LA to email GH, cc JM JM to email GH GH adding improvements via UM:#5751 |
|
101 | Errors in reaction rates and diagnostics | UM:#4466 | Luke Abraham / Scott Archer-Nicolls | There is a need to conserve science config for UKESM1, so changes need to be on a flag. Colin has been looking at DMSO, as this is not currently used in StratTrop but is in Offline. Might be able to merge the rate with CEJ changes, but best to leave the diagnostic issues as a separate change. LA, JM, CEJ, SAN to discuss offline and keep GM in loop. Use versioning system. Talk about further offline. Keep ongoing. LA to do with #103 Do for vn11.6. Quite a large ticket. Hope to progress for vn11.6 Bump for vn11.7 LA may be able to progress. Versioning has been done separately. Another issue was added to this recently, to do with formation of H2O - may have an impact Some CRI fixes needed too. May not be able to progress for vn11.7 Some of the fixes have been incorporated in UM:#5057 #5057 didn’t make vn11.8 |
|
109 | Bug in dry deposition of CH4 for 13/17/27 tiles | UM:#4913 | Gerd Folberth | To be Lodged after 10.6 Gerd noted unusual CH4 output in CMIP6 runs. May need a rose stem test for H2 and CH4 drdydep code lodged at UM:#4157 Makes a huge difference to CH4. Needs testing. AJH or GF will do during vn11.4, progress for vn11.6. May need UKESM1 approval. Code is difficult to understand. AJH thinks GarryH is aware of this problem. LA to ask GF to add more details GF has emailed. LA to update wiki pages (see below) Only impacts emission driven runs. LA to check with Gerd about this issue NLA to chase GF |
|
110 | UKCA solar-cycle routine does not pass NRUN-CRUN test | UM:#4828 | Luke Abraham | To be Lodged after 11.5 Have a solution. Solution is a bit tricky to check the logic correctly. Currently testing. May not make vn11.7. Will not make UMvn11.7 - need to analyse runs to understand what is going on NLA to investigate further |
ACTION: NLA to chat with OM offline |
119 | Incorrect calculation for gridbox area - only affects GLOMAP-mode diagnostics | um:#5506 | Calculation of volume appears correct, but area is not. This calculation is missing the factor fv_cos_theta_latitude so will not cause area to tend towards zero at the poles. It also seems to use an unnecessary approximation to the latitude difference (based on sinv_latitude ) which may break down at high latitudes. The calculation should perhaps be consistent with subsequent volume calculation but without the delta_r factor. This will likely only affect diagnostics calculated in ukca_emiss_diags_mode , where it is used to convert from kg/m2/s to mol/gridbox/s, although it is also printed out in ukca_aero_ctl . JH spotted. Just diagnostic affected. LA to open a ticket. cc JM, NC LA to update with MOSRS ticket No owner yet Needs testing. |
GM will have a look, but can’t take it on | |
124 | H2 and N2 cannot be set for RCP and WMOA1 scenario routines | UM:#5620 | Mohit Dalvi | Values are hard-wired in the code and cannot be changed without a branch Will be handled at UM11.8 At least for the next version |
|
128 | Bug in UKCA stratospheric chemistry diagnostics | UM:#5839 | Luke Abraham | pre-UMvn10.0. The indexing/values of diagnostics m01s50i101-i104 are incorrect. Fixes provided for vn11.1 and vn11.7. May affect CMIP6 diagnostics May be best to lump #101. |
ACTION: NLA to contact JM |
129 | Warning messages in ukca_radaer_lut_in | UM:#5918 | maybe Alan J Hewitt | All versions. Code does not crash or clean fail insetad runs until wallclock runs out Fix is obvious - change negative numbers to positive. AJH to take ticket. |
ACTION: AJH to take on, as has to change routine anyway for LFRic port |
130 | Bug in aerosol dry deposition when using 27 surface types in JULES | UM:#6088 Partially fixed using UM:#5622 |
The bug exists because of an erroneous dependency on the array size NTYPE (from JULES_SURFACE_TYPES_MOD) in subroutines UKCA_DDEPAER and UKCA_DDEPAER_INCL_SEDI. This is not appropriate because most of the processing is not written to be dependent on the surface types being used in JULES. A surface type at each location is instead inferred (in UKCA_AERO_CTL) from the roughness length and its index is held in the array ILSCAT. The index values are intended to be interpreted with reference to the JULES 9 type scheme. Only the values 1, 3, 7 and 9 are used, corresponding to forest, grassland, water and ice respectively. The indices are used to access arrays YR, CR and ALPHA that are set up in these aerosol dry deposition routines. The arrays should therefore always be of size 9 with values based on the 9 type scheme, irrespective of the JULES NTYPE value. However, alternative, incompatible arrays are set up if NTYPE is 13, 17 or 27. These were presumably intended to be accessed using the type indices consistent with the other JULES schemes but are actually accessed with the 9-type scheme indices available in ILSCAT. Inappropriate values will therefore be obtained when NTYPE is not 9. Note that Ticket UM:#5622 will fix the bug when interactive dry deposition is not used in the chemistry or when NTYPE in JULES is 13 or 17 but will not fix it when interactive dry deposition is selected (L_UKCA_INTDD = .TRUE.) and NTYPE is 27. Fixing it for NTYPE = 27 will change results so will require a separate bug fix ticket with appropriate testing and approval. Fix is simple (~1-line) but testing will take time. Does affect UKESM1 - anything that doesn’t have 9 surface types. CH to pass on to JM. |
Ticket opened. MD to take ownership | |
131 | Invalid surface roughness input to aerosol dry deposition over land | UM:#6088 | In the UM, UKCA uses the roughness length with stash code 00026 to infer the surface type for aerosol dry deposition (and for the old non-interactive gas dry deposition scheme). However, this value is only valid over the sea. It seems to be set to 0.5 over the land from the UM data I’ve looked at. This seems to be simply a filler value with no scientific significance. As a consequence of the 0.5 value, the UKCA aerosol dry deposition routines categorise all land as forest in the UM runs (routine ukca_aero_ctl uses z0>0.1m for forest). This will interact with bug #130 in aerosol dry deposition (see above) that occurs when using any number of surface types in JULES other than 9. In LFRic JH is using the tile values. Not sure how easy this would be to do in the UM. Could be out of scope for this group - potentially ask Adrian Lock or Ian Boutle. NLA to email. |
Found e-mail from Ian/Adrian (Feb:2016) suggesting 3-028 instead. Ticket opened. MD to take ownership. |
Redundant Code (LA)
Code for deletion
PLEASE REMEMBER TO EDIT THE ORIGINAL REDUNDANT CODE TABLE AND NOT THIS PAGE
- The NG-Composition Project Board has agreed to use the Redundant code table as the Technical Debt Log for UKCA as part of that project. This has highlighted the issue that some areas of code are not listed as being redundant, but are also not tested/used routinely and are not in any rose-stem tests, for example 2D photolysis and the TropIsop and Strat chemistry schemes. We have been asked to review this code and list untested code which is still desirable to be included for testing, and to ensure that other areas of code are listed for deletion.
2D photolysis
At the CMG, we previously agreed to retain the 2D photolysis on the trunk, but a number of issues have been encountered while working on the relevant LFRic/API ticket. This code does not work from the trunk due to missing input files for StratTrop, there is no a rose-stem test for it, and there is currently no way to run the Cambridge 2D model to generate new files when required (e.g., StratTrop, CRI etc. ).
ACTION: NLA to make new table for code not tested listing tickets for rose-stem tests that need to be added.
# | Description | Ticket | Remarks | Status | Timeframe | Owner | Comments from previous meetings | Comments from this meeting |
---|---|---|---|---|---|---|---|---|
4 | ukca_strat_update.F90 - relax_ozone | UM:#1505 | None | PROPOSED | Medium-term | Luke Abraham | ||
5 | Backward Euler Standard Tropospheric Chemistry | UM:#1507 | None | PROPOSED | Medium-term | Fiona O'Connor | ||
6 | Wet oxidation of SO2 in GLOMAP | UM:#4475 | None | PROPOSED | Medium-term | Colin Johnson | JM to look at this re: UKESM1 c.f. offline GLOMAP. JM to take ownership. Still needed for CTM. Relevance for standalone code and later coupling to CTM. |
|
7 | Updraft velocity - l_conv_tke | UM:#3668 | Proposed in December 2017. Approval being sought for GA series (likely GA8). | PROPOSED | Medium-term | Alan J Hewitt / Ian Boutle | ||
8 | All backward Euler chemistry schemes | None | Note that this is only for the “hard-coded” BE schemes. Quite a lot of work would be involved, but good as a longer-term aim. Note that problems for AQUM would need to be addressed first. | PROPOSED | Long-term | Luke Abraham | ||
9 | Deprecated UKCA netcdf attribute options | UM:#2513 | 'emission_type', string type attributes no longer recommeneded. Need some time for most new code to adopt these. Still used in AQUM files. | DEPRECATED | Medium-term | Nick Savage | AQUM are addressing the files. Have gone from operational jobs, have ticket to remove from rose-stem. LN to document. These have been removed from AQUM. AQUM work has been done. From AS: In the ancil ticket linked from #2513 (https://code.metoffice.gov.uk/trac/ancil/ticket/782) it says that GA7 depends on the old format files, and as far as I know this is still true. In the 3-step timeline suggested, we haven't reached step 2 unfortunately. That ticket was about ensuring GA8 did not use the old attributes. I believe that GA8 is due for release in the next few months (it has been a very long time coming), so UM vn11.5 may be the last version at which GA7 is supported. Do you mind waiting a bit longer? AJH to follow-up with LA |
|
10 | ASAD families chemistry | None | Families chemistry has not been used for many years (since at least UM6.1. It is likely that it no longer works if turned on, and it complicates the code making other developments harder | DEPRECATED | Short-term | Luke Abraham | ||
11 | Dry deposition 9 PFT option | None | The 13 PFT option has been adopted for UKESM1. The 9 PFT option is possibly not being used at all. Proposed in June 2018. |
PROPOSED | Long-term | Alan J Hewitt | ||
12 | Remove explicit C_SPECIES constants and calculate these from M_SPECIES and the mass of air |
None | Will have wide-ranging rose-stem KGO breakages | PROPOSED | Short-term | Luke Abraham | ||
13 | Remove temporary logical l_fix_nacl_density Introduced at UM:#2405 |
UM:#3490 | Requires AMIP testing and a GMED ticket Could be accepted to GA9 |
PROPOSED | Medium-term | Jane Mulcahy | See JM GMED action | |
15 | Remove temporary logical l_fix_ukca_impscav Introduced at UM:#1421 |
UM:#4145 | Approval being sought for GA series (likely GA9) GMED:#402 |
PROPOSED | Medium-term | Alan J Hewitt | See JM GMED action | |
16 | Remove temporary logical l_fix_improve_drydep Introduced at UM:#1250 |
UM:#3997 | Approval being sought for GA series (likely GA9) GMED:#393 |
PROPOSED | Medium-term | Alan J Hewitt | See JM GMED action | |
17 | Remove temporary logical l_fix_ukca_h2dd_x Introduced at UM:#4157 |
UM:#4245 | Approval being sought for GA series (likely GA9) GMED:#410 |
PROPOSED | Medium-term | Alan J Hewitt | See JM GMED action | |
18 | Remove temporary logical l_fix_neg_pvol_wat Introduced at UM:#4383 |
UM:#4416 | Approval being sought for GA series (likely GA9) GMED:#421 |
PROPOSED | Medium-term | Alan J Hewitt | ||
21 | Remove code behind L_use_stratclim from ukca_fastjx_mod Introduced prior to MOSRS |
UM:#4854 | This could either be properly included via a namelist or removed from the UM | PROPOSED | Medium-term | Alan J Hewitt | LA contact OW about progressing FastJX-RADAER work Ties in with API work. MD/NC to take a look, can be removed. Do as part of API. Relevant person take ownership as part of API work NC to take ownership as part of his work |
|
22 | Remove temporary logical l_fix_ukca_input_meoh_emiss Introduced at UM:#3031 |
UM:#4848 | Approval being sought for GA series (likely GA9) GMED:#461 |
PROPOSED | Medium-term | Mohit Dalvi | ||
24 | Remove ability to use of Cambridge 2D model data from ukca_stratf.F90 | PROPOSED | Neal Carr | Need to confirm from wider community for possible tropospheric configs and other choices for top-bdy |
Update on UKCA API redesign
Please see the document UKCA_API_design_proposal_v5.pdf for the latest API design.
Work progressing. When making changes around API, please refer to appendix as to which routines are UM and which are UKCA to ensure to keep the namespaces separate
Progress with Priorities (LA)
ACTION REQUEST FROM UKCA SCIENCE & MANAGEMENT BOARD: CMG to provide list of priorities & their status, with particular thought to UKESM2 changes that will be required.
The current priorities are:
# | Description | Owner | Desirable for UKESM2? | Status |
---|---|---|---|---|
1 | UKCA API | John Hemmings (main) | No | UKCA API work is progressing (AJH). Radaer API on hold until after GA7 port to LFRic. Functionality of RADAER in LFRic will work without an API although it will be preferable to update this in future. |
2 | Review UKCA tests to cover all requirements | Luke Abraham / Mohit Dalvi | ||
3 | Updates/fixes to GLOMAP routines | YES | ||
4 | Standalone UKCA configurations | Luke Abraham/Scott Archer-Nicholls | ExCALIBUR work progressing, currently working on build system. | |
5 | Emissions flexibility | Mohit Dalvi | ||
6 | Separare UKCA repository | Luke Abraham | Separate UKCA trac environment set-up on MOSRS here, will hold UKCA code on completion of current ExCALIBUR project. | |
7 | Non-API LFRic work, e.g. Porting science functionality etc. | Alan J Hewitt | No | Porting calculation of CDNC from GLOMAP Climatology aerosols complete. I will pass required elements of ukca_radaer to ukca_radaer_band_average (um:#5957) Porting calculation of fields in radaer_band_average to be passed to Socrates underway(lfric:#2488) |
8 | RADAER API | Alan Hewitt (main) | No | Radaer API on hold until after GA7 port to LFRic. Functionality of RADAER in LFRic will work without an API although it will be preferable to update this in future. Not clear if RADAER will work independently of UM in single column mode. AJH and AH and JH and Ben Johnson and NB need to work out if desirable and practical to run independent of UM. |
9 | RADAER optimisation | NO | ||
10 | Dry deposition and Jules tiles | YES | If Jules adds extra tiles for UKESM2 we may need to hard code options in dry deposition (like in um:#1761) | |
11 | FastJX-GLOMAP coupling | YES | UMvn10.8 branches exist - see ticket UM:#2689 Note: FastJX rates also need updating - potential for code as well as input-file changes | |
12 | Cloud-J | |||
13 | Fixes for chemistry rates and ASAD (bugs #101 & #103) | Luke Abraham | YES | #101 still open, #103 completed (necessary for GC5 ticket GMED:#551) |
14 | Rc fixes for water (bug #104) | UKESM Core Team | YES | Lodged at UMvn11.7 in ticket UM:#5167 |
15 | CH4 deposition fixes (bug #109) | Gerd Folberth | YES | |
16 | UKCA STASHmaster PP-codes | Luke Abraham | Verification system makes some use of existing PP-codes. New items should have a PP-code of 0 and verification team will use STASH items for these if required. | |
17 | Removal of redundant code | ALL | No | Parked whilst GA7.0 is a protected suite. |
18 | GLOMAP-CLIM API refresh | Alan J Hewitt | No | Parked until science functionality around GLOMAP aerosols ported to LFRic. |
19 | N48L85 configuration | |||
20 | ACTIVATE optimisation | TBC | No | A quick win here is to reduce nwbins in namelist ( see um:#4740 ). Parked for now. |
21 | Stratospheric Ozone improvements | James Keeble/Luke Abraham | YES | Work ongoing, covered in ticket UM:#5641. Note: FastJX rates also need updating - potential for code as well as input-file changes |
22 | Solar cycle in photolysis | Luke Abraham | YES | Some code lodged at UMvn10.9 (UM:#2587) but needs bugfixing. This work has stalled in ticket UM:#4828 (bug #110). |
23 | Deposition updates from CSIRO | YES | Code lodged in UM:#4020 (UMvn11.4, oceanic dry deposition) and UM:#5713 (UMvn11.8, Lightning NOx updates). |
UM11.9 Tickets
See here for UM vn11.9 tickets on MOSRS.
See here for UM vn12.0 tickets on MOSRS.
PLEASE REMEMBER TO USE THE UKCA KEYWORD FOR ALL UKCA CHANGES
NOTE: New requirement for LFric testing
of all changes under src/atmosphere
Round (virtual) Table of current activities, plans
Person | Activity |
---|---|
NLA | ARCHER2 eCSE project to work on UKCA flight-track simulator, to begin in June 2021 |
MR | involved with NGMS as external collab. - talking about training. Models not quite ready yet. May make someone in their team who would be a LFRic |
LFRic port taking shape - might have exciting results for next time
ACTION: MR to contact Leighton Regayre regarding flight-track work. Duncan Watson-Paris might have something too.
AOB/DNM (LA)
AJH - SUSSBCOCDU7mode is being ported to LFRic. There is currently no rose stem test for i_mode_setup==8. There is a parked branch (see um:#6032). Would someone be able to take over?
- this is with dust, might not happen this release. Get in touch with AJH if interested
AL asked about the UKCA evaluation suite and whether there was a mechanism to add new datasets. NLA suggested contacting MD. CH asked about ESMValTool use. There is a summer student coming to look at this. People should contact CH if they want to discuss this further.
Date of next meeting (LA)
ACTION: NLA to send round Doodle Poll for end of April.