anchor.diff:    subroutine S2PDOT(Bi,B,Zln,Colno,Ks,Ke)
       implicit none
       !------
  -    integer, intent(IN):: Ke
  -    integer, intent(IN):: Ks
  -    integer, intent(IN):: Colno(*)
  -    real(KIND=8), intent(IN):: Zln(4,*)
  -    real(KIND=8), intent(IN):: B(2,*)
  -    real(KIND=8), intent(OUT):: Bi(2)
  +    integer, intent(in):: Ke
  +    integer, intent(in):: Ks
  +    integer, intent(in):: Colno(*)
  +    real(kind=8), intent(in):: Zln(4,*)
  +    real(kind=8), intent(in):: B(2,*)
  +    real(kind=8), intent(out):: Bi(2)
       !------
       integer:: j
       integer:: jj
  
positive.code:   subroutine ADDR3(I,J,Aij,Invp,Xlnzr,Colno,Diag,Zln,Dsln,Nstop,Ir)
      implicit none
      !------
      integer, intent(IN):: I
      integer, intent(IN):: J
      integer, intent(IN):: Nstop
      integer, intent(IN):: Invp(*)
      integer, intent(IN):: Xlnzr(*)
      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:   subroutine QQSORT(Iw,Ik)
      implicit none
      !------
      integer, intent(in):: Ik
      integer, intent(out):: Iw(*)
      !------
      integer:: itemp
      integer:: l
  