anchor.diff:  
       end subroutine cumulative_population
   
  -    pure function get_pop_contrib_real(pops, real_factor) result(contrib)
  +    pure function get_pop_contrib(pops, real_factor, complx) result(contrib)
   
           ! Get contribution from a given position in a population list,
           ! given real or complex.
  
positive.code:             if (d0_pos > 1) cumulative_pops(d0_pos) = cumulative_pops(d0_pos-1)
              do i = d0_pos+1, nactive
                  cumulative_pops(i) = cumulative_pops(i-1) + &
                                          nint(get_pop_contrib_ptr(pops(:,i), real_factor))
              end do
          else
              ! V simple on other processors: no reference to get in the way!
              do i = 2, nactive
                  cumulative_pops(i) = cumulative_pops(i-1) + &
                                          nint(get_pop_contrib_ptr(pops(:,i), real_factor))
              end do
          end if
          if (nactive > 0) then
  
negative.code: use const, only: i0, p, int_p
  
  implicit none
  
  contains
  
  