anchor.diff:    subroutine NUFCTX(Xlnzr,Colno,Dsln,Zln,Diag,Indx,Temp,Neqns,Parent,Nch,Nstop,Ndeg,Ndegl,Ir)
       implicit none
       !------
  -    integer, intent(IN):: Ndeg
  -    integer, intent(IN):: Ndegl
  -    integer, intent(IN):: Neqns
  -    integer, intent(IN):: Nstop
  -    integer, intent(IN):: Xlnzr(*)
  -    integer, intent(IN):: Colno(*)
  -    integer, intent(IN):: Parent(*)
  -    integer, intent(OUT):: Indx(*)
  -    integer, intent(OUT):: Nch(*)
  -    real(KIND=8), intent(OUT):: Zln(Ndeg*Ndeg,*)
  -    real(KIND=8), intent(OUT):: Diag(Ndegl,*)
  -    real(KIND=8), intent(OUT):: Temp(Ndeg*Ndeg,*)
  -    real(KIND=8), intent(OUT):: Dsln(Ndeg*Ndeg,*)
  +    integer, intent(in):: Ndeg
  +    integer, intent(in):: Ndegl
  +    integer, intent(in):: Neqns
  +    integer, intent(in):: Nstop
  +    integer, intent(in):: Xlnzr(*)
  +    integer, intent(in):: Colno(*)
  +    integer, intent(in):: Parent(*)
  +    integer, intent(out):: Indx(*)
  +    integer, intent(out):: Nch(*)
  +    real(kind=8), intent(out):: Zln(Ndeg*Ndeg,*)
  +    real(kind=8), intent(out):: Diag(Ndegl,*)
  +    real(kind=8), intent(out):: Temp(Ndeg*Ndeg,*)
  +    real(kind=8), intent(out):: Dsln(Ndeg*Ndeg,*)
       !------
       integer:: ic
       integer:: Ir
       integer:: l
  -    real(KIND=8):: t1
  -    real(KIND=8):: t2
  -    real(KIND=8):: t3
  -    real(KIND=8):: t4
  -    real(KIND=8):: t5
  -    real(KIND=8):: tt
  -    real(KIND=8):: zz(100)
  -    real(KIND=8):: t(100)
  +    real(kind=8):: t1
  +    real(kind=8):: t2
  +    real(kind=8):: t3
  +    real(kind=8):: t4
  +    real(kind=8):: t5
  +    real(kind=8):: tt
  +    real(kind=8):: zz(100)
  +    real(kind=8):: t(100)
   
       !
       ! phase I
  
positive.code:   subroutine SUM1(Ic,Xlnzr,Colno,Zln,Temp,Indx)
      implicit none
      !------
      integer, intent(IN):: Ic
      integer, intent(IN):: Xlnzr(*)
      integer, intent(IN):: Colno(*)
      integer, intent(OUT):: Indx(*)
      real(KIND=8), intent(OUT):: Temp(*)
      real(KIND=8), intent(OUT):: Zln(*)
      !------
      integer:: j
      integer:: jc
  
negative.code:   subroutine INV3(Dsln,Ir)
      implicit none
      !------
      integer, intent(out):: Ir
      real(kind=8), intent(out):: Dsln(6)
      !------
      real(kind=8):: t(2)
  
      Ir = 0
      if ( DABS(Dsln(1))<RMIn ) goto 100
  