     vy_L = Q_Ry[2*n_cells + jmo] / Q_Ry[jmo];
     vz_R = Q_Lz[3*n_cells + id]  / Q_Lz[id]; 
     vz_L = Q_Rz[3*n_cells + kmo] / Q_Rz[kmo]; 
+    
     #endif
 
     // update the conserved variable array

  Real d, d_inv, vx, vy, vz;
  #endif
  #ifdef DE
  // Real vx_imo, vx_ipo, vy_jmo, vy_jpo, vz_kmo, vz_kpo;
  // int ipo, jpo, kpo;
  Real  P, E, E_kin, GE;
  Real vx_L, vx_R, vy_L, vy_R, vz_L, vz_R;
  #endif

    dev_conserved[(n_fields-1)*n_cells + id] += dtodx * (dev_F_x[(n_fields-1)*n_cells + imo] - dev_F_x[(n_fields-1)*n_cells + id])
                                  +  dtody * (dev_F_y[(n_fields-1)*n_cells + jmo] - dev_F_y[(n_fields-1)*n_cells + id])
                                  +  dtodz * (dev_F_z[(n_fields-1)*n_cells + kmo] - dev_F_z[(n_fields-1)*n_cells + id])
                                  +  0.5*P*(dtodx*(vx_imo-vx_ipo) + dtody*(vy_jmo-vy_jpo) + dtodz*(vz_kmo-vz_kpo));
                                  // -  P * ( dtodx * ( vx_R - vx_L ) + dtody * ( vy_R - vy_L ) + dtodz * ( vz_R - vz_L ) );
                                  
    #endif
    
