anchor.diff:       &       NBATCHL,WORK(KSLBT),WORK(KSLEBT),WORK(KSI1BT),
        &       WORK(KSIBT),
        &       NBATCHR,WORK(KCLBT),WORK(KCLEBT),WORK(KCI1BT),
  -     &       WORK(KCIBT),int_mb(KCONSPA),int_mb(KCONSPB),
  +     &       int_mb(KCIBT),int_mb(KCONSPA),int_mb(KCONSPB),
        &       dbl_mb(KLSCLFCL),dbl_mb(KLSCLFCR),S2_TERM1,IUSE_PH,IPHGAS,
        &       IDOSRHO1,SRHO1,IDOSRHO2,RHO2AA,RHO2AB,RHO2BB,
        &       NDACTORB,IDACTSPC,IDTFREORD,IFTDREORD,IOBPTS_SEL,
  
positive.code:         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)
  *. Final vector in now in KCOMVEC1_SD and we need sigma in the SD basis, so
  *. use 
             KKVEC1 = KCOMVEC1_SD
  
negative.code:      &       NBATCHL,WORK(KSLBT),WORK(KSLEBT),WORK(KSI1BT),
       &       WORK(KSIBT),
       &       NBATCHR,WORK(KCLBT),WORK(KCLEBT),WORK(KCI1BT),
       &       int_mb(KCIBT),int_mb(KCONSPA),int_mb(KCONSPB),
       &       dbl_mb(KLSCLFCL),dbl_mb(KLSCLFCR),S2_TERM1,IUSE_PH,IPHGAS,
       &       IDOSRHO1,SRHO1,IDOSRHO2,RHO2AA,RHO2AB,RHO2BB,
       &       NDACTORB,IDACTSPC,IDTFREORD,IFTDREORD,IOBPTS_SEL,NINOB,
  