anchor.diff:        CALL QENTER('HCONF')
   *. Fetch K matrix
   C     CALL GTJK(XDUM1       ,H,NTOOB,XDUM3       ,IREOTS)
  -      CALL GTJK(WORK(KLZSCR),H,NTOOB,WORK(KLZSCR),IREOTS)
  +      CALL GTJK(dbl_mb(KLZSCR),H,NTOOB,dbl_mb(KLZSCR),IREOTS)
   *
         CALL HCONFDIA_BBS(NAEL,NBEL,IJAGRP,IJBGRP,
        &           IASM,IATP,IAOC,NIA,IBSM,IBTP,IBOC,NIB,
        &           JASM,JATP,JAOC,NJA,JBSM,JBTP,JBOC,NJB,
  -     &           NTOOB,H,IOBPTS,NOBPT,IREOTS,WORK(KLZ(1)),
  -     &           WORK(KLZ(2)),WORK(KLZSCR),ISMFTO,NGAS,
  -     &           WORK(KLOCSTR(1)),WORK(KLOCSTR(2)),
  -     &           WORK(KLREO(1)),WORK(KLREO(2)),CB,SB) 
  +     &           NTOOB,H,IOBPTS,NOBPT,IREOTS,int_mb(KLZ(1)),
  +     &           int_mb(KLZ(2)),dbl_mb(KLZSCR),ISMFTO,NGAS,
  +     &           int_mb(KLOCSTR(1)),int_mb(KLOCSTR(2)),
  +     &           int_mb(KLREO(1)),int_mb(KLREO(2)),CB,SB) 
         CALL QEXIT('HCONF')
   *
         RETURN
  
positive.code: *. Generate occupation of K STRINGS
  
         CALL GETSTR2_TOTSM_SPGP(KGRP,NGAS,KSM,NELK,NSTRK,
       &                        WORK(KLLOC),NOCOB,
       &                        0,ISUM,IDUM)
  C     GETSTR2_TOTSM_SPGP(IGRP,NIGRP,ISPGRPSM,NEL,NSTR,ISTR,
  C    &                              NORBT,IDOREO,IZ,IREO)
  
negative.code: *
        CALL ADAADAS1_GAS(NK,I1,XI1S,LI1,
       &          IIOB,NIOB,IAC,JJOB,NJOB,JAC,
       &          int_mb(KLLOC),NELK,NSTRK,int_mb(KLLREO),int_mb(KLLZ),
       &          NOCOB,KMAX,KMIN,IEND,SCLFAC,NSTRIA(I12))
  *
         IF(NTEST.GE.1000) THEN
           WRITE(6,*) ' Reorder array after ADAADAS1'                  
           CALL IWRTMA(int_mb(KLLREO),1,NSTRI,1,NSTRI)
         END IF
   9999 CONTINUE
  *
  