anchor.diff:  C     COMMON/COMJEP/MXACJ,MXACIJ,MXAADST
           IF(NSTA*NSTB*NIAC*NKAC.NE.0)
        &  CALL ADAST_GASSM(NSTB,NSTA,IKAC,IIAC,IBSTRINI,KSTRBS,   
  -     &                 WORK(KSTSTM(KACGRP,1)),WORK(KSTSTM(KACGRP,2)),
  -     &                 IBORBSPS,IBORBSP,NORBTS,NKAC,NKACT,NIAC,
  -     &                 NKSTR,KBSTRIN,NELB,NACGSOB,I1,XI1S,SCLFAC,IAC,
  -     &                 LROW_IN,IEC)
  +     &            int_mb(KSTSTM(KACGRP,1)),int_mb(KSTSTM(KACGRP,2)),
  +     &            IBORBSPS,IBORBSP,NORBTS,NKAC,NKACT,NIAC,
  +     &            NKSTR,KBSTRIN,NELB,NACGSOB,I1,XI1S,SCLFAC,IAC,
  +     &            LROW_IN,IEC)
           KSTRBS = KSTRBS + NKSD     
   C       IF(NGASL-1.GT.0) GOTO 1000
           GOTO 1000
  
positive.code:       END IF
  *
  *. Batches  of C vector
        CALL PART_CIV2(IDC,dbl_mb(KPCBLTP),WORK(KNSTSO(IATP)),
       &              WORK(KNSTSO(IBTP)),
       &              NOCTPA,NOCTPB,NSMST,LBLOCK,dbl_mb(KLCIOIO),
       &              ISMOST(1,ISM),
       &              NBATCH,int_mb(KPCLBT),int_mb(KPCLEBT),
  
negative.code:        IF(NOCSF.EQ.1) THEN
          CALL GASDN2(I12,RHO1,RHO2,WORK(KVEC1P),WORK(KVEC2P),WORK(KC2),
       &       WORK(KCIOIO),int_mb(KSIOIO),ISMOST(1,ICSM),
       &       ISMOST(1,ISSM),int_mb(KCBLTP),int_mb(KSBLTP),NACOB,
       &       int_mb(KNSTSO(IATP)),int_mb(KISTSO(IATP)),
       &       int_mb(KNSTSO(IBTP)),int_mb(KISTSO(IBTP)),
       &       NAEL,IATP,NBEL,IBTP,IOCTPA,IOCTPB,NOCTPA,NOCTPB,
       &       NSMST,NSMOB,NSMSX,NSMDX,MXPNGAS,NOBPTS,IOBPTS,      
       &       MAXK,MAXI,LSCR1,LSCR1,WORK(KCSCR),WORK(KSSCR),
  