Thread: SAC interact with csh

Started: 2008-10-21 22:14:23
Last activity: 2008-10-21 22:14:23
Topics: SAC Help
Ching-Wen Chen
2008-10-21 22:14:23
Dear All,

I am writing an instrument response correction script. I set stations' name
and event seperately. It works for "read" in SAC, but not for trans
polezero. Any suggestion on this one? Thank your for your time and effort.

--

set STA = "GAC"
set EVENT = "200411150906"
set YYYY = "2004"
set DISN = `cat stadis | grep $STA | awk '{ print $4 }'`
set STANUM = `cat stadis | grep $STA | awk '{ print $2}'`
/bin/cp 20*.$STA..BHZ.D.SAC D$DISN.$STANUM.$STA.Z
if (-e D$DISN.$STANUM.$STA.Z ) then
sac2000 << !

### read works####
r D$DISN.$STANUM.$STA.Z
rmean
taper
p1
########## This doesn't work due to $YYYY and $EVENT
trans from polezero s
/data/cchen/Ontario_CN/CNSN_sac/$YYYY\_sac/$EVENT/Resp_files/$STA\_BHZ to
polezero s
/data/cchen/Ontario_CN/CNSN_sac/generic_resp_files/SAC_PZs_resp03.HHZ
##########

w over
!
endif

-----


Chingwen Chen

Page built 18:36:01 | v.b'4482ed