1#!/usr/bin/env bash
2# Testing gmt grdcut -Zmin/max
3
4ps=cut.ps
5
6gmt grdmath -R0/10/0/10 -I1 X Y MUL = t.nc
7cat << EOF > replace.txt
80 1 NaN
90 2 NaN
100 3 NaN
113 10 NaN
124 10 NaN
138 9 NaN
148 10 NaN
159 10 NaN
1610 10 NaN
177 4 NaN
18EOF
19gmt grdedit t.nc -Nreplace.txt
20gmt makecpt -Cjet -T0/100/5 > t.cpt
21gmt grdcut t.nc -Z20/70 -Gnew.nc
22gmt psscale -Ct.cpt -Dx0/4.2i+w8.4i/0.2i+jML+m -Ba10 -P -K > $ps
23gmt grdimage new.nc -Ct.cpt -JX4i -R -Baf -O -K -X1i >> $ps
24echo 10 5 -Z20/70 | gmt pstext -R -J -O -K -N -Dj0.5i/0 -F+f24p+jLM >> $ps
25gmt grdcut t.nc -Zr20/70 -Gnew.nc
26gmt grdimage new.nc -Ct.cpt -JX4i -R -Baf -BWsNe+t"Find region with 20 < z < 70" -O -K -Y4.4i >> $ps
27echo 10 5 -Zr20/70 | gmt pstext -R -J -O -N -Dj0.5i/0 -F+f24p+jLM >> $ps
28