*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
Thursday, August 26, 2010
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment