# --------------------------------------------------------------
# ------------------- MAIN PROGRAM -----------------------------
# --------------------------------------------------------------
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}
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] == []:
