Description:
This routine
deconvolves the 3 position 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 may happen if the
observation was aborted, then no further action is taken. If a 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_3POS_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_3POS_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
STATUS = INTEGER (Given and returned)
global status
Copyright
Copyright
©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.