Make MPAS to Antarctic Stereo Mapping

This example demonstrates how to create a mapping file from an MPAS mesh to an Antarctic stereographic grid.

Script Overview

The script uses pyremap to:

  • Read an MPAS mesh file.

  • Define a target Antarctic stereographic grid.

  • Generate a mapping file for regridding.

Example Code

#!/usr/bin/env python

"""
Creates a mapping file to remap MPAS files to an Antarctic stereographic grid.
"""

from pyremap import Remapper, get_polar_descriptor

# Input MPAS mesh file and name
src_mesh_name = 'oQU240'
src_mesh_filename = 'ocean.QU.240km.151209.nc'

# Configure the remapper
remapper = Remapper(ntasks=4, method='bilinear')
remapper.src_from_mpas(filename=src_mesh_filename, mesh_name=src_mesh_name)

# Define the target Antarctic stereographic grid descriptor
remapper.dst_descriptor = get_polar_descriptor(
    lx=6000.,  # Grid width in km
    ly=5000.,  # Grid height in km
    dx=10.,    # Grid cell size in x-direction (km)
    dy=10.,    # Grid cell size in y-direction (km)
    projection='antarctic'  # Antarctic stereographic projection
)

# Build the mapping file
remapper.build_map()

Expected Output

The script generates a mapping file named:

map_oQU240_to_6000.0x5000.0km_10.0km_Antarctic_stereo_esmfbilin.nc

that can be used for regridding MPAS data to an Antarctic stereographic grid.