Physics–dynamics–chemistry coupling across different meshes in LFRic‐Atmosphere: Formulation and idealised tests

The main components of an atmospheric model for numerical weather prediction are the “dynamical core,” which describes the resolved flow, and the “physical parametrisation,” which capture the effects of non‐fluid and non‐resolved fluid processes. Additionally, models used for air quality or climate applications may include a component that represents the evolution of chemicals and aerosols within the atmosphere. Though, traditionally, all these components use the same mesh with the same grid spacing, we present a formulation for the different components to use a series of nested meshes, with different horizontal grid spacings. This gives the model greater flexibility in the allocation of computational resources, so that resolution can be targeted to those parts that provide the greatest benefits in accuracy. The formulation presented here concerns the methods for mapping fields between meshes and is designed for the compatible finite‐element discretisation used by LFRic‐Atmosphere, the Met Office's next‐generation atmosphere model. Key properties of the formulation include the consistent and conservative transport of tracers on a mesh that is coarser than the dynamical core, and the handling of moisture to ensure mass conservation without generation of unphysical negative values. Having presented the formulation, it is then demonstrated through a series of idealised test cases that show the feasibility of this approach.