 
             if (ix < mx + 1 && iy < my + 1)   /* Is this needed? */
             {
-                /* Limit waves */
                 for(mw = 0; mw < mwaves; mw++)
                 {
                     /* X-faces */

                I_q = I + mq*zs;
                fp[I_q] = -apdq[mq]; 
                fm[I_q] = amdq[mq];
                if (order[0] > 0)
                {
                    apdq_trans[I_q] = apdq[mq];  
                    amdq_trans[I_q] = amdq[mq];                    
                }
            }

            if (order[0] == 2)
            {
                for (m = 0; m < meqn*mwaves; m++)
                {
                    I_waves = I + m*zs;
                    waves[I_waves] = wave[m];
                }                
            }

            for (mw = 0; mw < mwaves; mw++)
            {
                I_speeds = I + mw*zs;
                speeds[I_speeds] = s[mw];
                cfl = abs(s[mw]*dtdx);
                if (cfl > maxcfl)
                {

        }
        return;
    }


    /* ------------------------ Transverse Propagation : X-faces ---------------------- */
