SURFLIB_CLIP_GRID

Remove spikes from grid

Description:

Given the binned data calculated by SURFLIB_FILL_GRID and the limits calculated by SURFLIB_STATS_GRID, clip all points that lie outside the statistical limits provided by the stats array.

Invocation

CALL SURFLIB_CLIP_GRID(N_FILES, N_PTS, NX, NY, NMAX, NSIGMA, N_BOLS, BITNUM, QUALITY_PTR, BINS, BIN_POS, STATS, IPOS, JPOS, NSPIKES, STATUS)

Arguments

N_FILES = INTEGER (Given)
Number of input data sets
N_PTS(N_FILES) = INTEGER (Given)
Number of points in each data set
NX = INTEGER (Given)
Size of X dimension
NY = INTEGER (Given)
Size of Y dimension
NMAX = INTEGER (Given)
Maximum value allowed for third dimension of BINS
N_BOLS ( N_FILES ) = INTEGER (Given)
Number of bolometers per file (used for constructing message for users)
BITNUM = INTEGER (Given)
Bit number that is modified when set bad
DATA_PTR( N_FILES ) = INTEGER (Given & Returned)
The input data with each spike replaced with a median value
QUALITY_PTR(N_FILES) = INTEGER (Given & Returned)
Pointers to quality arrays. Note that quality arrays are modified. (The whole point of this routine) even though the pointers are not.
BINS(NX, NY, NMAX) = REAL (Given)
The data stored in relation to its position
BIN_POS (NX, NY, NMAX) = INTEGER (Given)
Position of each data point in quality arrays
STATS (NX, NY, 3) = REAL (Given)
Statistics of each bin. 1=Median, 2=High, 3=Low
PNTS(NMAX) = REAL (Given)
Scratch space for copying in the data from each I,J
IPOS(NX∗NY) = INTEGER (Given)
I position in array
JPOS(NX∗NY) = INTEGER (Given)
J position in array
NSPIKES = INTEGER (Returned)
Number of spikes removed for a given file
STATUS = INTEGER (Given & Returned)
Global Status

Copyright

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