             {
                 for(mw = 0; mw < mwaves; mw++)
                 {
-                    /* X-faces */
 
                     I_speeds = I + mw*zs;
                     s[mw] = speeds[I_speeds];

                        cqyy = fabs(s[mw])*(1.0 - fabs(s[mw])*dtdy)*wave[mq];
                        gm[I_q] += 0.5*cqyy;   
                        gp[I_q] += 0.5*cqyy;   
                        if (order[1] > 0)
                        {
                            /* Propagate second order corrections 
                               in transverse dir. */
                            bmdq_trans[I_q] += cqyy;     
                            bpdq_trans[I_q] -= cqyy;      
                        }
                    }                
                }  /* End of mwaves loop */
            } /* End of thread conditional */
        } /* End of thread loop */
        __syncthreads();
    } /* End of check on order[0] == 2 */



    if (order[1] == 0)


            for(m = 0; m < maux; m++)
            {
                I_aux = I + m*zs;
                auxl[m] = aux[I_aux - 1];
                auxr[m] = aux[I_aux];
