                 auxd[m] = aux[I_aux - ys];
             }                        
 
+            /* ------------------------ Normal solve in X direction ------------------- */
             rpn2(0, meqn, mwaves, maux, ql, qr, auxl, auxr, wave, s, amdq, apdq);
 
             for (mq = 0; mq < meqn; mq++) 

        }
        return;
    }



    /* ------------------------ Transverse Propagation : X-faces ---------------------- */

                }
            } 

            /* ------------------------ 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 */
