Traceback (most recent call last):
  File "fold.py", line 37, in <module>
    WF = sr.Fold(rh, dispersion_measure, frequency, sideband, polyco_file, polarization, fullpol, nthreads)
  File "/gpfs/fs0/scratch/p/pen/syedfard/B1133/gk049c/ar/scintelometry/scintillometry_reduction.py", line 79, in __init__
    self.integrator = self.initialize_pipeline(fh)
  File "/gpfs/fs0/scratch/p/pen/syedfard/B1133/gk049c/ar/scintelometry/scintillometry_reduction.py", line 95, in initialize_pipeline
    integrator = integration.Fold(power,n_phase=1024, phase=psr_polyco, step=self.interval, average=False)
  File "/home/p/pen/syedfard/scintillometry/scintillometry/integration.py", line 319, in __init__
    samples_per_frame=samples_per_frame)
  File "/home/p/pen/syedfard/scintillometry/scintillometry/integration.py", line 129, in __init__
    samples_per_frame=samples_per_frame, dtype=dtype)
  File "/home/p/pen/syedfard/scintillometry/scintillometry/base.py", line 379, in __init__
    polarization=polarization, dtype=dtype)
  File "/home/p/pen/syedfard/scintillometry/scintillometry/base.py", line 60, in __init__
    frequency = self._check_shape(frequency)
  File "/home/p/pen/syedfard/scintillometry/scintillometry/base.py", line 82, in _check_shape
    np.broadcast_to(value, self.sample_shape)
  File "/scinet/niagara/software/2018a/opt/base/anaconda3/2018.12/lib/python3.7/site-packages/numpy/lib/stride_tricks.py", line 176, in broadcast_to
    return _broadcast_to(array, shape, subok=subok, readonly=True)
  File "/scinet/niagara/software/2018a/opt/base/anaconda3/2018.12/lib/python3.7/site-packages/numpy/lib/stride_tricks.py", line 128, in _broadcast_to
    op_flags=[op_flag], itershape=shape, order='C')
ValueError: ('operands could not be broadcast together with remapped shapes [original->remapped]: (1,4,2) and requested shape (1,4,4)', 'value cannot be broadcast to sample shape')