anchor.diff:  *. Size of C(Ka,Jb,j),C(Ka,KB,ij)  resolution matrices
         IOCTPA = IBSPGPFTP(IATP)
         IOCTPB = IBSPGPFTP(IBTP)
  -      CALL MXRESCPH(WORK(KCIOIO),IOCTPA,IOCTPB,NOCTPA,NOCTPB,
  +      CALL MXRESCPH(dbl_mb(KCIOIO),IOCTPA,IOCTPB,NOCTPA,NOCTPB,
        &            NSMST,NSTFSMSPGP,MXPNSMST,
        &            NSMOB,MXPNGAS,NGAS,NOBPTS,IPRCIX,MAXK,
        &            NELFSPGP,
  
positive.code:            IF(NROOT.GT.1) THEN
              CALL REWINO(LUSC1)
              IF(NOCSF.EQ.1) THEN
                CALL GASANA(WORK(KVEC1),NBLOCK,WORK(KCIBT),WORK(KCBLTP),
       &                    LUSC1,ICISTR)
              ELSE
                IOUT = 6
                THRES = 0.1D0
                MAXTRM = 100
                CALL ANACSF2(LUC,NCOCCLS_ACT,WORK(KCIOCCLS_ACT),ISM,
       &             WORK(KPVEC1),WORK(KICONF_OCC(ISM)),
       &             NCONF_PER_OPEN(1,ISM),WORK(KCFTP),THRES,
       &             MAXTRM,IOUT)
              END IF
             ELSE 
              CALL REWINO(LUC)
              IF(NOCSF.EQ.1) THEN
                CALL GASANA(WORK(KVEC1),NBLOCK,WORK(KCIBT),WORK(KCBLTP),
       &                    LUC,ICISTR)
              ELSE
                IOUT = 6
                THRES = 0.1D0
                MAXTRM = 100
                CALL ANACSF2(LUC,NCOCCLS_ACT,WORK(KCIOCCLS_ACT),ISM,
       &             WORK(KPVEC1),WORK(KICONF_OCC(ISM)),
       &             NCONF_PER_OPEN(1,ISM),WORK(KCFTP),THRES,
       &             MAXTRM,IOUT)
              END IF !NOCSF switch
             END IF !Nroot switch
  
negative.code: *. And for the Sigma vector
        CALL Z_BLKFO_FOR_CISPACE(ISPC,ISM,LBLOCK,ICOMP,
       &     IPRNT,NSBLOCK,NSBATCH,
       &     dbl_mb(KSIOIO),dbl_mb(KSBLTP),NSOCCLS_ACT,
       &     dbl_mb(KSIOCCLS_ACT),
       &     int_mb(KSLBT),int_mb(KSLEBT),int_mb(KSLBLK),int_mb(KSI1BT),
       &     int_mb(KSIBT),
       &     int_mb(KSNOCCLS_BAT),int_mb(KSIBOCCLS_BAT),ILTEST)
  *. Number of BLOCKS
        NBLOCK = NCBLOCK
        NBATCH = NCBATCH
  