Note
Go to the end to download the full example code.
LFRic Orography#
This example demonstrates how to render a warped unstructured cubed-sphere mesh.
π Summary#
Creates a mesh from 1-D latitude and longitude unstructured cell points.
The resulting mesh contains quad cells and is constructed from CF UGRID unstructured cell points and connectivity.
It uses an unstructured Met Office LFRic C48 cubed-sphere of surface altitude data.
Note that, the data is located on the mesh nodes/points which results in mesh interpolation across the cell faces. The point surface altitudes are used to extrude the mesh to reveal the global surface topography. Also, Natural Earth coastlines are rendered.
from __future__ import annotations
import geovista as gv
from geovista.pantry.data import lfric_orog
import geovista.theme
def main() -> None:
"""Plot a warped LFRic unstructured mesh.
Notes
-----
.. versionadded:: 0.1.0
"""
# Load the sample data.
sample = lfric_orog()
# Create the mesh from the sample data.
mesh = gv.Transform.from_unstructured(
sample.lons,
sample.lats,
connectivity=sample.connectivity,
data=sample.data,
name=sample.name,
)
# Warp the mesh nodes by the surface altitude.
mesh.compute_normals(cell_normals=False, point_normals=True, inplace=True)
mesh.warp_by_scalar(scalars=sample.name, inplace=True, factor=2e-5)
# Plot the unstructured mesh.
plotter = gv.GeoPlotter()
sargs = {"title": f"{sample.name} / {sample.units}", "shadow": True}
plotter.add_mesh(mesh, scalar_bar_args=sargs)
plotter.add_axes()
plotter.add_text(
"LFRic C48 Unstructured Cube-Sphere",
position="upper_left",
font_size=10,
shadow=True,
)
plotter.camera.zoom(1.3)
plotter.show()
if __name__ == "__main__":
main()