anchor.diff:                icnt=icnt+1
                 mm=me+1
                 BTtKT%item(mm)=ll
  -              !if (i==1) write(*,*) 'u', mm, jj, kk, ll
  +              !if (i==1) write(0,*) 'u', mm, jj, kk, ll
               endif
               TtKTp=>BTtKT%A((mm-1)*ndof2+1:mm*ndof2)
               call blk_trimatmul_add(ndof, Ttp, Kp, Tp, TtKTp)
  
positive.code:             icnt=icnt+1
              mm=me+1
              BTtKT%item(mm)=ll
              !if (i==1) write(*,*) 'd', mm, jj, kk, ll
            endif
            TtKTp=>BTtKT%A((mm-1)*ndof2+1:mm*ndof2)
            call blk_trimatmul_add(ndof, Ttp, Kp, Tp, TtKTp)
  
negative.code:             enddo
              icnt=icnt+1
              iw(icnt)=ll
              !if (i==1) write(0,*) 'u', icnt, jj, kk, ll
            enddo ll3
          enddo
        enddo
        if (icnt == 0) icnt=1
        !if (i==1) write(0,*) iw(1:icnt)
        BTtKT%index(i)=BTtKT%index(i-1)+icnt
      enddo
      !write(0,*) BTtKT%index(1:n)-BTtKT%index(0:n-1)
  
      BTtKT%nnz=BTtKT%index(n)
      allocate(BTtKT%item(BTtKT%nnz))
  