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: