search

MPAS-Dev / MPAS-Tools

MPAS Tools Repository
Other
39 stars 65 forks source link

Error on vertex mask when creating transect masks with new tool #413

Closed milenaveneziani closed 3 years ago

milenaveneziani commented 3 years ago

I started using compute_mpas_region_masks and compute_mpas_transect_masks today, with my most recent mpas-analysis conda environment on anvil, and while the region_mask tool seems to work fine, I get the following error when running the transect_masks tool:

  Computing vertex masks:
Traceback (most recent call last):
  File "/home/ac.milena/.conda/envs/mpas-analysis/bin/compute_mpas_transect_masks", line 11, in <module>
    sys.exit(entry_point_compute_mpas_transect_masks())
  File "/home/ac.milena/.conda/envs/mpas-analysis/lib/python3.8/site-packages/mpas_tools/mesh/mask.py", line 370, in entry_point_compute_mpas_transect_masks
    dsMasks = compute_mpas_transect_masks(
  File "/home/ac.milena/.conda/envs/mpas-analysis/lib/python3.8/site-packages/mpas_tools/mesh/mask.py", line 258, in compute_mpas_transect_masks
    _get_polygons(dsMesh, maskType)
  File "/home/ac.milena/.conda/envs/mpas-analysis/lib/python3.8/site-packages/mpas_tools/mesh/mask.py", line 1071, in _get_polygons
    assert numpy.all(vertexIndices >= 0)
AssertionError

I am not sure if this has something to do with my geojson file or something else. Here is the full command I issued: compute_mpas_transect_masks -m EC30to60E2r2mesh -g arcticSections20210323.geojson -o EC30to60E2r2_arcticSections.20210323_newtool.nc -s 10e3 --add_edge_sign in this directory: /lcrc/group/e3sm/ac.milena/mpas-region_masks.

Note that the cell and edge masks completed without error, it only stopped at the vertex mask calculation.