subroutine contour(a,x,y,dim1,dim2,title)
c
c Make a contour plot of the array a, an array with dimensions (dim1,dim2),
c with an x-axis specified by x(dim1) and a y-axis specified by y(dim2).
c Place a title on the graph.
c 
c
c These arrays should be dimensioned such that lz > dim1, mz > dim2
c 
      real a(lz,mz),a2(lz*mz)
      real xdata(lz*mz),ydata(lz*mz)
      real work(13*lz*mz),scrarr(100**2)
      integer iwork(32*lz*mz),dim1,dim2,index
      character*80 title

      index=0
      do 30 l=1,dim1
      do 30 m=1,dim2
      index=index+1
      a2(index)=a(l,m)
      xdata(index)=x(l)
      ydata(index)=y(m)
 30   continue
c
      call set(0.0,1.0,0.0,1.0,0.0,1.0,0.0,1.0,1) 
      call plchlq(.5,.95,title,20.,0.,0.)
c
c Make the contour lines dashed if a < 0:
c
      call conop4('DAS=LSS','$$''''$$''''$''',0,0)
c
c Turn off the annoying message at the point if desired:
c
c      call conop1('MES=OFF')
c
c Label local maxima and minima if desired: 
c
c      call conop1('PMM=ON')
c
c Turn up the resolution (invisible array points) to 100.  
c The default is 40.
c
     call conop2('SSZ=ON',100)
c
c Make the plot. 
c
      call conran(xdata,ydata,a2,dim1*dim2,work,iwork,scrarr)
c
c Draw a box.
c
      call set(0.0,1.0,0.0,1.0,0.0,1.0,0.0,1.0,1) 
      CALL PLOTIT (    0,    0,0)
      CALL PLOTIT (32767,    0,1)
      CALL PLOTIT (32767,32767,1)
      CALL PLOTIT (    0,32767,1)
      CALL PLOTIT (    0,    0,1)
      call frame
c
      return
      end