anchor.diff:  *
   c      IMPLICIT REAL*8(A-H,O-Z)
   c      INCLUDE 'mxpdim.inc'
  +#include "mafdecls.fh"
  +#include "global.fh"
         INCLUDE 'wrkspc.inc'
    
   *.GLobal pointers
  
positive.code: *
        RETURN
        END
        FUNCTION IFRMR(WORK,IROFF,IELMNT)
  *
  * An integer array is stored in real array WORK,
  * starting from WORK(IROFF). Obtain element
  * IELMNT of this array
  *
        INTEGER WORK(*)
  *
        INCLUDE 'irat.inc'
  *. offset when work is integer array
        IIOFF = 1 + IRAT * (IROFF-1)
        IFRMR = WORK(IIOFF-1+IELMNT)
  *
        RETURN
        END
  
negative.code: CINA &            ISMFTO,IBSO,NACOB)
  COLD  CALL GT1DIS(H1DIA,IREOTS(1),WORK(KPINT1),WORK(KINT1O),
  COLD &            ISMFTO,IBSO,NTOOB)
  c..dongxia
  c..  work(kpint1) and work(kint1) are static arrays about integrals
  c..  may need to convert to GA. Leave here for the time being.
        CALL GT1DIS(H1DIA,IREOTS(1),WORK(KPINT1),WORK(KINT1),
       &            ISMFTO,IBSO,NTOOB)
  *
  