anchor.diff:          IF(NTEST.GE.1000)
        &  WRITE(6,*) ' NCBLK, NSBLK, NXBLOCK = ', NCBLK, NSBLK, NXBLOCK
           CALL GASDIAT(VECSCR,LUSC,SHIFT,ICISTR,I12,
  -     &               WORK(KXBLTP),NXBLOCK,WORK(KXIBT),IUSE_EXP)
  +     &               int_mb(KXBLTP),NXBLOCK,WORK(KXIBT),IUSE_EXP)
   *
   * Average diagonal in each configuration
   *
  
positive.code: *
  * 3: Generate diagonal contribitions to the density
  *
        CALL DIAG_DET_TO_DENSI(CB,LUQ1,WORK(KCBLTP),
       &     NQBLOCK,WORK(KCIBT),RHO1_QQ, RHO2_QQ,I12)
  C     DIAG_DET_TO_DENSI(CB,LUC,IBLTP,NBLOCK,IBLKFO,
  C    &                             RHO1D,RHO2D,I12)
  
negative.code:       IF(ISFIRST.EQ.1) THEN
          CALL Z_BLKFO_FOR_CISPACE(ISSPC,ISSM,LBLOCK,ICOMP,
       &       NTESTL,NSBLOCK,NSBATCH,
       &       int_mb(KSIOIO),int_mb(KSBLTP),NSOCCLS_ACT,
       &       WORK(KSIOCCLS_ACT),
       &       WORK(KSLBT),WORK(KSLEBT),WORK(KSLBLK),WORK(KSI1BT),
       &       WORK(KSIBT),
       &       WORK(KSNOCCLS_BAT),WORK(KSIBOCCLS_BAT),ILTEST)
  