anchor.diff:          IF(NOCSF.EQ.0) THEN
            IF(ICNFBAT.EQ.1) THEN
   *. In core
  -           CALL CSDTVCM(WORK(KKVEC1),
  -     &          WORK(KCOMVEC1_SD),WORK(KCOMVEC2_SD),1,0,ICSM,ICSPC,2)
  +           CALL CSDTVCM(KKVEC1,
  +     &          KCOMVEC1_SD,KCOMVEC2_SD,1,0,ICSM,ICSPC,2)
   *. Final vector in now in KCOMVEC1_SD and we need sigma in the SD basis, so
   *. use 
              KKVEC1 = KCOMVEC1_SD
  
positive.code: C?        WRITE(6,*) ' KCLBT, WORK(KCLBT) = ', KCLBT
  C?        CALL IWRTMA(WORK(KCLBT),1,1,1,1)
  *
            CALL CSDTVCMN(WORK(KKVEC1),WORK(KKVEC2),WORK(KVEC3),
       &         1,0,ICSM,ICSPC,2,2,LUC_SD,LUCEFF,NCOCCLS,
       &         WORK(KCIOCCLS_ACT),WORK(KCIBT),WORK(KCLBT))
            LUC2 = LUC_SD
  C?        WRITE(6,*) ' After CSDTVCMN '
           END IF
  
negative.code: *. KCIBT have been constructed
          CALL CSDTVCMN(L,R,WORK(KVEC3),
       &       1,0,ICSM,ICSPC,2,2,LU_RDET,LUR,NCOCCLS_ACT,
       &       dbl_mb(KCIOCCLS_ACT),int_mb(KCIBT),WORK(KCLBT))
  *. 
          CALL REWINO(LUL)
          CALL REWINO(LU_LDET)
          CALL CSDTVCMN(L,R,WORK(KVEC3),
       &       1,0,ISSM,ISSPC,2,2,LU_LDET,LUL,NSOCCLS_ACT,
       &       dbl_mb(KSIOCCLS_ACT),WORK(KSIBT),WORK(KSLBT))
         END IF ! Incore 
        END IF ! CSFs are in use
  
  