anchor.diff:          DO IGAS = 1, IGASL-1
             NELI = NELFGP(ITPFGS(IGAS))
   C         ICOPVE2(IIN,IOFF,NDIM,IOUT)
  -          CALL ICOPVE2(WORK(KOCSTR(ITPFGS(IGAS))),
  +          CALL ICOPVE2(int_mb(KOCSTR(ITPFGS(IGAS))),
        &                 (ISTRNM(IGAS)+IBSTFGS(IGAS)-2)*NELI+1,
        &                 NELI,ISTROC(IBEL,ISTR))
   C         DO JEL = 1, NELI
  
positive.code:      &             WRITE(6,*) '  (alpha spin)'
                IF (NSPIN.EQ.2.AND.IJSPIN.EQ.2)
       &             WRITE(6,*) '  (beta spin)'
                CALL WRTMAT(WORK(KLFBLK),NI,NJ,NI,NJ)
              END IF
  *. Block has been constructed , transfer to -complete- 
  *. symmetry blocked Fock matrix
  
negative.code:                 CALL WRTMAT(dbl_mb(KLINT),NI,NK,NI,NK)
                  WRITE(6,*) 
       &          ' 1-e densi for ISM JGAS KGAS ',ISM,JGAS,KGAS
                  CALL WRTMAT(dbl_mb(KLDEN),NJ,NK,NJ,NK)
                END IF
  *. And then a matrix multiply( they are pretty much in fashion 
  *. these days )
                CALL MATML7(dbl_mb(KLFBLK),dbl_mb(KLINT),dbl_mb(KLDEN),
       &                    NI,NJ,NI,NK,NJ,NK,ONE,ONE,2)
                 IF(NTEST.GE.1000) THEN
                   WRITE(6,*) ' Updated block '
                   CALL WRTMAT(dbl_mb(KLFBLK),NI,NJ,NI,NJ)
                 END IF
   
              END DO
              IF(NTEST.GE.1000) THEN
                WRITE(6,*) ' One-electron contributions'
                WRITE(6,*) ' =========================='
                CALL WRTMAT(dbl_mb(KLFBLK),NI,NJ,NI,NJ)
              END IF
              IF(I12L.EQ.2) THEN
  *. 2 : Two-electron part
  