%N grpixl %F $MIR/src/grpixl.f
%D Solid-fill multiple rectangular areas
%B
      SUBROUTINE GRPIXL (IA, IDIM, JDIM, I1, I2, J1, J2, 
     1                   X1, X2, Y1, Y2)
      INTEGER IDIM, JDIM, I1, I2, J1, J2
      INTEGER IA(IDIM,JDIM)
      REAL    X1, X2, Y1, Y2

 Determine the size of each rectangular element. If it is equal
 to the device pen width and the device supports pixel primitives,
 use pixel primitives. Otherwise, if the size is smaller than the
 device pen width emulate pixel output by plotting points. If the
 size is larger than the device pen width, emulate by outputting
 solid-filled rectangles.

 Arguments:
  IA     (input)  : the array to be plotted.
  IDIM   (input)  : the first dimension of array A.
  JDIM   (input)  : the second dimension of array A.
  I1, I2 (input)  : the inclusive range of the first index
                    (I) to be plotted.
  J1, J2 (input)  : the inclusive range of the second
                    index (J) to be plotted.
  X1, Y1 (input)  : world coordinates of one corner of the output
                    region
  X2, Y2 (input)  : world coordinates of the opposite corner of the
                    output region
