Riemann solvers for 2d advection equations:

Use one of the following pairs of Riemann solvers 
   (only rpn2 is needed if dimensional splitting is used, method(2) < 0)

rpn2ad1.f    Advective form, assumes edge velocities are specified in aux
rpt2ad1.f