Closed
Description
Hi!
I am trying to use pygimly to perform a gravimetric inversion. The gravimetric profile i want to invert has a length of about 35 kilometers, but I only need the mesh to be a few hundred meters deep, and this has brought me issues generating an adequate mesh. Here you can see my code, I also attach file with the boundaries of the mesh, and an image showing the mesh I have been able to generate so far. I would really appreciate if someone has any recomendation on how I can improve this.
import numpy as np
import pygimli as pg
import pygimli.meshtools as mt
mesh_bounds = np.loadtxt(r'./Files/mesh_bounds.txt', delimiter=',')
mesh_bounds = mesh_bounds.tolist()
topo = mt.createPolygon(mesh_bounds, isClosed=True, marker=1, area=1)
mesh = mt.createMesh(topo, quality=32)
ax, _ = pg.show(mesh)
ax.set_aspect(30)
Thank you very much!
Metadata
Metadata
Assignees
Labels
No labels