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.