anchor.diff:  use const, only: i0, lint, p
   
   implicit none
  +
  +! TEMPORARY ONLY
  +integer :: D0_normalisation
   
   contains
   
  
positive.code:                 call convert_excitor_to_determinant(cdet%f, excitation_level, cluster_population)
  
                  ! Normalisation factor for amplitudes...
                  amplitude = real(cluster_population,p)/(D0_population**(cluster_size-1)) ! TODO: don't accumulate D0 over report loops
  
                  ! We chose excitors uniformly.
                  !  -> prob. of each excitor is 1./(number of excitors on processor).
  
negative.code:             ! Must be the reference.
              cdet%f = f0
              excitation_level = 0
              amplitude = D0_normalisation ! TODO: don't accumulate D0 over report loops
              ! Only one cluster of this size to choose => p_clust = 1
          case default
              ! Select cluster from the excitors on the current processor in
  