                 for(mq = 0; mq < meqn; mq++)
                 {
                     I_q = I + mq*zs;
-                    cqxx += fabs(s[mw])*(1.0 - fabs(s[mw])*dtdx)*wlimitr*wave[mq];
-                    fm[I_q] += 0.5*cqxx;   /* amdq + cqxx */
-                    fp[I_q] += 0.5*cqxx;   /* apdq - cqxx */
+                    cqxx_x += fabs(s[mw])*(1.0 - fabs(s[mw])*dtdx)*wlimitr*wave[mq];
+                    fm[I_q] += 0.5*cqxx_x;   /* amdq + cqxx */
+                    fp[I_q] += 0.5*cqxx_x;   /* apdq - cqxx */
                 }
 
 

                apdq[mq] = fp[I_q];
            }
#endif            
            double cqxx = 0;
            for(mw = 0; mw < mwaves; mw++)
            {
                /* x-faces */

                for(mq = 0; mq < meqn; mq++)
                {
                    I_q = I + mq*zs;
                    cqxx_y += fabs(s[mw])*(1.0 - fabs(s[mw])*dtdy)*wlimitr*wave[mq];
                    gm[I_q] += 0.5*cqxx_y;   /* amdq + cqxx */
                    gp[I_q] += 0.5*cqxx_y;   /* apdq - cqxx */
                }


