bands = range(int(bands[0]), int(bands[-1])+1)
band_period = {'1': 30.0, '2': 21.2, '3': 15.0, '4': 10.6, '5': 7.5, '6': 5.3, '7': 3.7, '8': 2.7}
proc_ev_ls = glob.glob(os.path.join(remote_dir, '*.*.*.*'))
print '%s processed events found!' % len(proc_ev_ls)
for i in range(len(bands)):
all_passed_staev = []
for j in range(len(proc_ev_ls)):
# [bands[i]] is passed like this because reader gets list as an input
# reader(all_stations=False, just_high_cc=False, remove_GSN_median=False):
all_staev = uf.reader(proc_ev_ls[j], [bands[i]], band_period, all_stations=all_stations,
just_high_cc=just_high_cc, remove_GSN_median=remove_GSN_median)
if not all_staev:
continue
passed_staev = uf.filters(all_staev, [bands[i]], xcorr_limit=xcorr_limit)
if not passed_staev[0]:
continue
all_passed_staev.append(passed_staev)
t_shift_array = []
for j in range(len(all_passed_staev)):
# [0] in all_passed_staev[j][0] shows the current band
# in general we have a loop over bands and in each step there is just
# one band that we are working with which is accessible by [0]
for k in range(len(all_passed_staev[j][0])):
if not all_passed_staev[j][0][k] == []:
if not nr_cc:
t_shift_array.append(all_passed_staev[j][0][k][2])
else:
intro += '\npython dispersion_dT_dA.py bands[1-3] address-for-FFM-event.\n'
intro += 20*'-'
bands = sys.argv[1]
bands = range(int(bands[0]), int(bands[-1])+1)
band_period = {'1': 30.0, '2': 21.2, '3': 15.0, '4': 10.6, '5': 7.5, '6': 5.3, '7': 3.7, '8': 2.7}
if not all_processed:
evadd = sys.argv[2]
evname = os.path.baseanme(evadd)
if not evname:
evname = evadd.split('/')[-2]
all_staev = uf.reader(evadd, bands, band_period, all_stations=all_stations, just_high_cc=just_high_cc,
remove_GSN_median=remove_GSN_median)
passed_staev = uf.filters(all_staev, bands, xcorr_limit=xcorr_limit)
if plot_scatter:
uf.ffpscatter(passed_staev)
if plot_mean:
per, dt_mean, a_mean, flag = uf.stamean(passed_staev)
uf.ffplotmean(per, dt_mean)
else:
if not just_selected_events:
proc_ev_ls = glob.glob(os.path.join(remote_dir, '*.*.*.*'))
else:
proc_ev_ls = []
for ev in list_events:
proc_ev_ls.append(os.path.join(remote_dir, ev))
print '%s processed events found!' % len(proc_ev_ls)
if update_all: