Release Job RJ4.0

From UKCA
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

This page documents the vn8.4 GA4.0 CheST+GLOMAP-mode release job RJ4.0 (xlava on MONSooN and xlavb on ARCHER).

This job has been developed by Luke Abraham with help from many others.

It is a development of Release Candidate RC6.0 and Release Candidate RC6.2. Further details on the development of this job can be found on those pages.

Suitability for Release

Tropospheric Chemistry Users should note the tropospheric OH biases
Stratospheric Chemistry Users should note the low stratospheric NOy.
Tropospheric Aerosol Surface properties are green, while column properties are amber
Stratospheric Aerosol

Key

Release candidate still under evaluation
Release candidate not scientifically suitable to be released
Release candidate is suitable for development jobs,
and may be scientifically suitable to be released for some applications
Release candidate scientifically suitable to be released
Release candidate does not consider the required chemistry or aerosol processes

Overview

MONSooN job xlava ARCHER job xlavb
UKCA branch used fcm:um_br/pkg/Config/vn8.4_UKCA@18311
Original revision #: 16827
Current revision #: 18311
fcm:um_br/pkg/Config/vn8.4_UKCA@18311
Original revision #: 16827
Current revision #: 18311
Decomposition 12EW x 16NS on 6 nodes (see here) 12EW x 12NS on 12 nodes (see here)
Run-time (per model month) ~2 hours (see here) ~2 hours (see here)
Job-step (recommended) 1 model month in 10800 seconds (3 hours) 3 model months in 28800 seconds (8 hours)
Cost (per model year) 150 node-hours (see here) 115 kAU (see here)
Storage Requirements
(current STASH settings)
110GB per model year (32-bit pp-files & seasonal 64-bit dumps)
MOOSE costs: £12.00 per model year
110GB per model year (32-bit pp-files & seasonal 64-bit dumps) copied to the /nerc disk as the model runs,
using the branch fcm:um_br/dev/luke/vn8.4_hector_monsoon_archiving_ff2pp/src as a central script modification.


Functionality

Package Branch

The UKCA package branch fcm:um_br/pkg/Config/vn8.4_UKCA/src has been purposefully left without a revision number so that the most recent one will always be taken.

It is recommended to fix the revision number when performing any production runs, as this will ensure that there is a record of the state of the code when the run was done. To see which revision number is the most recent, please see the branch on PUMA here:

Bugfixes

2015-02-02 r18072

There were 2 updates required to the original release revision of 16827:

  • AA=aa_modes(IMODE) needed to be added to a loop in ukca_conden.F90.
  • all_tracers needed to be references with the levels defined as 1:model_levels to all for compatibility with ENDGame.

These were added at revision 18072. For further details, please see changeset 18702.

2015-02-13 r18311

There was 1 update required to the previous release revision of 18072:

  • This changeset adds a trap on SQD(:) to prevent floating point exceptions in UKCA_SOLVECOAGNUCL_V when H2SO4 evaporation is switched on.

This was added at revision 18311. For further details, please see changeset 18311.

2015-06-02 SST and Sea-Ice ancillaries

It was spotted that the dates in the SST and sea-ice ancillaries used by this job were incorrect. They were set to be daily from 0000-01-16 00:00:00, when in fact they should have been daily from 0000-01-01 12:00:00. These files have been regenerated, and updated files have been placed in the same directories as the originals. The UMUI jobs under ukca have been updated, but users should update their own jobs manually as follows:

  • reynolds.sst.avg2000 should be replaced by reynolds.sst.avg2000_fix
  • reynolds.seaice.avg2000 should be replaced by reynolds.sice.avg2000_fix (note shortening of seaice)

Using the old ancillaries would mean that the SSTs and sea-ice fields were actually incorrect by -15 days (i.e. the 1st January fields were actually from the 16th December etc.).

2015-07-31 RADAER pcalc_hadgem_v2.ukca file

Due to a bug in the program that creates the pcalc files, whereby input refractive indices did not cover the full range of wavelengths for organic carbon (The first LW waveband goes to 1 cm), this file has been updated to pcalc_hadgem_v2.ukca_fix, and can be found in the

