anchor.diff:         END DO
   *. Mapping of active orbitals from general symmetry to type order
          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))
   *. Reordering array from actual to supersymmetry-blocked order.
   C       GET_IACT_TO_GENSM_REO(IACT_TO_GENSM_REO,
   C    &  ISTA_TO_GENSM_REO, MO_STA_TO_ACT_REO, NTOOB)
  
positive.code:         END DO
          CALL ICOPVE(ISMFSO,WORK(KMO_GNSYM),NTOOB)
          CALL REO_ACT_ORB_TO_GNSM(
       &     WORK(KMO_GNSYM),WORK(KIREO_GNSYM_TO_TS_ACOB))
        END IF! Supersymmetry
  *
        IF(I_USE_SUPSYM.EQ.1) THEN
  
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)
  