                 }
             } 
 
+            /* ------------------------ Normal solve in Y direction ------------------- */
             rpn2(1, meqn, mwaves, maux, qd, qr, auxd, auxr, wave, s, bmdq, bpdq);
 
             /* Set value at bottom interface of cell I */

                I_q = I + mq*zs;
                fp[I_q] = -apdq[mq]; 
                fm[I_q] = amdq[mq];
                if (order[0] > 0)
                {
                    apdq_trans[I_q] = apdq[mq];  
                    amdq_trans[I_q] = amdq[mq];                    
                }
            }

            if (order[0] == 2)
            {
                for (m = 0; m < meqn*mwaves; m++)
                {
                    I_waves = I + m*zs;
                    waves[I_waves] = wave[m];
                }                
            }

            for (mw = 0; mw < mwaves; mw++)
            {
                I_speeds = I + mw*zs;
                speeds[I_speeds] = s[mw];
                cfl = abs(s[mw]*dtdx);
                if (cfl > maxcfl)
                {

                }  /* End of mwaves loop */
            } /* End of thread conditional */
        } /* End of thread loop */
        __syncthreads();
    } /* End of check on order[0] == 2 */



    if (order[1] == 0)
