anchor.diff:          if ( IERror/=0 ) stop "ALLOCATION ERROR, parent: SUB. matini.f"
           allocate (FCT%NCH(FCT%NEQns+1),stat=IERror)
           if ( IERror/=0 ) stop "ALLOCATION ERROR, nch: SUB. matini.f"
  -        call POSORD(FCT%PARent,BTRee,FCT%INVp,FCT%IPErm,FCT%NEQns,Quarent)
  +        call POSORD(FCT%PARent,BTRee,FCT%INVp,FCT%IPErm,FCT%NCH,FCT%NEQns,Quarent)
           !
           !   generate skelton graph
           !
  
positive.code:     integer(kind=kint), intent(out):: Parent(:)
      integer(kind=kint), intent(inout):: Invp(:)
      integer(kind=kint), intent(out):: Iperm(:)
      !------
      integer(kind=kint), allocatable:: Pordr(:)
      integer(kind=kint), allocatable:: Nch(:)
      integer(kind=kint), allocatable:: Iw(:)
      integer(kind=kint), allocatable:: Mch(:)
      integer(kind=kint):: i
  
negative.code:     enddo
  
      deallocate (Pordr)
      deallocate (Iw)
      deallocate (Mch)
    end subroutine POSORD
  