#!/bin/sh
DEST=../FINAL
FHIGH=0.002
#####
#	script for processing IRIS Digital Data
#####

if [ -d  ${DEST} ]
then
	echo ${DEST} exists
else
	mkdir ${DEST}
fi

if [ -d TEMP ]
then
	echo TEMP exists
else
	mkdir TEMP
fi

for i in *.sac
do
gsac > /dev/null 2>&1 << EOF
r $i
rtr
hp c ${FHIGH} np 2
w TEMP/$i
quit
EOF
done

cd TEMP
for i in *BHZ.sac
do
if [ -f $i ] 
then
echo $i

KSTNM=`saclhdr -KSTNM $i`
KCMPNM=`saclhdr -KCMPNM $i`
#####
#	safety if we do not have horizontals
#####

cp $i ../${DEST}/${KSTNM}${KCMPNM} 
gsac > /dev/null 2>&1 << EOF
r ${KSTNM}BH*.sac 
rotate3 to gc
w 
quit
EOF
rm ${KSTNM}BH?.sac
mv ${KSTNM}BH* ../${DEST}
else
	echo $i does not exist
fi
done
cd ..

cd TEMP
for i in *LHZ.sac
do
if [ -f $i ] 
then
echo $i

KSTNM=`saclhdr -KSTNM $i`
KCMPNM=`saclhdr -KCMPNM $i`
#####
#	safety if we do not have horizontals
#####

cp $i ../${DEST}/${KSTNM}${KCMPNM} 
gsac > /dev/null 2>&1 << EOF
r ${KSTNM}BH*.sac 
rotate3 to gc
w 
quit
EOF
rm ${KSTNM}LH?.sac
mv ${KSTNM}LH* ../${DEST}
else
	echo $i does not exist
fi
done
cd ..
#SCRIPT FOR IRIS DATA