anchor.diff:       & MINST,LUDIA,LUC,LUHC,LUSC1,LUSC2,LUSC3,LUSC34,LUSC35,LU8,
        & NVAR,NBLK,NROOT,MXCIV,MAXIT,LUCIVI,IPRNT,WORK(KLH0_EIGVEC),
        & NPRDET,WORK(KH0),WORK(KLH0_SUBDT),
  -     & NP1,NP2,NQ,WORK(KH0SCR),EADD,ICISTR,LBLK,
  +     & NP1,NP2,NQ,dbl_mb(KH0SCR),EADD,ICISTR,LBLK,
        & IDIAG,dbl_mb(KVEC3),THRES_E,NBATCH,
        & int_mb(KCLBT),int_mb(KCLEBT),int_mb(KCLBLK),int_mb(KCI1BT),
        & int_mb(KCIBT),
  
positive.code:       END IF ! ISBSPC_SEL .ne. 0
  *. Space for LZ and L2  for all roots of a CI
          MXLROOT = MAX(NROOT,INI_NROOT)
          CALL MEMMAN(KLZEXP,MXLROOT,'ADDL  ',2,'LZEXP ')
          CALL MEMMAN(KL2EXP,MXLROOT,'ADDL  ',2,'L2EXP ')
  *
  *. indeces for pair of orbitals symmetry ordered
  *. Lower half packed
  
negative.code:               KRHO2 = 1
              END IF
              CALL DENSI2(IDENSI,dbl_mb(KRHO1),dbl_mb(KRHO2),
       &      KPVEC1,KPVEC2,0,0,EXPS2,ISPNDEN,dbl_mb(KSRHO1),
       &      dbl_mb(KRHO2AA),dbl_mb(KRHO2AB),dbl_mb(KRHO2BB),1)
            ELSE
              CALL REWINO(LUSC1)
              CALL REWINO(LUC)
              CALL COPVCD(LUC,LUSC1,dbl_mb(KVEC1),0,LBLK)
              CALL DENSI2(IDENSI,dbl_mb(KRHO1),dbl_mb(KRHO2),
       &          dbl_mb(KVEC1),dbl_mb(KVEC2),LUC,LUSC1,EXPS2,
       &          ISPNDEN,dbl_mb(KSRHO1),dbl_mb(KRHO2AA),dbl_mb(KRHO2AB),
       &          dbl_mb(KRHO2BB),1 )
            END IF ! ICISTR switch
        END IF ! densities should be calculated
  *
  