anchor.diff:    character(len=string_size) :: mode         ! run mode: EVAL, SIM, LEARN
     integer                    :: rseed        ! random seed
     real(kflt)                 :: beta         ! temperature of the run
  +  integer                         :: iproc=0,nproc=1
     integer                         :: err
     character(long_string_size)     :: err_string
     integer                         :: dim1,dim2
  
positive.code:      call read_rst(urst,data_format,nvars,nclasses,iproc,nproc,seq,seqs_table,prm,fmodel,err)
       if (err /= 0) then
          if (iproc == 0) call dump_error(err,'')
          call mpi_wrapper_finalize(err)
          stop
       end if
       close(urst)
  
negative.code:   call units_open(trim(mode)//'.log','unknown',ulog,err)
    if(err /= 0) then
       if (iproc == 0) write(0,*) "error opening file ", trim(mode)//'.log'
       stop
    end if
  
  