anchor.diff:    call units_open_unf(rst_file,'old',unt,err)
     
     ! read binary file
  +  read(unt) data_type
     read(unt) data_format
     read(unt) nvars
     read(unt) nclasses
  
positive.code:   read(unt) prm
  
    ! dump file
    write(*,'(a)') '# <data_format> <nvars> <nclasses> <nseq>'
    write(*,'(a,3(1x,i4))') trim(data_format), nvars, nclasses, np
    if (np > 0) then
       do p = 1,np
          write(*,'(a,1000(1x,i2))') '# ',seq(:,p)
  
negative.code:   use arguments, only: read_opt,read_arg
    
    integer                    :: unt=11
    character(len=string_size) :: data_type
    character(len=string_size) :: data_format
    integer                    :: nvars,nclasses
    integer,       allocatable :: seq(:,:)
  