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

	
20

	
21

	
22

	
23

	
24

	
25

	
26

	
27

	
28

	
29

	
30

	
31

	
32

	
33

	
34

	
35

	
36

	
37

	
38

	
39

	
40

	
41

	
42

	
43

	
44

	
45

	
46

	
47

	
48

	
49

	
            #!/usr/bin/env bash
# Testing gmt grdcut -S for 3 different points

ps=origin.ps

# Create global grid by evaluating distances to 0,0
gmt set PROJ_ELLIPSOID Sphere
gmt grdmath -Rd -I1 0 0 SDIST KM2DEG = tmp.nc
gmt grdcontour tmp.nc -JQ0/7i -P -K -A10 -C5 -Baf -Y6.5i -Xc > $ps
# Plot three origins and circles
gmt psxy -R -J -O -K -SE << EOF >> $ps
> -W1p,red
30 30 0 3000 3000
> -W1p,green
-170 -60 0 3000 3000
> -W1p,blue
145 78 0 3000 3000
EOF
gmt psxy -R -J -O -K -Sx0.1i -W1p << EOF >> $ps
30 30
-170 -60
145 78
EOF
# 1st point
gmt grdcut tmp.nc -Sn30/30/1500k -Gt.nc
gmt psbasemap -Rt.nc -JM3i -O -K -Baf -Y-4i >> $ps
gmt grd2xyz t.nc -s | gmt psxy -R -J -O -K -Ss0.02i -Gcyan -N >> $ps
gmt grd2xyz t.nc -s+r | gmt psxy -R -J -O -K -Ss0.02i -Gorange -N >> $ps
gmt grdcontour t.nc -J -O -K -A10 -C5 -Gd2i >> $ps
echo 30 30 0 3000 3000 | gmt psxy -R -J -O -K -SE -W1p,red >> $ps
echo 30 30 | gmt psxy -R -J -O -K -Sx0.1i -W1p >> $ps
# 2nd point [set nodes to NaN outside]
gmt grdcut tmp.nc -Sn-170/-60/1500k -Gt.nc
gmt psbasemap -Rt.nc -JM3i -O -K -Baf -X4i >> $ps
gmt grd2xyz t.nc -s | gmt psxy -R -J -O -K -Ss0.02i -Gcyan -N >> $ps
gmt grd2xyz t.nc -s+r | gmt psxy -R -J -O -K -Ss0.02i -Gorange -N >> $ps
gmt grdcontour t.nc -J -O -K -A10 -C5 -Gd2i >> $ps
echo -170 -60 0 3000 3000 | gmt psxy -R -J -O -K -SE -W1p,green >> $ps
echo -170 -60 | gmt psxy -R -J -O -K -Sx0.1i -W1p >> $ps
# 3rd point
gmt grdcut tmp.nc -Sn145/78/1500k -Gt.nc
gmt psbasemap -Rt.nc -JQ0/7i -O -K -Baf -X-4i -Y-1.5i >> $ps
gmt grd2xyz t.nc -s | gmt psxy -R -J -O -K -Ss0.02i -Gcyan -N >> $ps
gmt grd2xyz t.nc -s+r | gmt psxy -R -J -O -K -Ss0.02i -Gorange -N >> $ps
gmt grdcontour t.nc -J -O -K -A10 -C5 -Gd2i >> $ps
echo 145 78 0 3000 3000 | gmt psxy -R -J -O -K -SE -W1p,blue >> $ps
echo 145 78 | gmt psxy -R -J -O -K -Sx0.1i -W1p >> $ps

gmt psxy -R -J -O -T >> $ps