subroutine wframe
c
      implicit none
      include 'itg.par'
      include 'itg.cmn'
      complex phi1(nzz,mz*lz)
      real phi2(mzz,nzz*lz),amin,amax
      integer i,m,n,istat,DFR8addimage
      character*1 image(malias,nalias)

      amin=0.
      amax=0.      
      call xspacec(potential,phi1)
      call yspace(phi1,phi2)
c
c Assume that we are interested in the l=1 point (theta=0):
c
      do n=1,nalias
         do m=1,malias
            amin=min(amin,phi2(m,n))
            amax=max(amin,phi2(m,n))
         enddo
      enddo
         
      if(amin.eq.amax) then
         image=char(0)
      else
         do n=1,nalias
            do m=1,malias
               image(m,n)=char(int(255*(phi2(m,n)-amin)/(amax-amin)))
            enddo
         enddo
      endif

      if(ihdf.eq.1) then
         istat=DFR8addimage('phi.hdf',image,malias,nalias,0)
      else
         write(20,20)((phi2(m,n),m=1,malias),n=1,nalias*ldb)
      endif

 20   format(1x,3e22.14)
      return 
      end