anchor.diff:  * total symmetry of the operator
   c      IMPLICIT REAL*8(A-H,O-Z)
   c      INCLUDE 'mxpdim.inc'
  +#include "errquit.fh"
  +#include "mafdecls.fh"
  +#include "global.fh"
         INCLUDE 'wrkspc.inc'
         INCLUDE 'glbbas.inc'
         INCLUDE 'lucinp.inc'
  
positive.code: *
  c      IMPLICIT REAL*8(A-H,O-Z)
  c      INCLUDE 'mxpdim.inc'
        INCLUDE 'wrkspc.inc'
        INCLUDE 'crun.inc'
        INCLUDE 'glbbas.inc'
  
negative.code: *.0: Pointers to one-integrals, all symmetries, Lower half matrices
        DO ISM = 1, NSMOB
          CALL PNT2DM(1,NSMOB,NSMSX,ADSXA,NTOOBS,NTOOBS,
       &       ISM  ,ISL1,dbl_mb(KPGINT1(ISM)),MXPOBS)
        END DO
  *.0.5: Pointers to one-electron integrals, all symmetries, complete form
        DO ISM = 1, NSMOB
          CALL PNT2DM(0,NSMOB,NSMSX,ADSXA,NTOOBS,NTOOBS,
       &       ISM  ,ISL1,dbl_mb(KPGINT1A(ISM)),MXPOBS)
        END DO
  *.1: Number of one-electron integrals
        CALL PNT2DM(1,NSMOB,NSMSX,ADSXA,NTOOBS,NTOOBS,
  