Seismometer response calculation

The instrument specifications show the poles and zeroes of the analog model.
Example: Guralp CMG T3 from Passcal

CMG
   T3
     5 
Real
Poles
Imag
-3.70E-02 +3.70E-02
-3.70E-02 -3.70E-02
-5.03E+02 0.00E+00
-1.01E+03 0.00E+00
-1.13E+03 0.00E+00
  2 
Zeros
0.00E+00 0.00E+00
0.00E+00 0.00E+00



Here's the basics of the calculation:


In the Guralp case there is a big conversion factor, order of 10-9. This could be m/s/cnt, where cnt is the count (the raw-data) unit.
 
The tsf-edit instructions were
TSF EDIT CMGT3
! Guralp CMG T3 120s
! http://www.passcal.nmt.edu/content/instrumentation/sensors/
! ... sensor-comparison-chart/poles-and-zeroes#cmg_3t
ADD 1.0 At #0
ZSPF P&Z:5,2 Hz UC=-1 +PC=0. +STEP SG=1.e9
ZP(1) (-3.7e-2 ,  3.7e-2)
ZP(2) (-3.7e-2 , -3.7e-2)
ZP(3) (-5.03e2 ,  0.0   )
ZP(4) (-1.01e3 ,  0.0   )
ZP(5) (-1.13e3 ,  0.0   )
REMDC
END
used in
tslist _24000 -rs0.01 -B2013,1,1 -Eseismometer-response,CMGT3 -N -n1 \
       -w resp.dat > ! resp.spec
The script producing the calculations and plot is found in  sas/seismometer-response.

.bye