anchor.diff:    subroutine S6UM1(Ic,Xlnzr,Colno,Zln,Temp,Indx)
       implicit none
       !------
  -    integer, intent(in):: Ic
  -    integer, intent(in):: Xlnzr(*)
  -    integer, intent(in):: Colno(*)
  -    integer, intent(inout):: Indx(*)
  -    real(kind=8), intent(inout):: Temp(9,*)
  -    real(kind=8), intent(inout):: Zln(9,*)
  -    !------
  -    integer:: j
  -    integer:: jc
  -    integer:: jj
  -    integer:: k
  -    integer:: ke
  -    integer:: ks
  -    integer:: l
  -    real(kind=8):: s(9)
  +    integer(kind=kint), intent(in):: Ic
  +    integer(kind=kint), intent(in):: Xlnzr(*)
  +    integer(kind=kint), intent(in):: Colno(*)
  +    integer(kind=kint), intent(inout):: Indx(*)
  +    real(kind=kreal), intent(inout):: Temp(9,*)
  +    real(kind=kreal), intent(inout):: Zln(9,*)
  +    !------
  +    integer(kind=kint):: j
  +    integer(kind=kint):: jc
  +    integer(kind=kint):: jj
  +    integer(kind=kint):: k
  +    integer(kind=kint):: ke
  +    integer(kind=kint):: ks
  +    integer(kind=kint):: l
  +    real(kind=kreal):: s(9)
   
       ks = Xlnzr(Ic)
       ke = Xlnzr(Ic+1)
  
positive.code:   subroutine QMDMRG(Xadj,Adjncy,Deg,Qsize,Qlink,Marker,Deg0,Nhdsze,Nbrhd,Rchset,Ovrlp)
      implicit none
      !------
      integer, intent(in):: Deg0
      integer, intent(in):: Nhdsze
      integer, intent(in):: Adjncy(*)
      integer, intent(in):: Nbrhd(*)
      integer, intent(in):: Xadj(*)
      integer, intent(inout):: Deg(*)
      integer, intent(inout):: Qsize(*)
      integer, intent(inout):: Qlink(*)
      integer, intent(inout):: Marker(*)
      integer, intent(out):: Rchset(*)
      integer, intent(out):: Ovrlp(*)
      !------
      integer:: deg1
      integer:: head
      integer:: inhd
      integer:: iov
      integer:: irch
      integer:: j
      integer:: jstrt
      integer:: jstop
      integer:: link
      integer:: lnode
      integer:: mark
      integer:: mrgsze
      integer:: nabor
      integer:: node
      integer:: novrlp
      integer:: rchsze
      integer:: root
  
      if ( Nhdsze<=0 ) return
      do inhd = 1, Nhdsze
  
negative.code:   !> @brief SETIJ
    !======================================================================!
    subroutine SETIJ(hecMESH,hecMAT,NEQns,NDEg,NTTbr,ISYm,JCOl,IROw)
      implicit none
      !------
      type (HECMWST_LOCAL_MESH), intent(in)::hecMESH
  