LFRic Orography

This example demonstrates how to render a warped unstructured cubed-sphere mesh.

📋 Summary

Uses an unstructured Met Office LFRic C48 cubed-sphere mesh of surface altitude data.

The mesh contains quad cells and is constructed from CF UGRID unstructured cell points and connectivity.

Note that the scalar elevation values are located on the mesh nodes/points which results in the rendered colours being interpolated across the cell faces. A pyvista “warp” operation extrudes the mesh, using the same node altitude values, to highlight the global surface topography.

The warp uses compute_normals() and warp_by_scalar(). See Computing Surface Normals for further details.

Tags: domain: orography filter: warp sample: unstructured

---

Static Scene

Interactive Scene

from __future__ import annotations

import geovista as gv
from geovista.pantry.meshes import lfric_orog
import geovista.theme


def main() -> None:
    """Plot a warped LFRic unstructured mesh.

    Notes
    -----
    .. versionadded:: 0.1.0

    """
    # Load the sample mesh.
    mesh = lfric_orog()

    # Warp the mesh nodes by the surface altitude.
    mesh.compute_normals(cell_normals=False, point_normals=True, inplace=True)
    mesh.warp_by_scalar(inplace=True, factor=2e-5)

    # Plot the unstructured mesh.
    p = gv.GeoPlotter()
    sargs = {"title": "Surface Altitude / m", "fmt": "%.0f"}
    p.add_mesh(mesh, scalar_bar_args=sargs)
    p.add_axes()
    p.add_text(
        "LFRic C48 Unstructured Cube-Sphere",
        position="upper_left",
        font_size=10,
    )
    p.camera.zoom(1.3)
    p.show()


if __name__ == "__main__":
    main()