     double ur = hur/hr;
     double vr = hvr/hr;
 
+#if 0
     double eta = hr  + aux2[0];
     double topo1 = aux1[0];
     double topo3 = aux3[0];
+#endif    
 
     /* Construct Roe variables */
     double hrs = sqrt(hr);

        szp[mw] = (1+z)/2*beta[mw];
    }

    for(mq = 0; mq < meqn; mq++)
    {
        bmasdq[mq]  = szm[0]*waveb[mq];
        bmasdq[mq] += szm[1]*waveb[meqn + mq];

    beta[2] = (delf3/(s3-s1))-(s1*delf1/(s3-s1));

    double wave[9];
    wave[0]  = 1;
    wave[mu] = s2;
    wave[mv] = s1;

    wave[meqn + 0]  = 0;
    wave[meqn + mu] = 1;
    wave[meqn + mv] = 0;

    wave[2*meqn + 0]  = 1;
    wave[2*meqn + mu] = s2;
    wave[2*meqn + mv] = s3;


    double szm[3], szp[3];
