#=========================================================
dicLog = {}#Diccionario que albergará listas de logA por estación
Gamma=[]
heads={}
sfiles=np.loadtxt('sfiles.txt',dtype=str)
for evento in sfiles:
print evento
#Lista que albergara las distancias de las estaciones a estudiar
logr=[]
r=[]
logA=[]
amp_t=[]
stations=[]
traces,head=cut_traces(evento,cut)#diccionario que alberga trazas cortadas
dist=distance(evento)#Diccionario que alberga la distancia de cada estacion al evento
for station in traces:
#print station
trN=traces[station][0]
trE=traces[station][1]
#print trN
seedrespN=respN(station,units)
#print seedrespN
seedrespE=respE(station,units)
if seedrespE != {} and seedrespN != {}:
# print station
cut=1
#file_name=raw_input('06-0503-41L.S201506')
#units=raw_input('Ingrese las unidades de salida de la repuesta instrumental (VEL, DIS, ACC), si NO desea filtrar 0: ')
units='DIS'
sfile = '13-1807-32L.S201405' #Magnitud 3.7
#11-0119-41L.S201305 27-1021-16L.S201403
#31-0253-08L.S201401 13-1807-32L.S201405
#01-2250-03L.S201403 04-2306-20L.S201406
#15-0655-42L.S201403 25-0727-39L.S201406
#17-0430-31L.S201403 20-2047-32L.S201407
#17-1026-05L.S201403 24-0215-50L.S201508
#24-2319-31L.S201403
traces,head=cut_traces(sfile,cut)#diccionario que alberga trazas cortadas
#Diccionario que alberga la distancia de cada estacion al evento
dist=distance(sfile)
#Lista que albergara las distancias de las estaciones a estudiar
logr=[]
logA=[]
#Diccionario que alberga los logaritmos de las amplitudes para encontrar diferencia entre cut_traces4 y 5
#30 35
#=========================================================
pre_filt = (0.02, 0.03, 20.0, 25.0)
<<<<<<< HEAD
=======
<<<<<<< HEAD
from cut_traces3 import cut_traces
from resp import respN
from obspy.signal.invsim import seisSim,simulate_seismometer
units = 'DIS'
station = 'GUA'
pre_filt = (0.2, 0.3, 20.0, 25.0)
dic,head=cut_traces('01-2250-03L.S201403',0)
tr=dic[station][0]
tr2 = tr.copy()
tr1 = tr.copy()
seedrespN=respN(station,units)
tr1.simulate(paz_remove=None, pre_filt=pre_filt, seedresp=seedrespN)
seedrespN['date']=tr2.stats.starttime
seedrespN['network']=tr2.stats.network
seedrespN['station']=tr2.stats.station
seedrespN['location']=tr2.stats.location
seedrespN['network']=tr2.stats.network
seedrespN['channel']=tr2.stats.channel
tr2.data=simulate_seismometer(tr2.data, 100, remove_sensitivity=False, simulate_sensitivity=False,
water_level=0, seedresp=seedrespN, pre_filt=None )
tr.plot()
#tr1.plot()
tr2.plot()