anchor.diff:        iw2(imax)=i
         iwS(i)=imax; wSL(i)=-1.d0/vmax
       enddo
  -!!$    write(*,*) 'iw2:'
  +!!$    write(0,*) 'iw2:'
   !!$    do i=1,hecMAT%N*ndof
  -!!$      if (iw2(i) > 0) write(*,*) i, iw2(i), iw1L(i), iw1U(i)
  +!!$      if (iw2(i) > 0) write(0,*) i, iw2(i), iw1L(i), iw1U(i)
   !!$    enddo
  -!!$    write(*,*) 'iwS:'
  -!!$    write(*,*) iwS(:)
  +!!$    write(0,*) 'iwS:'
  +!!$    write(0,*) iwS(:)
   
       deallocate(iw1L, iw1U)
     end subroutine choose_slaves
  
positive.code:       enddo
      enddo
      if (Ttmat%nnz /= Ttmat%index(Ttmat%n)) then
        write(*,*) Ttmat%nnz, Ttmat%index(Ttmat%n)
        stop 'ERROR: Ttmat%nnz wrong'
      endif
      ! item and A
  
negative.code:             icnt=icnt+1
              mm=me+1
              BTtKT%item(mm)=ll
              !if (i==1) write(0,*) 'd', mm, jj, kk, ll
            endif
            TtKTp=>BTtKT%A((mm-1)*ndof2+1:mm*ndof2)
            call blk_trimatmul_add(ndof, Ttp, Kp, Tp, TtKTp)
  