anchor.diff:          DO IGRP = 1, NSPGPFTP(ITP)
             IGRPABS = IGRP-1 + IBSPGPFTP(ITP)
             CALL NSTPTP_GAS_NEW(NGAS,ISPGPFTP(1,IGRPABS),
  -     &                    WORK(KNSTSGP(1)),NSMST,
  +     &                    int_mb(KNSTSGP(1)),NSMST,
        &                    WORK(KNSTSO(ITP)),IGRP,MXNSTRFSG,
        &                    NSMCLS,NSMCLSE,NSMCLSE1,NSTR_AS)
   *
  
positive.code:       KACGRP = ITPFGS(IOBTP)
  *. Number of strings per symmetry distribution      
        DO IGAS = 1, NGAS
          CALL ICOPVE2(WORK(KNSTSGP(1)),(ITPFGS(IGAS)-1)*NSMST+1,NSMST,
       &               NNSTSGP(1,IGAS))
          CALL ICOPVE2(WORK(KISTSGP(1)),(ITPFGS(IGAS)-1)*NSMST+1,NSMST,
       &               IISTSGP(1,IGAS))
        END DO
  *
  
negative.code:       NSTPSM(1) = 1
  *
        DO IGRP = 1, NGRPA
          CALL ICOPVE2(int_mb(KNSTSGP(1)),(IGRPA(IGRP)-1)*NSMST+1,NSMST,
       &               NSTFTP)
          CALL ICOPVE(NSTPSM,NSTFTP2,NSMST)
          CALL DIM_PROD_TWO_STGRPS(NSTPSM,NSTFTP2,NSTFTP,NSMST)
  