PGPLOT Module: PGVECT -- vector map of a 2D data array, with blanking


SUBROUTINE PGVECT (A, B, IDIM, JDIM, I1, I2, J1, J2, C, NC, TR,
1 BLANK)
INTEGER IDIM, JDIM, I1, I2, J1, J2, NC
REAL A(IDIM,JDIM), B(IDIM, JDIM), TR(6), BLANK, C

Draw a vector map of two arrays. This routine is similar to
PGCONB in that array elements that have the "magic value" defined by
the argument BLANK are ignored, making gaps in the vector map. The
routine may be useful for data measured on most but not all of the
points of a grid. Vectors are displayed as arrows; the style of the
arrowhead can be set with routine PGSAH, and the the size of the
arrowhead is determined by the current character size, set by PGSCH.


Arguments
A (input) : horizontal component data array.
B (input) : vertical component data array.
IDIM (input) : first dimension of A and B.
JDIM (input) : second dimension of A and B.
I1,I2 (input) : range of first index to be mapped (inclusive).
J1,J2 (input) : range of second index to be mapped (inclusive).
C (input) : scale factor for vector lengths, if 0.0, C will be
set so that the longest vector is equal to the
smaller of TR(2)+TR(3) and TR(5)+TR(6).
NC (input) : vector positioning code.
<0 vector head positioned on coordinates
>0 vector base positioned on coordinates
=0 vector centered on the coordinates
TR (input) : array defining a transformation between the I,J
grid of the array and the world coordinates. The
world coordinates of the array point A(I,J) are
given by:
X = TR(1) + TR(2)*I + TR(3)*J
Y = TR(4) + TR(5)*I + TR(6)*J
Usually TR(3) and TR(5) are zero - unless the
coordinate transformation involves a rotation
or shear.
BLANK (input) : elements of arrays A or B that are exactly equal to
this value are ignored (blanked).