/work/n02/n02/ukca/spectral/radv2/

directory on ARCHER, and the

/projects/ukca/inputs/spectral/radv2/

directory on MONSooN.

The impact of this bug is negligible because:

  1. refractive indices cannot effectively exceed the limits set in RADAER look-up tables, which are 3 - 1i in the LW;
  2. specific extinction in the LW is weak for Aitken and accumulation modes, where organic aerosols are located.

The shortwave spectrum, AOD and backscatter diagnostics are unaffected. James Mollard from the University of Reading has tested the impact, and found that the differences are within the inter-annual variability.

As well as producing this bugfixed file, Nicolas Bellouin has also updated to v2.1:

  • pcalc_hadgem_v2.1.ukca: using Bond and Bergstrom (2006) values
  • pcalc_hadgem_v2.1.ukca_bs: for use with backscatter diagnostics at CSIRO

It is recommend to upgrade to pcalc_hadgem_v2.1.ukca

The RJ4.0/xlav jobs have been upgraded to use the pcalc_hadgem_v2.1.ukca file. It is possible to use either the original pcalc_hadgem_v2.ukca file, or the pcalc_hadgem_v2.1.ukca_fix file with minimal UMUI changes.

MONSooN

UMUI changes made, compared to original release:

Difference in window subindep_HandEdit
 -> Model Selection
   -> Input/Output Control and Resources
     -> User hand edit files 
Differences in Table Hand edits
 5c5 
<  ~mdalvi/umui_jobs/hand_edits/vn8.4/raderv2_vn84.ed Y 
---
>  ~ukca/hand_edits/VN8.4/raderv2.1_vn84_MONSooN.ed Y

Difference in window atmos_Science_Section_UKCA_Rad
 -> Model Selection
   -> Atmosphere
     -> Scientific Parameters and Sections
       -> Section by section choices
         -> Section 34: UKCA Chemistry and Aerosols
           -> UKCA Chemistry Coupling
             -> MODE Aerosols in Radiation Scheme
Entry box: Directory path to UKCA_RADAER input files:
 Job xlaba: Entry is set to '/projects/ukca/inputs/spectral/radv2'
 Job xlava: Entry is set to '/projects/ukca/inputs/spectral/radv2.1'
Entry box: File of precomputed values:
 Job xlaba: Entry is set to 'pcalc_hadgem_v2.ukca'
 Job xlava: Entry is set to 'pcalc_hadgem_v2.1.ukca'

To use the pcalc_hadgem_v2.1.ukca_fix (or pcalc_hadgem_v2.1.ukca) file, change the hand-edit back to ~mdalvi/umui_jobs/hand_edits/vn8.4/raderv2_vn84.ed and make the required changes to the MODE Aerosols in Radiation Scheme UMUI panel.

ARCHER

UMUI changes made, compared to original release:

Difference in window subindep_HandEdit
 -> Model Selection
   -> Input/Output Control and Resources
     -> User hand edit files
Differences in Table Hand edits
 5c5 
<  ~ukca/hand_edits/VN8.4/raderv2_vn84_ARCHER.ed Y
---
>  ~ukca/hand_edits/VN8.4/raderv2.1_vn84_ARCHER.ed Y

Difference in window atmos_Science_Section_UKCA_Rad
 -> Model Selection
   -> Atmosphere
     -> Scientific Parameters and Sections
       -> Section by section choices
         -> Section 34: UKCA Chemistry and Aerosols
           -> UKCA Chemistry Coupling
             -> MODE Aerosols in Radiation Scheme
Entry box: Directory path to UKCA_RADAER input files:
 Job xlabb: Entry is set to '/projects/ukca/inputs/spectral/radv2'
 Job xlavb: Entry is set to '/work/n02/n02/ukca/spectral/radv2.1'
Entry box: Look-up table for accumulation-mode aerosols in the longwave:
 Job xlabb: Entry is set to 'nml_ac_lw'
 Job xlavb: Entry is set to 'nml_ac_lw_new'
Entry box: Look-up table for accumulation-mode aerosols in the shortwave:
 Job xlabb: Entry is set to 'nml_ac_sw'
 Job xlavb: Entry is set to 'nml_ac_sw_new'
