anchor.diff:      INTEGER, INTENT(OUT):: Iperm(*)
       INTEGER, INTENT(OUT):: Anc(*)
       INTEGER, INTENT(OUT):: Adjt(*)
  -
  +    !------
       INTEGER:: i
       INTEGER:: IDBg1
       INTEGER:: iy
       INTEGER:: izzz
       INTEGER:: joc
       INTEGER:: k
  -    INTEGER:: kk
       INTEGER:: l
       INTEGER:: ll
       INTEGER:: locc
  
positive.code:   !                20  building up matrix
    !                30  after LU decomposition
    !                40  after solving
    !
    !         #coded by t.arakawa of RIST on 040329
    !======================================================================!
    SUBROUTINE MATINI(Ir)
      IMPLICIT NONE
  
      INTEGER, INTENT(OUT):: Ir
  
      INTEGER:: IDBg
      INTEGER:: ir1
      INTEGER:: iv1
  
negative.code:     ELSE
         Ir = 0
      ENDIF
  
      ALLOCATE (TEMp(NDEg*NDEg*NEQns),STAT=IRR)
      IF ( IRR/=0 ) THEN
         WRITE (*,*) '##Error : Not enough memory'
  