anchor.diff:        NBEL = NELEC(IBTP)
   *. Arrays giving allowed type combinations 
         CALL MEMMAN(KSIOIO,NOCTPA*NOCTPB,'ADDL  ',2,'SIOIO ')
  -      CALL IAIBCM(ISSPC,WORK(KSIOIO))
  +      CALL IAIBCM(ISSPC,dbl_mb(KSIOIO))
   *. Arrays for additional symmetry operation
         IF(IDC.EQ.3.OR.IDC.EQ.4) THEN
           CALL MEMMAN(KSVST,NSMST,'ADDL  ',2,'SVST  ')
  
positive.code:         ELSE ! ICISTR switch
            IF(IDENSI.GE.1.OR.NROOT.GT.1) THEN
               CALL REWINO(LUSC1)
               CALL COPVCD(LUC,LUSC1,WORK(KVEC1),0,LBLK)
               CALL COPVCD(LUSC1,LUHC,WORK(KVEC1),1,LBLK)
            END IF
            IF(IDENSI.GE.1)
       &    CALL DENSI2(IDENSI,WORK(KRHO1),WORK(KRHO2),
       &          WORK(KVEC1),WORK(KVEC2),LUSC1,LUHC,EXPS2,
       &          ISPNDEN,WORK(KSRHO1),WORK(KRHO2AA),WORK(KRHO2AB),
       &          WORK(KRHO2BB),1 )
            IF(IPRNT.GT.1.AND.IDENSI.EQ.2)
  
negative.code: *. Largest number of strings of given symmetry and type
        MAXA = 0
        IF(NAEL.GE.1) THEN
          MAXA1 = IMNMX(int_mb(KNSTSO(IATPM1)),NSMST*NOCTYP(IATPM1),2)
          MAXA = MAX(MAXA,MAXA1)
        END IF
        IF(NAEL.GE.2) THEN
          MAXA1 = IMNMX(int_mb(KNSTSO(IATPM2)),NSMST*NOCTYP(IATPM2),2)
          MAXA = MAX(MAXA,MAXA1)
        END IF
        MAXB = 0
        IF(NBEL.GE.1) THEN
          MAXB1 = IMNMX(int_mb(KNSTSO(IBTPM1)),NSMST*NOCTYP(IBTPM1),2)
          MAXB = MAX(MAXB,MAXB1)
        END IF
        IF(NBEL.GE.2) THEN
          MAXB1 = IMNMX(int_mb(KNSTSO(IBTPM2)),NSMST*NOCTYP(IBTPM2),2)
          MAXB = MAX(MAXB,MAXB1)
        END IF
        MXSTBL = MAX(MAXA,MAXB)
  