anchor.diff:    !======================================================================!
     !> @brief NUFORM
     !======================================================================!
  -  subroutine NUFORM(hecMESH,hecMAT,Ir)
  +  subroutine NUFORM(hecMESH,hecMAT,NEQns,NDEg,LEN_colno,NSTop,LEN_dsln,INVp,XLNzr,COLno,DIAg,ZLN,DSLn,RALoc,STAge,Ir)
       use HECMW_UTIL
       implicit none
       !------
       type (HECMWST_LOCAL_MESH), intent(in)::hecMESH
       type (HECMWST_MATRIX), intent(in)::hecMAT
  +    integer(kind=kint), intent(in):: NEQns
  +    integer(kind=kint), intent(in):: NDEg
  +    integer(kind=kint), intent(in):: LEN_colno
  +    integer(kind=kint), intent(in):: NSTop
  +    integer(kind=kint), intent(in):: LEN_dsln
  +    integer(kind=kint), intent(in):: INVp(:)
  +    integer(kind=kint), intent(in):: XLNzr(:)
  +    integer(kind=kint), intent(in):: COLno(:)
  +    real(kind=8), pointer :: DIAg(:)
  +    real(kind=8), pointer :: ZLN(:)
  +    real(kind=8), pointer :: DSLn(:)
  +    integer(kind=kint), intent(inout):: RALoc
  +    integer(kind=kint), intent(inout):: STAge
       integer(kind=kint), intent(out):: Ir
       !------
       integer:: i
  
positive.code:         deallocate (QLInk)
          deallocate (ADJncy)
          deallocate (ZPIv)
          !*Nullify pointers
          nullify (IA)
          nullify (JA)
          nullify (DEG)
          nullify (MARker)
          nullify (RCHset)
          nullify (NBRhd)
          nullify (QSIze)
          nullify (QLInk)
          nullify (ADJncy)
          nullify (ZPIv)
          !
          !   build up xlnzr,colno  (this is the symbolic fct.)
          !
  
negative.code:     else
        Ir = 0
      endif
  
      allocate (wk(NDEg*NEQns),stat=IERror)
      if ( IERror/=0 ) then
  