             /* Start at first non-ghost interior cell */
             int I = (iy + mbc)*ys + ix + mbc - 1;
 
+            double *s = start;
+            double *wave = s + mwaves;
             for(int mw = 0; mw < mwaves; mw++)
             {
                 /* ------------------------------- Y-directions --------------------------- */

    double *const aux3   = aux2   + 2*maux;       /* 2*maux      */
    double *const bmasdq = aux3   + 2*maux;       /* meqn        */
    double *const bpasdq = bmasdq + meqn;         /* meqn        */


    /* --------------------------------- Start code ----------------------------------- */

        }      
        __syncthreads(); /* Needed to be sure all aux variables are available below */ 
    } 
#endif    

    /* -------------------------- Compute fluctuations -------------------------------- */

