#shell to plot RF wiggles from the "grid" files generated by 
#spectral coherence code of J. Park, arranged by source distance (epi files)
#Input - $1 - file name $2 file name $3 - scale

rm wigtab? tmp?

echo "EDIT PLOTTING PARAMETERS TO SUIT YOUR DATA"
set BOX = -R-2/15/20/160
set FRAME = -JX3/7.5

awk '{print $2}' $1 | sort -u -n > bazlist
foreach i ( `cat bazlist` )
awk '$2 == A && $1 >= -2' A=$i $1 >> wigtabR
awk '$2 == A && $1 >= -2' A=$i $2 >> wigtabT

echo ">" >> wigtabR
echo ">" >> wigtabT

end

set scale = $3

pswiggle wigtabR $FRAME $BOX -Z$scale -M -G0/0/255 -Ba4f1:"time, s":/f5a30:"epicentral distance, deg":\SWen -P -K > RFwig-epi.ps
pswiggle wigtabR $FRAME $BOX -Z$scale -W1/0 -M  -P -K -O >>  RFwig-epi.ps
echo 0 3 14 0 5 6 $1 scale $3 | pstext $FRAME -R-1/1/-2/2 -N -O -K >> RFwig-epi.ps
psxy $FRAME $BOX -W2/255/0/0 -O -K <<END >> RFwig-epi.ps
0 -10
0 380
END

#set scale = `echo $scale|awk '{print $1/2}' `
#SOMETIMES it is helpfull to boost the scale on T by a factor...

pswiggle wigtabT $FRAME $BOX -Z$scale -M -G0/0/255 -Ba4f1:"time, s":/f5a30\Swen -P -K -O -X3.3 >> RFwig-epi.ps
pswiggle wigtabT $FRAME $BOX -Z$scale -M -W1/0 -P -K -O >>  RFwig-epi.ps
psxy $FRAME $BOX -W2/255/0/0 -O -K <<END >> RFwig-epi.ps
0 -10
0 380
END

echo 0 3 14 0 5 6 $2 scale $scale | pstext $FRAME -R-1/1/-2/2 -N -O >> RFwig-epi.ps