anchor.diff:            DO JROOT = 1, IROOT
               CALL REWINO(LUSCR)
               IF(ICISTR.EQ.1) THEN
  -              CALL FRMDSC(WORK(KLBLI),NCOMBI,-1,LUIN,IMZERO,IAMPACK)
  -              CALL  TODSC(WORK(KLBLI),NCOMBI,-1,LUSCR)
  +              CALL FRMDSC(dbl_mb(KLBLI),NCOMBI,-1,LUIN,IMZERO,IAMPACK)
  +              CALL  TODSC(dbl_mb(KLBLI),NCOMBI,-1,LUSCR)
               ELSE
  -              CALL COPVCD(LUIN,LUSCR,WORK(KLBLI),0,-1)
  +              CALL COPVCD(LUIN,LUSCR,dbl_mb(KLBLI),0,-1)
               END IF
             END DO
             CALL REWINO(LUSCR)
  
positive.code: c      IMPLICIT REAL*8(A-H,O-Z)
  *
  c      INCLUDE 'mxpdim.inc'
        INCLUDE 'wrkspc.inc'
        INCLUDE 'orbinp.inc'
        INCLUDE 'strbas.inc'
  
negative.code: *. Input
  *     (and /LUCINP/ not occuring here )
  c      INCLUDE 'mxpdim.inc'
  #include "errquit.fh"
  #include "mafdecls.fh"
  #include "global.fh"
        INCLUDE 'wrkspc.inc'
        INCLUDE 'orbinp.inc'
        INCLUDE 'cgas.inc'
  