             fm[i] = amdq[mq];
         }
 
-#if 0
         for (m = 0; m < meqn*mwaves; m++)
         {
             int i = I_waves + m;

            }

        }
#endif        
    }

    if (ix < mx + 2*(mbc-1) && iy < my + 2*mbc-1) 

    double amdq[MEQN];
    double apdq[MEQN];

    int mq, mw,m;
    int x_stride_q, y_stride_q, I_q;
    int x_stride_aux, y_stride_aux, I_aux;
    int x_stride_waves, y_stride_waves, I_waves;
    int x_stride_s, y_stride_s, I_speeds;

    int ix = threadIdx.x + blockIdx.x * blockDim.x;
    int iy = threadIdx.y + blockIdx.y * blockDim.y;

    if (ix < mx + 2*mbc-1 && iy < my + 2*(mbc-1)) 
    {
        x_stride_q = meqn;
        y_stride_q = (2 * mbc + mx) * x_stride_q;
        I_q = (ix + mbc-1) * x_stride_q + (iy + mbc-1) * y_stride_q;

        x_stride_aux = maux;
        y_stride_aux = (2 * mbc + mx) * x_stride_aux;
        I_aux = (ix + mbc-1) * x_stride_aux + (iy + mbc-1) * y_stride_aux;

        x_stride_waves = mwaves*meqn;
        y_stride_waves = (2 * mbc + mx) * x_stride_waves;
        I_waves = (ix + mbc-1) * x_stride_waves + (iy + mbc-1) * y_stride_waves;

        x_stride_s = mwaves;
        y_stride_s = (2 * mbc + mx) * x_stride_s;
        I_speeds = (ix + mbc-1) * x_stride_s + (iy + mbc-1) * y_stride_s;

        ql[0] = qold[I_q - x_stride_q];
        qr[0] = qold[I_q];
