Using specfem3d to derive vertical layering problems


The following error occurs:
error: surface normal, dot=0 [0.0, -1.0, 0.0] (-1, 0, 0) 0.0 (49.0, 0.0, -19.0) (49.0, 0.0, -20.0) (50.0, 0.0, -20.0)
Traceback (most recent call last):
File “1.py”, line 134, in
exportlib.e2SEM(outdir=‘MESH/’)
File “/usr/local/lib/python2.7/dist-packages/geocubitlib/exportlib.py”, line 269, in e2SEM
top_absorbing=top_absorbing, hex27=hex27)
File “/usr/local/lib/python2.7/dist-packages/geocubitlib/cubit2specfem3d.py”, line 1553, in export2SPECFEM3D
sem_mesh.write(path=path_exporting_mesh_SPECFEM3D)
File “/usr/local/lib/python2.7/dist-packages/geocubitlib/cubit2specfem3d.py”, line 1284, in write
self._write_ascii(path=path)
File “/usr/local/lib/python2.7/dist-packages/geocubitlib/cubit2specfem3d.py”, line 1536, in _write_ascii
self.abs_write(path + self.absname)
File “/usr/local/lib/python2.7/dist-packages/geocubitlib/cubit2specfem3d.py”, line 1203, in abs_write
h, f, normal=normal, cknormal=cknormal)
File “/usr/local/lib/python2.7/dist-packages/geocubitlib/cubit2specfem3d.py”, line 874, in create_facenode_string
txt = str(hexa) + ’ ’ + ’ '.join(str(x) for x in nodes_ok)
TypeError: ‘NoneType’ object is not iterable

Seek Assist!!!

Hi @Xushuai,
we don’t support specfem3d.

When i look at specfem3d/CUBIT_GEOCUBIT at master · SPECFEM/specfem3d · GitHub
It seems like the code isn’t actively maintained for newer cubit versions.

I would suggest that you open a new issue Issues · SPECFEM/specfem3d · GitHub and seek help in the specfem community.