[obspy-users] gain/sensitivity in simulate

Stephen Bannister S.Bannister at gns.cri.nz
Mon Dec 19 04:50:50 CET 2011


Hi,

A problem with sensitivity in "simulate"  - any suggestions welcome 

I have imported details of a RESP response file, using  "attach_resp" from 
obspy.sac,   and then have the following 
in the trace (tr) stats  for that trace : 

         network: NZ
         station: KNZ
        location: 10
         channel: HHN
       starttime: 2011-12-10T00:00:00.968393Z
         endtime: 2011-12-10T23:58:59.998393Z
   sampling_rate: 100.0
           delta: 0.01
            npts: 8633904
           calib: 1.0
         _format: MSEED
             paz: AttribDict({'digitizer_gain': 1.0, 'seismometer_gain': 
838861000.0, 'zeros': [0j, 0j, (159+0j), 0j, 0j, (159+0j), 0j, 0j, 
(159+0j), 0j, 0j, 0j], 'gain': 1932953847860000.0, 't_shift': 
1.2483732000000001, 'poles': 
[(-0.023560000000000001+0.023560000000000001j), 
(-0.023560000000000001-0.023560000000000001j), (-50+0j), 
(-0.023560000000000001+0.023560000000000001j), 
(-0.023560000000000001-0.023560000000000001j), (-50+0j), 
(-0.023560000000000001+0.023560000000000001j), 
(-0.023560000000000001-0.023560000000000001j), (-50+0j), 
(-0.011780000000000001+0.011780000000000001j), 
(-0.011780000000000001-0.011780000000000001j), (-180+0j), (-160+0j), 
(-80+0j)]})
           mseed: AttribDict({'dataquality': 'D', 'record_length': 512, 
'byteorder': '>', 'encoding': 'STEIM2'})


If  I then try   tr.simulate (paz_remove='self', remove_sensitivty=True) 
then I get the following : 

Traceback (most recent call last):
  File "plot_2011_displacement.py", line 893, in <module>
    tr4 = scbresp.stream_remove_response(tr3)
  File "responsework.py", line 156, in stream_remove_response
    tr2.simulate(paz_remove='self',remove_sensitivity=True)
  File 
"/sp/prj/tvz95/mypython/lib/python2.6/site-packages/obspy.core-0.4.8.dev_r2862-py2.6.egg/obspy/core/trace.py", 
line 1055, in simulate
    simulate_sensitivity=simulate_sensitivity, **kwargs)
  File 
"/sp/prj/tvz95/mypython/lib/python2.6/site-packages/obspy.signal-0.4.9.dev_r2863-py2.6-linux-x86_64.egg/obspy/signal/invsim.py", 
line 448, in seisSim
    data /= paz_remove['sensitivity']
  File 
"/sp/prj/tvz95/mypython/lib/python2.6/site-packages/obspy.core-0.4.8.dev_r2862-py2.6.egg/obspy/core/util.py", 
line 134, in __getitem__
    return super(AttribDict, self).__getitem__(name)
KeyError: 'sensitivity'


Any suggestions anyone ??

Thanks
Stephen

Notice: This email and any attachments are confidential. If received in error please destroy and immediately notify us. Do not copy or disclose the contents.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.obspy.org/pipermail/obspy-users/attachments/20111219/b601a10c/attachment.html>


More information about the obspy-users mailing list