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


            for(mw = 0; mw < mwaves; mw++)
            {
                I_speeds = I + mw*zs;
                s[mw] = speeds[I_speeds];


                        apdq_trans[I_q] -= cqxx;                                 
                    }  
                }
            

                /* --------------------------- Y-directions --------------------------- */
                I_speeds = I + (mwaves + mw)*zs;
                s[mw] = speeds[I_speeds];

