anchor.diff:    double precision, allocatable :: val(:)
     integer(bit_kind), allocatable :: det(:,:,:)
     integer, allocatable :: task_id(:)
  -  integer :: done, Nindex
  +  integer :: Nindex
     integer, allocatable :: index(:)
     double precision, save :: time0 = -1.d0
     double precision :: time, timeLast, Nabove_old
  
positive.code:       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
  
  
negative.code:     end if
    end do
  
  end subroutine
  
  
  