anchor.diff:            CALL MEMCHK2('BE_NWM')
   C              GET_EXP_MKAPPA(EXPMK,KAPPAP,IOOEXC,NOOEXC)
             CALL GET_EXP_MKAPPA(WORK(KLMO1),WORK(KLKAPPA),
  -     &                        WORK(KLOOEXCC),NOOEXC)
  +     &                        int_mb(KLOOEXCC),NOOEXC)
             CALL MEMCHK2('AF_EMK')
   *. And new MO-coefficients
             CALL MULT_BLOC_MAT(WORK(KLMO2),WORK(KMOREF),WORK(KLMO1),
  
positive.code: *
        CALL FOCK_MAT_STANDARD(WORK(KF),2,WORK(KINT1),WORK(KFA))
        IF(IPRNT.GE.100) WRITE(6,*) ' F constructed '
        CALL E1_FROM_F(WORK(KLE1),WORK(KF),1,WORK(KLOOEXC),
       &               WORK(KLOOEXCC),
       &               NOOEXC,NTOOB,NTOOBS,NSMOB,IBSO,IREOST)
        E1NRM_ORB = SQRT(INPROD(WORK(KLE1),WORK(KLE1),NOOEXC))
        VNFINAL = E1NRM_ORB + ERROR_NORM_FINAL_CI
  
negative.code:       ISSPC = IPSPC
        ICSPC = IQSPC
        WRITE(6,*) ' DENSI2 will be called for PQ density '
        CALL DENSI2(I12,dbl_mb(KLRHO1),dbl_mb(KLRHO2),L,R,LUSC1,LUP,EXPS2,
       &     0,XDUM,XDUM,XDUM,XDUM,0)
        ONE = 1.0D0
        CALL VECSUM(RHO1,RHO1,dbl_mb(KLRHO1),ONE,ONE,LRHO1)
        IF(I12.EQ.2)
       &CALL VECSUM(RHO2,RHO2,dbl_mb(KLRHO2),ONE,ONE,LRHO2)
  *
  * Obtain Q-P densities - from  P-Q densities
  *
        CALL TRPMAT(dbl_mb(KLRHO1),NACOB,NACOB,dbl_mb(KLRHO1B))
        CALL VECSUM(RHO1,RHO1,dbl_mb(KLRHO1B),ONE,ONE,LRHO1)
        IF(I12.EQ.2) THEN
          CALL CONJ_RHO2(NACOB,dbl_mb(KLRHO2),dbl_mb(KLRHO2B))
          CALL VECSUM(RHO2,RHO2,dbl_mb(KLRHO2B),ONE,ONE,LRHO2)
        END IF
  *
  * And then to the interesting: the Q-Q density
  *
  C     GET_APPROX_DEN_AKBK_QQ(RHO1_QQ, RHO2_QQ,LUP,E,CB,HCB,I12) 
        CALL GET_APPROX_DEN_AKBK_QQ(dbl_mb(KLRHO1),dbl_mb(KLRHO2),LUP,E,
       &    L,R,I12)
  *
        IF(NTEST.GE.1000) THEN
  