#!/bin/sh
#
PROJ="-R-125/-114/31.0/42.0 -Jm1.7"
gmtset COLOR_MODEL hsv
makecpt -Cno_green -T1/8/1 -Z -N > step.cpt
makecpt -Cgray -T-1.4/.7/.1 -Z > grad.cpt
#
#  mask the deep ocean areas of the strain grid.
#
grdmath topo.grd -2000 GT 0. NAN = mask.grd
cp $1_strain.grd tmp.grd
grdmath tmp.grd LOG10 2 MUL NEG 7.0 ADD FLOOR 1 MAX mask.grd OR  = sample.grd
#
#   image the topo and sample grids
#
grdimage grad.grd $PROJ -Cgrad.cpt -X1.5 -Y1.5 -K -V >sample_rate.ps
grdimage sample.grd -Igrad.grd $PROJ -Cstep.cpt -Q -O -K -V >>sample_rate.ps
pscoast $PROJ -W6 -Di -A200 -B2WSen:."downsample": -O -V -N1 -N2 -K >> sample_rate.ps
psxy fault_data.gmt -m $PROJ -W2 -O -V -K >> sample_rate.ps
psscale -Cstep.cpt -D5/1.5/8/0.3h -L -O  -K >> sample_rate.ps
pstext $PROJ -O <<!>> sample_rate.ps
-124.2 32 12 0 5 1 downsample factor
!