anchor.diff:              ! Calculate and output all requested estimators based on the reduced
               ! density matrix. This is for ground-state RDMs only.
               if (dmqmc_in%rdm%calc_ground_rdm) call call_ground_rdm_procedures(dmqmc_estimates, beta_cycle, dmqmc_in%rdm)
  -            ! Calculate and output new weights based on the psip distirubtion in
  +            ! Calculate and output new weights based on the psip distriubtion in
               ! the previous loop.
               if (dmqmc_in%find_weights) call output_and_alter_weights(dmqmc_in, sys%max_number_excitations, &
                                                                        dmqmc_estimates%excit_dist, weighted_sampling)
  
positive.code: 
              do ireport = 1, qmc_in%nreport
  
                  call init_dmqmc_report_loop(bloom_stats, dmqmc_estimates, qs%spawn_store%rspawn)
                  tot_nparticles_old = qs%psip_list%tot_nparticles
  
                  do icycle = 1, qmc_in%ncycles
  
negative.code: 
      end subroutine init_dmqmc_beta_loop
  
      subroutine init_dmqmc_report_loop(calc_excit_dist, bloom_stats, dmqmc_estimates, rspawn)
  
          ! Initialise a report loop (basically zero quantities accumulated over
          ! a report loop).
  
          ! In:
          !    calc_excit_dist: true if the excitation distribution is being
          !        calculated at each report loop.
          ! In/Out:
          !    bloom_stats: type containing information regarding blooming events.
          !    dmqmc_estimates: type containing estimates of observables.
  