anchor.diff:      !
       !  build ia,ja
       !
  -    lia = last - neqns1
  -    lja = lia - NTTbr*2
  -    last = lja
  -    !rmem
       allocate (IA(NEQns+1),stat=IERror)
       if ( IERror/=0 ) stop "ALLOCATION ERROR, ia: SUB. matini"
       !WINDEBUG
  
positive.code:     integer(kind=4), allocatable :: leaf(:)
      integer:: ir1
      integer:: irr
      integer:: iv1
      integer:: izz
      integer:: izz0
      integer:: ladp
      integer:: last
      integer:: lbtree
      integer:: lcolno
      integer:: lcpt
      integer:: left
      integer:: lenv
      integer:: lenv2
      integer:: lia
      integer:: lja
      integer:: lleaf
      integer:: lncol
      integer:: lnleaf
      integer:: lpordr
      integer:: lwk1
      integer:: lwk2
      integer:: lwk3
      integer:: lwk4
      integer:: lwk5
      integer:: lwk6
      integer:: lwk7
      integer:: lwk8
      integer:: lxleaf
      integer:: maxl
      integer:: neqns1
      integer:: neqnsz
      integer:: ierror
  
      izz0 = 0
  
      Ir = 0
      lenv = 1000000000
      lenv2 = LRAtio*lenv
      neqns1 = NEQns + 2
      LEN_dsln = lenv2
      iv1 = 51
  
      !*Initialize allocation measure variables
      IALoc = 0
  
negative.code:         !
          !   post ordering
          !
          allocate (PARent(NEQns),stat=IERror)
          if ( IERror/=0 ) stop "ALLOCATION ERROR, parent: SUB. matini.f"
          allocate (NCH(NEQns+1),stat=IERror)
  