#!/bin/sh #set -x ##### # run three tests of quality # P , SH SV ##### ##### # create the DAT directory if it does not exist ##### if [ ! -d ../DAT ] then mkdir ../DAT fi ##### # move to the directory for QC # then examine Z traces for P, T traces for SH and R traves for SV # Note the theoretical arrival times have already been placed into # the trace headers ##### cd FINAL.QC gsac << EOF ##### # process P ##### fileid name r *Z rtr hp c 0.01 n 2 lp c 0.0833 n 2 int xlim a -30 a 140 ppk q relative perplot 4 wh ##### # process SH ##### r *T rtr hp c 0.01 n 2 lp c 0.05 n 2 int xlim t1 -100 t1 200 ppk q relative perplot 4 wh ##### # process SV ##### r *R rtr hp c 0.01 n 2 lp c 0.05 n 2 int xlim t1 -100 t1 200 ppk q relative perplot 4 wh q EOF ##### # now move the selected traces to the processing area ##### for i in *[ZRT] do IHDR20=`saclhdr -IHDR20 $i` ANS=`echo $IHDR20 | awk '{if( $1 < 1)print "NO";else if($1 >1)print "NO" ; else print "YES" }' ` if [ $ANS = "YES" ] then echo $i $IHDR20 $ANS cp $i ../../DAT fi done