#!/bin/sh
#
PROJ="-R-125/-114/31.0/42.0 -Jm1.7"
gmtset COLOR_MODEL hsv D_FORMAT %lg PLOT_DEGREE_FORMAT D
makecpt -Cno_green -Qo -T1/1000/1 -D -Z > strain.cpt
makecpt -Cgray -T-1.4/.7/.1 -Z > grad.cpt
#
#   image the topo and strain grids
#
grdimage grad.grd $PROJ -Cgrad.cpt -X1.5 -Y1.5 -K -V -P > $1.ps
grdimage $1_strain.grd -Igrad.grd $PROJ -Cstrain.cpt -Q -O -K -V >> $1.ps
pscoast $PROJ -W6 -Di -A200 -B2WSen:."$1": -O -V -N1 -N2 -K >> $1.ps
psxy fault_data.gmt -M $PROJ -W2 -O -V -K >> $1.ps
psxy model.geom -M $PROJ -K -O -Wthick/255 -A >> $1.ps
# psxy ../profiles/profiles.gmt -M $PROJ -W8,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)
!

open /Applications/Preview.app $1.ps