             gp[I_q - 1] -= gupdate;   
         }
 
-        __threadfence_block();
+        //__syncthreads();
 
         for(mq = 0; mq < meqn; mq++)
         {

        }


    /*     transverse-x

            |     |     | 
            |     |     | 
        ----|--O--|-----|-----
            |  ^--X     | 
            |     X     |
            |  v--X     |
        ----|--O--|-----|-----
            |     |     |
            |     |     |

    */              

        rpt2(0,meqn,mwaves,maux,ql,qr,aux1,aux2,aux3,0,0,amdq,bmasdq,bpasdq);


            fp[I_q] -= gupdate;
        }

        //__syncthreads();

        for(mq = 0; mq < meqn; mq++)
        {
