Compute gamma1 instead of loading it from a common block.

Fixes #79; see also https://github.com/clawpack/pyclaw/issues/447.

This will break tests since PyClaw examples need to be updated too.