anchor.diff:          ZERO = 0.0D0
           CALL SETVEC(dbl_mb(KRHO1),ZERO,NACOB**2)
           CALL REFORM_RHO1_TO_GNSM(dbl_mb(KRHO1),WORK(KLMAT2),2,
  -     &       WORK(KIREO_GNSYM_TO_TS_ACOB))
  +     &       int_mb(KIREO_GNSYM_TO_TS_ACOB))
           IF(NTEST.GE.0000) THEN
            WRITE(6,*) ' Supersymmetry averaged density '
            CALL  WRTMAT(dbl_mb(KRHO1),NACOB,NACOB,NACOB,NACOB)
  
positive.code: *. 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))
          IF(NTEST.GE.1000) THEN
            WRITE(6,*) ' Reformed density matrix '
            CALL APRBLM2(WORK(KLMAT2),NACOB_GS,NACOB_GS,NGENSMOB,0)
  
negative.code:         END DO
          CALL ICOPVE(ISMFSO,WORK(KMO_GNSYM),NTOOB)
          CALL REO_ACT_ORB_TO_GNSM(
       &     WORK(KMO_GNSYM),int_mb(KIREO_GNSYM_TO_TS_ACOB))
        END IF! Supersymmetry
  *
        IF(I_USE_SUPSYM.EQ.1) THEN
  