maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git


  
1

	
2

	
3

	
4

	
5

	
6

	
7

	
8

	
9

	
10

	
11

	
12

	
13

	
14

	
15

	
16

	
17

	
18

	
19

	
            #!/usr/bin/env bash
#
# Compute the gravity anomaly of the Flemish Cap Guyot

ps=fc_okb.ps

lim=-41:50/-41:20/47:30/47:50

# Get rough gmt surface description
gmt grd2xyz @earth_relief_10m_g -R$lim -fg > fc.dat

# Calculate the triangles
gmt triangulate fc.dat > fc_tri.dat

# Compute the grav anomaly using a contras density of 1700 kg/m^3
gmt gmtgravmag3d -C1700 -Gfc_okb.nc -R -I1m -Z-4300 -fg fc.dat -Tvfc_tri.dat

gmt grdcontour fc_okb.nc -C2.5 -A5 -JM14c -Ba -BWSen -P > $ps