!:DecomposeJump:! in q onto eigenbases
          DO j=2-mbc,mx+mbc
            ! Find coef(1) thru coef(4), the coefficients of the 4 eigenvectors
            delta(1) = q1D(j,1) - q1D(j-1,1)
            delta(2) = q1D(j,mu) - q1D(j-1,mu)
            delta(3) = q1D(j,mv) - q1D(j-1,mv)
            delta(4) = q1D(j,4) - q1D(j-1,4)
            coef(3) = g1c2(j) * (euv(j)*delta(1) + uR(j)*delta(2) + vR(j)*delta(3) - delta(4))
            coef(2) = delta(3) - vR(j)*delta(1)
            coef(4) = (delta(2) + (cR(j)-uR(j))*delta(1) - cR(j)*coef(3)) / (two*cR(j))
            coef(1) = delta(1) - coef(3) - coef(4)
            IF (NrVars==5) THEN
              ! Tracer variable
              coef(5) = q1D(j,5) - q1D(j-1,5)
            END IF
            DO mw=1,mwaves
              wave(j,1:NrVars,mw) = coef(mw)*wave(j,1:NrVars,mw)
            END DO
          END DO