Entry box: Look-up table for coarse-mode aerosols in the longwave:
 Job xlabb: Entry is set to 'nml_cr_lw'
 Job xlavb: Entry is set to 'nml_cr_lw_new'
Entry box: Look-up table for coarse-mode aerosols in the shortwave:
 Job xlabb: Entry is set to 'nml_cr_sw'
 Job xlavb: Entry is set to 'nml_cr_sw_new'
Entry box: File of precomputed values:
 Job xlabb: Entry is set to 'pcalc_hadgem_v2.ukca'
 Job xlavb: Entry is set to 'pcalc_hadgem_v2.1.ukca'

There were further changes made to the ARCHER job (xlavb) as compared to the MONSooN job, as when the job was ported from MONSooN, the RADAER files were set entirely by hand-edit, rather than changes being made to the UMUI. Additionally, the Cray Fortran compiler (cce/8.2.1, and also cce/8.3.7 as is used currently) requires fortran namelists to be reformatted. This was done to the RADAER files, hence these being nml_??_?w_new.

To use the pcalc_hadgem_v2.1.ukca_fix (or pcalc_hadgem_v2.1.ukca) file, change the hand-edit to ~ukca/hand_edits/VN8.4/raderv2.0_vn84_ARCHER.ed and make the required changes to the MODE Aerosols in Radiation Scheme UMUI panel. This new hand-edit only sets the required UKCAANSW and UKCAANLW namelist entries.

Pre-compiled builds

ARCHER

With the move to the Cray cce/8.3.7 compiler, pre-compiled builds do not currently work (as these were for the previous Cray compiler, cce/8.2.1). This has been updated in the jobs on the PUMA UMUI, but you should turn them off when compiling jobs afresh. While executables built previously will still work, new jobs should use the new compiler. If you build the reconfiguration executable this should be set to

  • Directory for the Reconfiguration executable: $DATAW/bin
  • Filename for the Reconfiguration executable: qxreconf

in the Compile options for the atmosphere and reconfiguration panel.

Tests comparing the RJ4.0 release job using cce/8.2.1 and cce/8.3.7 have shown that the change in the compiler should not change the results, i.e. executables built using the different compilers produce output that bit-compares.

MONSooN ibm02

These jobs make use of pre-compiled builds for the base GA4.0 model. This significantly increases compile time. However, a side-effect of this is that the reconfiguration executable must be taken from the pre-build.

When compiling the job, you MUST NOT ask for the reconfiguration executable to be built. Doing so will cause an error.

If you wish to compile the reconfiguration executable, you will need to turn off pre-compiled build in the FCM panel.

ARCHER: /nerc archiving

The ARCHER job is configured to convert the fieldsfile (64-bit) output to .pp (32-bit) output and then copy these across to the /nerc disk. To use this functionality, you will need to set up automatic archiving as detailed here:

Turn off archiving altogether

If you do not wish to archive at all, please turn this off in the

Model Selection
-> Post Processing
 -> Main Switch + General Questions

To turn off conversion to .pp (32-bit) format change the value of FF2PP_HECTOR from Y to N in

Model Selection
-> Input/Output Control and Resources
 -> Script Inserts and Modifications

For long runs, some form of archiving is encouraged, however for short testing runs it may be best to turn off archiving.

Archiving to /work

It is possible to archive to /work instead of /nerc. This avoids the additional steps required for automatic archiving. However archiving to /nerc is encouraged as this puts less pressure on the /work filesystem.

To change the directory that the data is archived to, go to

Model Selection
-> Post Processing
 -> Main Switch + General Questions

and change the name of the archive directory to

$DATADIR/$RUNID/archive

This will still allow for conversion to .pp format and the removal of dumps.

UKCA Evaluation Suite Diagnosics

Output which is needed for the UKCA Evaluation Suite has been automatically added to STASH via the inclusion of the hand-edit

~mdalvi/umui_jobs/hand_edits/vn8.4/add_ukca_eval1_diags_l85.ed

