!  drawcgm_prb2.f90  11 June 2000
!
program drawcgm_prb2
!
!*******************************************************************************
!
!! DRAWCGM_PRB2 carries out a test of some DRAWCGM routines.
!
  integer, parameter :: nxsmall=10
  integer, parameter :: nysmall=10
  integer, parameter :: nxbig=300
  integer, parameter :: nybig=300
  integer, parameter :: nxmask=3
  integer, parameter :: nymask=3
!
  character ( len = 8 ) date
  character ( len = 20 ) dev
  integer i
  real input(nxsmall,nysmall)
  integer ipixel(nxbig,nybig)
  integer j
  integer maxclr
  integer minclr
  integer mfield(nxmask,nymask)
  real pixel(nxbig,nybig)
  character ( len = 10 ) time
  real xmax
  real xmin
  real ymax
  real ymin
!
  data mfield/1,1,1,1,0,0,1,1,1/
!
  call date_and_time ( date, time )
  write ( *, * ) ' '
  write ( *, * ) 'DRAWCGM_PRB2'
  write ( *, * ) '  A sample problem for DRAWCGM.'
  write ( *, * ) ' '
  write ( *, * ) '  Today''s date: ', date
  write ( *, * ) '  Today''s time: ', time
!
!  Specify the output file type.
!
  dev = 'ps'

  call device ( dev )

  write ( *, * ) 'DEBUG - call OUTFIL'

  call outfil ( 'drawcgm_prb2.ps' )

  write ( *, * ) 'CALLED OUTFIL'
!
!  Set the data.
!
  do i = 1, nxsmall
    do j = 1, nysmall
      input(i,j) = real ( i + j ) / ( 2.0 * nxsmall )
    end do
  end do
!
!  Stretch the data.
!
  call rmat_expand_linear ( input, nxsmall, nysmall, pixel, nxbig, nybig )
!
!  Convert the real data to integer values.
!
  maxclr = 247
  minclr = 20

  call rtoint ( pixel, ipixel, nxbig, nybig, 0.0, 1.0, minclr, maxclr )
!
!  Draw the cells.
!
  xmin = 0.1
  ymin = 0.2
  xmax = 0.7
  ymax = 0.8

  call drwpix ( 2, xmin, ymin, xmax, ymax, ipixel, nxbig, nybig, 2 )

  write ( *, * ) ' '
  write ( *, * ) 'DRAWCGM_PRB2'
  write ( *, * ) '  Normal end of the sample problem for DRAWCGM.'

  stop
end