anchor.diff:        INCLUDE 'glbbas.inc'
         INCLUDE 'cgas.inc'
   *
  -      CALL EXCCLS2_S(NGAS,WORK(KIOCCLS),NCLS,IACTIN,IACTOUT,IEXC,
  +      CALL EXCCLS2_S(NGAS,int_mb(KIOCCLS),NCLS,IACTIN,IACTOUT,IEXC,
        &               IBASSPC_MX,IBASSPC)
   *
         RETURN 
  
positive.code:         CALL COPVCD(LUC,LUSC1,WORK(KVEC1),0,LBLK)
          CALL COPVCD(LUSC1,LUSC2,WORK(KVEC1),1,LBLK)
          XDUM = 0.0D0
          CALL DENSI2(1,WORK(KRHO1),WORK(KRHO2),
       &        WORK(KVEC1),WORK(KVEC2),LUSC1,LUSC2,EXPS2,
       &        0,XDUM,XDUM,XDUM,XDUM,1)
  *. Obtain MO-MO transformation matrix 
  
negative.code: *
        IRHO1SM = 1
  *. Extract symmetry blocks from complete one-electron density 
        CALL REORHO1(dbl_mb(KRHO1),dbl_mb(KLRHO1S),IRHO1SM,1)
  *. Number of elements in symmetry blocks of integrals and density
        LRHO1S = 0
        DO ISM = 1, NSMOB
  