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


  
1

	
2

	
3

	
4

	
5

	
6

	
7

	
8

	
9

	
10

	
11

	
12

	
13

	
            #!/usr/bin/env bash
ps=filtersample.ps

if [[ ${HAVE_GLIB_GTHREAD} =~ TRUE|ON ]]; then
  _thread_opt=-x+a
fi

echo "0.5 0.5 1" | gmt xyz2grd -di0 -I0.1 -An -Gtemp.grd -R0/1/0/1

gmt grdfilter temp.grd -Gtemp2.grd -R0/1/0/1 -I0.01 -D0 -Fg0.5 ${_thread_opt}
gmt grd2cpt temp2.grd -E100 -D -Cjet > sn.cpt
gmt grdimage temp2.grd -Baf -BWSne -P -JX4i -K -Csn.cpt -Xc  > $ps
gmt grdimage temp.grd  -Baf -BWSne -J -O -Csn.cpt -Y5i >> $ps