anchor.diff:        CALL T_TO_NK_VECS
        &  (T,KKORB,C,LUCIN,LUCOUT,int_mb(KNSTSO(IATP)),
        &      int_mb(KNSTSO(IBTP)),
  -     &      NBLOCK,WORK(KLCIBT),
  -     &      NAEL,NBEL,WORK(KLASTR),WORK(KLBSTR),
  -     &      WORK(KLCBLTP),NSMST,
  -     &      ICISTR,NTOOB,WORK(KLKAOC),WORK(KLKBOC))
  +     &      NBLOCK,int_mb(KLCIBT),
  +     &      NAEL,NBEL,int_mb(KLASTR),int_mb(KLBSTR),
  +     &      dbl_mb(KLCBLTP),NSMST,
  +     &      ICISTR,NTOOB,int_mb(KLKAOC),int_mb(KLKBOC))
      
         CALL MEMMAN(IDUM,IDUM,'FLUSM',IDUM,'T_TO_N')
   *
  
positive.code:       IF(ICISTR.GE.2) THEN
  *. Out of core version
          CALL GASDN2(I12,RHO1,RHO2,L,R,dbl_mb(KC2),
       &       WORK(KCIOIO),int_mb(KSIOIO),ISMOST(1,ICSM),
       &       ISMOST(1,ISSM),int_mb(KCBLTP),int_mb(KSBLTP),NACOB,
       &       int_mb(KNSTSO(IATP)),int_mb(KISTSO(IATP)),
       &       int_mb(KNSTSO(IBTP)),int_mb(KISTSO(IBTP)),
  
negative.code:       CALL TS_SYM_PNT2(IGRP,NIGRP,MXVLI,MNVLI,ISPGPSM,
       &                 IOFFI,LOFFI)
  *. Offset and dimension for active group in I strings
        CALL ICOPVE2(int_mb(KISTSGP(1)),(IGRP(IACGRP)-1)*NSMST+1,NSMST,
       &               IACIST)
  C?    WRITE(6,*) ' IACIST for IACGRP,IGRP = ', IACGRP,IGRP(IACGRP)
  C?    CALL IWRTMA(IACIST,1,NSMST,1,NSMST)
  