Following previous work on an inherently mass‐conserving semi‐implicit (SI) semi‐Lagrangian (SL) discretization of the two‐dimensional (2D) shallow‐water equations and 2D vertical slice equations, that approach is here extended to the 3D deep‐atmosphere, non‐hydrostatic global equations. As with the reduced‐dimension versions of this model, an advantage of the approach is that it preserves the same basic structure as a standard, non‐mass‐conserving, SISL version of the model. Additionally, the model is simply switchable to hydrostatic and/or shallow‐atmosphere forms. It is also designed to allow simple switching between various geometries (Cartesian, spherical, spheroidal). The resulting mass‐conserving model is applied to a standard set of test problems for such models in spherical geometry and compared with results from the standard SISL version of the model.

We describe Global Atmosphere 6.0 and Global Land 6.0 (GA6.0/GL6.0): the latest science configurations of the Met Office Unified Model and JULES (Joint UK Land Environment Simulator) land surface model developed for use across all timescales. Global Atmosphere 6.0 includes the ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment) dynamical core, which significantly increases mid-latitude variability improving a known model bias. Alongside developments of the model's physical parametrisations, ENDGame also increases variability in the tropics, which leads to an improved representation of tropical cyclones and other tropical phenomena. Further developments of the atmospheric and land surface parametrisations improve other aspects of model performance, including the forecasting of surface weather phenomena. We also describe GA6.1/GL6.1, which includes a small number of long-standing differences from our main trunk configurations that we continue to require for operational global weather prediction. Since July 2014, GA6.1/GL6.1 has been used by the Met Office for operational global numerical weather prediction, whilst GA6.0/GL6.0 was implemented in its remaining global prediction systems over the following year.

Numerical weather, climate, or Earth system models involve the coupling of components. At a broad level, these components can be classified as the resolved fluid dynamics, unresolved fluid dynamical aspects (i.e., those represented by physical parameterizations such as subgrid-scale mixing), and nonfluid dynamical aspects such as radiation and microphysical processes. Typically, each component is developed, at least initially, independently. Once development is mature, the components are coupled to deliver a model of the required complexity. The implementation of the coupling can have a significant impact on the model. As the error associated with each component decreases, the errors introduced by the coupling will eventually dominate. Hence, any improvement in one of the components is unlikely to improve the performance of the overall system. The challenges associated with combining the components to create a coherent model are here termed physics–dynamics coupling. The issue goes beyond the coupling between the parameterizations and the resolved fluid dynamics. This paper highlights recent progress and some of the current challenges. It focuses on three objectives: to illustrate the phenomenology of the coupling problem with references to examples in the literature, to show how the problem can be analyzed, and to create awareness of the issue across the disciplines and specializations. The topics addressed are different ways of advancing full models in time, approaches to understanding the role of the coupling and evaluation of approaches, coupling ocean and atmosphere models, thermodynamic compatibility between model components, and emerging issues such as those that arise as model resolutions increase and/or models use variable resolutions.