This routine corrects demodulated data for the effect of sky opacity. All data is corrected assuming the same airmass of observation, but the zenith sky opacity used will depend which SCUBA sub-instrument each bolometer belongs to.

After checking status on entry, the routine will loop through the bolometers measured. For each bolometer it will then ascertain the appropriate zenith sky opacity from the parent sub-instrument and wavelength of observation. The total sky optical depth will then be calculated and, if this is in the range 0 to 20, the flux correction factor derived. The routine will then loop through the positions measured by this bolometer, correcting the fluxes and variances. If the bolometer data quality was bad, or if sky optical depth lay outside the above range, the data quality will be set bad and no correction applied.

CALL SCULIB_REMOVE_OPACITY (N_BOLS, N_POS, BOL_SELECT_CHAN, BOL_SELECT_ADC, NUM_CHAN, NUM_ADC, BOL_TYPE, N_SUB, SUB_INSTRUMENT, TAUZ, AIRMASS, EXP_DATA, EXP_VARIANCE, EXP_QUALITY, STATUS)

N_BOLS = INTEGER (Given)

N_POS = INTEGER (Given)

BOL_SELECT_CHAN (N_BOLS) = INTEGER (Given)

BOL_SELECT_ADC (N_BOLS) = INTEGER (Given)

NUM_CHAN = INTEGER (Given)

NUM_ADC = INTEGER (Given)

BOL_TYPE (NUM_CHAN, NUM_ADC) = CHARACTER$\ast$($\ast$) (Given)

N_SUB = INTEGER (Given)

SUB_INSTRUMENT (N_SUB) = CHARACTER$\ast$($\ast$)

TAUZ (NSUB) = REAL (Given)

AIRMASS = REAL (Given)

EXP_DATA (N_BOLS, N_POS) = REAL (Given and returned)

EXP_VARIANCE (N_BOLS, N_POS) = REAL (Given and returned)

EXP_QUALITY (N_BOLS, N_POS) = REAL (Given and returned)

STATUS = INTEGER (Given and returned)

Copyright ©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.