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.