depthdat = config.get('FileIO', 'depthdat') stackfile = config.get('FileIO', 'stackfile') Velmod = config.get('FileIO', 'Velmod') stalist = config.get('FileIO', 'stalist') domperiod = float(config.get('para', 'domperiod')) Profile_width = float(config.get('para', 'Profile_width')) bin_radius = float(config.get('para', 'bin_radius')) Stack_range = np.arange(350, 751) azi = distaz.distaz(lat1, lon1, lat2, lon2).baz dis = distaz.distaz(lat1, lon1, lat2, lon2).delta Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width) Profile_lat = [] Profile_lon = [] for Profile_loca in Profile_range: (lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi, distaz.km2deg(Profile_loca)) Profile_lat = np.append(Profile_lat, [lat_loca], axis=1) Profile_lon = np.append(Profile_lon, [lon_loca], axis=1) # ----- Read depth .mat file -----# depthmat = sio.loadmat(depthdat) RFdepth = depthmat['YN_RFdepth'] # find stations beside the profile stalat_all = RFdepth[0, 0::]['stalat'] stalon_all = RFdepth[0, 0::]['stalon'] stalat_full = [] stalon_full = [] staidx = [] for i in range(stalon_all.shape[0]): stalat_full = np.append(stalat_full, stalat_all[i][0])
lat2 = float(line.split('/')[2]) lon2 = float(line.split('/')[3]) depthdat = config.get('FileIO', 'depthdat') stackfile = config.get('FileIO', 'stackfile') Velmod = config.get('FileIO', 'Velmod') stalist = config.get('FileIO', 'stalist') domperiod = float(config.get('para', 'domperiod')) Profile_width = float(config.get('para', 'Profile_width')) Stack_range = np.arange(1, 101) azi = distaz.distaz(lat1, lon1, lat2, lon2).baz dis = distaz.distaz(lat1, lon1, lat2, lon2).delta Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width) Profile_lat = [] Profile_lon = [] for Profile_loca in Profile_range: (lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi, distaz.km2deg(Profile_loca)) Profile_lat = np.append(Profile_lat, [lat_loca], axis=1) Profile_lon = np.append(Profile_lon, [lon_loca], axis=1) # read depth data depthmat = sio.loadmat(depthdat) RFdepth = depthmat['YN_RFdepth'] # find stations beside the profile stalat_all = RFdepth[0, 0::]['stalat'] stalon_all = RFdepth[0, 0::]['stalon'] stalat_full = [] stalon_full = [] staidx = [] for i in range(stalon_all.shape[0]): stalat_full = np.append(stalat_full, stalat_all[i][0])
lon2 = float(line.split('/')[3]) depthdat = config.get('FileIO', 'depthdat') stackfile = config.get('FileIO', 'stackfile') Velmod = config.get('FileIO', 'Velmod') stalist = config.get('FileIO', 'stalist') domperiod = float(config.get('para', 'domperiod')) Profile_width = float(config.get('para', 'Profile_width')) bin_radius = float(config.get('para', 'bin_radius')) Stack_range = np.arange(350, 751) azi = distaz.distaz(lat1, lon1, lat2, lon2).baz dis = distaz.distaz(lat1, lon1, lat2, lon2).delta Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width) Profile_lat = [] Profile_lon = [] for Profile_loca in Profile_range: (lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi, distaz.km2deg(Profile_loca)) Profile_lat = np.append(Profile_lat, [lat_loca], axis=1) Profile_lon = np.append(Profile_lon, [lon_loca], axis=1) # ----- Read depth .mat file -----# depthmat = sio.loadmat(depthdat) RFdepth = depthmat['YN_RFdepth'] # find stations beside the profile stalat_all = RFdepth[0, 0::]['stalat'] stalon_all = RFdepth[0, 0::]['stalon'] stalat_full = [] stalon_full = [] staidx = [] for i in range(stalon_all.shape[0]): stalat_full = np.append(stalat_full, stalat_all[i][0])
lon2 = float(line.split('/')[3]) depthdat = config.get('FileIO', 'depthdat') stackfile = config.get('FileIO', 'stackfile') Velmod = config.get('FileIO', 'Velmod') stalist = config.get('FileIO', 'stalist') domperiod = float(config.get('para', 'domperiod')) Profile_width = float(config.get('para', 'Profile_width')) Stack_range = np.arange(1, 101) azi = distaz.distaz(lat1, lon1, lat2, lon2).baz dis = distaz.distaz(lat1, lon1, lat2, lon2).delta Profile_range = np.arange(0, distaz.deg2km(dis), Profile_width) Profile_lat = [] Profile_lon = [] for Profile_loca in Profile_range: (lat_loca, lon_loca) = distaz.latlon_from(lat1, lon1, azi, distaz.km2deg(Profile_loca)) Profile_lat = np.append(Profile_lat, [lat_loca], axis=1) Profile_lon = np.append(Profile_lon, [lon_loca], axis=1) # read depth data depthmat = sio.loadmat(depthdat) RFdepth = depthmat['YN_RFdepth'] # find stations beside the profile stalat_all = RFdepth[0, 0::]['stalat'] stalon_all = RFdepth[0, 0::]['stalon'] stalat_full = [] stalon_full = [] staidx = [] for i in range(stalon_all.shape[0]): stalat_full = np.append(stalat_full, stalat_all[i][0])