         }  
         __syncthreads();
     }  
+
+    printf("here ...4\n");
 
     /* ------------------------- First order final update ----------------------------- */
 

        }
    }


    maxcflblocks[blockIdx.z] = BlockReduce(temp_storage).Reduce(maxcfl,cub::Max());

    //__syncthreads();  /* Does block reduce take care of this sync? */


                }
            }
        }

        printf("here ...3\n");

        ifaces_x = mx + 2;  
        ifaces_y = my + 1;
