ESCOMP / mizuRoute

Reach-based river routing model
http://escomp.github.io/mizuRoute/
GNU General Public License v3.0
42 stars 52 forks source link

HDMAmz_amazon grid fails in reading in mesh file #382

Open ekluzek opened 1 year ago

ekluzek commented 1 year ago

the HDMAmz_amazon grid is failing in reading the mesh file:rof/mizuRoute/meshes/HDMAmz_5x5_amazon_ctrcrd_cdf5_ESMFmesh_c20200731.nc

This is with using mapping files:

LND2ROF_FMAPNAME: /glade/p/cesmdata/cseg/inputdata/rof/mizuRoute/gridmaps/map_5x5_amazon_TO_HDMAmz_5x5_amazon_aave.201028.nc
ROF2LND_FMAPNAME: /glade/p/cesmdata/cseg/inputdata/rof/mizuRoute/gridmaps/map_HDMAmz_5x5_amazon_TO_5x5_amazon_aave.201028.nc

PET file gives:

20230517 124819.229 ERROR            PET71 /glade/p/cesmdata/cseg/PROGS/build/42450/esmf-8.4.0/src/Infrastructure/Mesh/src/ESMCI_Mesh_Glue.C:687 triangulate_warea() Arguments are incompatible  -  - there was a problem (e.g. repeated points, clockwise poly, etc.) with the triangulation of the element with id=4288
20230517 124819.230 ERROR            PET71 /glade/p/cesmdata/cseg/PROGS/build/42450/esmf-8.4.0/src/Infrastructure/Mesh/src/ESMCI_Mesh_Glue.C:1587 ESMCI_meshaddelements() Arguments are incompatible  - Internal subroutine call returned Error
20230517 124819.230 ERROR            PET71 /glade/p/cesmdata/cseg/PROGS/build/42450/esmf-8.4.0/src/Infrastructure/Mesh/src/ESMCI_Mesh_FileIO.C:634 ESMCI_mesh_create_from_ESMFMesh_ Arguments are incompatible  - Internal subroutine call returned Error
20230517 124819.230 ERROR            PET71 /glade/p/cesmdata/cseg/PROGS/build/42450/esmf-8.4.0/src/Infrastructure/Mesh/src/ESMCI_Mesh_FileIO.C:670 ESMCI_mesh_create_from_ESMFMesh_ Arguments are incompatible  - Internal subroutine call returned Error
20230517 124819.230 ERROR            PET71 /glade/p/cesmdata/cseg/PROGS/build/42450/esmf-8.4.0/src/Infrastructure/Mesh/src/ESMCI_Mesh_FileIO.C:284 ESMCI_mesh_create_from_file() Arguments are incompatible  - Internal subroutine call returned Error
20230517 124819.230 ERROR            PET71 ESMCI_MeshCap.C:2594 MeshCap::meshcreatefromfilenew() Arguments are incompatible  - Internal subroutine call returned Error
20230517 124819.230 ERROR            PET71 ESMF_Mesh.F90:1969 ESMF_MeshCreateFromFile() Arguments are incompatible  - Internal subroutine call returned Error
20230517 124819.230 ERROR            PET71 rof_comp_nuopc.F90:590 Arguments are incompatible  - Passing error in return code
20230517 124819.230 ERROR            PET71 ESM0001:src/addon/NUOPC/src/NUOPC_Driver.F90:2725 Arguments are incompatible  - Phase 'IPDv01p3' Initialize for modelComp 4: ROF did not return ESMF_SUCCESS
20230517 124819.230 ERROR            PET71 ESM0001:src/addon/NUOPC/src/NUOPC_Driver.F90:1891 Arguments are incompatible  - Passing error in return code
20230517 124819.230 ERROR            PET71 ensemble:src/addon/NUOPC/src/NUOPC_Driver.F90:2725 Arguments are incompatible  - Phase 'IPDv02p3' Initialize for modelComp 1: ESM0001 did not return ESMF_SUCCESS
20230517 124819.230 ERROR            PET71 ensemble:src/addon/NUOPC/src/NUOPC_Driver.F90:1895 Arguments are incompatible  - Passing error in return code
20230517 124819.230 ERROR            PET71 ensemble:src/addon/NUOPC/src/NUOPC_Driver.F90:481 Arguments are incompatible  - Passing error in return code
20230517 124819.230 ERROR            PET71 esmApp.F90:134 Arguments are incompatible  - Passing error in return code
20230517 124819.230 INFO             PET71 Finalizing ESMF
ekluzek commented 1 year ago

Maybe we need a polygon-centroid mesh file to read in? Is this grid important enough to keep?

nmizukami commented 1 year ago

If i remember, amazon HDMA used to use this unstructured mesh and do online remapping? so not using mapping files? Error seems to indicate some topology error in the mesh.

Also if I remember, this online mapping used to work, but not any more? if so, with some ESMF update (if any), it get pickier?

I guess a polygon centroid mesh with the mapping files could work.

ekluzek commented 1 month ago

This failed in the latest baseline:

/glade/campaign/cgd/tss/ctsm_baselines//mizu_c-cpl.n02_v2.1.1_ctsm5.1.dev159/SMS_D.5x5_amazon_rHDMA.I2000Clm50SpMizGs.derecho_intel.mizuroute-default