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:     case('int')
         
         ! get number of fields
         read(udata,'(a)',iostat=err) line
         call parser_nfields(line,newline,nfields)
         rewind(udata)
         
         ! set nvars
  
negative.code: 
      end if
  
      if (wid > 0.0_kflt) then 
  
         call data_reweight(seqs, wid, iproc)
  
      end if
  
  