 #endif    
 
     /* -------------------------- Compute fluctuations -------------------------------- */
-
-#if 0
-    ifaces_x = mx + 2*mbc-1;
-    ifaces_y = my + 2*mbc-1;
-    num_ifaces = ifaces_x*ifaces_y;
-#endif    
 
     for(int thread_index = threadIdx.x; thread_index < num_ifaces; thread_index += blockDim.x)
     {



    /* ---------------------- Second order corrections and limiters --------------------*/  
    
    if (order[0] == 2)
    {
        int ifaces_x = mx + 1;  
        int ifaces_y = my + 2;
        int num_ifaces = ifaces_x*ifaces_y;

        for(int thread_index = threadIdx.x; thread_index < num_ifaces; thread_index += blockDim.x)
        { 


    double maxcfl = 0;

    __shared__ int ifaces_x = mx + 2*mbc-1;
    __shared__ int ifaces_y = my + 2*mbc-1;
    __shared__ int num_ifaces = ifaces_x*ifaces_y;

#if 0
    if (b4step2 != NULL)
