anchor.diff:  *, In core version
          IF(NOCSF.EQ.1) THEN
           CALL GASDN2(I12,RHO1,RHO2,WORK(KVEC1P),WORK(KVEC2P),WORK(KC2),
  -     &       WORK(KCIOIO),WORK(KSIOIO),ISMOST(1,ICSM),
  +     &       WORK(KCIOIO),int_mb(KSIOIO),ISMOST(1,ICSM),
        &       ISMOST(1,ISSM),WORK(KCBLTP),WORK(KSBLTP),NACOB,
        &       WORK(KNSTSO(IATP)),WORK(KISTSO(IATP)),
        &       WORK(KNSTSO(IBTP)),WORK(KISTSO(IBTP)),
  
positive.code:      &WRITE(6,*) ' Density scratch space ',INTSCR
        CALL MEMMAN(KINSCR,INTSCR,'ADDL  ',2,'INSCR ') !done
  *
        CALL IAIBCM(ISSPC,WORK(KSIOIO))
        CALL IAIBCM(ISSPC,WORK(KCIOIO))
  *. Scratch space for CJKAIB resolution matrices
        CALL MXRESCPH(WORK(KCIOIO),IOCTPA,IOCTPB,NOCTPA,NOCTPB,
  
negative.code:       IF(ICISTR.GE.2) THEN
  *. Out of core version
          CALL GASDN2(I12,RHO1,RHO2,L,R,WORK(KC2),
       &       WORK(KCIOIO),int_mb(KSIOIO),ISMOST(1,ICSM),
       &       ISMOST(1,ISSM),WORK(KCBLTP),WORK(KSBLTP),NACOB,
       &       WORK(KNSTSO(IATP)),WORK(KISTSO(IATP)),
       &       WORK(KNSTSO(IBTP)),WORK(KISTSO(IBTP)),
  