 
     f_d_r  = mxr;
     f_mx_r = mxr*vxr + pr;
-    // f_mx_r = mxr*vxr;
     f_my_r = myr*vxr;
     f_mz_r = mzr*vxr;
     f_E_r  = (Er + pr)*vxr;
-    // f_E_r  = (Er)*vxr;
     #ifdef DE
     f_ge_r = dger*vxr;
     #endif

      scl[i] = dscl[i] / dl;
    }
    #endif
    #ifdef DE
    gel = dgel / dl;
    #endif
    vxr = mxr / dr;
    vyr = myr / dr;
    vzr = mzr / dr;

  Real f_d_r, f_mx_r, f_my_r, f_mz_r, f_E_r;
  // Real dls, drs, mxls, mxrs, myls, myrs, mzls, mzrs, Els, Ers;
  Real f_d, f_mx, f_my, f_mz, f_E;
  Real Sl, Sr, cfl, cfr;
  #ifdef DE
  Real dgel, dger, f_ge_l, f_ge_r, f_ge, E_kin;
  #endif
  #ifdef SCALAR
  Real dscl[NSCALARS], dscr[NSCALARS], scl[NSCALARS], scr[NSCALARS], scls[NSCALARS], scrs[NSCALARS], f_sc_l[NSCALARS], f_sc_r[NSCALARS], f_sc[NSCALARS];
