#!/bin/sh
#
PROJ="-R-125/-114/31.0/42.0 -Jm1.7"
#PROJ="-R-123/-114/31.5/37.5 -Jm1.7"
gmtset COLOR_MODEL rgb
makecpt -Cno_green -Qo -T1/1000/1 -D -Z > strain.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
mv $1_strain.grd tmp.grd
grdmath tmp.grd mask.grd OR = $1_strain.grd
#
#   image the topo and strain grids
#
grdimage grad.grd $PROJ -Cgrad.cpt -X1.5 -Y1.5 -K -V >$1.ps
grdimage $1_strain.grd -Igrad.grd $PROJ -Cstrain.cpt -Q -O -K -V >>$1.ps
pscoast $PROJ -W2 -Di -A200 -B2WSen:."$1": -O -V -N1 -N2 -K >> $1.ps
psxy fault_data.gmt  $PROJ -W.5 -O -V -K >> $1.ps
#psxy ../profiles/profiles.gmt $PROJ -W2p,white  -O -V -K >> $1.ps
psscale -Cstrain.cpt -D5/1.5/8/0.3h -L -O  -K >> $1.ps
pstext $PROJ -O <<!>> $1.ps
-124.2 32 12 0 5 1 2nd invariant (nanostrain/yr)
!