claw/applications/acoustics/1d/rp

Riemann solvers for one-dimensional acoustics

rp1ac.f   constant coefficients.   Parameters must be set and passed into
          rp1 via the common block 
	         common /cparam/ rho,bulk,cc,zz

rp1acv.f  variable coefficients.  
          aux(i,1) should contain the impedance Z in i'th cell
          aux(i,2) should contain the sound speed c in i'th cell