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

	
            #!/bin/bash
#	$Id: subset.sh 12115 2013-09-03 23:22:48Z fwobbe $
# Testing gmt grdcut -Zmin/max

ps=cut.ps

gmt grdmath -R0/10/0/10 -I1 X Y MUL = t.nc
cat << EOF > replace.txt
0 1 NaN
0 2 NaN
0 3 NaN
3 10 NaN
4 10 NaN
8 9 NaN
8 10 NaN
9 10 NaN
10 10 NaN
7 4 NaN
EOF
gmt grdedit t.nc -Nreplace.txt
gmt makecpt -Cjet -T0/100/5 > t.cpt
gmt grdcut t.nc -Z20/70 -Gnew.nc
gmt psscale -Ct.cpt -A -D0/4.2i/8.4i/0.2i -Ba10 -P -K > $ps
gmt grdimage new.nc -Ct.cpt -JX4i -R -Baf -O -K -X1i >> $ps
echo 10 5 -Z20/70 | gmt pstext -R -J -O -K -N -Dj0.5i/0 -F+f24p+jLM >> $ps
gmt grdcut t.nc -Zr20/70 -Gnew.nc
gmt grdimage new.nc -Ct.cpt -JX4i -R -Baf -BWsNe+t"Find region with 20 < z < 70" -O -K -Y4.4i >> $ps
echo 10 5 -Zr20/70 | gmt pstext -R -J -O -N -Dj0.5i/0 -F+f24p+jLM >> $ps