                         bpdq_trans[I_q] -= cqyy;      
                     } 
                 }   
-                //__syncthreads();                
             }  
         }  
         __syncthreads();

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

        /*  transverse-y
    

                for(mq = 0; mq < meqn; mq++)
                {
                    I_q = I + mq*zs;
                    cqxx = abs(s[mw])*(1.0 - abs(s[mw])*dtdx)*wave[mq];

                    fm[I_q] += 0.5*cqxx;   
                    fp[I_q] += 0.5*cqxx;  
