anchor.diff:        else
           if (Nabove(tooth) > Nabove_old) then
             print '(G10.3, 2X, F16.10, 2X, G16.3, 2X, F16.4, A20)', Nabove(tooth), avg+E, eqt, time-time0, ''
  -!print "(4(G23.13), 4(I9))", time - time0, avg, eqt, Nabove(tooth), tooth, first_det_of_teeth(tooth)-1, done, first_det_of_teeth(tooth+1)-first_det_of_teeth(tooth)
             Nabove_old = Nabove(tooth)
           endif
         endif
  -!print "(4(G22.13), 4(I9))", time - time0, avg, eqt, Nabove(tooth), tooth, first_det_of_teeth(tooth)-1, done, first_det_of_teeth(tooth+1)-first_det_of_teeth(tooth)
       end if
     end do pullLoop
   
  
positive.code:       if(Nabove(1) < 5d0) cycle
        call get_first_tooth(actually_computed, tooth)
       
        done = 0
        do i=first_det_of_teeth(tooth), first_det_of_teeth(tooth+1)-1
          if(actually_computed(i)) done = done + 1
        end do
  
        E0 = sum(pt2_detail(1,:first_det_of_teeth(tooth)-1))
        if (tooth <= comb_teeth) then
          prop = ((1d0 - dfloat(comb_teeth - tooth + 1) * comb_step) - pt2_cweight(first_det_of_teeth(tooth)-1))
  
negative.code: end subroutine
  
  
  
  subroutine get_comb(stato, dets, ct)
    implicit none
  