#!/bin/bash
# plot CMTs for Tonga
dfile=../cmt_012112.mdat	# data
dmin=100 			# minimum depth
ofile=tonga.ps
reg=-R172/190/-30/-14
proj=-JH`echo $reg | gawk -f reg2midlon.awk`/7
makecpt -T$dmin/700/10 -Cno_green -D > tmp.cpt
ann=-Ba5f1WeSn
# data selection 
gmtselect $reg -fg $dfile > regional_cmt.dat 
pscoast -Df -A5000 $reg $proj -K -P -G128 -S200 > $ofile
# for momement tensor 
sort -n +2 -r regional_cmt.dat | psmeca $reg $proj -Ztmp.cpt  -Sm0.15/-1 -D$dmin/1000 -K -O >> $ofile
# for actual double couple
#sort -n +2 -r regional_cmt.dat | psmeca $reg $proj -Ztmp.cpt  -Sd0.15/-1 -D$dmin/1000 -K -O >> $ofile
psxy ../nuvel.360.xy -O $reg $proj -K -m -fg -W2,darkorange >> $ofile
psscale -D6/1.5/2/.2 -Ctmp.cpt -Ba200/:"z [km]": -O -K >> $ofile
psbasemap $ann $reg $proj -O  >> $ofile
echo $0: output in $ofile
modifybb $ofile
rm tmp.*