anchor.diff:  *. 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)),
  +     &       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)),
        &       NAEL,IATP,NBEL,IBTP,IOCTPA,IOCTPB,NOCTPA,NOCTPB,
        &       NSMST,NSMOB,NSMSX,NSMDX,MXPNGAS,NOBPTS,IOBPTS,      
        &       MAXK,MAXI,LSCR1,LSCR1,WORK(KCSCR),WORK(KSSCR),
  
positive.code: *. Orbital K in type ordering
        KKORB = IREOST(KORB)
        CALL T_TO_NK_VECS
       &  (T,KKORB,C,LUCIN,LUCOUT,WORK(KNSTSO(IATP)),WORK(KNSTSO(IBTP)),
       &      NBLOCK,WORK(KLCIBT),
       &      NAEL,NBEL,WORK(KLASTR),WORK(KLBSTR),
       &      WORK(KLCBLTP),NSMST,
  
negative.code:        IF(NOCSF.EQ.1) THEN
          CALL GASDN2(I12,RHO1,RHO2,WORK(KVEC1P),WORK(KVEC2P),WORK(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)),
       &       NAEL,IATP,NBEL,IBTP,IOCTPA,IOCTPB,NOCTPA,NOCTPB,
       &       NSMST,NSMOB,NSMSX,NSMDX,MXPNGAS,NOBPTS,IOBPTS,      
       &       MAXK,MAXI,LSCR1,LSCR1,WORK(KCSCR),WORK(KSSCR),
  