 #include"global_cuda.h"
 #include"roe_cuda.h"
 
-
+#ifdef DE //PRESSURE_DE
+#include"hydro_cuda.h"
+#endif
 
 /*! \fn Calculate_Roe_Fluxes_CUDA(Real *dev_bounds_L, Real *dev_bounds_R, Real *dev_flux, int nx, int ny, int nz, int n_ghost, Real gamma, Real *dev_etah, int dir, int n_fields)
  *  \brief Roe Riemann solver based on the version described in Stone et al, 2008. */

    vxl = mxl / dl;
    vyl = myl / dl;
    vzl = mzl / dl;
    pl  = (El - 0.5*dl*(vxl*vxl + vyl*vyl + vzl*vzl)) * (gamma - 1.0);
    pl  = fmax(pl, (Real) TINY_NUMBER);
    #ifdef SCALAR
    for (int i=0; i<NSCALARS; i++) {

  Real test0, test1, test2, test3, test4;
  int hlle_flag = 0;
  #ifdef DE
  Real dgel, gel, dger, ger, f_ge_l, f_ge_r, E_kin;
  #endif
  #ifdef SCALAR
  Real dscalarl[NSCALARS], scalarl[NSCALARS], dscalarr[NSCALARS], scalarr[NSCALARS], f_scalar_l[NSCALARS], f_scalar_r[NSCALARS];
