SCULIB_3POS_DECONV

deconvolve 3-position chop from scan

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.