Incorrect Map Projection for Regional MPAS Mesh #328
Description
Describe the bug
As NCAR is transitioning away from developing WRF toward further developing MPAS, the user community will begin using the relatively recently-added limited-area capability of MPAS more to produce regional MPAS simulations.
A recent regional MPAS simulation was completed and we discovered that VERDI does not correctly display the regional map when plotting the grid mesh (or any of the output data from the simulation).
To Reproduce
This test was conducted using the 20231013 build of VERDI 2.1.4 under Linux:
- Open VERDI.
- Load the regional MPAS mesh file from the Atmos system:
/work/MOD3DEV/jherwehe/mpas/meshes/92-25km/conus_25km.grid.nc - Doubleclick on variable "meshStructure (unitless)[1]" to add it to Formulas.
- Click "Tile Plot" to see the incorrect map (however, based on the scale range, the data are being read correctly).
- The user cannot zoom in or out of the resulting incorrect regional map.
- And cell borders cannot be shown.
Expected behavior
The expected plot should initially be a colored contour plot illustrating the different MPAS cell sizes (nominally 25 km) centered over CONUS. The "Show Cell Borders" under "Controls" could then be selected to explicitly show the cell edges.
Screenshots
Here is the incorrect map for the regional MPAS mesh:
Operating system:
- Red Hat Enterprise Linux release 8.9 (Ootpa)
Additional context
VERDI must be expecting only global data from MPAS. So in some ways, this issue is about adding a new feature to VERDI.
Tried uploading here a zip of the test regional MPAS mesh from the path above, but it was too large (33MB).