count_uniq_bll = count_uniq_bll[sortind]
total_baselines = bl_length.size
baseline_bin_indices = range(0, int(NP.ceil(1.0*total_baselines/baseline_chunk_size)+1)*baseline_chunk_size, baseline_chunk_size)
if n_bl_chunks is None:
n_bl_chunks = int(NP.ceil(1.0*total_baselines/baseline_chunk_size))
bl_chunk = range(len(baseline_bin_indices)-1)
bl_chunk = bl_chunk[:n_bl_chunks]
bl = bl[:min(baseline_bin_indices[n_bl_chunks], total_baselines),:]
bl_length = bl_length[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
bl_id = bl_id[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
n_channels = nchan
frac_width = DSP.window_N2width(n_window=None, shape=bpass_shape, fftpow=fftpow, area_normalize=True, power_normalize=False)
# window = NP.sqrt(frac_width * n_channels) * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=False, power_normalize=True)
PDB.set_trace()
# window = n_channels * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
window = n_channels * DSP.window_fftpow(n_channels, shape=bpass_shape, fftpow=fftpow, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
bw = n_channels * freq_resolution
bandpass_str = '{0:0d}x{1:.1f}_kHz'.format(nchan, freq_resolution/1e3)
pfb_instr = ''
pfb_outstr = ''
if pfb_method is None:
pfb_instr = '_no_pfb'
pfb_outstr = 'no_pfb_'
snapshot_type_str = ''
if avg_drifts:
count_uniq_bll = count_uniq_bll[sortind]
total_baselines = bl_length.size
baseline_bin_indices = range(0, int(NP.ceil(1.0*total_baselines/baseline_chunk_size)+1)*baseline_chunk_size, baseline_chunk_size)
if n_bl_chunks is None:
n_bl_chunks = int(NP.ceil(1.0*total_baselines/baseline_chunk_size))
bl_chunk = range(len(baseline_bin_indices)-1)
bl_chunk = bl_chunk[:n_bl_chunks]
bl = bl[:min(baseline_bin_indices[n_bl_chunks], total_baselines),:]
bl_length = bl_length[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
bl_id = bl_id[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
n_channels = nchan
frac_width = DSP.window_N2width(shape=bpass_shape, area_normalize=False, power_normalize=True)
# window = NP.sqrt(frac_width * n_channels) * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=False, power_normalize=True)
window = n_channels * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
bw = n_channels * freq_resolution
bandpass_str = '{0:0d}x{1:.1f}_kHz'.format(nchan, freq_resolution/1e3)
pfb_instr = ''
pfb_outstr = ''
if pfb_method is None:
pfb_instr = '_no_pfb'
pfb_outstr = 'no_pfb_'
snapshot_type_str = ''
if avg_drifts:
snapshot_type_str = 'drift_averaged_'
if beam_switch:
int(NP.ceil(1.0 * total_baselines / baseline_chunk_size) + 1) *
baseline_chunk_size, baseline_chunk_size)
if n_bl_chunks is None:
n_bl_chunks = int(NP.ceil(1.0 * total_baselines / baseline_chunk_size))
bl_chunk = range(len(baseline_bin_indices) - 1)
bl_chunk = bl_chunk[:n_bl_chunks]
bl = bl[:min(baseline_bin_indices[n_bl_chunks], total_baselines), :]
bl_length = bl_length[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
bl_id = bl_id[:min(baseline_bin_indices[n_bl_chunks], total_baselines)]
n_channels = nchan
frac_width = DSP.window_N2width(n_window=None,
shape=bpass_shape,
fftpow=fftpow,
area_normalize=True,
power_normalize=False)
# window = NP.sqrt(frac_width * n_channels) * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=False, power_normalize=True)
PDB.set_trace()
# window = n_channels * DSP.windowing(n_channels, shape=bpass_shape, pad_width=0, centering=True, area_normalize=True, power_normalize=False)
window = n_channels * DSP.window_fftpow(n_channels,
shape=bpass_shape,
fftpow=fftpow,
pad_width=0,
centering=True,
area_normalize=True,
power_normalize=False)
bw = n_channels * freq_resolution
bandpass_str = '{0:0d}x{1:.1f}_kHz'.format(nchan, freq_resolution / 1e3)
