program main

!*****************************************************************************80
!
!! MAIN is the main program for DRAWCGM_PRB2.
!
!  Discussion:
!
!    DRAWCGM_PRB2 carries out a test of some DRAWCGM routines.
!
!  Modified:
!
!    25 September 2002
!
  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 = 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)
  real xmax
  real xmin
  real ymax
  real ymin

  data mfield/1,1,1,1,0,0,1,1,1/

  call timestamp ( )

  write ( *, * ) ' '
  write ( *, * ) 'DRAWCGM_PRB2'
  write ( *, '(a)' ) '  FORTRAN90 version'
  write ( *, * ) '  A sample problem for DRAWCGM.'
!
!  Specify the output file type.
!
  dev = 'cgmb'

  call device ( dev )

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

  call outfil ( 'drawcgm_prb2.cgmb' )

  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 ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'DRAWCGM_PRB2'
  write ( *, '(a)' ) '  Normal end of execution.'

  write ( *, '(a)' ) ' '
  call timestamp ( )

  stop
end