anchor.diff:  *. Obtain L and R in SD/Combination basis in vector KCOMVEC1_SD, KCOMVEC2_SD
   *. A scratch file is used to allow calc to be done with two complete vectors
   *L(SD)
  -         CALL CSDTVCM(L,WORK(KCOMVEC1_SD),WORK(KCOMVEC2_SD),
  +         CALL CSDTVCM(L,KCOMVEC1_SD,KCOMVEC2_SD,
        &        1,0,ISSM,ISSPC,2)
            CALL FILEMAN_MINI(LUSCX,'ASSIGN')
   *. Save result in LUSCX
  -         CALL REWINO(LUSCX)
  -         CALL TODSC(WORK(KCOMVEC1_SD),NSVAR,NSVAR,LUSCX)
  +cnw      CALL REWINO(LUSCX)
  +cnw      CALL TODSC(WORK(KCOMVEC1_SD),NSVAR,NSVAR,LUSCX)
  +         if(.not.dra_write_section(.false.,KCOMVEC1_SD,1,NSVAR,1,1,
  +     &      LUSCX,1,NSVAR,1,1)) call errquit('dra err',911)
   *R(SD)
  -         CALL CSDTVCM(R,WORK(KCOMVEC2_SD),WORK(KCOMVEC1_SD),
  +         CALL CSDTVCM(R,KCOMVEC2_SD,KCOMVEC1_SD,
        &        1,0,ICSM,ICSPC,2)
   *Retrieve L(SD) from scratch 
  -         CALL REWINO(LUSCX)
  -         CALL FRMDSCO(WORK(KCOMVEC1_SD),NSVAR,NSVAR,LUSCX,IMZERO)
  +cnw      CALL REWINO(LUSCX)
  +cnw      CALL FRMDSCO(WORK(KCOMVEC1_SD),NSVAR,NSVAR,LUSCX,IMZERO)
  +         if (.not.dra_read_section(.false.,KCOMVEC2_SD,1,NSVAR,1,1,
  +     &       LUSCX,1,NSVAR,1,1)) call errquit('dra err',911)
            CALL FILEMAN_MINI(LUSCX,'FREE  ')
          ELSE
   *. Not in core
  
positive.code:      &       dbl_mb(KINSCR),
       &       MXPOBS,IPRDEN,dbl_mb(KRHO1S),-1,-1,
       &       PSSIGN,PSSIGN,dbl_mb(KRHO1P),dbl_mb(KXNATO),
       &       NBATCHL,WORK(KSLBT),WORK(KSLEBT),WORK(KSI1BT),
       &       WORK(KSIBT),
       &       NBATCHR,WORK(KCLBT),WORK(KCLEBT),WORK(KCI1BT),
       &       int_mb(KCIBT),int_mb(KCONSPA),int_mb(KCONSPB),
       &       dbl_mb(KLSCLFCL),dbl_mb(KLSCLFCR),S2_TERM1,IUSE_PH,IPHGAS,
       &       IDOSRHO1,SRHO1,IDOSRHO2,RHO2AA,RHO2AB,RHO2BB,
  
negative.code:      &       dbl_mb(KINSCR),
       &       MXPOBS,IPRDEN,dbl_mb(KRHO1S),-1,-1,
       &       PSSIGN,PSSIGN,dbl_mb(KRHO1P),dbl_mb(KXNATO),
       &       NBATCHL,int_mb(KSLBT),WORK(KSLEBT),WORK(KSI1BT),
       &       int_mb(KSIBT),
       &       NBATCHR,int_mb(KCLBT),WORK(KCLEBT),WORK(KCI1BT),
       &       int_mb(KCIBT),int_mb(KCONSPA),int_mb(KCONSPB),
       &       dbl_mb(KLSCLFCL),dbl_mb(KLSCLFCR),S2_TERM1,IUSE_PH,IPHGAS,
       &       IDOSRHO1,SRHO1,IDOSRHO2,RHO2AA,RHO2AB,RHO2BB,
  