Well‐balanced high‐order solver for blood flow in networks of vessels with variable properties. (31st July 2013)