program gfm1
      implicit none
C
C     type functions referenced
      integer pgopen
C
C     local variables
C
      integer istat
      integer i,j
      real a(10,10)
      real c(5), tr(6)
C
C     ok - ready to start
C
      istat= pgopen('plot.ps/PS')
      if (istat.le.0) stop
C
C     initialise an array
C
      do i = 1,10
      do j =1,10
        a(i,j)=sqrt (float(abs((i-6)*j*(j-4))))
      enddo
      enddo
      write (6,10) a
10    format (10f6.2)
C
C     establish view port
C
      call pgenv(0.0,10.0, 0.0,10.0 , 1, 1)
C
C     draw a circle in middle of field
C
      call pgcirc (4.0,6.0,2.0)
C
C
C     do a contour plot
C
      c(1)=1.0
      c(2)=2.0
      c(3)=3.0
      c(4)=4.0
      c(5)=5.0
      tr(1)=0
      tr(2)=1
      tr(3)=0
      tr(4)=0
      tr(5)=0
      tr(6)=1
      call pgcons (a,10,10, 1,10, 1,10, C,5, TR )
C
C     close device
C
      call pgclos
      end