tfluct as input to solver, following rp

This PR

- adds the ability to provide a _tfluct solver_ to _SharpClaw_ directly in the same fashion as Riemann solvers are passed. When `solver.tfluct_solver=True`, a fortran python-wrapped routine can be given through `import _name-of-tfluct-solver_`, `solver.tfluct=_name-of-tfluct-solver_` 
- solver follows API convention of Riemann solvers
- works in n-d
- euler_1d has been modified to provide an example of usage

The correct name for this variable here is adq, not amdq2.  The function is supposed to compute the total fluctuation.