anchor.diff:         ! count data lines
          nseqs = 0
          do
  -          read(udata,'(a)',iostat=err) line
  +          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
  
positive.code:        rewind(udata)
  
         ! allocate memory for data
         allocate(seqs(nvars,nseqs),stat=err)
  
         ! read data
         do i = 1,nseqs
            read(udata,*,iostat=err) seqs(:,i)
            if(err > 0) then 
               error_code = 20
               return
  
negative.code:     case('int')
         
         ! get number of fields
         read(udata, '(a)', iostat=err) line
         call parser_nfields(line, newline, nfields)
         rewind(udata)
         
         ! set nvars
  