Conservative advection equation
   q_t + (u(x)*q)_x = 0
solved using the f-wave approach.  Cell-centered velocities are stored
in the aux array, see setaux.f.   The Riemann solver rp1adfw.f returns
fwaves, and the modified library routine step1fw.f routine must be used
in conjunction with this, as specified in the Makefile.