anchor.diff:                  DO M = 1, NM
                     IIOFF = KLINT + (M-1)*NKL*NI
                     IDOFF = KLDEN + (M-1)*NKL*NJ
  -                  CALL MATML7(dbl_mb(KLFBLK),WORK(IIOFF),WORK(IDOFF),
  -     &                        NI,NJ,NKL,NI,NKL,NJ,ONE,ONE,1)
  +                  CALL MATML7(dbl_mb(KLFBLK),dbl_mb(IIOFF),
  +     &                   dbl_mb(IDOFF),NI,NJ,NKL,NI,NKL,NJ,ONE,ONE,1)
                   END DO
                 END DO
                 END DO
  
positive.code:                 IXCHNG = 0
                  ICOUL  = 1
                  ONE = 1.0D0
                  CALL GETINT(WORK(KLINT),
       &               KGAS,KSM,LGAS,LSM,IGAS,ISM,MGAS,MSM,
       &               IXCHNG,0,0,ICOUL,ONE,ONE)
                  
  
negative.code: *
  * type of each symmetry block ( full, lower diagonal, absent )
  *
        CALL ZBLTP(ISMOST(1,ISM),NSMST,IDC,int_mb(KCBLTP),IDUMMY)
        CALL ZBLTP(ISMOST(1,ISM),NSMST,IDC,int_mb(KSBLTP),IDUMMY)
  
  *
  *. Number of dets etc per TTS block
  