diff --git a/src/dm/impls/plex/plexmed.c b/src/dm/impls/plex/plexmed.c index ad7f90c..7eaa39e 100644 --- a/src/dm/impls/plex/plexmed.c +++ b/src/dm/impls/plex/plexmed.c @@ -78,7 +78,7 @@ PetscErrorCode DMPlexCreateMedFromFile(MPI_Comm comm, const char filename[], Pet ierr = PetscLayoutSetUp(vLayout);CHKERRQ(ierr); ierr = PetscLayoutGetRanges(vLayout, &vrange);CHKERRQ(ierr); numVerticesLocal = vrange[rank+1]-vrange[rank]; - ierr = MEDfilterBlockOfEntityCr(fileID, numVertices, 1, spaceDim, MED_ALL_CONSTITUENT, MED_FULL_INTERLACE, MED_COMPACT_STMODE, + ierr = MEDfilterBlockOfEntityCr(fileID, numVertices, 1, spaceDim, MED_ALL_CONSTITUENT, MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE, MED_NO_PROFILE, vrange[rank]+1, 1, numVerticesLocal, 1, 1, &vfilter);CHKERRQ(ierr); /* Read mesh coordinates */ if (numVertices < 0) SETERRQ1(comm, PETSC_ERR_ARG_WRONG, "No nodes found in .med mesh file: %s", filename); @@ -106,7 +106,7 @@ PetscErrorCode DMPlexCreateMedFromFile(MPI_Comm comm, const char filename[], Pet ierr = PetscLayoutSetUp(cLayout);CHKERRQ(ierr); ierr = PetscLayoutGetRanges(cLayout, &crange);CHKERRQ(ierr); numCellsLocal = crange[rank+1]-crange[rank]; - ierr = MEDfilterBlockOfEntityCr(fileID, numCells, 1, numCorners, MED_ALL_CONSTITUENT, MED_FULL_INTERLACE, MED_COMPACT_STMODE, + ierr = MEDfilterBlockOfEntityCr(fileID, numCells, 1, numCorners, MED_ALL_CONSTITUENT, MED_FULL_INTERLACE, MED_GLOBAL_PFLMODE, MED_NO_PROFILE, crange[rank]+1, 1, numCellsLocal, 1, 1, &cfilter);CHKERRQ(ierr); /* Read cell connectivity */ if (numCells < 0) SETERRQ1(comm, PETSC_ERR_ARG_WRONG, "No cells found in .med mesh file: %s", filename);