SPECFEM3D_Globe with Cubit/Trelis

Good afternoon,

I am a fairly new user to the SPECFEM code family (and this site overall).

Under the Cartesian Package, I know we are able to input a mesh we generate with CUBIT/TRELIS through code included like run_cubit2specfem3d.py

However, is it possible (or has it already been attempted) to port a mesh generated in Cubit to SPECFEM3D_GLOBE? I don’t see the associated code or examples in the GLOBE package compared with Cartesian and am curious if there is some aspect that prevents this transfer.

Thank you for any advice or assistance.

well, the global version uses the in-house mesher based on a cubed-sphere mesh, which is perfectly balanced between different partitions, and works extensively with geographic/geocentric coordinate transformations. it deals with internal discontinuities and geochemical layering as well as higher sampling for the D’’ region. this mesher is meant to construct a planet.

the Cartesian version allows for much more mesh flexibility. mesh slices with different number of elements are common to load-balance between acoustic/elastic/poroelastic domains. coordinate transformations if needed are done between UTM and Cartesian xyz values. CUBIT/Trelis meshing is only the first step, before partitioning and creating a GLL mesh. this mesher ignores gravitational and rotational effects on waves, thus there is no need to account for radial directions and geographic position.

so in order to port from CUBIT/Trelis to the globe version, it would involve quite some work. it hasn’t been attempted yet - at least not to my knowledge. but, nothing is impossible…

Thank you for the response. That makes sense according to the way you explained it. I would still be interested in eventually trying it, but that may be a project for a different time.