 
         double *const qr     = start;          /* meqn   */
         double *const ql     = qr + meqn;      /* meqn   */
-        double *const amdq   = ql + meqn;      /* meqn   */
-        double *const apdq   = amdq + meqn;    /* meqn   */
+        double *const apdq   = ql + meqn;    /* meqn   */
         double *const aux1   = apdq + meqn;    /* 2*maux */
         double *const aux2   = aux1 + 2*maux;  /* 2*maux */
         double *const aux3   = aux2 + 2*maux;  /* 2*maux */

                int I_q = I + mq*zs;
                qr[mq] = qold[I_q];        /* Right */
                ql[mq] = qold[I_q - 1];    /* Left  */
                //qd[mq] = qold[I_q - ys];   /* Down  */  
            }

            for(int m = 0; m < maux; m++)

            ql[mq] = qold[I_q-1];
            qr[mq] = qold[I_q];

            //amdq[mq] = amdq_trans[I_q];
            apdq[mq] = apdq_trans[I_q];
        }            

