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

    // left and right fluxes 
    f_d_l  = mxl;
    f_mx_l = mxl*vxl + pl;
    // f_mx_l = mxl*vxl;
    f_my_l = myl*vxl;
    f_mz_l = mzl*vxl;
    f_E_l  = (El + pl)*vxl;
    // f_E_l  = (El)*vxl;
    #ifdef DE
    f_ge_l = dgel*vxl;
    #endif


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