Reference: GMT Chinese community
Here I plot the topography and add the locations of the seismometers on it.
1
2
3
4
5
6
7
8
9
10
11
12
| R=76.1167/76.6667/-69.6167/-69.35/-100/1000
gmt begin elevation_3d pdf
gmt grdtrack -R76.1167/76.6667/-69.6167/-69.35 -GREMA.nc staloc.txt > staloc.xyz # generate the corresponding elevations of the stations by track locations on topography file REMA.nc
gmt makecpt -Cgray -T0/900
gmt grdview REMA.nc -R$R -JM10c -JZ3c -N-100+ggray -C -Qi300 -I -Ba -BwsENZ -p30/25
cat staloc.xyz | awk '{print $1,$2,$5}' | gmt plot3d -R$R -JM10c -JZ3c -St0.05c -W1p,white -Gwhite -p
gmt colorbar -C -Ba -DJTC+o0/1c -p
gmt basemap -TdjRB+w3c+l+o0.5c/0.5c -p30/25/700 # plot the North arrow
gmt end show
rm staloc.xyz
|