The STASH that is added by this hand-edit is in the table below. It has been sent to the climate meaning stream (UPMEAN). Even if all the other diagnostics have been turned off in STASH, these fields will continue to be output unless the hand-edit is also turned off.

STASH Code Diagnostic
10 SPECIFIC HUMIDITY AFTER TIMESTEP
408 PRESSURE AT THETA LEVELS AFTER TS
30453 Height at Tropopause Level
34001 O3 MASS MIXING RATIO AFTER TIMESTEP
34002 NO MASS MIXING RATIO AFTER TIMESTEP
34007 HONO2 MASS MIXING RATIO AFTER TSTEP
34009 CH4 MASS MIXING RATIO AFTER TSTEP
34010 CO MASS MIXING RATIO AFTER TSTEP
34042 ClO MASS MIXING RATIO AFTER TSTEP
34049 N2O MASS MIXING RATIO AFTER TIMESTEP
34081 OH MASS MIXING RATIO AFTER TIMESTEP
34150 AGE OF AIR IN SECONDS
34153 NO2 MASS MIXING RATIO AFTER TSTEP
34159 OZONE COLUMN
34301 Ox PROD: HO2+NO
34302 Ox PROD: MeOO+NO
34303 Ox PROD: NO+RO2
34304 Ox PROD: OH+INORGANIC ACID
34305 Ox PROD: OH+ORGANIC NITRATE
34306 Ox PROD: ORGANIC NITRATE PHOTOLYSIS
34307 Ox PROD: OH + PAN-TYPE REACTIONS
34311 Ox LOSS: O(1D)+H2O
34312 Ox LOSS: MINOR LOSS REACTIONS
34313 Ox LOSS: HO2+O3
34314 Ox LOSS: OH+O3
34315 Ox LOSS: O3+ALKENE
34316 Ox LOSS: N2O5+H2O
34317 Ox LOSS: NO3 CHEMICAL LOSS
34321 Ox BUDGET: O3 DRY DEPOSITION (3D)
34322 Ox BUDGET: NOy DRY DEPOSITION (3D)
34331 Ox BUDGET: NOy WET DEPOSITION (3D)
34341 RXN FLUX: OH+CH4 (CH4 LIFETIME) TROP
34351 STE: O3
34361 AIR MASS DIAGNOSTIC (TROP ONLY)
34362 TROPOSPHERIC MASK

Start dump

This job is initialised to the file xkjgja.da19991201_00, which has been spun-up for 15 years in the RC6.2 job xkjgj.

TEM Diagnostics

Due to the issues documented in CMS Helpdesk ticket #1357, the TEM diagnostics (STASH m01s30i310-316) are only output to T6H UPA. Otherwise, the STASH is the same as for RC6.0. It is the same on both ARCHER and MONSooN.

Initial conditions and forcing

  • xkawa was initialised from the final dump of xjcin (dated 2015-12-01), then re-dated using change_dump_date to 1999-12-01.
    • See /projects/ukca/inputs/initial/vn84GA4_UKCA.19991201_00.txt
  • SSTs and sea-ice use daily values from the Reynolds dataset, produced by meaning over the transient values from 1995-01-01 to 2005-01-01 using cdo ydaymean.
    • See /projects/ukca/inputs/ancil/surf/reynolds.qrclim.sst.avg2000.txt for SSTs.
    • See /projects/ukca/inputs/ancil/surf/reynolds.qrclim.seaice.avg2000.txt for the sea-ice.
    • Note: there was a minor error with these files, in that the time is set to midnight and not noon. However, this is unlikely to cause any major problems.
  • Time-slice conditions for the year 2000 (actual date 2000-12-01), known as TS2000, were used using the RCP scenario conditions for GHGs and WMO2011 values for ODSs.
    • These are the forcings specified for the CCMI transient REF-C1/2 simulations
    • These values are set in the Spec of trace gases section of the UMUI
      • Note that CFC-114 is only used in UKCA and not in the radiation scheme (as the spectral file does not consider CFC-114). This is done with the hand-edit ~ukca/hand_edits/VN8.4/CFC-114_not_in_Rad.ed.
    • These were calculated using the scenario function on PUMA:
