anchor.diff:      INTEGER, INTENT(IN):: Par(*)
       INTEGER, INTENT(OUT):: Indx(*)
       INTEGER, INTENT(OUT):: Nch(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Diag(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Temp(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Zln(*)
  -    !------
  -    INTEGER:: ISEm
  +    REAL(KIND=8), INTENT(OUT):: Diag(*)
  +    REAL(KIND=8), INTENT(OUT):: Temp(*)
  +    REAL(KIND=8), INTENT(OUT):: Zln(*)
  +    !------
       INTEGER:: j
       INTEGER:: jc
       INTEGER:: jj
  
positive.code:     INTEGER, INTENT(IN):: Ik
      INTEGER, INTENT(OUT):: Iw(*)
      !------
      INTEGER:: IDBg1
      INTEGER:: itemp
      INTEGER:: l
      INTEGER:: m
      COMMON /DEBUG / IDBg1
  
      IF ( Ik<=1 ) RETURN
      DO l = 1, Ik - 1
  
negative.code:   !                              =-1   non positive index
    !                              =1    too big index
    !                              =10   insufficient storage
    !        contents of iv
    !               pointers 1 &zpiv(1)  2 &iperm(1)  3 &invp(1)
    !                        4 &parent(1)5 &nch(1)    6 &xlnzr(1)
  