#!/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 -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 >$1.ps
grdimage $1_cut.grd -Igrad.grd $PROJ -Cstrain.cpt -Q -O -K -V >>$1.ps
pscoast $PROJ -W3 -Di -A200 -B2WSen:."$1": -O -V -N1 -N2 -K >> $1.ps
psxy fault_data.gmt -m $PROJ -W1 -O -V -K >> $1.ps
psxy ../04_profiles/profiles.gmt -m $PROJ -W3,white  -O -V -K >> $1.ps
pstext labels.txt $PROJ -F+f+a+j -Ggray -O -K >> $1.ps
psscale -Cstrain.cpt -D5/1.5/8/0.3h -L -O  -K >> $1.ps
pstext $PROJ -O <<!>> $1.ps
-121.2 32 12 0 5 1 2nd invariant (nanostrain/yr)
!