I am trying to calculate a simple model of free air anomaly using the gmt talwani3d by placing cylinders 0.05km high at 0.15km and 0.1km below the surface of a 10km x 10km flat surface. In the description of the 3D model table, I have written the z level and density contrast in the header for the two segments, but how do I specify the height of cylinders?
Can someone provide a suitable sample model description?
One thing that is recommended when examples in docs are not enough is to try look at the tests. I found that these thee tests use talwani3d (I never used it myself)
You should look also to gmtgravmag3d. It must be able to do what you want.
Dear Joaquim,
Thank you very much. These sample scripts are very much useful for me.
Model cylinder should be described a series of circle with different vertical levels. Model cylinder of Cyltest.sh is described with 10 same circle with same depth interval, 375.1 m. I found reducing the mid-level circles did not change the free air anomaly much, only top and bottom circles are crucial.
I would like to check gmt gravmag3d, too.