         szp[mw] = (1+z)/2*beta[mw];
     }
 
-    for(mq = 0; mq < meqn; mq++)
+    for(int mq = 0; mq < meqn; mq++)
     {
         bmasdq[mq]  = szm[0]*waveb[mq];
         bmasdq[mq] += szm[1]*waveb[meqn + mq];

        bpasdq[mq] += max(sb[1],0.)*waveb[meqn + mq];
        bpasdq[mq] += max(sb[2],0.)*waveb[2*meqn + mq];
    }
}
#endif


__device__ cudaclaw_cuda_rpt2_t geoclaw_rpt2 = geoclaw_rpt2shallow;

void geoclaw_assign_rpt2(cudaclaw_cuda_rpt2_t *rpt2)
{

    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];
