\* Source Pyomo model name=Model *\

min 
objective:
+25 flow(line2_C2_b_gas3__0)
+25 flow(line2_C2_b_gas3__1)
+25 flow(line3_C1_b_gas4__0)
+25 flow(line3_C1_b_gas4__1)
+25 flow(line4_C1_b_gas4__0)
+25 flow(line4_C1_b_gas4__1)
+25 flow(source_1_C1_b_gas1__0)
+25 flow(source_1_C1_b_gas1__1)

s.t.

c_e_Bus_balance(C1_b_gas1__0)_:
-1 flow(C1_b_gas1__line1_0)
-1 flow(C1_b_gas1__line2_0)
+1 flow(source_1_C1_b_gas1__0)
= 0

c_e_Bus_balance(C1_b_gas1__1)_:
-1 flow(C1_b_gas1__line1_1)
-1 flow(C1_b_gas1__line2_1)
+1 flow(source_1_C1_b_gas1__1)
= 0

c_e_Bus_balance(C1_b_gas4__0)_:
+1 flow(line3_C1_b_gas4__0)
+1 flow(line4_C1_b_gas4__0)
= 25

c_e_Bus_balance(C1_b_gas4__1)_:
+1 flow(line3_C1_b_gas4__1)
+1 flow(line4_C1_b_gas4__1)
= 25

c_e_Bus_balance(C2_b_gas2__0)_:
-1 flow(C2_b_gas2__line3_0)
+1 flow(line1_C2_b_gas2__0)
= 0

c_e_Bus_balance(C2_b_gas2__1)_:
-1 flow(C2_b_gas2__line3_1)
+1 flow(line1_C2_b_gas2__1)
= 0

c_e_Bus_balance(C2_b_gas3__0)_:
-1 flow(C2_b_gas3__line4_0)
+1 flow(line2_C2_b_gas3__0)
= 0

c_e_Bus_balance(C2_b_gas3__1)_:
-1 flow(C2_b_gas3__line4_1)
+1 flow(line2_C2_b_gas3__1)
= 0

c_e_GasLineBlock_pressure_flow(line1_0)_:
-50 GasLineBlock_pressure(C1_b_gas1__0)
+1 flow(line1_C2_b_gas2__0)
= 0

c_e_GasLineBlock_pressure_flow(line1_1)_:
-50 GasLineBlock_pressure(C1_b_gas1__1)
+1 flow(line1_C2_b_gas2__1)
= 0

c_e_GasLineBlock_pressure_flow(line2_0)_:
-100 GasLineBlock_pressure(C1_b_gas1__0)
+100 GasLineBlock_pressure(C2_b_gas3__0)
+1 flow(line2_C2_b_gas3__0)
= 0

c_e_GasLineBlock_pressure_flow(line2_1)_:
-100 GasLineBlock_pressure(C1_b_gas1__1)
+100 GasLineBlock_pressure(C2_b_gas3__1)
+1 flow(line2_C2_b_gas3__1)
= 0

c_e_GasLineBlock_pressure_flow(line3_0)_:
+100 GasLineBlock_pressure(C1_b_gas4__0)
+1 flow(line3_C1_b_gas4__0)
= 0

c_e_GasLineBlock_pressure_flow(line3_1)_:
+100 GasLineBlock_pressure(C1_b_gas4__1)
+1 flow(line3_C1_b_gas4__1)
= 0

c_e_GasLineBlock_pressure_flow(line4_0)_:
+100 GasLineBlock_pressure(C1_b_gas4__0)
-100 GasLineBlock_pressure(C2_b_gas3__0)
+1 flow(line4_C1_b_gas4__0)
= 0

c_e_GasLineBlock_pressure_flow(line4_1)_:
+100 GasLineBlock_pressure(C1_b_gas4__1)
-100 GasLineBlock_pressure(C2_b_gas3__1)
+1 flow(line4_C1_b_gas4__1)
= 0

c_e_GasLineBlock__equate_pressure_flows(line1_0)_:
-1 flow(C1_b_gas1__line1_0)
+1 flow(line1_C2_b_gas2__0)
= 0

c_e_GasLineBlock__equate_pressure_flows(line1_1)_:
-1 flow(C1_b_gas1__line1_1)
+1 flow(line1_C2_b_gas2__1)
= 0

c_e_GasLineBlock__equate_pressure_flows(line2_0)_:
-1 flow(C1_b_gas1__line2_0)
+1 flow(line2_C2_b_gas3__0)
= 0

c_e_GasLineBlock__equate_pressure_flows(line2_1)_:
-1 flow(C1_b_gas1__line2_1)
+1 flow(line2_C2_b_gas3__1)
= 0

c_e_GasLineBlock__equate_pressure_flows(line3_0)_:
-1 flow(C2_b_gas2__line3_0)
+1 flow(line3_C1_b_gas4__0)
= 0

c_e_GasLineBlock__equate_pressure_flows(line3_1)_:
-1 flow(C2_b_gas2__line3_1)
+1 flow(line3_C1_b_gas4__1)
= 0

c_e_GasLineBlock__equate_pressure_flows(line4_0)_:
-1 flow(C2_b_gas3__line4_0)
+1 flow(line4_C1_b_gas4__0)
= 0

c_e_GasLineBlock__equate_pressure_flows(line4_1)_:
-1 flow(C2_b_gas3__line4_1)
+1 flow(line4_C1_b_gas4__1)
= 0

c_e_ONE_VAR_CONSTANT: 
ONE_VAR_CONSTANT = 1.0

bounds
    -inf <= flow(C1_b_gas1__line1_0) <= +inf
    -inf <= flow(C1_b_gas1__line1_1) <= +inf
    -inf <= flow(C1_b_gas1__line2_0) <= +inf
    -inf <= flow(C1_b_gas1__line2_1) <= +inf
    -inf <= flow(C2_b_gas2__line3_0) <= +inf
    -inf <= flow(C2_b_gas2__line3_1) <= +inf
    -inf <= flow(C2_b_gas3__line4_0) <= +inf
    -inf <= flow(C2_b_gas3__line4_1) <= +inf
   0 <= flow(line1_C2_b_gas2__0) <= 100
   0 <= flow(line1_C2_b_gas2__1) <= 100
   0 <= flow(line2_C2_b_gas3__0) <= 100
   0 <= flow(line2_C2_b_gas3__1) <= 100
   0 <= flow(line3_C1_b_gas4__0) <= 100
   0 <= flow(line3_C1_b_gas4__1) <= 100
   0 <= flow(line4_C1_b_gas4__0) <= 100
   0 <= flow(line4_C1_b_gas4__1) <= 100
   0 <= flow(source_1_C1_b_gas1__0) <= 100
   0 <= flow(source_1_C1_b_gas1__1) <= 100
   0 <= GasLineBlock_pressure(C1_b_gas1__0) <= 1
   0 <= GasLineBlock_pressure(C1_b_gas1__1) <= 1
   0 <= GasLineBlock_pressure(C1_b_gas4__0) <= 1
   0 <= GasLineBlock_pressure(C1_b_gas4__1) <= 1
   0 <= GasLineBlock_pressure(C2_b_gas3__0) <= 1
   0 <= GasLineBlock_pressure(C2_b_gas3__1) <= 1
end