anchor.diff:            NKL = NK*NL
   *. Obtain block B(KSM,KGAS,LSM,LGAS)
             CALL EXTR_SYMGAS_BLK_FROM_ORBMAT
  -     &    (B,WORK(KLBBLK),KSM,KGAS,LSM,LGAS)
  +     &    (B,dbl_mb(KLBBLK),KSM,KGAS,LSM,LGAS)
   *. Fetch RHO2(ISM,JSM,KSM,LSM)
   C         GETD2(RHO2B,ISM,IGAS,JSM,JGAS,KSM,KGAS,LSM,LGAS,ISPC)
             CALL GETD2
  -     &    (WORK(KLRHO2B),ISM,IGAS,JSM,JGAS,KSM,KGAS,LSM,LGAS,1)
  +     &    (dbl_mb(KLRHO2B),ISM,IGAS,JSM,JGAS,KSM,KGAS,LSM,LGAS,1)
   * sum(kl) RHO2(ij,kl) B(kl)
             FACTORC = 0.0D0
             FACTORAB = 1.0D0
             ZERO = 0.0D0
  -          CALL SETVEC(WORK(KLVEC),ZERO,NIJ)
  -          CALL MATML7(WORK(KLVEC),WORK(KLRHO2B),WORK(KLBBLK),NIJ,1,
  -     &                NIJ,NKL,NKL,1,FACTORC,FACTORAB,0)
  +          CALL SETVEC(dbl_mb(KLVEC),ZERO,NIJ)
  +          CALL MATML7(dbl_mb(KLVEC),dbl_mb(KLRHO2B),dbl_mb(KLBBLK),NIJ,
  +     &                1,NIJ,NKL,NKL,1,FACTORC,FACTORAB,0)
             IF(NTEST.GE.1000) THEN
               WRITE(6,*) ' KLVEC: '
  -            CALL WRTMAT(WORK(KLVEC),1,NIJ,1,NIJ)
  +            CALL WRTMAT(dbl_mb(KLVEC),1,NIJ,1,NIJ)
             END IF
   *. sum(ij) (A(ij) (sum(kl) RHO2(ij,kl)B(kl) )
  -          AB2 = AB2 + INPROD(WORK(KLABLK),WORK(KLVEC),NIJ)
  +          AB2 = AB2 + INPROD(dbl_mb(KLABLK),dbl_mb(KLVEC),NIJ)
             IF(NTEST.GE.100) THEN
               WRITE(6,'(A,4I4,E20.12)') 
        &      ' Two-electron contribution after ISM, JSM, KSM, LSM = ',
  
positive.code: *. Obtain block A(ISM,IGAS,JSM,JGAS)
  C       EXTR_SYMGAS_BLK_FROM_ORBMAT(A,ABLK,ISM,IGAS,JSM,JGAS)
          CALL EXTR_SYMGAS_BLK_FROM_ORBMAT
       &  (A,WORK(KLABLK),ISM,IGAS,JSM,JGAS)
          DO KSM = 1, NSMOB
           LSM = MULTD2H(KSM,IBSM)
           DO KGAS = 1, NGAS
  
negative.code:       END DO
  *
        NTTS = MXNTTS
        CALL MEMMAN(KLSCLFCL,NTTS, 'ADDL  ',2,'SCLF_L') !done
        CALL MEMMAN(KLSCLFCR,NTTS, 'ADDL  ',2,'SCLF_R') !done
  
        NVARL = NSVAR
        NVARR = NCVAR
  