anchor.diff:  *
           IF(ICISTR.EQ.1) THEN
             IF(ICOPY.EQ.0) THEN
  -            CALL FRMDSC(WORK(KLBLU),NCOMBU,-1,LUUT,IMZERO,IAMPACK)
  -            CALL WRTMAT(WORK(KLBLU),1,NCOMBU,1,NCOMBUT)
  +            CALL FRMDSC(dbl_mb(KLBLU),NCOMBU,-1,LUUT,IMZERO,IAMPACK)
  +            CALL WRTMAT(dbl_mb(KLBLU),1,NCOMBU,1,NCOMBUT)
             ELSE
  -            CALL FRMDSC(WORK(KLBLU),NCOMBU,-1,LUIN,IMZERO,IAMPACK)
  -            CALL WRTMAT(WORK(KLBLU),1,NCOMBU,1,NCOMBUT)
  +            CALL FRMDSC(dbl_mb(KLBLU),NCOMBU,-1,LUIN,IMZERO,IAMPACK)
  +            CALL WRTMAT(dbl_mb(KLBLU),1,NCOMBU,1,NCOMBUT)
             END IF
           ELSE
  -          CALL WRTVCD(WORK(KLBLU),LUUT,0,-1)
  +          CALL WRTVCD(dbl_mb(KLBLU),LUUT,0,-1)
             END IF
           END DO
         END IF
  
positive.code: *. input  
  c      INCLUDE 'mxpdim.inc'
        INCLUDE 'wrkspc.inc'
        INCLUDE 'strbas.inc'
        INCLUDE 'cgas.inc'
  *. Input and Output ( NELFSPGP(MXPNGAS,MXPSTT) )
  
negative.code: *
  c      IMPLICIT REAL*8(A-H,O-Z)
  c      INCLUDE 'mxpdim.inc'
  #include "errquit.fh"
  #include "mafdecls.fh"
  #include "global.fh"
        INCLUDE 'wrkspc.inc'
  *./ORBINP/
        INCLUDE 'orbinp.inc'
  