anchor.diff:  *. Bring first density to form with supersymmetry blocks
   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))
  +     &       int_mb(KIREO_GNSYM_TO_TS_ACOB))
           IF(NTEST.GE.1000) THEN
             WRITE(6,*) ' Reformed density matrix '
             CALL APRBLM2(WORK(KLMAT2),NACOB_GS,NACOB_GS,NGENSMOB,0)
  
positive.code: *. 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 '
            CALL APRBLM2(WORK(KLMAT2),NACOB_GS,NACOB_GS,NGENSMOB,0)
  
negative.code: CBERT VECTOR DOUBLES!!! LOCAL
        CALL MEMMAN(KVEC1,LBLOCK,'ADDL  ',2,'VEC1  ')
        CALL MEMMAN(KVEC2,LBLOCK,'ADDL  ',2,'VEC2  ')
  C...dongxia  so KVEC1P and KVEC2P also local
        KVEC1P = KVEC1
        KVEC2P = KVEC2
  * 
  