anchor.diff:         rewind(udata)
   
          ! allocate memory for data
  -       allocate(seqs(nvars,nseqs),stat=err)
  +       allocate(seqs(nvars, nseqs), stat=err)
   
          ! read data
  -       do i = 1,nseqs
  -          read(udata,*,iostat=err) seqs(:,i)
  +       do i = 1, nseqs
  +          read(udata, *, iostat=err) seqs(:, i)
             if(err > 0) then 
                error_code = 20
                return
  
positive.code:        ! count data lines
         nseqs = 0
         do
            read(udata,'(a)',iostat=err) line
            ! exit at end of file or empty line
            if( err < 0 .or. len_trim(line) == 0) exit
            nseqs = nseqs + 1
  
negative.code: 
      end if
  
      if (wid > 0.0_kflt) then 
  
         call data_reweight(seqs, wid, iproc)
  
      end if
  
  