#!/bin/sh
#
PROJ="-R-125/-114/31.0/42.0 -Jm1.7"
#PROJ="-R-122/-114/31.0/38.0 -Jm2.0"
gmtset COLOR_MODEL rgb
makecpt -Cno_green -Qo -T.0001/1./1 -D -Z > prob.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.grd tmp.grd
grdmath tmp.grd mask.grd OR = $1.grd
#
#   image the topo and strain grids
#
grdimage grad.grd $PROJ -Cgrad.cpt -X1.5 -Y1.5 -K -V >$1.ps
grdimage $1.grd -Igrad.grd $PROJ -Cprob.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 -m $PROJ -W.5 -O -V -K >> $1.ps
#psxy ../profiles/profiles.gmt -m $PROJ -W2,white  -O -V -K >> $1.ps
psscale -Cprob.cpt -D5/1.5/8/0.3h -L -O  -K >> $1.ps
pstext $PROJ -O <<!>> $1.ps
-124.2 32 12 0 5 1 30 yr probability, M > 6.7
!