./toolbox
dcfold.pro

Last modification date:

Thu Jun 2 23:56:08 2005

dcfold

function dcfold(sig, ref, [ftol=double])

Average two parts of in-band frequency switching (signal and reference phases, with the cal-switching phases already calibrated). Typically this happens after getfs. It does not matter which data container containers which part since their relative distance in frequency is used to determine how one is shifted to align with the other. The returned result is a new data container that the user must eventually free using data_free in order to avoid any memory leaks. sig and ref must have the same number of channels and the same spacing between channels. If sig and ref are already aligned, then this is a simple average (a warning message will be printed). If the shift necessary to align ref with sig is more than the total number of channels, there is no overlap and this procedure will print an error message and return without altering sig and ref.

The "ref" data container is shifted to align in sky frequency with the "sig" data container using dcshift and the two data containers are averaged - weighting each by the inverse of square of their system temperatures. The system temperature in the result is the weighted average of the two system temperatures.

Returns

data container. The user is responsible for freeing this. returns -1 on error.

Uses

dcshift data_free data_valid

Version

$Id: dcfold.pro,v 1.5 2005/05/19 20:43:56 bgarwood Exp $

Parameters
sig in/out, required variable	The data container to use as the signal part in the average.
ref in/out, required variable	The data container to use as the reference part in the average. This data is shifted using dcshift to align with "sig" before averaging.

Keywords

ftol in, optional double (def. 0.005)

The fractional channel shift tolerance. If the fractional part of the channel shift necessary to align the two parts is less than this value, no fractional shift as described in the documentatio for dcshift will be done. It might be useful to turn off the fractional shift because it can cause aliases and ringing in the case of very strong lines or other sharp features. If ftol > 0.5 no fractional shifts will be done.