anchor.diff:    use constants
     use fasta,       only : fasta_read
     use errors
  -  use mpi_wrapper
     use dump,        only: read_rst,read_prm_unit,dump_energies,dump_rst
     use eval_command_line
     use units,       only: units_initialize,units_open
  
positive.code:      stop
    end if
  
    !================================================ init. mpi
  
    call mpi_wrapper_initialize(err)
  
    !================================================ init. the random number generator
  
    call random_initialize(rseed,iproc)
  
negative.code: 
       if (err /= 0) then
          if (iproc == 0) call dump_error(err,err_string)
          stop
       end if
  
  