Inconsistency gravfft between 6.3 and 6.4

Q&A
1

Hi there,
I used gravfft with 6.3 version trying to get the computation of an airy isostatic response from a topography with :
GMT gravfft topo.grd -Gairy.grd -T0/2700/3300/1035+m -Z30000
And it seems to work fine (even though I still struggle with the real impact of the density…)
But when I use version 6.4, the same command has a segmentation fault due to the +m option. When I remove it, the result is totally different, with a short wavelength related to the topography. I think it is coming from the adding effect of topography and isostatic root ?
Can you please explain me ? And how to isolate the signal coming only from the isostatic root ?
Thanks

2

Man, we are short on personnel to address all incoming issues. Could this be the same as reported in gravfft problems computing moho geopotential · Issue #7255 · GenericMappingTools/gmt · GitHub ?

First, please use only last version when reporting new problems. If 6.5 still crashes please open an issue with all necessary data to reproduce it.

3

Hi, yes this is the same: segmentation fault + difference in the results. Thanks. I confirm that -Q option solves the segmentation (which is systematic for me with 6.5). The difference in results, I can’t find why. Still looking.