Cubit mesh for crustal-strikeslip-3d example

PyLith
1

I’m attempting to run the crustal-strikeslip-3d example on pylith 4.2.1 using the binary. I can run the gmsh version without problems, but I get an error saying there is an inverted cell in the mesh when I run the CUBIT version (full error code attached). I created the CUBIT mesh using the provided generate_cubit.py file, just had to update ids_from_entities to just: return " ".join(str(e.id()) for e in entities), because cubit.get_id_string doesn’t exist in my CUBIT version (Cubit 2021.5, I am worried that my version might just be rather old). No errors generating the CUBIT mesh, and I double-checked the surface and curve signatures, and they seem correct.

crustal_3d_cubit_error_code.txt (6.2 KB)

2

The error message indicates that the Jacobian of a cell is 0, so the cell is perfectly flat (zero volume). My guess is that there are collapsed cells as a result of differences in the ids of geometric or mesh entities or differences in Cubit. We can only support versions of Cubit that we have access to, so we generally support the current release.