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:   call units_open(trim(mode)//'.log','unknown',ulog,err)
    if(err /= 0) then
       if (iproc == 0) write(0,*) "error opening file ", trim(mode)//'.log'
       call mpi_wrapper_finalize(err)
       stop
    end if
  
  
negative.code:           prm,fmodel,data_format,err)
       if (err /= 0) then
          if (iproc == 0) call dump_error(err,'')
          stop
       end if
       close(uprm)
  