Good day! I have a problem replicating a result. Let us consider the top plot here:
On the plot, velocity is in mm/s. I think the problem is to remove instrument response correctly. I downloaded seismogram and pole-zero instrument response. I want to obtain velocity data. The instrument response is already given in velocity. I was able to get synthetic seismograms as shown in the figure, but I have a problem working with real data.
I need a signal to be filtered in 0.001–0.004 Hz
Using IRIS SAC software, I do:
SAC> r KYK_BHZ SAC> rmean SAC> rtrend SAC> taper SAC> trans from polezero subtype KYK_LHZ.zp SAC> bp p 2 n 4 c 0.001 0.004
- I do not set
FREQlimitsfor transfer function at all here, but I tried to do it like
0.0009 0.001 0.004 0.005or
0.0009 0.001 1 2, but it also doesn't work.
- My synthetic seismogram is 1 Hz data, so I don't do any resampling.
- I tried to do
trans from polezero subtype KYK_LHZ.zp to vel, but it also doesn't work.