Sedimentary corrections with Tesseroids

[jmdagert]

Hi Leo, I'm currently working on gravity corrections with Tesseroids using the old version of Fatiando. The topographic correction came out fine using _height=7000*np.ones(topomodel.size), but when I run the same script for the sedimentary correction a keep getting the following error:

If I run the script using a negative height (-7000 mts as in _height=-7000*np.ones(sedimentarymodel.size)), the outcome seems odd. Can you please help me with this issue?

Cheers, Juan Manuel Dagert

[leouieda]

Hi @jmdagert, I'll be a bit more information about what you're running and if possible the full code that generates this error.

Also, there is a new implementation of this in Harmonica: https://www.fatiando.org/harmonica/dev/gallery/forward/tesseroid_layer.html#sphx-glr-gallery-forward-tesseroid-layer-py It's not been officially released yet but there should be a new version with that out in the next month or so. It's a faster and better implemented code that runs on Python 3.

[jmdagert]

Hi Leo,

Thanks for your answer. Find attached a picture of the script I'm trying to run. The input data is actually a .txt that combines sedimentary information downloaded from https://ngdc.noaa.gov/mgg/sedthick/ and isopach polygons I drew from some maps.

Cheers, Juan Manuel Dagert