It looks like the seamount example for the documentation of grdflexure gives the densities in the wrong order.
Garrett
Hi Garrett,
Could you please make this more explicit (a one liner command and a bit of the output)?
Hi Joaquim,
The 1st example given (like the subsequent examples) is
“gmt grdflexure smt.nc -Gflex.nc -E10k -D2700/3300/1035”
2700/3300/1035 is the wrong order for " -Drm/rl[/ri]/rw". The correct order should be -D3300/2700/1035 (mantle/crust/water). Using the current order (2700/3300/1035) produces a bogus solution.
Thanks,
Garrett
OK, I see it now. Unfortunately grdflexure -D does not follow the order of gravfft -T and probably that’s what was in Paul’s mind when he wrote those examples but changing that would be a breaking change.
Ah, ok. But can the documentation pages (for all versions) for grdflexure be fixed so that the correct order is given? I was stuck on this issue for some time and almost abandoned using grdflexure until Chatgpt helped me find the problem. Thanks. Garrett
Yes, I’ll fix it but in these cases the best is to open an issue so it won’t get lost/forgotten.