anchor.diff:      ! Currently only implemented for DMQMC.
       logical :: replica_tricks = .false.
   
  -    ! [review] - NSB: This is used both for excit_distribution and ground-state
  -    ! [review] - NSB: I should create a new variable for the ground-state
  -    ! [review] - NSB: RDMs and add it to that derived type.
  -    ! When calculating certain DMQMC properties, we only want to start
  -    ! averaging once the ground state is reached. The below integer is input
  -    ! by the user, and gives the iteration at which data should start being
  -    ! accumulated for the quantity. This is currently only used for the
  -    ! reduced density matrix and calculating importance sampling weights.
  -    integer :: start_averaging = 0
  +    ! When performing a ground-state RDM calculation, on what iteration do we
  +    ! start accumulating the ground-state RDM?
  +    integer :: start_av_rdm = 0
  +
  +    ! When calculating the distributon of particles across excitation levels,
  +    ! on what iteration do we start performing this calculation?
  +    integer :: start_av_excit_dist = 0
   
       ! If this logical is true then the program runs the DMQMC algorithm with
       ! importance sampling.
  
positive.code:     ! to be calculated. Each element of this array corresponds to one of these RDMs.
      type(rdm_t), allocatable :: rdms(:) ! nrdms
  
      ! rdm_traces(i,j) holds the trace of replica i of the rdm with label j.
      real(p), allocatable :: rdm_traces(:,:) ! (sampling_size, nrdms)
  
      type(rdm_spawn_t), allocatable :: rdm_spawn(:) ! nrdms
  
      ! The length of the spawning array for RDMs. Each RDM calculated has the same
      ! length array.
      integer :: spawned_rdm_length
  end type dmqmc_inst_rdms_t
  
  
  
negative.code:     ! are reset on each processor to start the next report loop.
      real(p), allocatable :: estimator_numerators(:) !(number_dmqmc_estimators)
  
      ! This array is used to hold the number of particles on each excitation
      ! level of the density matrix.
      real(p), allocatable :: excit_distribution(:) ! (0:max_number_excitations)
  
      ! [todo] - (NSB) remove these two - I doubt they'll ever be used.
  