                         wave[mq] = waves[I_waves];
                     }                        
 
-                    if (mthlim[mw] > 0)
+                    // if (mthlim[mw] > 0)
                     {
                         wnorm2 = dotl = dotr = 0;
                         for(mq = 0; mq < meqn; mq++)

                gm[I_q - 1] -= gupdate;       
                gp[I_q - 1] -= gupdate;
            }
        } /* thread conditional */
    } /* Thread loop */

    __syncthreads();


                    waves[I_waves] = wave[mw*meqn + mq];
                }
            }

            /* ------------------------ Normal solve in Y direction ------------------- */
            rpn2(1, meqn, mwaves, maux, qd, qr, auxd, auxr, wave, s, bmdq, bpdq);
