anchor.diff:      INTEGER, INTENT(OUT):: Ir
       INTEGER, INTENT(OUT):: Indx(*)
       INTEGER, INTENT(OUT):: Nch(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Zln(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Diag(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Temp(*)
  -    DOUBLE PRECISION, INTENT(OUT):: Dsln(*)
  +    REAL(KIND=8), INTENT(OUT):: Zln(*)
  +    REAL(KIND=8), INTENT(OUT):: Diag(*)
  +    REAL(KIND=8), INTENT(OUT):: Temp(*)
  +    REAL(KIND=8), INTENT(OUT):: Dsln(*)
       !------
       INTEGER:: ic
  -    INTEGER:: ISEm
       INTEGER:: l
  -    INTEGER:: LRAtio
  -    DOUBLE PRECISION:: EPSm
  -    DOUBLE PRECISION:: RMAx
  -    DOUBLE PRECISION:: RMIn
  -    DOUBLE PRECISION:: t1
  -    DOUBLE PRECISION:: t2
  -    DOUBLE PRECISION:: t3
  -    DOUBLE PRECISION:: t4
  -    DOUBLE PRECISION:: t5
  -    DOUBLE PRECISION:: tt
  -    COMMON /MCHDPN/ RMAx, RMIn, EPSm, LRAtio
  -    COMMON ISEm
  +    REAL(KIND=8):: t1
  +    REAL(KIND=8):: t2
  +    REAL(KIND=8):: t3
  +    REAL(KIND=8):: t4
  +    REAL(KIND=8):: t5
  +    REAL(KIND=8):: tt
   
       ISEm = 1
       !
  
positive.code:     INTEGER:: ke
      INTEGER:: ks
      INTEGER:: l
      INTEGER:: ndeg2
      INTEGER:: ndeg2l
      DATA idbg, ndeg2, ndeg2l/0, 9, 6/
  
      Ir = 0
      ii = Invp(I)
  
negative.code:     INTEGER:: kk
      INTEGER:: ks
      INTEGER:: l
      REAL(KIND=8):: t(21)
      REAL(KIND=8):: zz(36)
  
      ks = Xlnzr(Ic)
      ke = Xlnzr(Ic+1)
  