anchor.diff:          END DO
           CALL ICOPVE(ISMFSO,WORK(KMO_GNSYM),NTOOB)
           CALL REO_ACT_ORB_TO_GNSM(
  -     &     WORK(KMO_GNSYM),WORK(KIREO_GNSYM_TO_TS_ACOB))
  +     &     WORK(KMO_GNSYM),int_mb(KIREO_GNSYM_TO_TS_ACOB))
         END IF! Supersymmetry
   *
         IF(I_USE_SUPSYM.EQ.1) THEN
  
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: *. 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,
       &       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)
  