#! /usr/bin/env python

import obspy
import obspy.xseed

# fix data to read here!
fname = 'test/1434692_data.fseed'

# set pre-filter corners ...
set_filter = ( 0.1, 1, 30, 40 )

# read the dataset (trace & response) from file ...
st = obspy.read(fname)
p  = obspy.xseed.Parser(fname)

# (2) evalresp method, calls external executable;

# RESP come from fullseed and are written by Obspy
if True:
    # set temporal filename of RESP file
    respdir = 'test/RESP-obspy'
    # write set of RESP files from parser
    p.writeRESP(folder=respdir)

# HERE: just referece an existing file ...
if False:
    respdir = 'RESP-1'

st_c2 = st.copy()
for tr in st_c2:
    respfile = '{:s}/RESP.{:s}'.format(respdir, tr.id)
    seedresp = {'filename': respfile,
                'date':     tr.stats.starttime,
                'units':    'VEL'
                }

    tr.simulate(paz_remove=None, pre_filt=set_filter, seedresp=seedresp)

    # invalidate respose object
    del seedresp