SCULIB_CALC_SKYDIP_TEMPS - Calculate all the skydip temps from raw data

Description:

This routine accepts the SKYDIP (3 dimensional) raw data and works out the temperature associated with each integration. The 2 dimensional array (bolometers, integrations) is then returned.

Invocation:

CALL SCULIB_CALC_SKYDIP_TEMPS( N_TEMPS, N_FITS, FITS, PARAM, N_INTEGRATIONS, N_MEASUREMENTS, N_POS, N_BOLS, DEM_PNTR, IN_DATA, OUT_DATA, OUT_VAR, OUT_QUAL, OUT_DEM_PNTR, OUT_DATA_AV, OUT_VAR_AV, OUT_QUAL_AV, SCRATCH, STATUS)

Arguments:

N_TEMPS = INTEGER (Given)

Number of SKYDIP temperatures.

SUB_REQUIRED = INTEGER (Given)

The number of the sub-instrument that is requested. If this number is less than 1 or greater than N_BOLS then all subinstrumentes are returned. This should be 0 when called from REDUCE_SWITCH and the required number if called from SKYDIP. Exists purely to prevent multiple requests for T_COLD if we are only interested in one wavelength

N_FITS = INTEGER (Given)

Number of FITS items

FITS = CHARACTER$*$80 (Given)

The FITS items

PARAM = CHARACTER$*$(PAR__SZNAM) (Given)

Name of parameter from which T_COLD can be read

N_INTEGRATIONS = INTEGER (Given)

Number of integrations

N_MEASUREMENTS = INTEGER (Given)

Number of measurements

N_POS = INTEGER (Given)

Number of data points

N_BOLS = INTEGER (Given)

Number of bolometers (sub instruments)

NUM_CHAN = INTEGER (Given)

number of A/D channels per A/D card

NUM_ADC = INTEGER (Given)

number of A/D cards

BOL_TYPE (NUM_CHAN, NUM_ADC)

= CHARACTER$*$($*$) (Given) the types of the bolometers

BOL_CHAN (N_BOLS) = INTEGER (Given)

channel numbers of selected bolometers

BOL_ADC (N_BOLS) = INTEGER (Given)

ADC numbers of selected bolometers

DEM_PNTR(1, 1, N_INTEGRATIONS, N_MEASUREMENTS) = INTEGER (Given)

Pointer to positions in the input data

IN_DATA(N_TEMPS, N_BOLS, N_POS) = REAL (Given)

Input data [NTEMPS, NBOLS, NPOS]

OUT_DATA(N_BOLS, N_POS) = REAL (Returned)

Processed sky brightness temperatures.

OUT_VAR(N_BOLS, N_POS) = REAL (Returned)

Variance on each point in OUT_DATA

OUT_QUAL(N_BOLS, N_POS) = BYTE (Returned)

Quality of OUT_DATA

OUT_DEM_PNTR(1, 1, N_INTEGRATIONS, N_MEASUREMENTS) = INTEGER (Returned)

Modified DEM_PNTR to reflect changes in data array

OUT_DATA_AV(N_MEASUREMENTS, N_BOLS) = REAL (Returned)

Average sky temp for each measurement

OUT_VAR_AV(N_MEASUREMENTS, N_BOLS) = REAL (Returned)

Variance of the Average sky temp for each measurement

OUT_QUAL_AV(N_MEASUREMENTS, N_BOLS) = BYTE (Returned)

Quality of the Average sky temp for each measurement

SCRATCH(N_TEMPS, N_BOLS, N_INTEGRATIONS) = REAL (Given)

Work space to store the data for each measurement (since I cant simply pass a pointer in as the data is non-contiguous)

STATUS = INTEGER (Given)

Global Status

Copyright

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