Thursday, August 26, 2010

sac remove instrument response

*InstRm: remove instrument signal from seismic signal
[lhsu@thought leslie]$ more InstRm

* read signal
r 2008.005.00.00.16.1150.XQ.ME30.01.BHZ.M.SAC

*possibly transfer first, but remove mean and trend
rmean
rtrend

*taper ends
taper w 0.1

*transfer function to remove instrument noise
transfer from evalresp to vel freqlimit 0.0009 0.001 15 17

*convert to m from nm (sac default)
mul 1e-9

*write out the resultant velocity vector (binary)
w Day005_vel

*fft
fft

*write out the power spectrum (binary)
wsp Day005_sac

No comments:

Post a Comment