program profil
      character filein*10
      integer i,j,nx,ny,k,kunit,kx(4)
      double precision a(258,256)
      kunit = 20
      kx(1) = 33
      kx(2) = 65
      kx(3) = 97
      kx(4) = 129
      nx = 258
      ny = 256
      print *,'input data filename'
      read(5,'(a)') filein
      open(11,file=filein,status='old')
      read(11,*) ((a(i,j),i=1,nx),j=1,ny)
      do 5 k=1,4
         kunit = kunit + 1
         write(kunit,101) (dfloat(j),a(kx(k),j),j=1,ny)
 5    continue
 101  format(2(f16.10,2x))
      end