anchor.diff:    integer,    allocatable :: seqs0(:,:)
     real(kflt), allocatable :: prm(:)       ! parameters array
     real(kflt), allocatable :: fmodel(:)    ! model frequencies
  -  real(kflt), allocatable :: fdata(:)     ! empirical frequencies
  -  real(kflt), allocatable :: scores(:,:)  ! final scores
     integer, allocatable    :: seqs_table(:,:)
     ! command line parameters
     integer                    :: udata        ! data unit
  
positive.code: 
  
    !================================================ print a header
  
    if ( iproc == 0 ) then
       write(ulog,'(a)')         '#'
       write(ulog,'(a)')         '#==========================================='
  
negative.code:      write(ulog,'(a,1x,i8)')   '#  n. of variables              = ', nvars
       write(ulog,'(a,1x,i8)')   '#  n. of classes                = ', nclasses
       write(ulog,'(a,1x,i8)')   '#  n. of seqs                   = ', nseqs
       write(ulog,'(a)')         '#'
       write(ulog,'(a)')         '#==========================================='
       write(ulog,'(a)')         '#'
  