                 fm[I_q - ys + 1] -= gupdate;
                 fp[I_q - ys + 1] -= gupdate;
             }
-
             __syncthreads();
+
 
             /* idir = 1; imp = 1;  Re-use amdp, apdq */
             rpt2(1,meqn,mwaves,maux,qd,qr,aux1,aux2,aux3,1,bpdq,bmasdq,bpasdq);

                gm[I_q - 1 + ys] -= gupdate;
                gp[I_q - 1 + ys] -= gupdate;
            }

            __syncthreads();



                gm[I_q + ys] -= gupdate;
                gp[I_q + ys] -= gupdate;
            }
            __syncthreads();
            

            /* ----------------------- Transverse : Y-faces --------------------------- */
