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.