anchor.diff:            ELSE
               CALL REWINO(LUSC1)
               CALL REWINO(LUC)
  -            CALL COPVCD(LUC,LUSC1,WORK(KVEC1),0,LBLK)
  +            CALL COPVCD(LUC,LUSC1,dbl_mb(KVEC1),0,LBLK)
               CALL DENSI2(IDENSI,WORK(KRHO1),WORK(KRHO2),
  -     &          WORK(KVEC1),WORK(KVEC2),LUC,LUSC1,EXPS2,
  +     &          dbl_mb(KVEC1),dbl_mb(KVEC2),LUC,LUSC1,EXPS2,
        &          ISPNDEN,WORK(KSRHO1),WORK(KRHO2AA),WORK(KRHO2AB),
        &          WORK(KRHO2BB),1 )
             END IF ! ICISTR switch
  
positive.code:              WRITE(6,*) ' ======================'
               WRITE(6,*)
               IF(NROOT.EQ.1) THEN
                 CALL WRTVCD(WORK(KVEC1),LUC  ,1,LBLK)
               ELSE
                 CALL WRTVCD(WORK(KVEC1),LUSC1,1,LBLK)
               END IF
  *         ^ End if print of CI vector
          END IF
  
negative.code:          ECOREL = 0.0D0
           IF(NOCSF.EQ.1) THEN
  *. Determinants: TTSS blocks
             CALL COMHAM(WORK(KLHMAT),NVAR,NBLOCK,int_mb(KCLBLK),
       &                 KPVEC1,KPVEC2,ECOREL)
           ELSE
  *. CSFs: occlass blocks
  *. Number of CSFs per occlass
             CALL MEMMAN(KLNCS_FOR_OCCLS,NCOCCLS_ACT,'ADDL  ',1,'NCS_OC')
             CALL GET_NCSF_PER_OCCLS_FOR_CISPACE(ISM,dbl_mb(KCIOCCLS_ACT),
       &          NCOCCLS_ACT,int_mb(KNCS_FOR_OCCLS),
       &          int_mb(KLNCS_FOR_OCCLS))
  C          GET_NCSF_PER_OCCLS_FOR_CISPACE(ISYM,IOCCLS_ACT,
  C    &           NOCCLS_ACT,NCS_FOR_OCCLS,NCS_FOR_OCCLS_ACT)
             CALL COMHAM(WORK(KLHMAT),NVAR,NCOCCLS_ACT,
       &          int_mb(KLNCS_FOR_OCCLS),KPVEC1,KPVEC2,ECOREL)
           END IF
           STOP ' Enforced stop after COMHAM'
        END IF
  