anchor.diff:          write_restart_shift = restart_in%write_restart_shift
           call init_restart_info_t(ri, write_id=restart_in%write_id)
   
  +        ! If averaging over spin spin polarisation will change in general. Copy
  +        ! system to prevent modification to original system type.
           sys_copy = sys
   
           outer_loop: do beta_cycle = 1, dmqmc_in%beta_loops
  
positive.code: 
              ! Distribute psips uniformly along the diagonal of the density
              ! matrix.
              call create_initial_density_matrix(rng, sys, qmc_in, dmqmc_in, qs%ref, annihilation_flags, &
                                                 init_tot_nparticles, qs%psip_list, qs%spawn_store%spawn)
  
              ! Allow the shift to vary from the very start of the beta loop, if
  
negative.code:         ! times, to accumulate statistics for each value for beta.
  
          ! In:
          !    sys: system being studied.
          !    restart_in: input options for HDF5 restart files.
          !    load_bal_in: input options for load balancing.
          !    reference_in: current reference determinant.  If not set (ie
          !       components allocated) then a best guess is made based upon the
          !       desired spin/symmetry.
          ! In/Out:
          !    qmc_in: input options relating to QMC methods.
          !    dmqmc_in: input options relating to DMQMC.
          !    dmqmc_estimates: type containing all DMQMC estimates.
  