Description:
This routine unpacks data from a jiggle observation into a rectangular 2-d
map.
If status is good on entry the data and variance of the output map will be initialised to ‘bad’ values,
the quality to 1. If there are any data to unpack the routine will then attempt to do so. The method
used depends on whether the switch covered part/all of the entire jiggle pattern or contains data for
several repeats of the jiggle pattern.
In the first case, the routine will check that the indices of the section of jiggle pattern covered by this
switch lie within the bounds of the full jiggle pattern. If not, an error will be reported and bad status
returned. Otherwise, the data will be unpacked into the map as specified by the I_JIGGLE, J_JIGGLE
arrays.
When the data cover several repeats of the jiggle pattern, the routine will check that the jiggle index of
the first point in the datablock is 1 and that the number of data points is an integer multiple of the size
of the jiggle. If not, an error will be reported and bad status returned. Otherwise, the data will be
unpacked into the map as specified by I_JIGGLE, J_JIGGLE. The unpacking differs from
the first case in that points in the unpacked map will be set to the average of the good
quality values contributed by the separate jiggles, variances will be calculated from the
dispersion of the values about the mean if there was more than 1, or set to the input variance
otherwise.
Invocation
CALL SCULIB_UNPACK_JIGGLE_SEPARATES (N_JIGGLES, N_BOLS,
J_DATA, J_VARIANCE, J_QUALITY, J_START, JIGGLE_COUNT, I_JIGGLE, J_JIGGLE, IDIM,
JDIM, MAP_DATA, MAP_VARIANCE, MAP_QUALITY, MAP_NUMPTS, J_END, BADBIT, STATUS)
Arguments
N_JIGGLES = INTEGER (Given)
number of jiggles in datablock
N_BOLS =
INTEGER (Given)
number of bolometers measured
J_DATA (N_BOLS, N_JIGGLES) = REAL
(Given)
data for each jiggle measured
J_VARIANCE (N_BOLS, N_JIGGLES) = REAL
(Given)
variance on J_DATA
J_QUALITY (N_BOLS, N_JIGGLES) = BYTE (Given)
quality on J_DATA
J_START = INTEGER (Given)
index of first jiggle in datablock in the
overall jiggle pattern
JIGGLE_COUNT = INTEGER (Given)
number of jiggles in jiggle
pattern
I_JIGGLE (JIGGLE_COUNT) = INTEGER (Given)
i index in output map of each
jiggle position in pattern
J_JIGGLE (JIGGLE_COUNT) = INTEGER (Given)
j index in
output map of each jiggle position in pattern
IDIM = INTEGER (Given)
i dimension
of output map
JDIM = INTEGER (Given)
j dimension of output map
MAP_DATA
(IDIM, JDIM, N_BOLS) = REAL (Returned)
output map data
MAP_VARIANCE (IDIM,
JDIM, N_BOLS) = REAL (Returned)
output map variance
MAP_NUMPTS (IDIM, JDIM,
N_BOLS) = INTEGER (Returned)
output map quality
MAP_QUALITY (IDIM, JDIM,
N_BOLS) = BYTE (Returned)
output map quality
J_END = INTEGER (Returned)
index of
last jiggle in datablock in the overall jiggle pattern
BADBIT = BYTE (Given)
bad bit
mask
STATUS = INTEGER (Given and returned)
global status
Copyright
Copyright
©1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.