Hi all,
I apologise if this isn’t the right place to post this.
I am trying to estimate the elastic thickness around Elysium Mons on Mars using the gravitational admittance method. However, my estimates, especially my observed admittance, look very off and don’t match the pattern i’d expect it to (from literature at least). I’m wondering if I have missed something in my code or if there is something wrong with my data.
Any help/tips would be greatly appreciated, thanks!
code:
gmt set PROJ_ELLIPSOID=Mars
#resample topography so that it is at 0.05 degree increments to match FAA
gmt grdconvert mola32.nc -Gtopo1.grd -R-180/180/-90/90
gmt grdsample topo1.grd -R-180/180/-90/90 -I0.05/0.05 -Gmola32.grd#trim the area to surround elysium mons
gmt grdcut faa_test.grd -R143/151/21/29 -Gfaa_cut.grd
gmt grdcut mola32.grd -R143/151/21/29 -Gtopo_cut.grdgmt gravfft topo_cut.grd faa_cut.grd -Iwkt -T60000/2900/3500/0.01 -Z48000 -N+a -Gadmittance_60km.txt