$ /home/ukca/bin/scenario 2000/12/01 USER 
USER-SUPPLIED SCENARIO SELECTED. PLEASE INPUT FILENAME
/home/ukca/tools/scenario/RCP6_MIDYR_CONC_WMO2011_CCMI.DAT
-----------------------------------------------
| 2000/12/01   USER SCENARIO:                 |
-----------------------------------------------
| CFCl3       =   1.23542E-09     CFC11/F11   |
| CF2Cl2      =   2.26638E-09     CFC12/F12   |
| CF2ClCFCl2  =   5.29666E-10     CFC113/F113 |
| CF2ClCF2Cl  =   9.73663E-11     CFC114/F114 |
| CF2ClCF3    =   4.25934E-11     CFC115/F115 |
| CCl4        =   5.19557E-10                 |
| MeCCl3      =   1.95889E-10     CH3CCl3     |
| CHF2Cl      =   4.31321E-10     HCFC22      |
| MeCFCl2     =   5.35084E-11     HCFC141b    |
| MeCF2Cl     =   4.24961E-11     HCFC142b    |
| CF2ClBr     =   2.34011E-11     H1211       |
| CF2Br2      =   2.74854E-13     H1202       |
| CF3Br       =   1.46645E-11     H1301       |
| CF2BrCF2Br  =   4.48740E-12     H2402       |
| MeCl        =   9.58777E-10     CH3Cl       |
| MeBr        =   2.80100E-11     CH3Br       |
| CH2Br2      =   1.80186E-11                 |
| N2O         =   4.80116E-07                 |
| CH4         =   9.67017E-07                 |
| CF3CHF2     =   6.18271E-12     HFC125      |
| CH2FCF3     =   5.43855E-11     HFC134a     |
| H2          =   3.45280E-08                 |
| N2          =   7.54682E-01                 |
| CO2         =   5.61246E-04                 |
-----------------------------------------------
UM/UKCA LBC MMRs for: 2000/12/01, using the USER scenario
VALUES FOR USE IN THE UMUI (ZERO VALUES CAN BE TREATED AS "Excluded"):
 CH4         =  9.670E-07
 N2O         =  4.801E-07
 CFC11       =  1.235E-09
 CFC12       =  2.266E-09
 CFC113      =  5.297E-10
 HCFC22      =  4.313E-10
 HFC125      =  6.183E-12
 HFC134a     =  5.439E-11
 CO2         =  5.61246E-04
VALUES FOR USE IN THE UKCA HAND-EDIT:
 MeBrMMR=2.80100E-11,
 MeClMMR=9.58777E-10,
 CH2Br2MMR=1.80186E-11,
 H2MMR=3.45280E-08,
 N2MMR=0.75468    ,
 CFC114MMR=9.73663E-11,
 CFC115MMR=4.25934E-11,
 CCl4MMR=5.19557E-10,
 MeCCl3MMR=1.95889E-10,
 HCFC141bMMR=5.35084E-11,
 HCFC142bMMR=4.24961E-11,
 H1211MMR=2.34011E-11,
 H1202MMR=2.74854E-13,
 H1301MMR=1.46645E-11,
 H2402MMR=4.48740E-12,
  • The following sources are used for the chemistry and aerosol emissions:
    • The 2D Sulphur-Cycle Emissions are the year 2000 values extracted from the standard CMIP5 dataset used for CLASSIC, which can be found at
      • /projects/um1/ancil/atmos/n96/classic_aerosol/cmip5/1970_2010/v0/qrclim.sulpsurf.
    • Year 2000 AR5 emissions are used for NO (s0i301), CO (s0i303), HCHO (s0i304), C2H6 (s0i305), C3H8 (s0i306), (CH3)2CO (s0i307; Me2CO), CH3CHO (s0i308) MeCHO), Black Carbon (BC) fossil fuel surface emissions (s0i310), Black Carbon (BC) biofuel surface emissions (s0i311), Organic Carbon (OC) fossil fuel surface emissions (s0i312), Organic Carbon (OC) biofuel surface emissions (s0i313), and 3D NO aircraft emissions (s0i340).
    • Year 2000 GEIA emissions are used for C5H8 (s0i309; isoprene) and Monoterpenes (s0i314).
    • Year 2005 MEGAN emissions are used for CH3OH, labelled as NVOC in the STASHmaster file (s0i315; MeOH).
    • Year 2000 GFED2 emissions are used for 3D Black Carbon (s0i322; BC) and 3D Organic Carbon (s0i323; OC).

