anchor.diff:  *
         IF(IWAY.EQ.1) THEN
   C            MATVCC2(A,VIN,VOUT,NROW,NCOL,ITRNS,FACIN)
  -        CALL MATVCC2(WORK(KDTOC + IOFFCD -1),C_CSF,C_SD,NPSD,NPCSF,
  +        CALL MATVCC2(dbl_mb(KDTOC + IOFFCD -1),C_CSF,C_SD,NPSD,NPCSF,
        &                0,FACTORC)
   *
           IF(NTEST.GE.10000) THEN
  
positive.code:         WRITE(6,*) 
          WRITE(6,*) ' IOPEN, IOFFCD = ',IOPEN, IOFFCD
          WRITE(6,*) ' The DTOC array '
          CALL WRTMAT(WORK(KDTOC + IOFFCD -1),NPSD,NPCSF,
       &  NPSD, NPCSF)
        END IF
  *
  
negative.code:         IALPHA = (IOPEN+MS2)/2
          IDET = NPCMCNF(IOPEN+1)
  C?      WRITE(6,*) ' IOPEN, IDET = ', IOPEN, IDET
          CALL REO_PTDET(IOPEN,IALPHA,int_mb(KZ_PTDT(ITP)),
       &                 int_mb(KREO_PTDT(ITP)),IPDTCNF(IDTBS),
       &                 IDET,int_mb(KLSCR2) )
  C?     WRITE(6,*) ' Z array as delivered by REO_PTDET'
  C?     CALL IWRTMA(WORK(KZ_PTDT(ITP)),IOPEN,IALPHA,IOPEN,IALPHA)
  