anchor.diff:  *. Set up the basespace for the various occupation classes
   C          Z_BASSPC_FOR_ALL_OCCLS(IOCCLS,NOCCLS,IBASSPC)
         CALL Z_BASSPC_FOR_ALL_OCCLS(int_mb(KIOCCLS),NOCCLS_MAX,
  -     &      WORK(KBASSPC))
  +     &      int_mb(KBASSPC))
   *. Dimensions for all occupation classes
   
   *
  
positive.code: C IFRMR(WORK,IROFF,IELMNT)
           DO JOCSBCLS = 1, NOCSBCLST
  C?         WRITE(6,*) ' JOCSBCLS = ', JOCSBCLS
             LOCC = IFRMR(dbl_mb(KLSBCNF),1,JOCSBCLS)
  C?         WRITE(6,*) ' LOCC = ', LOCC
             CALL MEMMAN(LPOINT,LOCC,'ADDL  ',1,'OCSBCN')
  C?         WRITE(6,*) ' LPOINT = ', LPOINT
  
negative.code:          CALL MEMMAN(KINT_2EINI,NINT2,'ADDS  ',2,'INT2_I')
        END IF
  *. Pointers to symmetry block of integrals
        CALL MEMMAN(KPINT1,NBINT1,'ADDS  ',2,'PINT1 ')  !done
        CALL MEMMAN(KPINT2,NBINT2,'ADDS  ',2,'PINT2 ')  !done
  *. Pointers to nonsymmetric one-electron integrals
        DO ISM = 1, NSMOB
  *. triangular packed
          CALL MEMMAN(KPGINT1(ISM),NSMOB,'ADDS  ',2,'PGINT1')  !done
  *. no packing
          CALL MEMMAN(KPGINT1A(ISM),NSMOB,'ADDS  ',2,'PGIN1A')  !done
        END DO
  *. Symmetry of last index as a function of initial index
        CALL MEMMAN(KLSM1,NBINT1,'ADDS  ',2,'LSM1   ')  !done
        CALL MEMMAN(KLSM2,NBINT2,'ADDS  ',2,'LSM2   ')  !done
  *.3 One-body density   
        CALL MEMMAN(KRHO1,NACOB ** 2,'ADDS  ',2,'RHO1  ')  !done
  *.3.1: One-body spin density
        IF(ISPNDEN.GE.1) THEN
          CALL MEMMAN(KSRHO1,NACOB **2, 'ADDS  ',2,'SRHO1 ')
  