Generate 1-3-12 km MPAS variable-resolution mesh

[bruceyoung01]

I plan to generate 1-3-12 km MPAS variable-resolution mesh for our MPAS simulation over the Baltimore-Washington area. I was wondering whether it is feasible based on the current capacity of MPAS. Thanks.

 

[Ming Chen]

Hi, bruceyoung01,

The MPAS Limited-Area Tool can be used to subset a variable-resolution regional mesh. The most closest mesh I can find is the 60-12km mesh posted in MPAS home page. You need to:

(1) run grid_rotate to relocate the fine-mesh area to Baltimore-Washington area,

(see grid_rotate at MPAS-Tools/mesh_tools/grid_rotate at master · MPAS-Dev/MPAS-Tools)

(2) run MPAS Limited-Area Tool to subset a regional mesh

(see GitHub - MPAS-Dev/MPAS-Limited-Area: Python tool to create a regional subset of a global MPAS Mesh)

(3) scale the regional mesh to higher resolution using scale_region utility.

(GitHub - mgduda/scale_region: A Python tool for scaling regional MPAS meshes)

This will give you a 12-3km mesh that roughly meet your need. However, this mesh cannot go up to 1-km resolution.

If you really need the 1-km mesh, I am thinking that you may need to create a 12-3-1km mesh yourself. The mesh-generation tool is probably able to do so. Honestly, I haven't used this tool to produce high-resolution nested mesh myself. So I am not 100% sure how it works.

More details can be found at 71_duda.pdf.

Michale gave a talk on how to work with the mesh-generation package. Please see

Hope this is helpful for you.