subroutine output(pf,u,v,nx,ny,kunit)
      integer i,j,nx,ny,kunit
      double precision pf(0:nx+1,0:ny+1),u(nx,ny),v(nx,ny)
c      write(21,101) ((u(i,j),i=1,nx),j=1,ny)
c      write(21,101) ((v(i,j),i=1,nx),j=1,ny)
      print *,'kunit=',kunit
      write(kunit,*) pf
 101  format(e16.10)
 102  format(2(e16.10,2x))
 103  format(/)
      end