     cudaDeviceSynchronize();
     CHECK(cudaPeekAtLastError());
 
-
-    // collect max cfl numbers
+    /* -------------------------------- Compute CFL ------------------------------------*/ 
     int n = 2*mwaves*(2*mbc+mx)*(2*mbc+my);
     int maxidx;
 

#include "../fc2d_cudaclaw.h"
#include "cudaclaw_allocate.h"

#include <fclaw2d_patch.h>
#include <fclaw2d_global.h>
#include <fclaw2d_clawpatch.h>

#include "../fc2d_cudaclaw_options.h"
#include "src/patches/clawpatch/fclaw2d_clawpatch_options.h"
#include <fclaw2d_global.h>

#include "cudaclaw_update_q.h"
#include "cudaclaw_flux2.h"

#include "../fc2d_cudaclaw_check.cu"  /* CHECK defined here */
#include <cublas_v2.h>
    
double cudaclaw_step2_batch(fclaw2d_global_t *glob,
        cudaclaw_fluxes_t* array_fluxes_struct, 

    return maxcfl;
}

