anchor.diff:  **. Diagonal of one-body integrals and coulomb and exchange integrals
   *   One-body integrals stored in KINT1 are used
   *
  -      CALL GT1DIA(WORK(KLH1D))
  -      CALL GTJK(WORK(KLJ),WORK(KLK),NTOOB,WORK(KLSCR2),IREOTS)
  +      CALL GT1DIA(dbl_mb(KLH1D))
  +      CALL GTJK(dbl_mb(KLJ),dbl_mb(KLK),NTOOB,dbl_mb(KLSCR2),IREOTS)
   *
         IF( LUDIA .GT. 0 ) CALL REWINO(LUDIA)
         IF(IEXP_PH.EQ.0) THEN
   *. Good old version where PH is not used for diagonal
  -        CALL GASDIAS(NAEL,WORK(KLASTR),NBEL,WORK(KLBSTR),
  -     &       NACOB,DIAG,NSMST,WORK(KLH1D),
  -     &       WORK(KLXA),WORK(KLXB),WORK(KLSCR),WORK(KLJ),
  -     &       WORK(KLK),WORK(KNSTSO(IATP)),WORK(KNSTSO(IBTP)),
  +        CALL GASDIAS(NAEL,int_mb(KLASTR),NBEL,int_mb(KLBSTR),
  +     &       NACOB,DIAG,NSMST,dbl_mb(KLH1D),
  +     &       dbl_mb(KLXA),dbl_mb(KLXB),dbl_mb(KLSCR),dbl_mb(KLJ),
  +     &       dbl_mb(KLK),int_mb(KNSTSO(IATP)),int_mb(KNSTSO(IBTP)),
        &       LUDIA,ECORE,PLSIGN,PSSIGN,IPRDIA,NTOOB,ICISTR,
  -     &       WORK(KLRJKA),I12,IBLTP,NBLOCK,IBLKFO)
  +     &       dbl_mb(KLRJKA),I12,IBLTP,NBLOCK,IBLKFO)
         ELSE
   *. use also ph in diagonal construction
  -        CALL GASDIAS_PH(NAEL,WORK(KLASTR),NBEL,WORK(KLBSTR),
  -     &       NACOB,DIAG,NSMST,WORK(KLH1D),
  -     &       WORK(KLXA),WORK(KLXB),WORK(KLSCR),WORK(KLJ),
  -     &       WORK(KLK),WORK(KNSTSO(IATP)),WORK(KNSTSO(IBTP)),
  +        CALL GASDIAS_PH(NAEL,int_mb(KLASTR),NBEL,int_mb(KLBSTR),
  +     &       NACOB,DIAG,NSMST,dbl_mb(KLH1D),
  +     &       dbl_mb(KLXA),dbl_mb(KLXB),dbl_mb(KLSCR),dbl_mb(KLJ),
  +     &       dbl_mb(KLK),int_mb(KNSTSO(IATP)),int_mb(KNSTSO(IBTP)),
        &       LUDIA,ECORE,PLSIGN,PSSIGN,IPRDIA,NTOOB,ICISTR,
  -     &       WORK(KLRJKA),I12,IBLTP,NBLOCK,IBLKFO,IPHGAS,ITPFTO,
  -     &       WORK(KLPHSTR),NPHELFSPGP,IOCTPA,IOCTPB,NINOB,WORK(KLISCR))
  +     &       dbl_mb(KLRJKA),I12,IBLTP,NBLOCK,IBLKFO,IPHGAS,ITPFTO,
  +     &       int_mb(KLPHSTR),NPHELFSPGP,IOCTPA,IOCTPB,NINOB,
  +     &       int_mb(KLISCR))
         END IF
   
   *.Flush local memory
  
positive.code: CINA &            ISMFTO,IBSO,NACOB)
  COLD  CALL GT1DIS(H1DIA,IREOTS(1),WORK(KPINT1),WORK(KINT1O),
  COLD &            ISMFTO,IBSO,NTOOB)
        CALL GT1DIS(H1DIA,IREOTS(1),WORK(KPINT1),WORK(KINT1),
       &            ISMFTO,IBSO,NTOOB)
  *
  
negative.code: *
  c      IMPLICIT REAL*8(A-H,O-Z)
  c      INCLUDE 'mxpdim.inc'
  #include "mafdecls.fh"
  #include "global.fh"
        INCLUDE 'wrkspc.inc'
   
  *.GLobal pointers
  