Description:
This
routine deconvolves the chopped beam response from the individual scans of a SCUBA
raster map. To achieve this it cycles through the scans making up the observation, calling
SCULIB_FIND_SWITCH to locate the start and finish indices of each scan in the demodulated data
array. If SCULIB_FIND_SWITCH indicates that there is no data for this scan, which might happen if
the observation was aborted, then no further action is taken. If the scan is too long to be handled by
the routine an error message will be output and the routine will return with bad status. Otherwise,
SCULIB_GENSYCONFN and SCULIB_2POS_CONFN will be called to generate the convolution
functions needed to deconvolve the chop. The routine then cycles through the bolometers, calling
SCULIB_CONVOLVE to do the required convolutions with the scan data for each.
Invocation
CALL
SCULIB_2POS_DECONV (N_EXPOSURES, N_INTEGRATIONS, N_MEASUREMENTS, DEMOD_POINTER, N_BOL,
N_POS, IN_DATA, IN_VARIANCE, IN_QUALITY, SAMPLE_DX, BEAM_SEP, OUT_DATA, OUT_VARIANCE,
OUT_QUALITY, STATUS)
Arguments
N_EXPOSURES = INTEGER (Given)
maximum
number of exposures per integration
N_INTEGRATIONS = INTEGER (Given)
number of
integrations in the observation
N_MEASUREMENTS = INTEGER (Given)
number of
measurements in the observation
DEMOD_POINTER (N_EXPOSURES, N_INTEGRATIONS,
N_MEASUREMENTS)
array pointing to start and finish of scans in IN_DATA
N_BOL = INTEGER
(Given)
the number of bolometers for which data was taken
N_POS = INTEGER (Given)
the number of positions measured in the scan
IN_DATA (N_BOL, N_POS) = REAL
(Given)
the measured data
IN_VARIANCE (N_BOL, N_POS) = REAL (Given)
the
variance on IN_DATA
IN_QUALITY (N_BOL, N_POS) = BYTE (Given)
the quality
on IN_DATA
SAMPLE_DX = REAL (Given)
the measurement spacing along the scan
(arcseconds)
BEAM_SEP = REAL (Given)
the beam separation (arcseconds)
OUT_DATA
(N_BOL, N_POS) = REAL (Returned)
the deconvolved data
OUT_VARIANCE (N_BOL,
N_POS) = REAL (Returned)
the variance on OUT_DATA
OUT_QUALITY (N_BOL,
N_POS) = BYTE (Returned)
the quality on OUT_DATA
BADBIT = BYTE (Given)
bad bit
mask
STATUS = INTEGER (Given and returned)
global status
Copyright
Copyright
©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.