SCULIB_2POS_DECONV - deconvolve square chop from scan

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.