 #include"global_cuda.h"
 #include"hllc_cuda.h"
 
+#ifdef DE //PRESSURE_DE
+#include"hydro_cuda.h"
+#endif
 
 
 /*! \fn Calculate_HLLC_Fluxes_CUDA(Real *dev_bounds_L, Real *dev_bounds_R, Real *dev_flux, int nx, int ny, int nz, int n_ghost, Real gamma, int dir, int n_fields)

  Real f_d, f_mx, f_my, f_mz, f_E;
  Real Sl, Sr, Sm, cfl, cfr, ps;
  #ifdef DE
  Real dgel, dger, gel, ger, gels, gers, f_ge_l, f_ge_r, f_ge;
  #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];

    vxr = mxr / dr;
    vyr = myr / dr;
    vzr = mzr / dr;
    #ifdef DE //PRESSURE_DE
    E_kin = 0.5 * dr * ( vxr*vxr + vyr*vyr + vzr*vzr );
    pr = Get_Pressure_From_DE( Er, Er - E_kin, dger, gamma );
    #else
    pr  = (Er - 0.5*dr*(vxr*vxr + vyr*vyr + vzr*vzr)) * (gamma - 1.0);
    #endif //PRESSURE_DE
    pr  = fmax(pr, (Real) TINY_NUMBER);    
    #ifdef SCALAR
    for (int i=0; i<NSCALARS; i++) {
