anchor.diff:          WRITE(6,*) 
           WRITE(6,*) ' IOPEN, IOFFCD = ',IOPEN, IOFFCD
           WRITE(6,*) ' The DTOC array '
  -        CALL WRTMAT(WORK(KDTOC + IOFFCD -1),NPSD,NPCSF,
  +        CALL WRTMAT(dbl_mb(KDTOC + IOFFCD -1),NPSD,NPCSF,
        &  NPSD, NPCSF)
         END IF
   *
  
positive.code:         IALPHA = (IOPEN+MS2)/2
          IDET = NPCMCNF(IOPEN+1)
  C?      WRITE(6,*) ' IOPEN, IDET = ', IOPEN, IDET
          CALL REO_PTDET(IOPEN,IALPHA,WORK(KZ_PTDT(ITP)),
       &                 WORK(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)
  
negative.code: C?      CALL IWRTMA(ISIGN,1,NPSD,1,NPSD)
  C?      WRITE(6,*) ' Result after ISIGN '
  C?      CALL WRTMAT(C_SD,1,NPSD,1,NPSD)
          CALL MATVCC2(dbl_mb(KDTOC + IOFFCD -1),C_SD,C_CSF,NPSD,NPCSF,
       &               1,FACTORC)
        END IF
  *
  