anchor.diff:          character(255) :: restart_file, junk
           integer :: io, i
           logical :: exists
  +        integer :: restart_version
   
           if (parent) then
               io = get_free_unit()
  
positive.code:             open(io, file=restart_file)
  
              read (io,*) junk
              read (io,*) mc_cycles_done
              read (io,*) junk
              read (io,*) nparticles_old_restart, shift, vary_shift
  
negative.code: 
              open(io, file=restart_file)
  
              write (io,*) 'Restart version'
              write (io,*) restart_version
              write (io,*) '# number of cycles'
              write (io,*) nmc_cycles
              write (io,*) '# shift'
  