Auxiliary conditions (6, 6)
(1): h2 = cp * T2 + b * P2
(1): h1 = cp * T1 + b * P1
(1): h3 = cp * T3 + b * P3
(1): P2 = 0.5 * (P2_l + P2_r)
(1): P1 = 0.5 * (P1_l + P1_r)
(1): 0.0 = Q1 + (-Q2) - Q3