!:RightGoingFluxes:!
            DO j=2-mbc,mx+mbc
              df=zero
              DO mwav=1,mwaves
                df(1:NrVars) = df(1:NrVars) + speed(j,mwav) * wave(j,1:NrVars,mwav)
              END DO
              Apdq(j,1:NrVars) = df(1:NrVars) - Amdq(j,1:NrVars)
            END DO
          !===
          ELSE
          !===