Making a transient run from a timeslice run

As noted above, this job is configured to run as a year 2000 timeslice. It is relatively straight-forward to change the emissions & forcings to make this job to run as a timeslice at a different time (e.g. pre-industrial or 2100 etc.) or as a transient run, although you will need to change the initial conditions of the long-lived chemical tracers as well as the SSTs etc.

You will need to:

  • Obtain or create a set of SST and sea-ice ancillaries for the time-period you are interested in.
  • Obtain or create a set of emissions ancillaries for the time-period you are interested in.

It is possible to use a start-dump from a timeslice job to initialise a transient run. However, care must be taken over the values for the long-lived chemical tracers, which will have lower boundary conditions set for their surface concentrations. As described above, the UMUI is used to set the values for these gases. It is possible to use the UKCA routine ukca_rcp_scenario to read the RCP forcing files produced for CMIP5, and an example hand-edit is

$ cat /home/ukca/hand_edits/VN8.4/UKCA_RCP6.0.ed
# UKCA_RCP6.0
# (vn8.5)
# control variables to tell the model to use the RCP6.0 scenario, and
# where the data for this is located.
#
# NOTE: FOR THIS HAND-EDIT TO WORK, THE OPTION 'OVERRIDE DEFAULTS' MUST 
#       BE SET IN THE UMUI 

ed CNTLATM <<\EOF1
/L_UKCA_USEUMUIVALS/
d
i
 I_UKCA_SCENARIO=2,
 UKCA_RCPdir='/projects/ukca/nlabra/scenario'
 UKCA_RCPfile='RCP6_MIDYR_CONC.DAT'
.
wq
EOF1

However, this only affects the values in UKCA - the values seen by the radiation will not be affected. You must therefore either make equivalent changes to the UMUI panels which specify these trace gases, or you will need to include additional Fortran code in atmos_physics1.F90 (please contact Luke Abraham for more information on this option).

As well as changing how the lower boundary condition for these gases is specified, it is advisable to re-scale the gases according to their new surface concentrations. As the CheST chemistry scheme lumps the Cl contribution into CFC11 and CFC12, and the Br contribution into CH3Br (MeBr), it is not easy to use the scenario routine described above to calculate these concentrations. However, there is a lumped scenario function, which can be found at ~ukca/bin/scenario.lumped which will produce these values, and is used in the same way, e.g.

$ ~ukca/bin/scenario.lumped 2000/12/01 USER
USER-SUPPLIED SCENARIO SELECTED. PLEASE INPUT FILENAME
/home/ukca/tools/scenario/RCP6_MIDYR_CONC_WMO2011_CCMI.DAT
-------------------------------------------------------
| 2000/12/01   USER SCENARIO (LUMPED MMR):            |
-------------------------------------------------------
| CFCl3       =   2.99423E-09     CFC11/F11 (LUMPED)  |
| CF2Cl2      =   3.16674E-09     CFC12/F12 (LUMPED)  |
| MeBr        =   7.07330E-11     CH3Br     (LUMPED)  |
| N2O         =   4.80116E-07                         |
| CH4         =   9.67017E-07                         |
| H2          =   3.45280E-08                         |
| COS         =   1.07640E-09                         |
-------------------------------------------------------

Note: the value for COS has been updated.

