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

                fm[I_q + 1] -= gupdate;
                fp[I_q + 1] -= gupdate;
            }   
            
            __syncthreads();

        } /* Thread conditional */

                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);
