     pr  = (Er - 0.5*dr*(vxr*vxr + vyr*vyr + vzr*vzr)) * (gamma - 1.0);
     #endif//DE
     pr  = fmax(pr, (Real) TINY_NUMBER);    
-    #ifdef SCALAR
-    for (int i=0; i<NSCALARS; i++) {
-      scr[i] = dscr[i] / dr;
-    }
-    #endif
+    // #ifdef SCALAR
+    // for (int i=0; i<NSCALARS; i++) {
+    //   scr[i] = dscr[i] / dr;
+    // }
+    // #endif
     // #ifdef DE
     // ger = dger / dr;
     // #endif

  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];
  #endif
  
  // Real etah = 0;

    pl  = (El - 0.5*dl*(vxl*vxl + vyl*vyl + vzl*vzl)) * (gamma - 1.0);
    #endif//DE
    pl  = fmax(pl, (Real) TINY_NUMBER);
    // #ifdef SCALAR
    // for (int i=0; i<NSCALARS; i++) {
    //   scl[i] = dscl[i] / dl;
    // }
    // #endif
    // #ifdef DE
    // gel = dgel / dl;
    // #endif
