Interpolates from an array of data on a "sparse" rectangular grid which is regularly spaced in x and y to an array of data on a "dense" rectangular grid and initializes contouring from the array on the dense grid. (by a "sparse" grid is meant one whose dimensions are smaller than one would like, so that contour lines constructed directly on it are composed of long straight segments.) cpsps1 may be viewed as a data smoothing routine. cpsprs is an alternate name for the routine cpsps1.
CALL CPSPS1 (ZSPS, KSPS, MSPS, NSPS, RWRK, LRWK, IWRK,
+ LIWK, ZDAT, LZDT)
#include <ncarg/ncargC.h>
void c_cpsps1 (float *zsps, int ksps, int msps, int nsps,
float *rwrk, int lrwk, int *iwrk, int liwk, float *zdat,
int lzdt)
(REAL array, dimensioned KSPS x n, where "n" is greater than or equal to NSPS, input) is the "sparse" array of data, from which the "dense" array is to be generated.
(INTEGER, input) is the first dimension of the array ZSPS.
(INTEGER, input) is the first dimension of the "sparse" array of data in ZSPS. MSPS must be less than or equal to KSPS.
(INTEGER, input) is the second dimension of the "sparse" array of data in ZSPS. NSPS must be less than or equal to the declared second dimension of the array ZSPS.
(REAL array, dimensioned LRWK, input/output) is the real work array.
(INTEGER, input) is the length of RWRK.
(INTEGER array, dimensioned LIWK, input/output) is the integer work array.
(INTEGER, input) is the length of IWRK.
(REAL array, dimensioned LZDT, output) is the array in which the interpolated "dense" array of data is to be returned. The dimensions of the interpolated array may be supplied by the user or determined by Conpack, depending on the value of the parameter 'ZDS'. Note that, if the size of the dense array is not a product of the size of the sparse array and some perfect square, the aspect ratio of the dense grid may be slightly different from that of the sparse grid.
(INTEGER, input) is the length of ZDAT.
The C-binding argument descriptions are the same as the FORTRAN argument descriptions with the following exceptions:
Dimensioned l by ksps, where l \(>= nsps.
The second dimension of the array zsps.
The second dimension of the sparse array of data in zsps. msps \(<= ksps.
The first dimension of the sparse array of data in zsps. nsps \(<= l, where l is the declared first dimension of the array zsps.
CPSPS1 performs the same functions as CPRECT, but, in addition, it interpolates from a sparse array of data to a dense array of data. CPSPS1 does this by using the routines BSURF1 and BSURF2, from the package Fitpack, by Alan K. Cline, to fit bicubic splines under tension to the sparse array of data and to compute the dense grid of data that is returned to you. The tension on the spline surfaces is specified by the parameter 'T3D'. By default, CPSPS1 selects the dimensions of the dense array of data; if desired, you can specify these dimensions by setting the parameter 'ZDS' non-zero and the parameters 'ZD1', 'ZDM', and 'ZDN' to the desired values. In either case, once 'ZD1', 'ZDM', and 'ZDN' are set, they should not be reset by you until the contour plot is complete and a different contour plot is to be drawn.
Because the routines BSURF1 and BSURF2 do not have a built-in special value feature, if the special value parameter 'SPV' is set non-zero and the sparse array contains occurrences of that value, special action must be taken. The indices of the special values in the sparse array are saved in a part of the integer workspace array; the special values are then replaced by values interpolated from adjacent grid points and the resulting array is used to obtain the dense array; then, the special values in the sparse array are restored and the corresponding elements of the dense array are also given the special value.
To use CPSPS1 or c_cpsps1, load the NCAR Graphics libraries ncarg, ncarg_gks, and ncarg_c, preferably in that order.
See the conpack man page for a description of all Conpack error messages and/or informational messages.
Online: conpack, cpback, cpchcf, cpchcl, cpchhl, cpchil, cpchll, cpcica, cpclam, cpcldm, cpcldr, cpcltr, cpcnrc, cpdrpl, cpezct, cpgetc, cpgeti, cpgetr, cplbam, cplbdr, cpmpxy, cpmviw, cpmvrw, cppkcl, cppklb, cprect, cprset, cpscae, cpsetc, cpseti, cpsetr, cpsps2, ncarg_cbind
Hardcopy: NCAR Graphics Contouring and Mapping Tutorial
Copyright (C) 1987-2009
University Corporation for Atmospheric Research
The use of this Software is governed by a License Agreement.