 
             for(mw = 0; mw < mwaves; mw++)
             {
+                maxcfl = max(maxcfl,fabs(s[mw])*dtdy);
+
                 I_speeds = I + (mwaves + mw)*zs;
                 speeds[I_speeds] = s[mw];
                 for(mq = 0; mq < meqn; mq++)

                        cqxx = fabs(s[mw])*(1.0 - fabs(s[mw])*dtdx)*wave[mq];
                        fm[I_q] += 0.5*cqxx;   
                        fp[I_q] += 0.5*cqxx;                               
                        if (order[1] > 0)
                        {
                            /* Propagate second order corrections 
                               in transverse dir. */
                            amdq_trans[I_q] += cqxx;   
                            apdq_trans[I_q] -= cqxx;      
                        }
                    }

                    /* Y-faces */

                fm[I_q - ys] -= gupdate;        
                fp[I_q - ys] -= gupdate;
            }
        } 
    } 

    __syncthreads();

