#!/bin/sh set -x ##### # script for processing IRIS Digital Data ##### ###### # shell script to make entries for the various command files ###### YEAR=2006 MO=05 DY=10 HR=02 MN=42 SEC=56 MSEC=300 LAT=52.71 LON=-169.24 DEP=56.0 MAG=6.3 ##### # No changes below here ##### if [ -d ../GOOD ] then echo GOOD exists else mkdir ../GOOD fi for i in *.SAC do NET=`saclhdr -KNETWK $i` KSTNM=`saclhdr -KSTNM $i` KCMPNM=`saclhdr -KCMPNM $i` LOC=`saclhdr -KHOLE $i` DOY=`saclhdr -NZJDAY $i` DELTA=`saclhdr -DELTA $i` FHH=`echo $DELTA | awk '{print 0.45/$1}' ` FHL=`echo $DELTA | awk '{print 0.25/$1}' ` cp RESP.${NET}.${KSTNM}.${LOC}.${KCMPNM} resp evalresp ${KSTNM} ${KCMPNM} ${YEAR} ${DOY} 0.001 ${FHH} 2049 -u 'vel' -f resp gsac << EOF r $i ch EVLA $LAT EVLO $LON EVDP $DEP ch O CAL $YEAR $MO $DY $HR $MN $SEC $MSEC ch lovrok true ch lcalda true wh rtr w ../GOOD/${KSTNM}${KCMPNM}.S transfer from eval subtype AMP.${NET}.${KSTNM}.${LOC}.${KCMPNM} PHASE.${NET}.${KSTNM}.${LOC}.${KCMPNM} TO NONE FREQLIMITS 0.005 0.01 ${FHL} ${FHH} w ../GOOD/${KSTNM}${KCMPNM}.sac quit EOF done #SCRIPT FOR IRIS DATA