You should run this routine for the date you wish to start on, and compare the values to those above (or, if applying this method to a different timeslice run, to the values generated by the scenario.lumped program for the date of that timeslice) - any differences will mean that the fields will need to be rescaled (by default, H2 and COS are actually constant). This can be easily done using the climate data operators (by first using Xconv to extract these to netCDF from the dump you wish to use). You will also need to subtract the CFC11/CFC12 and CH3Br contributions from the LUMPED Cl (s34i100) and LUMPED Br (s34i099) tracers. To do this, do:

  1. Extract your CFC11 (s34i055), CFC12 (s34i056), CH3Br (s34i057), LUMPED Br (s34i099; assumed to be as BrO), and LUMPED Cl (s34i100; assumed to be as HCl) tracers from the dump file.
  2. Rescale your CFC11, CFC12, and CH3Br fields to the correct values as calculated from the ~ukca/bin/scenario.lumped program (i.e. multiply by the factor new/original).
  3. Extract the CH4 (s34i009) and N2O (s34i049) fields from the dump file and rescale these in the same way.
  4. Convert all Cl and Br fields to vmr. The required conversion factors can be found in the src/atmosphere/UKCA/ukca_constants.F90. You should divide the tracer concentration by the c_species value for that species. Note that you should use the values for BrO for LUMPED Br and HCl for LUMPED Cl.
  5. Subtract the original CFC11 and CFC12 fields (in vmr) from the LUMPED Cl vmr field, and then add-in the new CFC11 and CFC12 fields (also in vmr).
  6. Subtract the original CH3Br field (in vmr) from the LUMPED Br vmr field, and then add-in the new CH3Br field (also in vmr).
  7. Convert your new LUMPED Cl to mmr (using the conversion factor for HCl, as before).
  8. Convert your new LUMPED Br to mmr (using the conversion factor for BrO, as before).
  9. Use Xancil to convert these netCDF fields to ancillary file format (use the generalised ancillary file option).
  10. In the UMUI, turn on reconfiguration, and in the Initialisation of user prognostics panel use option 7 to point to the ancillary file containing these fields.

Note: once you have made these changes, you should spin the model up again. We suggest at least a 10-year spin-up. This can either be done as another timeslice, or as a transient run started earlier than needed, e.g. for a 1960-2010 simulation (starting at 1959-12-01) you could either

  • Start the run in 1949-12-01 and run for 10 years to produce the 1959-12-01 start dump.
  • Make a TS1960 timeslice (with forcing values set to e.g. 1959-12-01) and run this for 10-years to get a 1959-12-01 start dump. In this case you could also start the model from 1949-12-01, but would have the SSTs, sea-ice, and emissions etc. fixed at 1959 values.

Results

MOOSE on MONSooN

Further information on how to use MOOSE can be found on the collaboration twiki.

Note: This job is configured to automatically archive data to MOOSE when running on MONSooN. This setting should not be changed, except for very short runs. If you take a copy of this job and run it, you must first manually make the MOOSE set to hold the data in the archive. This is done by

 moo mkset --project-owner=project-YOUR_MONSooN_PROJECT -v moose:/crum/jobid

For instance, if you were in the UKCA project you would have --project-owner=project-ukca. This can be done after the job has started running. The archiving intelligently knows which files need archiving through the use of files named archive_XXXXXX.do. The files also control the deleting of files and dumps once they have been archived or are no longer needed. If, for some reason, the files cannot be archived (e.g. MOOSE is down or the set has not yet been made) then the files will not be deleted. They will continue being generated and existing on the /projects disk until they can be archived.

You should not need to do anything with the fieldsfiles until the whole simulation has been completed (in this case, the whole 15-years). When it does you will find that, while the climate mean files and dumps have been archived, the last files in the e.g. *.pa* or *.pb* streams etc will not have been. You will need to archive these manually by, e.g.

moo put -f -vv -c=umpp jobida.pzYYYYmmm moose:/crum/jobid/apz.pp/jobida.pzYYYYmmm.pp

Where the -c=umpp converts the files from 64-bit fieldsfiles to 32-bit pp-files. Remember to put the .pp at the end of the name of the file in the set on MOOSE.

Evaluation Suite Output (MONSooN)

These were generated from the UKCA Evaluation Suite available on the MONSooN post-processor.

Lightning NOx

The lightning NOx produced in these years was

  • ARCHER = 3.98946 Tg(N)/year
  • MONSooN = 3.93295 Tg(N)/year

Additional Results from RC6.2

Please see the RC6.2 page for results for that job showing Aerosol Optical Depth and its Aerosol Evaluation.