anchor.diff:      real(kind=8), intent(in):: Zln(Ndeg,Ndeg,*)
       real(kind=8), intent(in):: Diag(Ndegl,*)
       real(kind=8), intent(in):: Dsln(Ndeg,Ndeg,*)
  -    real(kind=8), intent(out):: B(Ndeg,*)
  +    real(kind=8), intent(inout):: B(Ndeg,*)
       real(kind=8), intent(out):: Wk(Ndeg,*)
       !------
       integer:: i
  
positive.code:     integer, intent(in):: Colno(*)
      real(kind=8), intent(in):: Aij(9)
      integer, intent(out):: Ir
      real(kind=8), intent(out):: Zln(9,*)
      real(kind=8), intent(out):: Diag(6,*)
      real(kind=8), intent(out):: Dsln(9,*)
      !------
      integer:: i0
      integer:: ii
  
negative.code:     integer, intent(in):: Colno(*)
      integer, intent(in):: Parent(*)
      integer, intent(out):: Indx(*)
      integer, intent(inout):: Nch(*)
      real(kind=8), intent(inout):: Zln(Ndeg*Ndeg,*)
      real(kind=8), intent(inout):: Diag(Ndegl,*)
      real(kind=8), intent(out):: Temp(Ndeg*Ndeg,*)
      real(kind=8), intent(inout):: Dsln(Ndeg*Ndeg,*)
      !------
      integer:: ic
      integer:: Ir
  