anchor.diff:              IPRDEN_SAVE = IPRDEN
               IPRDEN = 0
               XDUM = 0
  -            CALL DENSI2(1     ,dbl_mb(KRHO1),WORK(KRHO2),
  +            CALL DENSI2(1     ,dbl_mb(KRHO1),dbl_mb(KRHO2),
        &           dbl_mb(KVEC1),dbl_mb(KVEC2),LUSC1,LUHC,EXPS2,
        &           0,XDUM,XDUM,XDUM,XDUM,1)
               IPRDEN = IPRDEN_SAVE
  
positive.code:          write(fn,'(E22.15)')enrep
           end if
           close(fn)
  c.. Now that the integrals are in g_lu1e and g_lu2e, fully packed. 
  
         ELSE IF (INTIMP .EQ. 5) THEN
  *
  
negative.code:             KRHO2 = 1
            END IF
            IF(IDENSI.GE.1)
       &    CALL DENSI2(IDENSI,dbl_mb(KRHO1),dbl_mb(KRHO2),
       &         KPVEC1,KPVEC2,0,0,EXPS2,ISPNDEN,WORK(KSRHO1),
       &         WORK(KRHO2AA),WORK(KRHO2AB),WORK(KRHO2BB),1)
            IF(IDENSI.EQ.2) THEN
  