Initialized HLLC flux computation

Implements method for characteristic velocity computation of the
intermediate wave, according to P. Batten et. al.