anchor.diff:    ci = 0
     do i=1,N_det_generators
       if(ltbc(i)) then
  -      ci += 1
  +      ci = ci+1
         tbc(ci) = i
       end if
     end do
  
positive.code:       pt2_detail(:, index(i)) += pt2_mwen(:,i)
        parts_to_get(index(i)) -= 1
        if(parts_to_get(index(i)) < 0) then 
          stop "PARTS ??"
        end if
        if(parts_to_get(index(i)) == 0) actually_computed(index(i)) = .true.
  
negative.code:   h = N_det-1
  
    do while(h >= l)
      i = ishft(h+l,-1)
      if(w(i+1) > v) then
        h = i-1
      else
  