I’d like to create a mesh for a simple model representing a 3-layer geological section using Gmsh, to be imported into the Specfem software. For the interface between the first and second layers, I plan to use points imported from a file. This will allow for variations in elevation, with a total of 31 points defining the interface. The interface between the second and third layers will be identical to the previous one, with a constant value subtracted. I’m using the Python programming language, specifically the ‘pygmsh’ module.