#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
tpf=0.05 #taper fraction
dt=0.01 #Define df
wt=17 #water level
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
if units =='0':#No se aplica respuesta instrumental