     }
     #endif
     #ifdef DE
-    ge_i = dev_conserved[(n_fields-1)*n_cells + id] / d_i;
+    ge_i = dge / d_i;
     #endif
     // cell i-1
     if (dir == 0) id = xid-1 + yid*nx + zid*nx*ny;

#include"global.h"
#include"global_cuda.h"
#include"plmp_cuda.h"


/*! \fn __global__ void PLMP_cuda(Real *dev_conserved, Real *dev_bounds_L, Real *dev_bounds_R, int nx, int ny, int nz, int n_ghost, Real dx, Real dt, Real gamma, int dir, int n_fields)

    }
    #endif
    #ifdef DE
    ge_ipo =  dge / d_ipo;
    #endif


