cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c                                                                      c
c      Drawcgm_test2.For v. 1.0 (12/87)                                c
c      Copyright 1987 Pittsburgh Supercomputing Center                 c
c      Authors Joel Welling and Jonathan Goldick                       c
c      This is one of two test routines for Drawcgm.for                c
c                                                                      c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C       IMPLICIT NONE                                         
        INTEGER Nxsmall,Nysmall,Nxbig,Nybig,Nxmed,Nymed,Nxmask,Nymask
        PARAMETER (Nxsmall=10,Nysmall=10,Nxbig=300,Nybig=300)
        PARAMETER (Nxmed=150,Nymed=150)
        PARAMETER (Nxmask=3,Nymask=3)
        REAL Xmin,Ymin,Xmax,Ymax,Xmincb,Xmaxcb,Xlabel,Ylabel
        PARAMETER (Xmin= 0.1,Ymin=0.2,Xmax= 0.7,Ymax=0.8)
        PARAMETER (Xmincb=0.78, Xmaxcb=0.88, Xlabel=0.1, Ylabel=0.93)
        REAL Xminpt,Yminpt,Xmaxpt,Ymaxpt,Xminhb,Yminhb,Xmaxhb,Ymaxhb
        PARAMETER (Xminpt= 0.1,Yminpt=0.5,Xmaxpt= 0.8,Ymaxpt=0.9)
        PARAMETER (Xminhb= 0.1,Yminhb=0.2,Xmaxhb= 0.8,Ymaxhb=0.45)
        REAL Xmnmed1,Ymnmed1,Xmxmed1,Ymxmed1
        PARAMETER (Xmnmed1= 0.1,Ymnmed1=0.5,Xmxmed1=0.4,Ymxmed1=0.8)
        REAL Xmnmed2,Ymnmed2,Xmxmed2,Ymxmed2
        PARAMETER (Xmnmed2= 0.4,Ymnmed2=0.5,Xmxmed2=0.7,Ymxmed2=0.8)
        REAL Xmnmed3,Ymnmed3,Xmxmed3,Ymxmed3
        PARAMETER (Xmnmed3= 0.1,Ymnmed3=0.2,Xmxmed3=0.4,Ymxmed3=0.5)
        REAL Xmnmed4,Ymnmed4,Xmxmed4,Ymxmed4
        PARAMETER (Xmnmed4= 0.45,Ymnmed4=0.15,Xmxmed4=0.75,Ymxmed4=0.45)
        INTEGER Maxclr,Minclr,Nclr
        PARAMETER (Maxclr=247,Minclr=20,Nclr=248)
        INTEGER Ipixel(Nxbig,Nybig),Mfield(Nxmask,Nymask)
        REAL Input(Nxsmall,Nysmall),Pixel(Nxbig,Nybig)
        EXTERNAL Getdat,Strlin,Rtoint,Drwpix
        DATA Mfield/1,1,1,1,0,0,1,1,1/
*
*       Get the data.  (This routine is in this module)
*
        CALL Device('tek')
        CALL Getdat(Input,Nxsmall,Nysmall)
*
*       Stretch and quantize it, then draw a single image with a single 
*       call to Drwpix.
*
        CALL Strlin(Input,Nxsmall,Nysmall,Pixel,Nxbig,Nybig)
        CALL Rtoint(Pixel,Ipixel,Nxbig,Nybig,0.0,1.0,Minclr,Maxclr)
        CALL Drwpix(2,Xmin,Ymin,Xmax,Ymax,Ipixel,Nxbig,Nybig,2)

        STOP
        END

        SUBROUTINE Getdat(Input,Nxdim,Nydim)
C       IMPLICIT NONE 
        INTEGER Nxdim,Nydim,I,J                     
        REAL Input(Nxdim,Nydim)
        INTRINSIC Float,Mod
        DO 10 I= 1,Nxdim
                DO 10 J= 1,Nydim
                        Input(i,j)= Float(I+J)/(2.0*Nxdim)   
*                        Input(i,j)= float(Mod((i*j),Nxdim))/Nxdim
10                      CONTINUE
        RETURN
        END