anchor.diff:  *. Obtain densities etc for chosen reference root
   *
         IF(ISROOT.NE.NROOT.AND.IDENSI.GE.1) THEN
  -          CALL REWINO(LUC)
  -          CALL FRMDSC(KPVEC1,NVAR,-1,LUC,IMZERO,IAMPACK)
  -          CALL COPVEC(KPVEC1,KPVEC2,NVAR)
  -          IF(IDENSI.EQ.1) THEN
  -            KRHO2 = 1
  -          END IF
  -          CALL DENSI2(IDENSI,WORK(KRHO1),WORK(KRHO2),
  -     &    KPVEC1,KPVEC2,0,0,EXPS2,ISPNDEN,WORK(KSRHO1),
  -     &    WORK(KRHO2AA),WORK(KRHO2AB),WORK(KRHO2BB),1)
  +          IF(ICISTR.EQ.1) THEN
  +            CALL REWINO(LUC)
  +            CALL FRMDSC(KPVEC1,NVAR,-1,LUC,IMZERO,IAMPACK)
  +            CALL COPVEC(KPVEC1,KPVEC2,NVAR)
  +            IF(IDENSI.EQ.1) THEN
  +              KRHO2 = 1
  +            END IF
  +            CALL DENSI2(IDENSI,WORK(KRHO1),WORK(KRHO2),
  +     &      KPVEC1,KPVEC2,0,0,EXPS2,ISPNDEN,WORK(KSRHO1),
  +     &      WORK(KRHO2AA),WORK(KRHO2AB),WORK(KRHO2BB),1)
  +          ELSE
  +            CALL REWINO(LUSC1)
  +            CALL REWINO(LUC)
  +            CALL COPVCD(LUC,LUSC1,WORK(KVEC1),0,LBLK)
  +            CALL DENSI2(IDENSI,WORK(KRHO1),WORK(KRHO2),
  +     &          WORK(KVEC1),WORK(KVEC2),LUC,LUSC1,EXPS2,
  +     &          ISPNDEN,WORK(KSRHO1),WORK(KRHO2AA),WORK(KRHO2AB),
  +     &          WORK(KRHO2BB),1 )
  +          END IF ! ICISTR switch
         END IF ! densities should be calculated
   *
   *. Supersymmetry
  
positive.code: *
         END IF
  * 
        IF(IPRNT.GT.1) WRITE(6,*)
        NDET = XISPSM(ISM,ISPC)
  *. Largest number of dets, any symmetry
  
negative.code:         INICI = -1
        END IF
  *
        IF(ICISTR.EQ.1) THEN
          LBLK = NVAR
        ELSE
          LBLK = - 1
        END IF
  *. This is a CI and not a perturbation calculation 
  *. ( what about B-W perturbation theory ???? ) 
        IPERTOP = 0
  