Error in plex.c when generating Green's functions on a large mesh

[willic3]

PETSc exits with error when generating Green's functions on a mesh with ~1.35 M cells

PETSC ERROR: PetscAbortErrorHandler: DMPlexGetTransitiveClosure() at /scratch/charlesw/cig/pylith3/source/petsc/src/dm/impls/plex/plex.c:4085 Point -253701943 is not in [0, 5924403)

To Reproduce

PETSc branch: [knepley/pylith]
PyLith fork: [https://github.com/willic3/pylith.git]
PyLith branch: [willic3/fix-subduction-3d-example]

I run the problem using a very simple .cfg file based on the step07 examples in examples/subduction-3d (based on my current branch).

Expected behavior

I expect the example to work, as similar examples are working.

Screenshots

If applicable, add screenshots to help explain your problem.

Environment

Please provide the following information:

• OS: [Fedora 39]
• PyLith version [Branch used is synced with current main]

Additional context

Same mesh is being used successfully to generate GF using PyLith v 2.2.2 binary.