anchor.diff:        CALL APRBLM2(WORK(KFI),NTOOBS,NTOOBS,NMSOB,1)
         END IF
         CALL FA_FROM_INIINT
  -     &(WORK(KFA),WORK(KMOMO),WORK(KMOMO),WORK(KRHO1),1)
  +     &(WORK(KFA),WORK(KMOMO),WORK(KMOMO),dbl_mb(KRHO1),1)
         KINT2_A(IE2ARR_F)  = KINT2_FSAVE
   *
         ONE = 1.0D0
  
positive.code:         CALL APRBLM2(WORK(KFI),NTOOBS,NTOOBS,NSMOB,1)
        END IF
        CALL FA_FROM_INIINT
       &(WORK(KFA),WORK(KMOMO),WORK(KMOMO),WORK(KRHO1),1)
        KINT2_A(IE2ARR_F)  = KINT2_FSAVE
  *
        ONE = 1.0D0
  
negative.code: *. Obtain and diagonalize symmetry ordered density matrix 
  *. over active orbitals
  C            REFORM_RHO1_TO_GNSM(RHO1_ST,RHO1_GNSM_ST,IWAY,IREO_GNSYM_TO_TS)
          CALL REFORM_RHO1_TO_GNSM(dbl_mb(KRHO1),WORK(KLMAT2),1,
       &       WORK(KIREO_GNSYM_TO_TS_ACOB))
          IF(NTEST.GE.1000) THEN
            WRITE(6,*) ' Reformed density matrix '
  