MPAS-A tool to create custom regional mesh

B

bgrasic

New member
Hello everyone!

I am familiar with creating regional meshes from existing meshes available on the "MPAS-Atmosphere Meshes" website.

I wonder if there exists an open source tool for creating a custom regional mesh for MPAS-A model or maybe a tutorial?

Is there a thread on this forum that I missed?

I'm mainly interested in whether it would be possible to create a mesh based on a set of points given in WGS coordinates?

Regards, Boštjan
 
Hello @bgrasic, this depends on what level of customization you're looking for. The most supported workflow is to cut out a smaller area from some existing (almost always global) mesh. In general this is documented in Section 3 of the MPAS Tutorial which you can access at MPAS Tutorial — Practice Session Guide.

General procedure:
  1. Download a (global) grid for whatever resolution you want (including variable resolution) from mpas-dev.github.io > MPAS-Atmosphere Download > MPAS-Atmosphere meshes.
  2. Sub-set the mesh with the MPAS Limited-Area tool from MPAS-Dev/MPAS-Limited-Area GitHub repo. The README explains the process, it supports a few specifications including using lat/lon points to define your new mesh's shape.
  3. If needed, use the grid_rotate tool to center your new grid over the region of interest. Tool comes from MPAS-Dev/MPAS-Tools GitHub repo / mesh_tools / grid_rotate.
Other tools may exist, but this is what I've used as a MPAS developer.
 
Thank you for your answer. I am already familiar with proposed procedure and I am using it already.
I had in mind more customized level, so I guess I will wait for a while until some other tools appear.
 
Is there any success @gdicker and @bgrasic creating a new mesh with the Jigsaw tools? I created a new mesh successfully, subsetted successfully with MPAS Limited area, but encountered a block with the init_atmosphere trying to generate the static file.
 
Not yet, but I also managed to make some minor progress. I recently managed to create a mesh with jigsaw, MpasMeshConverter.x and MPAS-Limited.Area. Now I am dealing with the problem of init_atmosphere crashing when creating a static NetCDF.

If anyone can take a look at my mesh or suggest how to check what else is wrong with the mesh, I am happy to share my mesh.
 
bgrasic said:
Not yet, but I also managed to make some minor progress. I recently managed to create a mesh with jigsaw, MpasMeshConverter.x and MPAS-Limited.Area. Now I am dealing with the problem of init_atmosphere crashing when creating a static NetCDF.

If anyone can take a look at my mesh or suggest how to check what else is wrong with the mesh, I am happy to share my mesh.
Click to expand...
By chance, would it happen to be segmentation faulting/segmentation aborting? I'm stuck at generating static NetCDF grid as well (just recently made a post for this as well here) and I've been having seg faults and seg aborts after I used those same tools (a Python script that uses MPAS-Tools which uses MpasMeshConverter and Jigsaw) then used MPAS-Limited-Area to try and take a subset of that area out and while visualizing it on Python shows it worked the init_atmosphere wont get past the static file generation part
 
Yes, it is crashing with segmentation faulting/segmentation aborting after entering subroutine init_atm_static.

And there is not output to .err file, this is the tail of .log file:

WARNING: Variable nominalMinDc not in input file.
WARNING: File Krsko3kmR500km.region.nc does not contain a seekable xtime variable. Forcing a read of the first time record.
real-data GFS test case

--- enter subroutine init_atm_static:
 
bgrasic said:
Yes, it is crashing with segmentation faulting/segmentation aborting after entering subroutine init_atm_static.

And there is not output to .err file, this is the tail of .log file:

WARNING: Variable nominalMinDc not in input file.
WARNING: File Krsko3kmR500km.region.nc does not contain a seekable xtime variable. Forcing a read of the first time record.
real-data GFS test case

--- enter subroutine init_atm_static:
Click to expand...
Practically the exact same error I've gotten as well:

Allocating fields ...
485 MB allocated for fields on this task
----- done allocating fields -----

WARNING: Variable nominalMinDc not in input file.
WARNING: File conus.grid.nc does not contain a seekable xtime variable. Forcing a read of the first time record.
real-data GFS test case

--- enter subroutine init_atm_static:
 
Has anyone successfully used "hex_projection" tool from MPAS-Tools to create regional domain?
Is there maybe some working example or tutorial available somewhere?
 
@ bgrasic

Would you please create a new thread for your question? We recommend this because it is easier for other users to search and access the discussion. thanks.
 
You must log in or register to reply here.
Back
Top