Running: ./testmodel.py --libraries=/home/hudson/saved_omc/libraries/.openmodelica/libraries --ompython_omhome=/usr ThermalSeparation_ThermalSeparation.Examples.ProPen_EQ_SteadyState.conf.json
loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/ModelicaServices 4.1.0+maint.om/package.mo", uses=false) [Timeout 180]
"Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/ModelicaServices 4.1.0+maint.om/package.mo): time 0.001351/0.001351, allocations: 88.33 kB / 20.99 MB, free: 292 kB / 14.72 MB
"
[Timeout remaining time 180]
loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/Complex 4.1.0+maint.om/package.mo", uses=false) [Timeout 180]
"Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/Complex 4.1.0+maint.om/package.mo): time 0.001418/0.001418, allocations: 163 kB / 24.29 MB, free: 3.25 MB / 14.72 MB
"
[Timeout remaining time 180]
loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/Modelica 4.1.0+maint.om/package.mo", uses=false) [Timeout 180]
"Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/Modelica 4.1.0+maint.om/package.mo): time 1.431/1.431, allocations: 177.1 MB / 204.6 MB, free: 9.043 MB / 190.1 MB
"
[Timeout remaining time 178]
loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/package.mo", uses=false) [Timeout 180]
"Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/package.mo): time 0.4827/0.4827, allocations: 50.75 MB / 311.8 MB, free: 2.242 MB / 286.1 MB
[/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/Components/Compressors/Compressor.mo:146:3-158:16:writable] Warning: External function annotation should occur on the external-clause, not on the function.
[/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/Components/Compressors/Compressor.mo:160:3-168:15:writable] Warning: External function annotation should occur on the external-clause, not on the function.
[/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/Components/Compressors/OldCompressorModels.mo:145:4-157:17:writable] Warning: External function annotation should occur on the external-clause, not on the function.
[/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/Components/Compressors/OldCompressorModels.mo:159:5-167:17:writable] Warning: External function annotation should occur on the external-clause, not on the function.
[/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/Components/Compressors/OldCompressorModels.mo:320:4-332:17:writable] Warning: External function annotation should occur on the external-clause, not on the function.
[/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/Components/Compressors/OldCompressorModels.mo:334:5-342:17:writable] Warning: External function annotation should occur on the external-clause, not on the function.
"
[Timeout remaining time 179]
Using package ThermalSeparation with version 0.2 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/ThermalSeparation master/package.mo)
Using package Modelica with version 4.1.0 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/Modelica 4.1.0+maint.om/package.mo)
Using package Complex with version 4.1.0 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/Complex 4.1.0+maint.om/package.mo)
Using package ModelicaServices with version 4.1.0 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/ModelicaServices 4.1.0+maint.om/package.mo)
Running command: translateModel(ThermalSeparation.Examples.ProPen_EQ_SteadyState,tolerance=1e-06,outputFormat="empty",numberOfIntervals=2500,variableFilter="",fileNamePrefix="ThermalSeparation_ThermalSeparation.Examples.ProPen_EQ_SteadyState")
translateModel(ThermalSeparation.Examples.ProPen_EQ_SteadyState,tolerance=1e-06,outputFormat="empty",numberOfIntervals=2500,variableFilter="",fileNamePrefix="ThermalSeparation_ThermalSeparation.Examples.ProPen_EQ_SteadyState") [Timeout 660]
"Warning: Ignoring unknown experiment annotation option: Algorithm = \"Dassl\"
Notification: Performance of FrontEnd - Absyn->SCode: time 1.788e-05/1.788e-05, allocations: 2.312 kB / 395.8 MB, free: 5.082 MB / 366.1 MB
Notification: Performance of NFInst.instantiate(ThermalSeparation.Examples.ProPen_EQ_SteadyState): time 0.0236/0.02362, allocations: 19.98 MB / 415.8 MB, free: 1.023 MB / 382.1 MB
Notification: Performance of NFInst.instExpressions: time 0.02286/0.04649, allocations: 15.42 MB / 431.2 MB, free: 1.555 MB / 398.1 MB
Notification: Performance of NFInst.updateImplicitVariability: time 0.003098/0.04958, allocations: 75.5 kB / 431.3 MB, free: 1.48 MB / 398.1 MB
Notification: Performance of NFTyping.typeComponents: time 0.006582/0.05617, allocations: 2.563 MB / 433.8 MB, free: 14.9 MB / 414.1 MB
Notification: Performance of NFTyping.typeBindings: time 0.00792/0.06409, allocations: 3.359 MB / 437.2 MB, free: 11.52 MB / 414.1 MB
Notification: Performance of NFTyping.typeClassSections: time 0.003724/0.06781, allocations: 1.54 MB / 438.7 MB, free: 9.98 MB / 414.1 MB
Notification: Performance of NFFlatten.flatten: time 0.01956/0.08737, allocations: 14.67 MB / 453.4 MB, free: 11.28 MB / 430.1 MB
Notification: Performance of NFFlatten.resolveConnections: time 0.004845/0.09222, allocations: 2.477 MB / 455.9 MB, free: 8.766 MB / 430.1 MB
Notification: Performance of NFEvalConstants.evaluate: time 0.01254/0.1048, allocations: 6.163 MB / 462 MB, free: 2.59 MB / 430.1 MB
Notification: Performance of NFSimplifyModel.simplify: time 0.00873/0.1135, allocations: 4.915 MB / 466.9 MB, free: 13.66 MB / 446.1 MB
Notification: Performance of NFPackage.collectConstants: time 0.001686/0.1152, allocations: 412 kB / 467.3 MB, free: 13.26 MB / 446.1 MB
Notification: Performance of NFFlatten.collectFunctions: time 0.00636/0.1215, allocations: 1.769 MB / 469.1 MB, free: 11.49 MB / 446.1 MB
Notification: Performance of combineBinaries: time 0.1961/0.3176, allocations: 11.68 MB / 480.8 MB, free: 62.89 MB / 446.1 MB
Notification: Performance of replaceArrayConstructors: time 0.007215/0.3248, allocations: 7.672 MB / 488.5 MB, free: 62.46 MB / 446.1 MB
Notification: Performance of NFVerifyModel.verify: time 0.002495/0.3273, allocations: 385.1 kB / 488.8 MB, free: 62.46 MB / 446.1 MB
Notification: Performance of FrontEnd: time 0.0006662/0.328, allocations: 156.8 kB / 489 MB, free: 62.46 MB / 446.1 MB
Notification: Model statistics after passing the front-end and creating the data structures used by the back-end:
 * Number of equations: 1366 (897)
 * Number of variables: 1366 (1065)
Notification: Performance of [SIM] Bindings: time 0.03784/0.3658, allocations: 40.03 MB / 0.5166 GB, free: 49 MB / 446.1 MB
Notification: Performance of [SIM] FunctionAlias: time 0.004676/0.3705, allocations: 3.878 MB / 0.5204 GB, free: 47.12 MB / 446.1 MB
Notification: Performance of [SIM] Early Inline: time 0.01616/0.3866, allocations: 17.56 MB / 0.5376 GB, free: 37.23 MB / 446.1 MB
Notification: Performance of [SIM] Simplify 1: time 0.003935/0.3906, allocations: 1.763 MB / 0.5393 GB, free: 36.64 MB / 446.1 MB
Warning: NBAlias.setStartFixed: Alias set with conflicting unfixed start values detected. Use -d=dumprepl for more information.
Notification: Performance of [SIM] Alias: time 0.01993/0.4105, allocations: 16.22 MB / 0.5551 GB, free: 27.96 MB / 446.1 MB
Notification: Performance of [SIM] Simplify 2: time 0.002362/0.4129, allocations: 1.541 MB / 0.5566 GB, free: 27.38 MB / 446.1 MB
Notification: Performance of [SIM] Remove Stream: time 0.001474/0.4143, allocations: 1.394 MB / 0.558 GB, free: 26.6 MB / 446.1 MB
Notification: Performance of [SIM] Detect States: time 0.002742/0.4171, allocations: 3.422 MB / 0.5613 GB, free: 24.71 MB / 446.1 MB
Notification: Performance of [SIM] Events: time 0.001335/0.4184, allocations: 1.634 MB / 0.5629 GB, free: 23.97 MB / 446.1 MB
Notification: Performance of [SIM] Partitioning: time 0.006248/0.4247, allocations: 6.412 MB / 0.5692 GB, free: 19.79 MB / 446.1 MB
Error: Internal error NBSorting.tarjan failed to sort system:
System Variables (424/558)
****************************
(1|1)          [ALGB] (1) Real column1.propsVapIn.MM = column1.mediumVapourIn.properties.MM (min = 0.001, max = 0.25, nominal = 0.032)
(2|2)          [ALGB] (1) protected Real[1] column1.pressureLoss.rho_l_calc = if column1.pressureLoss.rho_const then {column1.pressureLoss.rho_l_nom for $f1 in 1:1} else column1.pressureLoss.rho_l (min = {0.0 for $f2 in 1:1})
(3|3)          [ALGB] (1) protected Real[1] column1.holdup.eta = column1.propsLiq.eta (min = {0.0 for $f1 in 1:1})
(4|4)          [ALGB] (1) Real sourceLiquid.medium.u (nominal = 2e4)
(5|5)          [ALGB] (1) Real[1] column1.propsLiq.lambda = {column1.mediumLiquid[1].properties.lambda}
(6|6)          [ALGB] (1) Real column1.propsVapIn.cp = column1.mediumVapourIn.properties.cp
(7|7)          [ALGB] (1) Real sourceLiquid.medium.state.h (min = -1e8, max = 1e8, nominal = 1e6)
(8|8)          [DISC] (1) Boolean[1] column1.bool_eps
(9|9)          [ALGB] (1) Real column1.V_liq = sum(column1.eps_liq * column1.H * column1.A * column1.eps)
(10|10)        [ALGB] (1) Real sourceLiquid.medium.MM (min = 0.001, max = 0.25, nominal = 0.032)
(11|11)        [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.x_l = {column1.balanceEquations.x_l_star[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(12|13)        [ALGB] (1) Real column1.Ndot_v_in
(13|14)        [ALGB] (2) Real[2] column1.balanceEquations.homotopyMethod.K = {1.0 for $f1 in 1:2}
(14|16)        [ALGB] (2) Real[1, 2] column1.bubblePressure.p_partial
(15|18)        [ALGB] (2) Real[1, 2] column1.x_l_star (start = column1.x_l_start, min = max({0.0 for $f3 in 1:2, $f4 in 1:1}, {0.0 for $f3 in 1:2, $f4 in 1:1}), max = max({1.0 for $f1 in 1:2, $f2 in 1:1}, {1.0 for $f1 in 1:2, $f2 in 1:1}))
(16|20)        [ALGB] (1) Real[1] column1.h_v = if column1.homotopyMethod.bool_h and false then {homotopy(column1.mediumVapour[1].h, ({15000.0 for $f1 in 1:1})[1])} else column1.mediumVapour.h
(17|21)        [ALGB] (1) protected Real[1] column1.pressureLoss.rho_l = column1.propsLiq.rho (min = {0.0 for $f1 in 1:1})
(18|22)        [ALGB] (2) Real[1, 2] column1.bubblePressure.x_vap_liq = {column1.x_vap_liq[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(19|24)        [ALGB] (2) Real[2] sourceLiquid.c
(20|26)        [ALGB] (2) Real[1, 2] column1.mediumVapour.x_star = {column1.x_v_star[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(21|28)        [ALGB] (1) Real[1] column1.h_l_feed
(22|29)        [ALGB] (1) Real[1] column1.sum_xv = sum(column1.x_v_star[:, i] for i in 1:2)
(23|30)        [ALGB] (1) Real[1] column1.mediumVapourLink.mediumConIn.MM (min = {0.0 for $f1 in 1:1})
(24|31)        [ALGB] (2) Real[2] column1.x_l_in (min = max({0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}), max = max({1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}))
(25|33)        [ALGB] (1) protected Real[1] column1.mediumVapour.R
(26|34)        [ALGB] (1) flow Real sinkGas.gasPortIn.Ndot
(27|35)        [ALGB] (1) input Real[1] column1.reaction.propsLiq.u = {column1.mediumLiquid[1].properties.u} (nominal = {2e4 for $f54 in 1:1})
(28|36)        [ALGB] (1) Real[1] column1.n_mol_V
(29|37)        [ALGB] (1) Real[1] column1.balanceEquations.thermoEquilibrium.T = {column1.balanceEquations.T_star[1]} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(30|38)        [ALGB] (1) Real column1.propsLiqIn.MM = column1.mediumLiquidIn.properties.MM (min = 0.001, max = 0.25, nominal = 0.032)
(31|39)        [ALGB] (2) stream Real[1, 2] column1.internalFeedPort.feedVapourInternal.x_outflow (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(32|41)        [DISC] (1) Boolean column1.balanceEquations.homotopyMethod.bool_rho = column1.balanceEquations.homotopyMethod.whichHomotopy[5]
(33|42)        [ALGB] (1) Real[1] column1.propsVap.rho = {column1.mediumVapour[1].properties.rho} (start = {1.0 for $f12 in 1:1}, min = {0.0 for $f14 in 1:1}, max = {1e5 for $f13 in 1:1}, nominal = {1.0 for $f11 in 1:1})
(34|43)        [ALGB] (1) Real[1] column1.Vdot_l (start = {0.03 for $f2 in 1:1}, nominal = {1e-4 for $f1 in 1:1})
(35|44)        [ALGB] (1) Real column1.propsLiqIn.cp = column1.mediumLiquidIn.properties.cp
(36|45)        [DISC] (1) Boolean column1.balanceEquations.homotopyMethod.bool_Ndot_inter = column1.balanceEquations.homotopyMethod.whichHomotopy[1]
(37|46)        [ALGB] (1) Real[1] column1.bubblePressure.T = {column1.T_l[1]} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(38|47)        [ALGB] (1) Real[1] column1.bubblePressure.henryCoefficient.T = {column1.bubblePressure[1].T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(39|48)        [ALGB] (1) Real sourceGas_Vdot.medium.p_start (start = 1.6e5)
(40|49)        [ALGB] (2) Real[2] column1.propsLiqIn.x = column1.mediumLiquidIn.properties.x (min = {0.0 for $f5 in 1:2}, max = {1.0 for $f4 in 1:2})
(41|51)        [ALGB] (2) Real[2] column1.results.p_v = column1.p_hyd
(42|53)        [ALGB] (1) Real column1.propsVapIn.d = column1.mediumVapourIn.properties.d (min = 0.0)
(43|54)        [ALGB] (1) Real sourceLiquid.medium.properties.MM (min = 0.001, max = 0.25, nominal = 0.032)
(44|55)        [ALGB] (1) Real[1] column1.balanceEquations.thermoEquilibrium.v_v = column1.balanceEquations.MM_v / column1.balanceEquations.rho_v (min = {0.0 for $f1 in 1:1})
(45|56)        [ALGB] (2) output Real[2] sourceGas_Vdot.medium.properties.c
(46|58)        [ALGB] (2) Real[2] column1.balanceEquations.c_v_in = column1.propsVapIn.c
(47|60)        [ALGB] (1) Real[1] column1.mediumLiquid.mediumConOut.rho = {column1.mediumLiquid[1].d} (min = {0.0 for $f1 in 1:1})
(48|61)        [ALGB] (2) output Real[2] sourceGas_Vdot.medium.state.X (min = {0.0 for $f3 in 1:2}, max = {1.0 for $f2 in 1:2}, nominal = {0.1 for $f1 in 1:2})
(49|63)        [ALGB] (1) Real[1] column1.mediumVapour.T = {column1.T_v[1]} (start = {300.0 for $f2 in 1:1}, min = {1.0 for $f4 in 1:1}, max = {1e4 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(50|64)        [ALGB] (2) Real[1, 2] column1.n_mol_V_i
(51|66)        [ALGB] (2) Real[1, 2] column1.results.x_l_star = column1.x_l_star (min = {0.0 for $f25 in 1:2, $f26 in 1:1}, max = {1.0 for $f23 in 1:2, $f24 in 1:1})
(52|68)        [ALGB] (2) Real[2] column1.x_downStreamOut_act (min = max({0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}), max = max({1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}))
(53|70)        [ALGB] (1) protected Real[1] column1.pressureLoss.rho_v_calc = if column1.pressureLoss.rho_const then {column1.pressureLoss.rho_v_nom for $f1 in 1:1} else column1.pressureLoss.rho_v (min = {0.0 for $f2 in 1:1})
(54|71)        [ALGB] (2) protected Real[2] sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat = {231.1, 309.2} (start = {288.15 for $f2 in 1:2}, min = {0.0 for $f3 in 1:2}, nominal = {300.0 for $f1 in 1:2})
(55|73)        [ALGB] (1) Real[1] column1.propsVap.h = {column1.mediumVapour[1].properties.h}
(56|74)        [ALGB] (1) protected Real[1] column1.pressureLoss.sigma_l = column1.propsLiq.sigma
(57|75)        [ALGB] (1) Real[1] column1.mediumLiquid.eta (min = {0.0 for $f1 in 1:1})
(58|76)        [ALGB] (2) Real[1, 2] column1.c_v_feed_used
(59|78)        [ALGB] (1) Real column1.rho_l_in = column1.rho_l_in (min = 0.0, max = 1e5)
(60|79)        [ALGB] (2) Real[1, 2] column1.n_i_liq
(61|81)        [ALGB] (1) Real[1] column1.u_v = column1.u_v (nominal = {2e4 for $f2 in 1:1}, StateSelect = default)
(62|82)        [ALGB] (1) Real sourceLiquid.medium.properties.v (start = 1e-5, min = 0.0, StateSelect = never)
(63|83)        [ALGB] (1) Real[1] column1.mdot_v = column1.Vdot_v .* column1.rho_v
(64|84)        [ALGB] (1) Real sourceLiquid.medium.properties.rho (start = 1.0, min = 0.0, max = 1e5, nominal = 1.0)
(65|85)        [ALGB] (2) Real[1, 2] column1.balanceEquations.c_l_star (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(66|87)        [ALGB] (1) Real[1] column1.n_total = column1.n_liq + column1.n_vap
(67|88)        [ALGB] (1) Real sourceLiquid.p (min = 0.0, max = 1e8)
(68|89)        [ALGB] (1) Real[1] column1.T (start = {288.15 for $f2 in 1:1}, min = max({0.0}, {0.0}, {0.0 for $f3 in 1:1}, {0.0}), nominal = {300.0 for $f1 in 1:1})
(69|90)        [ALGB] (1) Real sourceGas_Vdot.rho (min = 0.0, max = 1e5)
(70|91)        [ALGB] (1) Real column1.propsLiqIn.lambda = column1.mediumLiquidIn.properties.lambda
(71|92)        [ALGB] (2) Real[2] column1.propsVapIn.x = column1.mediumVapourIn.properties.x (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(72|94)        [ALGB] (1) protected Real[1] column1.initOption.v = column1.propsLiq.v
(73|95)        [ALGB] (2) Real[1, 2] column1.mediumVapour.calcSpecificEnthalpy.x = {column1.mediumVapour[1].x} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(74|97)        [ALGB] (2) Real[1, 2] column1.X_l
(75|99)        [ALGB] (1) input Real[1] column1.balanceEquations.stateLiq.h = {column1.mediumLiquid[1].state.h} (min = {-1e8 for $f3 in 1:1}, max = {1e8 for $f2 in 1:1}, nominal = {1e6 for $f1 in 1:1})
(76|100)       [ALGB] (1) output Real sourceGas_Vdot.medium.state.p (start = 1e5, min = 0.0, max = 1e8, nominal = 1e5)
(77|101)       [ALGB] (1) Real[1] column1.propsLiq.MM = {column1.mediumLiquid[1].properties.MM} (min = {0.001 for $f28 in 1:1}, max = {0.25 for $f27 in 1:1}, nominal = {0.032 for $f26 in 1:1})
(78|102)       [ALGB] (1) Real column1.propsLiqIn.d = column1.mediumLiquidIn.properties.d (min = 0.0)
(79|103)       [ALGB] (1) Real[1] column1.T_star (start = column1.T_l_start, min = max({0.0 for $f2 in 1:1}, {0.0}), nominal = {300.0 for $f1 in 1:1})
(80|104)       [ALGB] (1) Real[1] column1.MM_l_feed (min = max({0.0}, {0.0 for $f1 in 1:1}))
(81|105)       [ALGB] (1) Real sourceLiquid.medium.state.p (start = 1e5, min = 0.0, max = 1e8, nominal = 1e5)
(82|106)       [ALGB] (1) Real[1] column1.propsVap.T = {column1.mediumVapour[1].properties.T} (start = {300.0 for $f2 in 1:1}, min = {1.0 for $f4 in 1:1}, max = {1e4 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(83|107)       [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.henryCoefficient.henryCoeff.x_l = {column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.x_l} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(84|109)       [ALGB] (1) Real[1] column1.mediumVapour.calcSpecificEnthalpy.p = {column1.mediumVapour[1].p}
(85|110)       [ALGB] (2) Real[2] column1.x_upStreamIn_act (min = max({0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}), max = max({1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}))
(86|112)       [ALGB] (1) Real sourceLiquid.medium.properties.p
(87|113)       [ALGB] (2) Real[1, 2] column1.results.x_v = column1.x_v_star (min = {0.0 for $f17 in 1:2, $f18 in 1:1}, max = {1.0 for $f15 in 1:2, $f16 in 1:1})
(88|115)       [ALGB] (1) Real[1] column1.propsLiq.eta = {column1.mediumLiquid[1].properties.eta} (min = {0.0 for $f10 in 1:1})
(89|116)       [ALGB] (1) Real[1] column1.balanceEquations.rho_l = column1.propsLiq.d (min = {0.0 for $f1 in 1:1})
(90|117)       [ALGB] (1) Real column1.mediumLiquidIn.state.p (start = 1e5, min = 0.0, max = 1e8, nominal = 1e5)
(91|118)       [ALGB] (1) Real column1.mediumVapourIn.cp
(92|119)       [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.x_vap_liq = {column1.balanceEquations.x_vap_liq[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(93|121)       [ALGB] (1) Real[1] column1.h_l
(94|122)       [ALGB] (1) Real column1.propsVapIn.lambda = column1.mediumVapourIn.properties.lambda
(95|123)       [ALGB] (1) Real[1] column1.p_bub = column1.p_bub
(96|124)       [ALGB] (1) Real[1] column1.sum_xl = sum(column1.x_l[:, i] for i in 1:2)
(97|125)       [ALGB] (1) Real sourceLiquid.rho (min = 0.0, max = 1e5)
(98|126)       [ALGB] (2) Real[2] column1.x_upStreamOut_act (min = max({0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}), max = max({1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}))
(99|128)       [ALGB] (2) protected Real[1, 2] column1.mediumVapour.calcSpecificEnthalpy.T_sat = {{231.1, 309.2} for $f1 in 1:1} (start = {288.15 for $f4 in 1:2, $f5 in 1:1}, min = {0.0 for $f6 in 1:2, $f7 in 1:1}, nominal = {300.0 for $f2 in 1:2, $f3 in 1:1})
(100|130)      [ALGB] (1) Real column1.propsLiqIn.p = column1.mediumLiquidIn.properties.p
(101|131)      [DISC] (1) Boolean $TEV_1
(102|132)      [ALGB] (2) Real[2] column1.propsLiqIn.eta_comp = column1.mediumLiquidIn.properties.eta_comp (min = {0.0 for $f1 in 1:2})
(103|134)      [DISC] (1) Boolean $TEV_0
(104|135)      [ALGB] (2) Real[1, 2] column1.mediumVapour.X (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(105|137)      [ALGB] (2) Real[2] column1.c_v_in
(106|139)      [ALGB] (1) Real[1] column1.n_mol_L
(107|140)      [ALGB] (1) Real column1.balanceEquations.h_v_in = column1.propsVapIn.h
(108|141)      [ALGB] (1) Real[1] column1.mediumVapour.p_start (start = {1.6e5 for $f1 in 1:1})
(109|142)      [ALGB] (1) Real[1] column1.balanceEquations.thermoEquilibrium.henryCoefficient.henryCoeff.T = {column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(110|143)      [ALGB] (1) Real[1] column1.Vdot_v (start = {0.01 for $f2 in 1:1}, nominal = {0.01 for $f1 in 1:1})
(111|144)      [ALGB] (1) Real[1] column1.propsLiq.rho = {column1.mediumLiquid[1].properties.rho} (start = {1.0 for $f20 in 1:1}, min = {0.0 for $f22 in 1:1}, max = {1e5 for $f21 in 1:1}, nominal = {1.0 for $f19 in 1:1})
(112|145)      [ALGB] (1) Real[1] column1.mediumLiquid.p = {column1.p_hyd[1]} (start = {1e5 for $f2 in 1:1}, min = {0.0 for $f4 in 1:1}, max = {1e8 for $f3 in 1:1}, nominal = {1e5 for $f1 in 1:1})
(113|146)      [ALGB] (1) Real[1] column1.balanceEquations.rho_v = column1.propsVap.d (min = {0.0 for $f1 in 1:1})
(114|147)      [ALGB] (2) Real[1, 2] column1.x_l (start = column1.x_l_start, min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(115|149)      [ALGB] (1) Real[1] column1.h_v_feed
(116|150)      [ALGB] (1) Real[1] column1.Vdot_l_feed
(117|151)      [ALGB] (1) Real[1] column1.mediumLiquid.R
(118|152)      [ALGB] (1) Real[1] column1.mediumVapour.lambda
(119|153)      [DISC] (1) Boolean column1.balanceEquations.homotopyMethod.bool_h = column1.balanceEquations.homotopyMethod.whichHomotopy[6]
(120|154)      [ALGB] (1) Real[1] column1.mediumLiquid.calcSpecificEnthalpy.T = {column1.mediumLiquid[1].T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(121|155)      [ALGB] (1) input Real[1] column1.balanceEquations.stateLiq.p = {column1.mediumLiquid[1].state.p} (start = {1e5 for $f21 in 1:1}, min = {0.0 for $f23 in 1:1}, max = {1e8 for $f22 in 1:1}, nominal = {1e5 for $f20 in 1:1})
(122|156)      [ALGB] (1) Real[1] column1.Ndot_l
(123|157)      [ALGB] (2) protected Real[2] column1.mediumLiquidIn.A = {-6.71791, -7.27889}
(124|159)      [ALGB] (1) Real[1] column1.evapEnthalpy.p = {column1.p_hyd[1]}
(125|160)      [ALGB] (2) Real[2] column1.balanceEquations.x_l_in = column1.propsLiqIn.x (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(126|162)      [ALGB] (1) Real column1.mediumLiquidIn.v (start = 1e-5, min = 0.0, StateSelect = never)
(127|163)      [ALGB] (2) protected Real[2] column1.mediumLiquidIn.D = {-1.2499, -2.22146}
(128|165)      [ALGB] (1) Real[1] column1.mediumLiquid.cp
(129|166)      [ALGB] (1) stream Real column1.upStreamIn.h_outflow
(130|167)      [ALGB] (1) Real column1.h_v_in = column1.h_v_in
(131|168)      [ALGB] (1) Real[1] column1.mediumVapour.d (start = {1.0 for $f2 in 1:1}, min = {0.0 for $f4 in 1:1}, max = {1e5 for $f3 in 1:1}, nominal = {1.0 for $f1 in 1:1})
(132|169)      [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.p_sat = {column1.balanceEquations.p_sat[1]}
(133|171)      [ALGB] (2) Real[1, 2] column1.mediumVapour.c = {column1.c_v[1]}
(134|173)      [ALGB] (1) Real[1] column1.u_l = column1.u_l (nominal = {2e4 for $f2 in 1:1}, StateSelect = default)
(135|174)      [ALGB] (1) Real[1] column1.results.Vdot_l = column1.Vdot_l
(136|175)      [ALGB] (2) Real[1, 2] column1.n_i_vap
(137|177)      [ALGB] (1) Real[1] column1.propsLiq.sigma = {column1.mediumLiquid[1].properties.sigma}
(138|178)      [ALGB] (1) output Real sourceGas_Vdot.medium.properties.MM (min = 0.001, max = 0.25, nominal = 0.032)
(139|179)      [ALGB] (1) output Real column1.mediumVapourIn.state.T (start = 300.0, min = 1.0, max = 1e4, nominal = 300.0)
(140|180)      [ALGB] (2) Real[1, 2] column1.bubblePressure.p_sat = {column1.p_sat_bulk[1]}
(141|182)      [ALGB] (1) Real column1.mediumVapourIn.p_start (start = 1.6e5)
(142|183)      [ALGB] (1) protected Real column1.mediumVapourIn.R
(143|184)      [ALGB] (1) Real[1] column1.balanceEquations.h_v = column1.propsVap.h
(144|185)      [ALGB] (1) Real[1] column1.bubblePressure.v_v = column1.MM_v / column1.rho_v (min = {0.0 for $f1 in 1:1})
(145|186)      [ALGB] (1) Real[1] column1.results.eps_liq = column1.eps_liq
(146|187)      [ALGB] (1) Real[1] column1.mediumLiquid.d (start = {1.0 for $f2 in 1:1}, min = {0.0 for $f4 in 1:1}, max = {1e5 for $f3 in 1:1}, nominal = {1.0 for $f1 in 1:1})
(147|188)      [ALGB] (1) Real[1] column1.mediumLiquid.T = {column1.T_l[1]} (start = {300.0 for $f2 in 1:1}, min = {1.0 for $f4 in 1:1}, max = {1e4 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(148|189)      [ALGB] (1) Real[1] column1.mediumLiquid.h = {if column1.homotopyMethod.bool_h and column1.homotopyMethod.useHomotopy then homotopy(column1.h_l[1], ({column1.homotopyMethod.h_liq for $f1 in 1:1})[1]) else column1.h_l[1]}
(149|190)      [ALGB] (2) stream Real[2] column1.upStreamOut.x_outflow (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(150|192)      [DISC] (1) Boolean column1.balanceEquations.homotopyMethod.bool_Edot_inter = column1.balanceEquations.homotopyMethod.whichHomotopy[2]
(151|193)      [ALGB] (1) stream Real[1] column1.internalFeedPort.feedVapourInternal.h_outflow
(152|194)      [ALGB] (2) Real[1, 2] column1.mediumLiquid.p_sat
(153|196)      [ALGB] (1) stream Real column1.downStreamIn.h_outflow
(154|197)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.h = {column1.mediumLiquid[1].properties.h}
(155|198)      [DISC] (1) Boolean column1.balanceEquations.homotopyMethod.bool_dp = column1.balanceEquations.homotopyMethod.whichHomotopy[3]
(156|199)      [ALGB] (1) Real[1] column1.mediumVapour.calcSpecificEnthalpy.T = {column1.mediumVapour[1].T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(157|200)      [ALGB] (1) Real[1] column1.liquidLoad = (column1.Vdot_l .* 3600.0) / column1.geometry.A
(158|201)      [ALGB] (1) input Real[1] column1.balanceEquations.stateLiq.T = {column1.mediumLiquid[1].state.T} (start = {300.0 for $f14 in 1:1}, min = {1.0 for $f16 in 1:1}, max = {1e4 for $f15 in 1:1}, nominal = {300.0 for $f13 in 1:1})
(159|202)      [ALGB] (1) Real[1] column1.propsVap.cp = {column1.mediumVapour[1].properties.cp}
(160|203)      [ALGB] (1) Real[1] column1.mediumLink.mediumConIn.MM (min = {0.0 for $f1 in 1:1})
(161|204)      [ALGB] (1) Real column1.mediumVapourIn.u (nominal = 2e4)
(162|205)      [ALGB] (1) Real[1] column1.mediumLiquid.mediumConOut.MM = {column1.mediumLiquid[1].MM} (min = {0.0 for $f1 in 1:1})
(163|206)      [ALGB] (1) Real column1.Vdot_l_in (start = 0.03, nominal = 1e-4)
(164|207)      [ALGB] (2) protected Real[2] column1.mediumVapourIn.calcSpecificEnthalpy.T_sat = {231.1, 309.2} (start = {288.15 for $f2 in 1:2}, min = {0.0 for $f3 in 1:2}, nominal = {300.0 for $f1 in 1:2})
(165|209)      [ALGB] (1) output Real sourceGas_Vdot.medium.properties.h
(166|210)      [ALGB] (1) flow Real[1] column1.internalFeedPort.feedLiquidInternal.Ndot (min = {-1e12 for $f2 in 1:1}, max = {0.0 for $f1 in 1:1})
(167|211)      [ALGB] (2) Real[2] column1.balanceEquations.c_l_in = column1.propsLiqIn.c
(168|213)      [ALGB] (1) Real[1] column1.balanceEquations.hu_stat = column1.balanceEquations.hu_stat
(169|214)      [ALGB] (1) protected Real[1] column1.pressureLoss.w_sup_l = column1.Vdot_l / column1.pressureLoss.geometry.A
(170|215)      [ALGB] (2) Real[1, 2] column1.bubblePressure.henryCoefficient.x_l = {column1.bubblePressure[1].x_l} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(171|217)      [ALGB] (1) Real column1.T_v_in (start = column1.T_v_start[1], min = 1.0, max = 1e4, nominal = 300.0)
(172|218)      [ALGB] (1) Real[1] column1.mediumLiquid.v (start = {1e-5 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, StateSelect = never)
(173|219)      [ALGB] (1) input Real[1] column1.balanceEquations.stateVap.T = {column1.mediumVapour[1].state.T} (start = {300.0 for $f9 in 1:1}, min = {1.0 for $f11 in 1:1}, max = {1e4 for $f10 in 1:1}, nominal = {300.0 for $f8 in 1:1})
(174|220)      [ALGB] (1) protected Real[1] column1.pressureLoss.w_sup_v = column1.Vdot_v / column1.pressureLoss.geometry.A
(175|221)      [ALGB] (2) Real[1, 2] column1.c_v (start = {10.0 for $f3 in 1:2, $f4 in 1:1}, StateSelect = default)
(176|223)      [ALGB] (2) protected Real[1, 2] column1.mediumLiquid.C = {{-2.23017, -2.87593} for $f1 in 1:1}
(177|225)      [ALGB] (2) Real[1, 2] column1.c_l_feed
(178|227)      [ALGB] (1) Real[1] column1.MM_v_feed (min = max({0.0}, {0.0 for $f1 in 1:1}))
(179|228)      [ALGB] (1) Real[1] column1.initOption.p_v = column1.p_hyd[1:1]
(180|229)      [ALGB] (1) Real[1] column1.results.Vdot_v = column1.Vdot_v
(181|230)      [ALGB] (1) Real[1] column1.heatTransferWall.T_wall (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(182|231)      [ALGB] (1) Real column1.balanceEquations.rho_l_in = column1.propsLiqIn.d (min = 0.0)
(183|232)      [ALGB] (2) input Real[1, 2] column1.reaction.propsLiq.eta_comp = {column1.mediumLiquid[1].properties.eta_comp} (min = {0.0 for $f13 in 1:2, $f14 in 1:1})
(184|234)      [ALGB] (2) protected Real[1, 2] column1.mediumLiquid.A = {{-6.71791, -7.27889} for $f1 in 1:1}
(185|236)      [ALGB] (1) Real[1] column1.internalFeedPort.feedVapourInternal.p
(186|237)      [ALGB] (2) Real[1, 2] column1.balanceEquations.c_v_star (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(187|239)      [ALGB] (1) Real column1.mediumLiquidIn.state.d (start = 1.0, min = 0.0, max = 1e5, nominal = 1.0)
(188|240)      [ALGB] (1) Real column1.h_upStreamOut_act
(189|241)      [ALGB] (1) Real sourceLiquid.medium.R
(190|242)      [ALGB] (1) output Real sourceGas_Vdot.medium.properties.u (nominal = 2e4)
(191|243)      [ALGB] (2) Real[1, 2] column1.results.x_v_star = column1.x_v_star (min = {0.0 for $f33 in 1:2, $f34 in 1:1}, max = {1.0 for $f31 in 1:2, $f32 in 1:1})
(192|245)      [ALGB] (2) Real[1, 2] column1.propsLiq.x = {column1.mediumLiquid[1].properties.x} (min = {0.0 for $f38 in 1:2, $f39 in 1:1}, max = {1.0 for $f36 in 1:2, $f37 in 1:1})
(193|247)      [ALGB] (1) Real column1.MM_v_in = column1.MM_v_in (start = 0.03, min = 0.001, max = 0.25)
(194|248)      [DISC] (1) Boolean $SEV_6
(195|249)      [DISC] (1) Boolean $SEV_5
(196|250)      [ALGB] (2) Real[2] sourceLiquid.X = ({0.0441, 0.07215} .* sourceLiquid.x[:]) / sourceLiquid.medium.MM
(197|252)      [DISC] (1) Boolean $SEV_4
(198|253)      [ALGB] (2) Real[1, 2] column1.balanceEquations.K
(199|255)      [ALGB] (2) protected Real[1, 2] column1.reaction.x = {column1.reaction[1].propsLiq.x} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(200|257)      [ALGB] (2) Real[1, 2] column1.balanceEquations.n_mol (min = {0.0 for $f1 in 1:2, $f2 in 1:1})
(201|259)      [ALGB] (1) Real[1] column1.Ndot_v
(202|260)      [ALGB] (2) protected Real[2] column1.mediumLiquidIn.B = {1.33932, 1.48754}
(203|262)      [ALGB] (1) Real column1.Ndot_l_in (min = 1e12, max = 0.0)
(204|263)      [ALGB] (2) Real[2] column1.x_v_in (min = max({0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}), max = max({1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}))
(205|265)      [ALGB] (1) Real column1.mediumVapourIn.lambda
(206|266)      [ALGB] (1) Real[1] column1.balanceEquations.h_l = column1.propsLiq.h
(207|267)      [ALGB] (2) Real[1, 2] column1.balanceEquations.p_sat
(208|269)      [ALGB] (2) Real[2] sourceLiquid.medium.c = (sourceLiquid.rho .* sourceLiquid.x) / sourceLiquid.medium.MM
(209|271)      [ALGB] (2) Real[1, 2] column1.mediumLiquid.calcSpecificEnthalpy.x = {column1.mediumLiquid[1].x} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(210|273)      [ALGB] (2) input Real[1, 2] column1.reaction.propsLiq.c = {column1.mediumLiquid[1].properties.c}
(211|275)      [ALGB] (2) Real[1, 2] column1.n_tot
(212|277)      [ALGB] (2) Real[2] column1.p_hyd
(213|279)      [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.henryCoefficient.henryCoeff.He
(214|281)      [ALGB] (1) Real[1] column1.propsLiq.d = {column1.mediumLiquid[1].properties.d} (min = {0.0 for $f44 in 1:1})
(215|282)      [ALGB] (1) Real[1] column1.MM_l = column1.MM_l (start = {0.018 for $f1 in 1:1}, min = max({0.0 for $f2 in 1:1}, {0.001}), max = {0.25})
(216|283)      [ALGB] (1) protected Real sourceGas_Vdot.medium.R
(217|284)      [DISC] (6) Boolean[6] column1.balanceEquations.homotopyMethod.whichHomotopy = {false for $f1 in 1:6}
(218|290)      [ALGB] (1) Real sourceLiquid.medium.state.d (start = 1.0, min = 0.0, max = 1e5, nominal = 1.0)
(219|291)      [ALGB] (1) Real sourceLiquid.medium.properties.d (min = 0.0)
(220|292)      [ALGB] (1) input Real[1] column1.balanceEquations.stateLiq.d = {column1.mediumLiquid[1].state.d} (start = {1.0 for $f7 in 1:1}, min = {0.0 for $f9 in 1:1}, max = {1e5 for $f8 in 1:1}, nominal = {1.0 for $f6 in 1:1})
(221|293)      [ALGB] (2) Real[2] sourceGas_Vdot.medium.calcSpecificEnthalpy.hX
(222|295)      [ALGB] (1) Real[1] column1.n_vap = sum(column1.n_i_vap[:, i] for i in 1:2)
(223|296)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.d = {column1.mediumLiquid[1].properties.d} (min = {0.0 for $f44 in 1:1})
(224|297)      [ALGB] (1) Real sourceLiquid.medium.properties.u (nominal = 2e4)
(225|298)      [ALGB] (1) Real sourceLiquid.medium.v (start = 1e-5, min = 0.0, StateSelect = never)
(226|299)      [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.p_partial
(227|301)      [ALGB] (2) Real[1, 2] column1.results.c_v = column1.c_v
(228|303)      [ALGB] (1) protected Real[1] column1.holdup.rho = column1.propsLiq.rho (min = {0.0 for $f1 in 1:1})
(229|304)      [ALGB] (2) Real[1, 2] column1.results.c_l = column1.c_l
(230|306)      [ALGB] (2) protected Real[2] column1.mediumLiquidIn.C = {-2.23017, -2.87593}
(231|308)      [ALGB] (1) Real sourceLiquid.h (start = -6e4)
(232|309)      [ALGB] (2) Real[1, 2] column1.propsVap.c = {column1.mediumVapour[1].properties.c}
(233|311)      [ALGB] (1) flow Real column1.heatPort.Qdot = -sum(column1.heatTransferWall.Qdot_out)
(234|312)      [ALGB] (2) Real[1, 2] column1.mediumVapour.calcSpecificEnthalpy.hX
(235|314)      [ALGB] (2) stream Real[2] column1.downStreamOut.x_outflow (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(236|316)      [ALGB] (1) Real[1] column1.mediumVapourLink.mediumConIn.rho (min = {0.0 for $f1 in 1:1})
(237|317)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.T = {column1.mediumLiquid[1].properties.T} (start = {300.0 for $f2 in 1:1}, min = {1.0 for $f4 in 1:1}, max = {1e4 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(238|318)      [ALGB] (1) Real[1] column1.eps_liq (start = {column1.eps_liq_start for $f2 in 1:1}, StateSelect = default)
(239|319)      [ALGB] (1) Real[1] column1.heatTransferWall.Qdot_out
(240|320)      [ALGB] (1) Real[1] column1.mediumVapour.cp
(241|321)      [ALGB] (1) protected Real[1] column1.pressureLoss.eta_l = column1.propsLiq.eta (min = {0.0 for $f1 in 1:1})
(242|322)      [ALGB] (1) Real[1] column1.evapEnthalpy.T = {column1.T_v[1]} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(243|323)      [ALGB] (2) input Real[1, 2] column1.reaction.propsLiq.x = {column1.mediumLiquid[1].properties.x} (min = {0.0 for $f38 in 1:2, $f39 in 1:1}, max = {1.0 for $f36 in 1:2, $f37 in 1:1})
(244|325)      [ALGB] (2) Real[1, 2] column1.mediumLiquid.x = {column1.x_l[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(245|327)      [ALGB] (1) Real[1] column1.rho_v_feed (min = max({0.0}, {0.0 for $f1 in 1:1}))
(246|328)      [ALGB] (1) Real column1.h_downStreamOut_act
(247|329)      [ALGB] (2) Real[1, 2] column1.bubblePressure.x_l = {column1.x_l[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(248|331)      [ALGB] (1) protected Real[1] column1.reaction.p = {column1.reaction[1].propsLiq.p}
(249|332)      [ALGB] (1) Real[1] column1.propsLiq.cp = {column1.mediumLiquid[1].properties.cp}
(250|333)      [ALGB] (1) Real sourceLiquid.medium.properties.h
(251|334)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.MM = {column1.mediumLiquid[1].properties.MM} (min = {0.001 for $f28 in 1:1}, max = {0.25 for $f27 in 1:1}, nominal = {0.032 for $f26 in 1:1})
(252|335)      [ALGB] (1) Real[1] column1.pressureLoss.hu_dyn = column1.pressureLoss.hu_dyn
(253|336)      [ALGB] (2) stream Real[2] column1.upStreamIn.x_outflow (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(254|338)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.rho = {column1.mediumLiquid[1].properties.rho} (start = {1.0 for $f20 in 1:1}, min = {0.0 for $f22 in 1:1}, max = {1e5 for $f21 in 1:1}, nominal = {1.0 for $f19 in 1:1})
(255|339)      [ALGB] (1) Real[1] column1.rho_l_feed (min = max({0.0}, {0.0 for $f1 in 1:1}))
(256|340)      [ALGB] (1) Real column1.mediumLiquidIn.cp
(257|341)      [ALGB] (2) Real[2] column1.mediumVapourIn.calcSpecificEnthalpy.hX
(258|343)      [ALGB] (1) Real $FUN_34
(259|344)      [ALGB] (1) Real $FUN_33
(260|345)      [ALGB] (1) Real $FUN_32
(261|346)      [ALGB] (1) Real column1.rho_v_in = column1.rho_v_in (min = 0.0, max = 1e5)
(262|347)      [ALGB] (2) Real[2] column1.mediumLiquidIn.eta_comp (min = {0.0 for $f1 in 1:2})
(263|349)      [ALGB] (1) Real $FUN_31
(264|350)      [ALGB] (2) Real[1, 2] column1.c_l (start = {1.0 for $f3 in 1:2, $f4 in 1:1}, StateSelect = default)
(265|352)      [ALGB] (1) Real $FUN_30
(266|353)      [ALGB] (1) Real sourceGas_Vdot.MM = sourceGas_Vdot.MM (start = 0.02, min = 0.001, max = 0.25)
(267|354)      [ALGB] (2) Real[2] column1.mediumVapourIn.X (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(268|356)      [ALGB] (1) stream Real[1] column1.internalFeedPort.feedLiquidInternal.h_outflow
(269|357)      [ALGB] (2) stream Real[2] sinkGas.gasPortIn.x_outflow (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(270|359)      [ALGB] (1) Real[1] column1.reaction.V = column1.H * column1.A * column1.eps_liq
(271|360)      [ALGB] (2) Real[1, 2] column1.bubblePressure.henryCoefficient.henryCoeff.x_l = {column1.bubblePressure[1].henryCoefficient.x_l} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(272|362)      [ALGB] (1) protected Real[1] column1.pressureLoss.rho_v = column1.propsVap.rho (min = {0.0 for $f1 in 1:1})
(273|363)      [ALGB] (1) Real column1.propsLiqIn.v = column1.mediumLiquidIn.properties.v (start = 1e-5, min = 0.0, StateSelect = never)
(274|364)      [ALGB] (2) Real[1, 2] column1.results.x_l = column1.x_l (min = {0.0 for $f9 in 1:2, $f10 in 1:1}, max = {1.0 for $f7 in 1:2, $f8 in 1:1})
(275|366)      [ALGB] (2) Real[2] column1.propsVapIn.c = column1.mediumVapourIn.properties.c
(276|368)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.cp = {column1.mediumLiquid[1].properties.cp}
(277|369)      [ALGB] (2) Real[1, 2] column1.Ndot_reac
(278|371)      [ALGB] (1) Real[1] column1.Vdot_v_feed
(279|372)      [ALGB] (1) Real[1] column1.Qdot_reac
(280|373)      [ALGB] (1) Real[1] column1.mediumLiquid.lambda
(281|374)      [ALGB] (1) Real[1] column1.mediumVapour.mediumConOut.MM = {column1.mediumVapour[1].MM} (min = {0.0 for $f1 in 1:1})
(282|375)      [ALGB] (1) Real[1] column1.balanceEquations.thermoEquilibrium.henryCoefficient.T = {column1.balanceEquations.thermoEquilibrium[1].T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(283|376)      [ALGB] (1) Real column1.propsVapIn.rho = column1.mediumVapourIn.properties.rho (start = 1.0, min = 0.0, max = 1e5, nominal = 1.0)
(284|377)      [ALGB] (1) Real[1] column1.Qdot_wall
(285|378)      [ALGB] (1) Real column1.mediumLiquidIn.p = column1.p_hyd[2] (start = 1e5, min = 0.0, max = 1e8, nominal = 1e5)
(286|379)      [ALGB] (2) Real[1, 2] column1.activityCoeff.activityCoeff.x_l = {column1.activityCoeff[1].x_l} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(287|381)      [ALGB] (1) Real[1] column1.T_v (start = column1.T_v_start, min = max({0.0}, {0.0 for $f2 in 1:1}), nominal = {300.0 for $f1 in 1:1})
(288|382)      [ALGB] (1) Real[1] column1.balanceEquations.n_mol_V (min = {0.0 for $f2 in 1:1}, StateSelect = default)
(289|383)      [ALGB] (1) Real $FUN_29
(290|384)      [ALGB] (2) Real[1, 2] column1.propsVap.x = {column1.mediumVapour[1].properties.x} (min = {0.0 for $f25 in 1:2, $f26 in 1:1}, max = {1.0 for $f23 in 1:2, $f24 in 1:1})
(291|386)      [ALGB] (1) Real[1] column1.Vdot_le
(292|387)      [ALGB] (2) Real[2, 1] $FUN_9
(293|389)      [DISC] (1) Integer $FUN_28
(294|390)      [ALGB] (1) Real $FUN_8
(295|391)      [ALGB] (1) Real $FUN_27
(296|392)      [ALGB] (1) Real $FUN_26
(297|393)      [ALGB] (1) Real column1.propsLiqIn.T = column1.mediumLiquidIn.properties.T (start = 300.0, min = 1.0, max = 1e4, nominal = 300.0)
(298|394)      [ALGB] (1) Real $FUN_6
(299|395)      [ALGB] (1) Real $FUN_25
(300|396)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.v = {column1.mediumLiquid[1].properties.v} (start = {1e-5 for $f32 in 1:1}, min = {0.0 for $f33 in 1:1}, StateSelect = never)
(301|397)      [ALGB] (1) Real $FUN_5
(302|398)      [ALGB] (1) Real column1.balanceEquations.h_l_in = column1.propsLiqIn.h
(303|399)      [ALGB] (1) Real $FUN_4
(304|400)      [ALGB] (1) Real $FUN_22
(305|401)      [ALGB] (1) Real $FUN_2
(306|402)      [ALGB] (1) Real $FUN_21
(307|403)      [ALGB] (1) Real[1] column1.rho_l = if column1.homotopyMethod.bool_rho and false then {homotopy(column1.mediumLiquid[1].d, ({1000.0 for $f1 in 1:1})[1])} else column1.mediumLiquid.d (start = {1.0 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1})
(308|404)      [ALGB] (1) output Real sourceGas_Vdot.medium.properties.d (min = 0.0)
(309|405)      [ALGB] (1) Real $FUN_20
(310|406)      [ALGB] (2) Real[1, 2] column1.c_v_feed
(311|408)      [ALGB] (1) Real[1] column1.balanceEquations.thermoEquilibrium.p = {column1.balanceEquations.p_v[1]}
(312|409)      [ALGB] (1) Real column1.balanceEquations.rho_v_in = column1.propsVapIn.d (min = 0.0)
(313|410)      [ALGB] (2) Real[1, 2] column1.mediumLiquid.c = {(column1.mediumLiquid[1].x .* column1.mediumLiquid[1].d) / column1.mediumLiquid[1].MM}
(314|412)      [ALGB] (2) Real[2] column1.mediumLiquidIn.c = (column1.rho_l_in .* column1.x_l_in) / column1.MM_l_in
(315|414)      [ALGB] (1) Real[1] column1.holdup.Vdot (start = {0.03 for $f1 in 1:1})
(316|415)      [ALGB] (1) Real[1] column1.mediumVapour.mediumConOut.rho = {column1.mediumVapour[1].d} (min = {0.0 for $f1 in 1:1})
(317|416)      [ALGB] (1) Real[1] column1.mediumVapour.mediumConOut.h = {column1.mediumVapour[1].h}
(318|417)      [ALGB] (2) stream Real[2] column1.downStreamIn.x_outflow (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(319|419)      [ALGB] (2) Real[2] column1.propsLiqIn.c = column1.mediumLiquidIn.properties.c
(320|421)      [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.henryCoefficient.x_l = {column1.balanceEquations.thermoEquilibrium[1].x_l} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(321|423)      [ALGB] (2) Real[1, 2] column1.x_v_star (start = column1.x_v_start, min = max({0.0 for $f3 in 1:2, $f4 in 1:1}, {0.0 for $f3 in 1:2, $f4 in 1:1}, {0.0 for $f3 in 1:2, $f4 in 1:1}, {0.0 for $f3 in 1:2, $f4 in 1:1}), max = max({1.0 for $f1 in 1:2, $f2 in 1:1}, {1.0 for $f1 in 1:2, $f2 in 1:1}, {1.0 for $f1 in 1:2, $f2 in 1:1}, {1.0 for $f1 in 1:2, $f2 in 1:1}))
(322|425)      [ALGB] (2) Real[2] column1.balanceEquations.x_v_in = column1.propsVapIn.x (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(323|427)      [ALGB] (1) Real[1] column1.propsLiq.u = {column1.mediumLiquid[1].properties.u} (nominal = {2e4 for $f54 in 1:1})
(324|428)      [ALGB] (1) output Real sourceGas_Vdot.medium.properties.rho (start = 1.0, min = 0.0, max = 1e5, nominal = 1.0)
(325|429)      [ALGB] (1) Real[1] column1.MM_v = column1.MM_v (start = 0.028 * {1.0 for $f1 in 1:1}, min = max({0.0 for $f2 in 1:1}, {0.001}), max = {0.25})
(326|430)      [ALGB] (1) Real[1] column1.mediumVapour.h
(327|431)      [ALGB] (1) Real[1] column1.propsVap.d = {column1.mediumVapour[1].properties.d} (min = {0.0 for $f31 in 1:1})
(328|432)      [ALGB] (1) Real column1.balanceEquations.MM_l_in = column1.propsLiqIn.MM (min = 0.0)
(329|433)      [ALGB] (1) Real[1] column1.activityCoeff.T = {column1.T_star[1]} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(330|434)      [ALGB] (1) Real column1.propsLiqIn.rho = column1.mediumLiquidIn.properties.rho (start = 1.0, min = 0.0, max = 1e5, nominal = 1.0)
(331|435)      [ALGB] (2) Real[2, 1] $FUN_19
(332|437)      [ALGB] (2) Real[1, 2] column1.x_vap_liq (min = max({0.0 for $f3 in 1:2, $f4 in 1:1}, {0.0 for $f3 in 1:2, $f4 in 1:1}), max = max({1.0 for $f1 in 1:2, $f2 in 1:1}, {1.0 for $f1 in 1:2, $f2 in 1:1}))
(333|439)      [ALGB] (1) Real column1.p_v_in (start = 1.5e5, min = 0.0, max = 1e8, StateSelect = default)
(334|440)      [ALGB] (1) Real $FUN_14
(335|441)      [ALGB] (1) Real[1] column1.n_liq = sum(column1.n_i_liq[:, i] for i in 1:2)
(336|442)      [ALGB] (2) Real[2, 1] $FUN_13
(337|444)      [ALGB] (1) protected Real column1.mediumLiquidIn.Theta = min(1.0, column1.T_l_in / 647.14)
(338|445)      [ALGB] (1) Real[1] column1.balanceEquations.thermoEquilibrium.p_bubble
(339|446)      [ALGB] (1) Real $FUN_12
(340|447)      [ALGB] (2) Real[2] $FUN_11
(341|449)      [ALGB] (1) Real[1] column1.balanceEquations.n_mol_L (min = {0.0 for $f2 in 1:1}, StateSelect = default)
(342|450)      [ALGB] (1) Real $FUN_10
(343|451)      [ALGB] (1) stream Real column1.upStreamOut.h_outflow
(344|452)      [ALGB] (1) Real column1.propsVapIn.T = column1.mediumVapourIn.properties.T (start = 300.0, min = 1.0, max = 1e4, nominal = 300.0)
(345|453)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.p = {column1.mediumLiquid[1].properties.p}
(346|454)      [ALGB] (2) Real[2] column1.c_l_in
(347|456)      [ALGB] (1) Real column1.h_upStreamIn_act
(348|457)      [ALGB] (2) Real[1, 2] column1.balanceEquations.x_l = column1.propsLiq.x (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(349|459)      [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.K
(350|461)      [ALGB] (1) protected Real[1] column1.mediumLiquid.Theta = {min(1.0, column1.mediumLiquid[1].T / 647.14)}
(351|462)      [ALGB] (1) Real[1] column1.propsVap.lambda = {column1.mediumVapour[1].properties.lambda}
(352|463)      [ALGB] (1) Real column1.mediumLiquidIn.u (nominal = 2e4)
(353|464)      [ALGB] (1) Real column1.balanceEquations.MM_v_in = column1.propsVapIn.MM (min = 0.0)
(354|465)      [ALGB] (1) Real column1.propsVapIn.u = column1.mediumVapourIn.properties.u (nominal = 2e4)
(355|466)      [ALGB] (1) Real column1.mediumLiquidIn.state.h (min = -1e8, max = 1e8, nominal = 1e6)
(356|467)      [ALGB] (2) Real[1, 2] column1.propsLiq.c = {column1.mediumLiquid[1].properties.c}
(357|469)      [ALGB] (1) output Real column1.mediumVapourIn.state.p (start = 1e5, min = 0.0, max = 1e8, nominal = 1e5)
(358|470)      [ALGB] (1) Real[1] column1.bubblePressure.henryCoefficient.henryCoeff.T = {column1.bubblePressure[1].henryCoefficient.T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(359|471)      [ALGB] (1) Real[1] column1.results.T_l = column1.T (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(360|472)      [ALGB] (1) Real[1] column1.eps_vap (start = {0.99 for $f1 in 1:1})
(361|473)      [ALGB] (2) Real[1, 2] column1.bubblePressure.x_v = {column1.x_v[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(362|475)      [ALGB] (1) Real sourceGas_Vdot.medium.u (nominal = 2e4)
(363|476)      [ALGB] (2) Real[2] column1.mediumLiquidIn.p_sat
(364|478)      [ALGB] (2) Real[1, 2] column1.p_sat_bulk (start = {column1.p_v_start[1] for $f1 in 1:2, $f2 in 1:1})
(365|480)      [ALGB] (2) Real[1, 2] column1.reaction.c = {column1.c_l[1]}
(366|482)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.lambda = {column1.mediumLiquid[1].properties.lambda}
(367|483)      [ALGB] (2) protected Real[1, 2] column1.mediumLiquid.B = {{1.33932, 1.48754} for $f1 in 1:1}
(368|485)      [ALGB] (1) protected Real[1] column1.holdup.sigma = column1.propsLiq.sigma
(369|486)      [ALGB] (1) flow Real sinkLiquid.liquidPortIn.Ndot (min = 1e12, max = 0.0)
(370|487)      [ALGB] (1) Real[1] column1.mediumLink.mediumConIn.rho (min = {0.0 for $f1 in 1:1})
(371|488)      [ALGB] (1) Real column1.propsLiqIn.eta = column1.mediumLiquidIn.properties.eta (min = 0.0)
(372|489)      [ALGB] (2) Real[2] sourceLiquid.medium.properties.c
(373|491)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.sigma = {column1.mediumLiquid[1].properties.sigma}
(374|492)      [ALGB] (1) Real[1] column1.balanceEquations.MM_v = column1.propsVap.MM (start = 0.028 * {1.0 for $f1 in 1:1}, min = {0.0 for $f2 in 1:1})
(375|493)      [ALGB] (1) Real[1] column1.mediumLiquid.sigma
(376|494)      [ALGB] (2) Real[1, 2] column1.balanceEquations.thermoEquilibrium.x_v = {column1.balanceEquations.x_v_star[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(377|496)      [ALGB] (2) Real[2] column1.x_downStreamIn_act (min = max({0.0 for $f2 in 1:2}, {0.0 for $f2 in 1:2}), max = max({1.0 for $f1 in 1:2}, {1.0 for $f1 in 1:2}))
(378|498)      [ALGB] (1) Real[1] column1.F = (column1.Vdot_v / column1.geometry.A) .* column1.rho_v .^ 0.5
(379|499)      [ALGB] (1) Real column1.mediumLiquidIn.lambda
(380|500)      [ALGB] (1) Real[1] column1.propsLiq.v = {column1.mediumLiquid[1].properties.v} (start = {1e-5 for $f32 in 1:1}, min = {0.0 for $f33 in 1:1}, StateSelect = never)
(381|501)      [ALGB] (1) Real[1] column1.internalFeedPort.feedLiquidInternal.p
(382|502)      [ALGB] (1) Real column1.propsLiqIn.u = column1.mediumLiquidIn.properties.u (nominal = 2e4)
(383|503)      [ALGB] (1) Real[1] column1.activityCoeff.activityCoeff.T = {column1.activityCoeff[1].T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(384|504)      [ALGB] (1) flow Real[1] column1.internalFeedPort.feedVapourInternal.Ndot
(385|505)      [ALGB] (1) Real[1] column1.bubblePressure.p = {column1.p_v[1]}
(386|506)      [ALGB] (2) protected Real[1, 2] column1.mediumLiquid.D = {{-1.2499, -2.22146} for $f1 in 1:1}
(387|508)      [ALGB] (2) Real[1, 2] column1.bubblePressure.henryCoefficient.henryCoeff.He
(388|510)      [ALGB] (2) Real[1, 2] column1.mediumLiquid.eta_comp (min = {0.0 for $f1 in 1:2, $f2 in 1:1})
(389|512)      [ALGB] (1) input Real[1] column1.balanceEquations.stateVap.p = {column1.mediumVapour[1].state.p} (start = {1e5 for $f2 in 1:1}, min = {0.0 for $f4 in 1:1}, max = {1e8 for $f3 in 1:1}, nominal = {1e5 for $f1 in 1:1})
(390|513)      [DISC] (1) Boolean column1.balanceEquations.homotopyMethod.bool_K = column1.balanceEquations.homotopyMethod.whichHomotopy[4]
(391|514)      [ALGB] (1) Real column1.Vdot_v_in (nominal = 0.01)
(392|515)      [ALGB] (1) Real column1.propsVapIn.h = column1.mediumVapourIn.properties.h
(393|516)      [ALGB] (1) Real[1] column1.propsLiq.h = {column1.mediumLiquid[1].properties.h}
(394|517)      [ALGB] (2) Real[2] sourceGas_Vdot.c
(395|519)      [ALGB] (1) Real[1] column1.mediumVapour.p = {column1.p_v[1]} (start = {1e5 for $f2 in 1:1}, min = {0.0 for $f4 in 1:1}, max = {1e8 for $f3 in 1:1}, nominal = {1e5 for $f1 in 1:1})
(396|520)      [ALGB] (1) Real[1] column1.balanceEquations.MM_l = column1.propsLiq.MM (start = {0.018 for $f1 in 1:1}, min = {0.0 for $f2 in 1:1})
(397|521)      [ALGB] (1) Real[1] column1.rho_v = if column1.homotopyMethod.bool_rho and false then {homotopy(column1.mediumVapour[1].d, ({1.0 for $f1 in 1:1})[1])} else column1.mediumVapour.d (min = {0.0 for $f2 in 1:1})
(398|522)      [ALGB] (2) Real[1, 2] column1.propsLiq.eta_comp = {column1.mediumLiquid[1].properties.eta_comp} (min = {0.0 for $f13 in 1:2, $f14 in 1:1})
(399|524)      [ALGB] (1) input Real[1] column1.reaction.propsLiq.eta = {column1.mediumLiquid[1].properties.eta} (min = {0.0 for $f10 in 1:1})
(400|525)      [ALGB] (1) Real column1.T_l_in (start = column1.T_l_start[1], min = 1.0, max = 1e4, nominal = 300.0)
(401|526)      [ALGB] (2) input Real[1, 2] column1.balanceEquations.stateVap.X = {column1.mediumVapour[1].state.X} (min = {0.0 for $f19 in 1:2, $f20 in 1:1}, max = {1.0 for $f17 in 1:2, $f18 in 1:1}, nominal = {0.1 for $f15 in 1:2, $f16 in 1:1})
(402|528)      [ALGB] (1) Real[1] column1.mediumLiquid.calcSpecificEnthalpy.p = {column1.mediumLiquid[1].p}
(403|529)      [ALGB] (1) Real[1] column1.mdot_l = column1.Vdot_l .* column1.rho_l
(404|530)      [ALGB] (2) stream Real[1, 2] column1.internalFeedPort.feedLiquidInternal.x_outflow (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(405|532)      [ALGB] (2) output Real[2] column1.mediumVapourIn.state.X (min = {0.0 for $f3 in 1:2}, max = {1.0 for $f2 in 1:2}, nominal = {0.1 for $f1 in 1:2})
(406|534)      [ALGB] (1) Real column1.mediumLiquidIn.state.T (start = 300.0, min = 1.0, max = 1e4, nominal = 300.0)
(407|535)      [ALGB] (2) Real[2] sourceGas_Vdot.medium.X (min = {0.0 for $f2 in 1:2}, max = {1.0 for $f1 in 1:2})
(408|537)      [ALGB] (1) Real column1.mediumLiquidIn.R
(409|538)      [ALGB] (1) Real[1] column1.balanceEquations.T = column1.propsLiq.T (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(410|539)      [ALGB] (2) Real[1, 2] column1.bubblePressure.K
(411|541)      [ALGB] (2) Real[1, 2] column1.n_mol_L_i
(412|543)      [ALGB] (1) protected Real[1] column1.reaction.T = {column1.reaction[1].propsLiq.T} (start = {288.15 for $f2 in 1:1}, min = {0.0 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(413|544)      [ALGB] (2) Real[1, 2] column1.X_v
(414|546)      [ALGB] (1) Real[1] column1.propsLiq.p = {column1.mediumLiquid[1].properties.p}
(415|547)      [ALGB] (2) Real[1, 2] column1.activityCoeff.x_l = {column1.x_l_star[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})
(416|549)      [ALGB] (1) Real column1.mediumLiquidIn.eta (min = 0.0)
(417|550)      [ALGB] (1) Real[1] column1.propsVap.MM = {column1.mediumVapour[1].properties.MM} (min = {0.001 for $f20 in 1:1}, max = {0.25 for $f19 in 1:1}, nominal = {0.032 for $f18 in 1:1})
(418|551)      [ALGB] (1) Real[1] column1.propsVap.u = {column1.mediumVapour[1].properties.u} (nominal = {2e4 for $f41 in 1:1})
(419|552)      [ALGB] (1) Real[1] column1.mediumLiquid.mediumConOut.h = {column1.mediumLiquid[1].h}
(420|553)      [ALGB] (1) Real[1] column1.propsLiq.T = {column1.mediumLiquid[1].properties.T} (start = {300.0 for $f2 in 1:1}, min = {1.0 for $f4 in 1:1}, max = {1e4 for $f3 in 1:1}, nominal = {300.0 for $f1 in 1:1})
(421|554)      [ALGB] (1) Real column1.MM_l_in = column1.MM_l_in (min = 0.001, max = 0.25)
(422|555)      [ALGB] (1) Real[1] column1.Ndot_source_startUp
(423|556)      [ALGB] (1) Real column1.propsLiqIn.h = column1.mediumLiquidIn.properties.h
(424|557)      [ALGB] (2) Real[1, 2] column1.mediumVapour.x = {column1.x_v[1]} (min = {0.0 for $f3 in 1:2, $f4 in 1:1}, max = {1.0 for $f1 in 1:2, $f2 in 1:1})


System Equations (431/557)
****************************
(1|1)          [FOR-] (2) ($RES_SIM_254)
(1|1)          [----] for {$i1 in 1:1, $i2 in 1:2} loop
(1|1)          [----]   [SCAL] (1) column1.bubblePressure[$i1].K[$i2] = (column1.bubblePressure[$i1].p_sat[({{1, 1}, {2, 2}})[$i2, 2]] * column1.bubblePressure[$i1].factor_K[({{1, 1}, {2, 2}})[$i2, 2]]) / column1.bubblePressure[$i1].p; ($RES_SIM_255)
(1|1)          [----] end for;
(2|3)          [SCAL] (1) column1.reaction[1].p = column1.reaction[1].propsLiq.p; ($RES_SIM_1043)
(3|4)          [SCAL] (1) column1.x_downStreamOut_act[2] = column1.downStreamOut.x_outflow[2]; ($RES_SIM_1042)
(4|5)          [FOR-] (2) ($RES_SIM_256)
(4|5)          [----] for {$i1 in 1:1, $i2 in 1:2} loop
(4|5)          [----]   [SCAL] (1) column1.bubblePressure[$i1].p_partial[$i2] = column1.bubblePressure[$i1].factor_K[({{1, 1}, {2, 2}})[$i2, 2]] * column1.bubblePressure[$i1].x_l[({{1, 1}, {2, 2}})[$i2, 2]] * column1.bubblePressure[$i1].p_sat[({{1, 1}, {2, 2}})[$i2, 2]]; ($RES_SIM_257)
(4|5)          [----] end for;
(5|7)          [SCAL] (1) column1.x_downStreamOut_act[1] = column1.downStreamOut.x_outflow[1]; ($RES_SIM_1041)
(6|8)          [SCAL] (1) column1.x_downStreamIn_act[2] = sourceLiquid.x[2]; ($RES_SIM_1040)
(7|9)          [SCAL] (1) $FUN_12 = sum({0.0441, 0.07215} .* column1.x_l_in); ($RES_AUX_1009)
(8|10)         [FOR-] (2) ($RES_AUX_1007)
(8|10)         [----] for {$i2 in 1:2, $i1 in 1:1} loop
(8|10)         [----]   [SCAL] (1) $FUN_13[$i2, $i1] = exp(column1.bubblePressure[$i1].henryCoefficient.henryCoeff.henry_C[$i2] * (1/column1.bubblePressure[$i1].henryCoefficient.henryCoeff.T - 1/column1.bubblePressure[$i1].henryCoefficient.henryCoeff.henry_T_norm)); ($RES_AUX_1008)
(8|10)         [----] end for;
(9|12)         [SCAL] (1) $FUN_14 = sum(column1.bubblePressure[1].p_partial); ($RES_AUX_1006)
(10|13)        [SCAL] (1) column1.propsVapIn.rho = column1.rho_v_in; ($RES_SIM_340)
(11|14)        [SCAL] (1) column1.propsVapIn.T = column1.T_v_in; ($RES_SIM_342)
(12|15)        [ARRY] (2) column1.mediumVapourIn.state.X = column1.mediumVapourIn.X; ($RES_SIM_343)
(13|17)        [SCAL] (1) column1.mediumLiquid[1].D[1] = -1.2499; ($RES_SIM_1129)
(14|18)        [SCAL] (1) column1.mediumVapourIn.state.p = column1.p_v_in; ($RES_SIM_344)
(15|19)        [SCAL] (1) column1.mediumVapourIn.state.T = column1.T_v_in; ($RES_SIM_345)
(16|20)        [SCAL] (1) column1.mediumVapourIn.lambda = 0.018 * column1.x_v_in[1] + 0.0144 * column1.x_v_in[2]; ($RES_SIM_346)
(17|21)        [SCAL] (1) column1.mediumLiquidIn.p = column1.p_hyd[2]; ($RES_BND_647)
(18|22)        [SCAL] (1) column1.mediumVapourIn.cp = 1549.0 * column1.x_v_in[1] + 2456.0 * column1.x_v_in[2]; ($RES_SIM_348)
(19|23)        [SCAL] (1) column1.mediumLiquidIn.A[2] = -7.27889; ($RES_SIM_1124)
(20|24)        [SCAL] (1) column1.rho_v_in = column1.p_v_in / (column1.T_v_in * column1.mediumVapourIn.R); ($RES_SIM_349)
(21|25)        [SCAL] (1) column1.mediumLiquidIn.A[1] = -6.71791; ($RES_SIM_1123)
(22|26)        [SCAL] (1) column1.mediumLiquidIn.B[2] = 1.48754; ($RES_SIM_1122)
(23|27)        [SCAL] (1) column1.mediumLiquidIn.B[1] = 1.33932; ($RES_SIM_1121)
(24|28)        [SCAL] (1) column1.mediumLiquidIn.C[2] = -2.87593; ($RES_SIM_1120)
(25|29)        [ARRY] (1) column1.pressureLoss.rho_l = column1.propsLiq.rho; ($RES_BND_908)
(26|30)        [SCAL] (1) column1.x_downStreamIn_act[1] = sourceLiquid.x[1]; ($RES_SIM_1039)
(27|31)        [ARRY] (1) column1.pressureLoss.sigma_l = column1.propsLiq.sigma; ($RES_BND_909)
(28|32)        [SCAL] (1) column1.x_upStreamOut_act[2] = if $SEV_6 then sinkGas.gasPortIn.x_outflow[2] else column1.upStreamOut.x_outflow[2]; ($RES_SIM_1038)
(29|33)        [SCAL] (1) -column1.Ndot_l_in = -(sourceLiquid.Flow / sourceLiquid.medium.MM) * sourceLiquid.rho; ($RES_SIM_435)
(30|34)        [FOR-] (2) ($RES_SIM_260)
(30|34)        [----] for {$i1 in 1:1, $i2 in 1:2} loop
(30|34)        [----]   [SCAL] (1) column1.bubblePressure[$i1].henryCoefficient.henryCoeff.He[$i2] = if column1.bubblePressure[$i1].henryCoefficient.henryCoeff.henry_temp then column1.bubblePressure[$i1].henryCoefficient.henryCoeff.henry_H[$i2] * $FUN_13[$i2, $i1] else column1.bubblePressure[$i1].henryCoefficient.henryCoeff.henry_H[$i2]; ($RES_SIM_261)
(30|34)        [----] end for;
(31|36)        [SCAL] (1) column1.x_upStreamOut_act[1] = if $SEV_6 then sinkGas.gasPortIn.x_outflow[1] else column1.upStreamOut.x_outflow[1]; ($RES_SIM_1037)
(32|37)        [SCAL] (1) column1.x_upStreamIn_act[2] = if $SEV_5 then sourceGas_Vdot.x[2] else column1.upStreamIn.x_outflow[2]; ($RES_SIM_1036)
(33|38)        [SCAL] (1) column1.x_upStreamIn_act[1] = if $SEV_5 then sourceGas_Vdot.x[1] else column1.upStreamIn.x_outflow[1]; ($RES_SIM_1035)
(34|39)        [SCAL] (1) column1.x_l_in[2] = sourceLiquid.x[2]; ($RES_SIM_1034)
(35|40)        [ARRY] (2) sourceLiquid.c = (sourceLiquid.x / sourceLiquid.medium.MM) .* sourceLiquid.rho; ($RES_SIM_439)
(36|42)        [SCAL] (1) column1.x_l_in[1] = sourceLiquid.x[1]; ($RES_SIM_1033)
(37|43)        [SCAL] (1) column1.x_v_in[2] = sourceGas_Vdot.x[2]; ($RES_SIM_1032)
(38|44)        [SCAL] (1) column1.propsLiqIn.h = sourceLiquid.h; ($RES_SIM_266)
(39|45)        [SCAL] (1) column1.x_v_in[1] = sourceGas_Vdot.x[1]; ($RES_SIM_1031)
(40|46)        [SCAL] (1) column1.propsLiqIn.d = column1.rho_l_in; ($RES_SIM_267)
(41|47)        [ARRY] (2) column1.propsLiqIn.c = column1.mediumLiquidIn.c; ($RES_SIM_268)
(42|49)        [ARRY] (2) column1.propsLiqIn.eta_comp = column1.mediumLiquidIn.eta_comp; ($RES_SIM_269)
(43|51)        [SCAL] (1) $SEV_5 = column1.Ndot_v_in > 0.0; ($RES_EVT_1399)
(44|52)        [SCAL] (1) $SEV_4 = column1.eps_liq[1] < 1e-5 and column1.Vdot_l_in < 1e-8; ($RES_EVT_1398)
(45|53)        [SCAL] (1) column1.mediumVapourIn.u = column1.h_v_in - column1.T_v_in * column1.mediumVapourIn.R * column1.MM_v_in; ($RES_SIM_350)
(46|54)        [ARRY] (1) column1.balanceEquations.thermoEquilibrium.v_v = column1.balanceEquations.MM_v / column1.balanceEquations.rho_v; ($RES_BND_825)
(47|55)        [ARRY] (2) column1.mediumLiquidIn.c = column1.rho_l_in * column1.x_l_in / column1.MM_l_in; ($RES_BND_650)
(48|57)        [SCAL] (1) column1.mediumVapourIn.R = 8.31446261815324 / column1.MM_v_in; ($RES_SIM_351)
(49|58)        [SCAL] (1) column1.MM_v_in = max(1e-6, $FUN_8); ($RES_SIM_353)
(50|59)        [SCAL] (1) column1.mediumLiquidIn.C[1] = -2.23017; ($RES_SIM_1119)
(51|60)        [SCAL] (1) column1.mediumVapourIn.p_start = column1.mediumVapourIn.state.p; ($RES_SIM_354)
(52|61)        [SCAL] (1) column1.mediumLiquidIn.D[2] = -2.22146; ($RES_SIM_1118)
(53|62)        [SCAL] (1) $TEV_1 = time < (ramp.startTime + ramp.duration); ($RES_EVT_1393)
(54|63)        [FOR-] (2) ($RES_SIM_355)
(54|63)        [----] for $i1 in 1:2 loop
(54|63)        [----]   [SCAL] (1) column1.mediumVapourIn.X[$i1] = max(1e-5, min(1.0, (({0.0441, 0.07215})[$i1] * column1.x_v_in[$i1]) / column1.MM_v_in)); ($RES_SIM_356)
(54|63)        [----] end for;
(55|65)        [SCAL] (1) column1.mediumLiquidIn.D[1] = -1.2499; ($RES_SIM_1117)
(56|66)        [SCAL] (1) $TEV_0 = time < ramp.startTime; ($RES_EVT_1392)
(57|67)        [SCAL] (1) column1.activityCoeff[1].T = column1.T_star[1]; ($RES_SIM_1116)
(58|68)        [SCAL] (1) column1.mediumVapourIn.calcSpecificEnthalpy.hX[2] = 0.07215 * (358000.0 + 2206.0 * ((-1406.0) + column1.mediumVapourIn.calcSpecificEnthalpy.T_sat[2]) + 1599.0 * (column1.T_v_in - column1.mediumVapourIn.calcSpecificEnthalpy.T_sat[2])); ($RES_SIM_357)
(59|69)        [ARRY] (2) column1.activityCoeff[1].x_l = column1.x_l_star[1, :]; ($RES_SIM_1115)
(60|71)        [SCAL] (1) column1.mediumVapourIn.calcSpecificEnthalpy.hX[1] = 0.0441 * (426000.0 + 2456.0 * ((-1406.0) + column1.mediumVapourIn.calcSpecificEnthalpy.T_sat[1]) + 1549.0 * (column1.T_v_in - column1.mediumVapourIn.calcSpecificEnthalpy.T_sat[1])); ($RES_SIM_358)
(61|72)        [SCAL] (1) column1.activityCoeff[1].activityCoeff.T = column1.activityCoeff[1].T; ($RES_SIM_1114)
(62|73)        [ARRY] (1) column1.pressureLoss.rho_v = column1.propsVap.rho; ($RES_BND_910)
(63|74)        [ARRY] (2) column1.activityCoeff[1].activityCoeff.x_l = column1.activityCoeff[1].x_l; ($RES_SIM_1113)
(64|76)        [ARRY] (1) column1.pressureLoss.rho_l_calc = if column1.pressureLoss.rho_const then {column1.pressureLoss.rho_l_nom} else column1.pressureLoss.rho_l; ($RES_BND_911)
(65|77)        [SCAL] (1) column1.evapEnthalpy[1].p = column1.p_hyd[1]; ($RES_SIM_1112)
(66|78)        [ARRY] (1) column1.pressureLoss.rho_v_calc = if column1.pressureLoss.rho_const then {column1.pressureLoss.rho_v_nom} else column1.pressureLoss.rho_v; ($RES_BND_912)
(67|79)        [SCAL] (1) column1.evapEnthalpy[1].T = column1.T_v[1]; ($RES_SIM_1111)
(68|80)        [SCAL] (1) column1.bubblePressure[1].T = column1.T[1]; ($RES_SIM_1110)
(69|81)        [ARRY] (1) column1.pressureLoss.w_sup_l = column1.Vdot_l / column1.pressureLoss.geometry.A; ($RES_BND_914)
(70|82)        [ARRY] (1) column1.pressureLoss.w_sup_v = column1.Vdot_v / column1.pressureLoss.geometry.A; ($RES_BND_915)
(71|83)        [SCAL] (1) sourceLiquid.medium.properties.h = sourceLiquid.h; ($RES_SIM_441)
(72|84)        [SCAL] (1) sourceLiquid.medium.properties.d = sourceLiquid.rho; ($RES_SIM_442)
(73|85)        [ARRY] (1) column1.pressureLoss.eta_l = column1.propsLiq.eta; ($RES_BND_917)
(74|86)        [ARRY] (2) sourceLiquid.medium.properties.c = sourceLiquid.medium.c; ($RES_SIM_443)
(75|88)        [SCAL] (1) sourceLiquid.medium.properties.p = sourceLiquid.p; ($RES_SIM_445)
(76|89)        [SCAL] (1) column1.propsLiqIn.p = column1.mediumLiquidIn.p; ($RES_SIM_270)
(77|90)        [SCAL] (1) sourceLiquid.medium.properties.u = sourceLiquid.medium.u; ($RES_SIM_446)
(78|91)        [SCAL] (1) column1.propsLiqIn.u = column1.mediumLiquidIn.u; ($RES_SIM_271)
(79|92)        [SCAL] (1) column1.propsLiqIn.cp = column1.mediumLiquidIn.cp; ($RES_SIM_272)
(80|93)        [SCAL] (1) column1.propsLiqIn.lambda = column1.mediumLiquidIn.lambda; ($RES_SIM_273)
(81|94)        [ARRY] (2) column1.propsLiqIn.x = column1.x_l_in; ($RES_SIM_274)
(82|96)        [SCAL] (1) column1.propsLiqIn.v = column1.mediumLiquidIn.v; ($RES_SIM_275)
(83|97)        [SCAL] (1) column1.propsLiqIn.MM = column1.MM_l_in; ($RES_SIM_276)
(84|98)        [SCAL] (1) column1.propsLiqIn.rho = column1.rho_l_in; ($RES_SIM_277)
(85|99)        [SCAL] (1) column1.propsLiqIn.eta = column1.mediumLiquidIn.eta; ($RES_SIM_279)
(86|100)       [ARRY] (1) column1.propsVap.h = column1.mediumVapour.h; ($RES_SIM_360)
(87|101)       [ARRY] (1) column1.propsVap.d = column1.mediumVapour.d; ($RES_SIM_361)
(88|102)       [SCAL] (1) column1.mediumLiquidIn.Theta = min(1.0, 0.00154526068547764 * column1.T_l_in); ($RES_BND_661)
(89|103)       [ARRY] (2) column1.propsVap[1].c[:] = column1.mediumVapour[1].c; ($RES_SIM_362)
(90|105)       [ARRY] (1) column1.propsVap.u = column1.u_v; ($RES_SIM_363)
(91|106)       [ARRY] (2) column1.bubblePressure[1].x_l = column1.x_l[1, :]; ($RES_SIM_1109)
(92|108)       [ARRY] (1) column1.propsVap.cp = column1.mediumVapour.cp; ($RES_SIM_364)
(93|109)       [ARRY] (2) column1.bubblePressure[1].x_v = column1.x_v_star[1, :]; ($RES_SIM_1108)
(94|111)       [ARRY] (1) column1.propsVap.lambda = column1.mediumVapour.lambda; ($RES_SIM_365)
(95|112)       [ARRY] (2) column1.bubblePressure[1].x_vap_liq = column1.x_vap_liq[1, :]; ($RES_SIM_1107)
(96|114)       [ARRY] (2) column1.propsVap[1].x[:] = column1.mediumVapour[1].x; ($RES_SIM_366)
(97|116)       [SCAL] (1) column1.bubblePressure[1].p = column1.p_hyd[1]; ($RES_SIM_1106)
(98|117)       [ARRY] (1) column1.propsVap.MM = column1.MM_v; ($RES_SIM_367)
(99|118)       [ARRY] (2) column1.bubblePressure[1].p_sat = column1.p_sat_bulk[1, :]; ($RES_SIM_1105)
(100|120)      [ARRY] (1) column1.propsVap.rho = column1.mediumVapour.d; ($RES_SIM_368)
(101|121)      [ARRY] (2) column1.bubblePressure[1].henryCoefficient.x_l = column1.bubblePressure[1].x_l; ($RES_SIM_1104)
(102|123)      [SCAL] (1) column1.bubblePressure[1].henryCoefficient.T = column1.bubblePressure[1].T; ($RES_SIM_1103)
(103|124)      [ARRY] (2) column1.bubblePressure[1].henryCoefficient.henryCoeff.x_l = column1.bubblePressure[1].henryCoefficient.x_l; ($RES_SIM_1102)
(104|126)      [ARRY] (1) column1.holdup.rho = column1.propsLiq.rho; ($RES_BND_922)
(105|127)      [SCAL] (1) column1.bubblePressure[1].henryCoefficient.henryCoeff.T = column1.bubblePressure[1].henryCoefficient.T; ($RES_SIM_1101)
(106|128)      [ARRY] (1) column1.holdup.sigma = column1.propsLiq.sigma; ($RES_BND_923)
(107|129)      [ARRY] (1) column1.holdup.eta = column1.propsLiq.eta; ($RES_BND_924)
(108|130)      [SCAL] (1) sourceLiquid.medium.properties.v = sourceLiquid.medium.v; ($RES_SIM_450)
(109|131)      [ARRY] (1) column1.balanceEquations.rho_v = column1.propsVap.d; ($RES_BND_750)
(110|132)      [SCAL] (1) sourceLiquid.medium.properties.MM = sourceLiquid.medium.MM; ($RES_SIM_451)
(111|133)      [SCAL] (1) column1.balanceEquations.rho_v_in = column1.propsVapIn.d; ($RES_BND_751)
(112|134)      [SCAL] (1) sourceLiquid.medium.properties.rho = sourceLiquid.rho; ($RES_SIM_452)
(113|135)      [ARRY] (1) column1.balanceEquations.MM_v = column1.propsVap.MM; ($RES_BND_752)
(114|136)      [ARRY] (1) column1.F = column1.Vdot_v / column1.geometry.A .* column1.rho_v .^ 0.5; ($RES_BND_928)
(115|137)      [SCAL] (1) column1.balanceEquations.MM_v_in = column1.propsVapIn.MM; ($RES_BND_753)
(116|138)      [ARRY] (1) column1.liquidLoad = column1.Vdot_l .* 3600.0 / column1.geometry.A; ($RES_BND_929)
(117|139)      [ARRY] (1) column1.balanceEquations.h_v = column1.propsVap.h; ($RES_BND_754)
(118|140)      [SCAL] (1) column1.propsLiqIn.T = column1.T_l_in; ($RES_SIM_280)
(119|141)      [SCAL] (1) column1.balanceEquations.h_v_in = column1.propsVapIn.h; ($RES_BND_755)
(120|142)      [SCAL] (1) sourceLiquid.medium.state.h = 55.55555555555556 * sourceLiquid.h; ($RES_SIM_456)
(121|143)      [SCAL] (1) column1.mediumLiquidIn.state.h = 55.55555555555556 * sourceLiquid.h; ($RES_SIM_281)
(122|144)      [SCAL] (1) sourceLiquid.medium.state.d = sourceLiquid.rho; ($RES_SIM_457)
(123|145)      [SCAL] (1) column1.mediumLiquidIn.state.d = column1.rho_l_in; ($RES_SIM_282)
(124|146)      [SCAL] (1) column1.mediumLiquidIn.state.T = column1.T_l_in; ($RES_SIM_283)
(125|147)      [ARRY] (1) column1.balanceEquations.rho_l = column1.propsLiq.d; ($RES_BND_758)
(126|148)      [SCAL] (1) sourceLiquid.medium.state.p = sourceLiquid.p; ($RES_SIM_459)
(127|149)      [SCAL] (1) column1.mediumLiquidIn.state.p = column1.mediumLiquidIn.p; ($RES_SIM_284)
(128|150)      [SCAL] (1) column1.balanceEquations.rho_l_in = column1.propsLiqIn.d; ($RES_BND_759)
(129|151)      [SCAL] (1) column1.mediumLiquidIn.cp = 2456.0 * column1.x_l_in[1] + 2206.0 * column1.x_l_in[2]; ($RES_SIM_285)
(130|152)      [SCAL] (1) column1.mediumLiquidIn.lambda = 0.108 * column1.x_l_in[1] + 0.122 * column1.x_l_in[2]; ($RES_SIM_286)
(131|153)      [SCAL] (1) column1.mediumLiquidIn.v = column1.MM_l_in / column1.rho_l_in; ($RES_SIM_287)
(132|154)      [FOR-] (2) ($RES_BND_840)
(132|154)      [----] for $i1 in 1:2 loop
(132|154)      [----]   [SCAL] (1) column1.balanceEquations.homotopyMethod.K[$i1] = 1.0; ($RES_BND_841)
(132|154)      [----] end for;
(133|156)      [SCAL] (1) column1.mediumLiquidIn.u = sourceLiquid.h - (column1.mediumLiquidIn.p / column1.rho_l_in) * column1.MM_l_in; ($RES_SIM_289)
(134|157)      [SCAL] (1) column1.eps_vap[1] = 1.0 - column1.eps_liq[1]; ($RES_SIM_9)
(135|158)      [ARRY] (1) column1.propsVap.T = column1.mediumVapour.T; ($RES_SIM_370)
(136|159)      [SCAL] (1) column1.n_mol_L[1] = column1.eps_liq[1] * column1.eps * column1.H * column1.A * $FUN_33; ($RES_SIM_6)
(137|160)      [ARRY] (2) column1.balanceEquations.stateVap[1].X[:] = column1.mediumVapour[1].X; ($RES_SIM_371)
(138|162)      [FOR-] (6) ($RES_BND_846)
(138|162)      [----] for $i1 in 1:6 loop
(138|162)      [----]   [SCAL] (1) column1.balanceEquations.homotopyMethod.whichHomotopy[$i1] = false; ($RES_BND_847)
(138|162)      [----] end for;
(139|168)      [SCAL] (1) column1.n_mol_V[1] = column1.eps_vap[1] * column1.eps * column1.H * column1.A * $FUN_34; ($RES_SIM_5)
(140|169)      [ARRY] (1) column1.balanceEquations.stateVap.p = column1.mediumVapour.p; ($RES_SIM_372)
(141|170)      [ARRY] (1) column1.balanceEquations.stateVap.T = column1.mediumVapour.T; ($RES_SIM_373)
(142|171)      [SCAL] (1) column1.balanceEquations.homotopyMethod.bool_Ndot_inter = column1.balanceEquations.homotopyMethod.whichHomotopy[1]; ($RES_BND_848)
(143|172)      [SCAL] (1) column1.mediumVapour[1].lambda = 0.018 * column1.mediumVapour[1].x[1] + 0.0144 * column1.mediumVapour[1].x[2]; ($RES_SIM_374)
(144|173)      [ARRY] (1) column1.rho_v = column1.mediumVapour.d; ($RES_BND_674)
(145|174)      [SCAL] (1) column1.balanceEquations.homotopyMethod.bool_Edot_inter = column1.balanceEquations.homotopyMethod.whichHomotopy[2]; ($RES_BND_849)
(146|175)      [SCAL] (1) column1.mediumVapour[1].cp = 1549.0 * column1.mediumVapour[1].x[1] + 2456.0 * column1.mediumVapour[1].x[2]; ($RES_SIM_376)
(147|176)      [SCAL] (1) column1.mediumVapour[1].d = column1.mediumVapour[1].p / (column1.mediumVapour[1].T * column1.mediumVapour[1].R); ($RES_SIM_377)
(148|177)      [SCAL] (1) column1.u_v[1] = column1.mediumVapour[1].h - column1.mediumVapour[1].T * column1.mediumVapour[1].R * column1.MM_v[1]; ($RES_SIM_378)
(149|178)      [ARRY] (1) column1.h_v = column1.mediumVapour.h; ($RES_BND_678)
(150|179)      [SCAL] (1) column1.mediumVapour[1].R = 8.31446261815324 / column1.MM_v[1]; ($RES_SIM_379)
(151|180)      [ARRY] (1) column1.balanceEquations.MM_l = column1.propsLiq.MM; ($RES_BND_760)
(152|181)      [SCAL] (1) column1.balanceEquations.MM_l_in = column1.propsLiqIn.MM; ($RES_BND_761)
(153|182)      [SCAL] (1) sourceLiquid.medium.v = sourceLiquid.medium.MM / sourceLiquid.rho; ($RES_SIM_462)
(154|183)      [ARRY] (1) column1.balanceEquations.h_l = column1.propsLiq.h; ($RES_BND_762)
(155|184)      [SCAL] (1) column1.balanceEquations.h_l_in = column1.propsLiqIn.h; ($RES_BND_763)
(156|185)      [SCAL] (1) sourceLiquid.medium.u = sourceLiquid.h - (sourceLiquid.p / sourceLiquid.rho) * sourceLiquid.medium.MM; ($RES_SIM_464)
(157|186)      [ARRY] (2) column1.balanceEquations.c_v_in = column1.propsVapIn.c; ($RES_BND_765)
(158|188)      [ARRY] (1) column1.results.T_l = column1.T; ($RES_BND_590)
(159|189)      [SCAL] (1) sourceLiquid.medium.R = 8.31446261815324 / sourceLiquid.medium.MM; ($RES_SIM_466)
(160|190)      [SCAL] (1) column1.mediumLiquidIn.R = 8.31446261815324 / column1.MM_l_in; ($RES_SIM_291)
(161|191)      [ARRY] (2) column1.balanceEquations.x_v_in = column1.propsVapIn.x; ($RES_BND_766)
(162|193)      [ARRY] (2) column1.results.x_l = column1.x_l; ($RES_BND_591)
(163|195)      [SCAL] (1) column1.mediumLiquidIn.eta = 1.1e-4 * column1.x_l_in[1] + 2.4e-4 * column1.x_l_in[2]; ($RES_SIM_292)
(164|196)      [ARRY] (2) column1.results.x_v = column1.x_v_star; ($RES_BND_592)
(165|198)      [SCAL] (1) sourceLiquid.rho = 528.0 * sourceLiquid.x[1] + 646.0 * sourceLiquid.x[2]; ($RES_SIM_468)
(166|199)      [SCAL] (1) column1.rho_l_in = 528.0 * column1.x_l_in[1] + 646.0 * column1.x_l_in[2]; ($RES_SIM_293)
(167|200)      [ARRY] (2) column1.results.x_l_star = column1.x_l_star; ($RES_BND_593)
(168|202)      [SCAL] (1) sourceLiquid.medium.MM = max(1e-5, $FUN_2); ($RES_SIM_469)
(169|203)      [SCAL] (1) column1.MM_l_in = max(1e-5, $FUN_12); ($RES_SIM_294)
(170|204)      [ARRY] (2) column1.results.x_v_star = column1.x_v_star; ($RES_BND_594)
(171|206)      [ARRY] (2) column1.results.c_l = column1.c_l; ($RES_BND_595)
(172|208)      [FOR-] (2) ($RES_SIM_296)
(172|208)      [----] for $i1 in 1:2 loop
(172|208)      [----]   [SCAL] (1) $FUN_11[$i1] = (({369.8, 469.6})[$i1] / column1.T_l_in) * (column1.mediumLiquidIn.C[$i1] * (1.0 - column1.T_l_in / ({369.8, 469.6})[$i1]) ^ 3.0 + column1.mediumLiquidIn.B[$i1] * (1.0 - column1.T_l_in / ({369.8, 469.6})[$i1]) ^ 1.5 + column1.mediumLiquidIn.A[$i1] * (1.0 - column1.T_l_in / ({369.8, 469.6})[$i1]) + column1.mediumLiquidIn.D[$i1] * (1.0 - column1.T_l_in / ({369.8, 469.6})[$i1]) ^ 6.0); ($RES_SIM_297)
(172|208)      [----] end for;
(173|210)      [ARRY] (2) column1.results.c_v = column1.c_v; ($RES_BND_596)
(174|212)      [ARRY] (2) column1.results.p_v = column1.p_hyd; ($RES_BND_597)
(175|214)      [SCAL] (1) sourceLiquid.h = (108.3096 * column1.x_l_in[1] + 159.1629 * column1.x_l_in[2]) * ((-1392.0) + column1.T_l_in); ($RES_SIM_298)
(176|215)      [ARRY] (1) column1.results.Vdot_l = column1.Vdot_l; ($RES_BND_598)
(177|216)      [ARRY] (1) column1.propsLiq.h = column1.mediumLiquid.h; ($RES_SIM_299)
(178|217)      [ARRY] (1) column1.results.Vdot_v = column1.Vdot_v; ($RES_BND_599)
(179|218)      [SCAL] (1) column1.balanceEquations.homotopyMethod.bool_dp = column1.balanceEquations.homotopyMethod.whichHomotopy[3]; ($RES_BND_850)
(180|219)      [SCAL] (1) column1.balanceEquations.homotopyMethod.bool_K = column1.balanceEquations.homotopyMethod.whichHomotopy[4]; ($RES_BND_851)
(181|220)      [SCAL] (1) column1.balanceEquations.homotopyMethod.bool_rho = column1.balanceEquations.homotopyMethod.whichHomotopy[5]; ($RES_BND_852)
(182|221)      [SCAL] (1) column1.balanceEquations.homotopyMethod.bool_h = column1.balanceEquations.homotopyMethod.whichHomotopy[6]; ($RES_BND_853)
(183|222)      [ARRY] (1) column1.initOption.p_v = column1.p_hyd[1:1]; ($RES_BND_854)
(184|223)      [SCAL] (1) column1.MM_v[1] = max(1e-6, $FUN_6); ($RES_SIM_381)
(185|224)      [ARRY] (1) column1.rho_l = column1.mediumLiquid.d; ($RES_BND_681)
(186|225)      [ARRY] (1) column1.mediumVapour.p_start = column1.balanceEquations.stateVap.p; ($RES_SIM_382)
(187|226)      [FOR-] (2) ($RES_SIM_383)
(187|226)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(187|226)      [----]   [SCAL] (1) column1.mediumVapour[$i1].X[$i2] = max(1e-5, min(1.0, (({0.0441, 0.07215})[$i2] * column1.mediumVapour[$i1].x[$i2]) / column1.MM_v[$i1])); ($RES_SIM_384)
(187|226)      [----] end for;
(188|228)      [SCAL] (1) column1.mediumVapour[1].calcSpecificEnthalpy.hX[2] = 0.07215 * (358000.0 + 2206.0 * ((-1406.0) + column1.mediumVapour[1].calcSpecificEnthalpy.T_sat[2]) + 1599.0 * (column1.mediumVapour[1].calcSpecificEnthalpy.T - column1.mediumVapour[1].calcSpecificEnthalpy.T_sat[2])); ($RES_SIM_385)
(189|229)      [SCAL] (1) column1.mediumVapour[1].calcSpecificEnthalpy.hX[1] = 0.0441 * (426000.0 + 2456.0 * ((-1406.0) + column1.mediumVapour[1].calcSpecificEnthalpy.T_sat[1]) + 1549.0 * (column1.mediumVapour[1].calcSpecificEnthalpy.T - column1.mediumVapour[1].calcSpecificEnthalpy.T_sat[1])); ($RES_SIM_386)
(190|230)      [SCAL] (1) column1.mediumVapour[1].h = $FUN_5; ($RES_SIM_387)
(191|231)      [SCAL] (1) column1.Ndot_v_in = ramp.offset + (if $TEV_0 then 0.0 else if $TEV_1 then (ramp.height * (time - ramp.startTime)) / ramp.duration else ramp.height); ($RES_SIM_388)
(192|232)      [FOR-] (2) ($RES_SIM_10)
(192|232)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(192|232)      [----]   [SCAL] (1) column1.x_v_star[$i1, $i2] = (column1.MM_v[$i1] * column1.c_v[$i1, $i2]) / column1.rho_v[$i1]; ($RES_SIM_11)
(192|232)      [----] end for;
(193|234)      [FOR-] (2) ($RES_SIM_12)
(193|234)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(193|234)      [----]   [SCAL] (1) column1.x_l[$i1, $i2] = column1.c_l[$i1, $i2] * column1.mediumLiquid[$i1].v; ($RES_SIM_13)
(193|234)      [----] end for;
(194|236)      [SCAL] (1) sourceLiquid.h = (108.3096 * sourceLiquid.x[1] + 159.1629 * sourceLiquid.x[2]) * ((-1392.0) + sourceLiquid.T); ($RES_SIM_473)
(195|237)      [ARRY] (2) column1.n_mol_L_i[1, :] = column1.n_mol_L[1] * column1.x_l[1, :]; ($RES_SIM_18)
(196|239)      [ARRY] (2) column1.n_mol_V_i[1, :] = column1.n_mol_V[1] * column1.x_v_star[1, :]; ($RES_SIM_19)
(197|241)      [ARRY] (2) column1.balanceEquations.c_l_in = column1.propsLiqIn.c; ($RES_BND_775)
(198|243)      [ARRY] (2) column1.balanceEquations.x_l_in = column1.propsLiqIn.x; ($RES_BND_776)
(199|245)      [ARRY] (2) column1.balanceEquations.x_l = column1.propsLiq.x; ($RES_BND_777)
(200|247)      [SCAL] (1) column1.internalFeedPort[1].feedVapourInternal.Ndot = 0.0; ($RES_SIM_479)
(201|248)      [ARRY] (1) column1.bubblePressure.v_v = column1.MM_v / column1.rho_v; ($RES_BND_690)
(202|249)      [ARRY] (2) sourceGas_Vdot.c = (sourceGas_Vdot.x / sourceGas_Vdot.MM) .* sourceGas_Vdot.rho; ($RES_SIM_392)
(203|251)      [SCAL] (1) sourceGas_Vdot.medium.properties.h = column1.h_v_in; ($RES_SIM_397)
(204|252)      [SCAL] (1) sourceGas_Vdot.medium.properties.d = sourceGas_Vdot.rho; ($RES_SIM_398)
(205|253)      [ARRY] (2) sourceGas_Vdot.medium.properties.c = sourceGas_Vdot.c; ($RES_SIM_399)
(206|255)      [SCAL] (1) column1.Vdot_l_in = (column1.Ndot_l_in / column1.rho_l_in) * column1.MM_l_in; ($RES_SIM_20)
(207|256)      [ARRY] (1) column1.Vdot_l = column1.Ndot_l / column1.rho_l .* column1.MM_l; ($RES_SIM_21)
(208|257)      [SCAL] (1) column1.Vdot_v_in = (column1.Ndot_v_in / column1.rho_v_in) * column1.MM_v_in; ($RES_SIM_22)
(209|258)      [ARRY] (1) column1.Vdot_v = column1.Ndot_v / column1.rho_v .* column1.MM_v; ($RES_SIM_23)
(210|259)      [SCAL] (1) column1.internalFeedPort[1].feedLiquidInternal.Ndot = 0.0; ($RES_SIM_480)
(211|260)      [SCAL] (1) -sinkLiquid.liquidPortIn.Ndot = -column1.Ndot_l[1]; ($RES_SIM_24)
(212|261)      [SCAL] (1) sourceLiquid.p = column1.p_hyd[2]; ($RES_SIM_26)
(213|262)      [SCAL] (1) -sinkGas.gasPortIn.Ndot = -column1.Ndot_v[1]; ($RES_SIM_27)
(214|263)      [SCAL] (1) sinkGas.p = column1.p_hyd[2]; ($RES_SIM_29)
(215|264)      [ARRY] (1) column1.initOption.v = column1.propsLiq.v; ($RES_BND_874)
(216|265)      [ARRY] (1) column1.reaction.V = column1.H * column1.A * column1.eps_liq; ($RES_BND_876)
(217|266)      [SCAL] (1) column1.downStreamIn.h_outflow = column1.h_l[1]; ($RES_SIM_31)
(218|267)      [ARRY] (2) column1.downStreamIn.x_outflow = column1.x_l[1, :]; ($RES_SIM_32)
(219|269)      [SCAL] (1) column1.h_downStreamOut_act = column1.h_l[1]; ($RES_SIM_33)
(220|270)      [ARRY] (2) column1.downStreamOut.x_outflow = column1.x_l[1, :]; ($RES_SIM_34)
(221|272)      [SCAL] (1) column1.upStreamIn.h_outflow = column1.h_v[1]; ($RES_SIM_35)
(222|273)      [ARRY] (2) column1.upStreamIn.x_outflow = column1.x_v_star[1, :]; ($RES_SIM_36)
(223|275)      [SCAL] (1) column1.upStreamOut.h_outflow = column1.h_v[1]; ($RES_SIM_37)
(224|276)      [ARRY] (2) column1.upStreamOut.x_outflow = column1.x_v_star[1, :]; ($RES_SIM_38)
(225|278)      [ARRY] (1) column1.sum_xl = sum(column1.x_l[:, i] for i in 1:2); ($RES_AUX_975)
(226|279)      [ARRY] (1) column1.sum_xv = sum(column1.x_v_star[:, i] for i in 1:2); ($RES_AUX_976)
(227|280)      [ARRY] (1) column1.balanceEquations.T = column1.propsLiq.T; ($RES_BND_798)
(228|281)      [SCAL] (1) column1.V_liq = sum(column1.eps_liq .* column1.H * column1.A * column1.eps); ($RES_AUX_979)
(229|282)      [SCAL] (1) column1.h_upStreamOut_act = if $SEV_6 then 1.0 else column1.upStreamOut.h_outflow; ($RES_SIM_41)
(230|283)      [SCAL] (1) column1.h_upStreamIn_act = if $SEV_5 then column1.h_v_in else column1.upStreamIn.h_outflow; ($RES_SIM_42)
(231|284)      [ARRY] (1) column1.n_liq = sum(column1.n_i_liq[:, i] for i in 1:2); ($RES_AUX_980)
(232|285)      [ARRY] (1) column1.n_vap = sum(column1.n_i_vap[:, i] for i in 1:2); ($RES_AUX_981)
(233|286)      [SCAL] (1) -column1.heatPort.Qdot = sum(column1.heatTransferWall.Qdot_out); ($RES_AUX_983)
(234|287)      [ARRY] (2) column1.c_l_in = (column1.x_l_in / column1.MM_l_in) .* column1.rho_l_in; ($RES_SIM_47)
(235|289)      [SCAL] (1) $FUN_34 = sum(column1.c_v[1, :]); ($RES_AUX_984)
(236|290)      [SCAL] (1) $FUN_33 = sum(column1.c_l[1, :]); ($RES_AUX_985)
(237|291)      [SCAL] (1) $FUN_32 = sum(column1.n_tot[1, :]); ($RES_AUX_986)
(238|292)      [SCAL] (1) $FUN_31 = sum(column1.c_v_feed[1, :]); ($RES_AUX_987)
(239|293)      [SCAL] (1) $FUN_30 = sum(column1.c_l_feed[1, :]); ($RES_AUX_988)
(240|294)      [SCAL] (1) $FUN_29 = exp(-(2.157463 * column1.holdup.rho[1]) / max(1e-7, column1.holdup.geometry.a * column1.holdup.sigma[1] * column1.holdup.geometry.a)); ($RES_AUX_989)
(241|295)      [ARRY] (2) column1.c_v_in = (column1.x_v_in / column1.MM_v_in) .* column1.rho_v_in; ($RES_SIM_50)
(242|297)      [ARRY] (2) column1.X_v[1, :] = {0.0441, 0.07215} .* (column1.x_v_star[1, :] / column1.MM_v[1]); ($RES_SIM_53)
(243|299)      [SCAL] (1) $FUN_28 = sign(column1.pressureLoss.hu_dyn[1]); ($RES_AUX_990)
(244|300)      [ARRY] (2) column1.X_l[1, :] = {0.0441, 0.07215} .* (column1.x_l[1, :] / column1.MM_l[1]); ($RES_SIM_54)
(245|302)      [SCAL] (1) $FUN_27 = abs(column1.pressureLoss.hu_dyn[1]); ($RES_AUX_991)
(246|303)      [ARRY] (2) column1.n_i_liq[1, :] = (column1.H .* column1.A .* column1.eps_liq[1] .* column1.c_l[1, :]) .* column1.eps; ($RES_SIM_55)
(247|305)      [SCAL] (1) $FUN_26 = sum(column1.c_l_feed[1, :]); ($RES_AUX_992)
(248|306)      [ARRY] (2) column1.n_i_vap[1, :] = (column1.H .* column1.A .* column1.eps_vap[1] .* column1.c_v[1, :]) .* column1.eps; ($RES_SIM_56)
(249|308)      [SCAL] (1) $FUN_25 = sum(column1.c_v_feed[1, :]); ($RES_AUX_993)
(250|309)      [SCAL] (1) 1.0 = sum(column1.balanceEquations.x_l[1, :]); ($RES_AUX_994)
(251|310)      [ARRY] (2) column1.x_vap_liq[1, :] = column1.n_tot[1, :] / $FUN_32; ($RES_SIM_58)
(252|312)      [SCAL] (1) 1.0 = sum(column1.x_v_star[1, :]); ($RES_AUX_995)
(253|313)      [FOR-] (2) ($RES_SIM_59)
(253|313)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(253|313)      [----]   [SCAL] (1) column1.n_tot[$i1, $i2] = column1.H * (column1.c_v[$i1, ({{1, 1}, {2, 2}})[$i2, 1]] * (1.0 - column1.eps_liq[$i1]) + column1.c_l[$i1, ({{1, 1}, {2, 2}})[$i2, 2]] * column1.eps_liq[$i1]) * column1.A; ($RES_SIM_60)
(253|313)      [----] end for;
(254|315)      [SCAL] (1) $FUN_22 = sum(column1.c_v[1, :]); ($RES_AUX_996)
(255|316)      [SCAL] (1) $FUN_21 = sum(column1.c_l[1, :]); ($RES_AUX_997)
(256|317)      [SCAL] (1) $FUN_20 = sum(column1.balanceEquations.thermoEquilibrium[1].p_partial); ($RES_AUX_998)
(257|318)      [FOR-] (2) ($RES_AUX_999)
(257|318)      [----] for {$i2 in 1:2, $i1 in 1:1} loop
(257|318)      [----]   [SCAL] (1) $FUN_19[$i2, $i1] = exp(column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.henry_C[$i2] * (1/column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.T - 1/column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.henry_T_norm)); ($RES_AUX_1000)
(257|318)      [----] end for;
(258|320)      [SCAL] (1) $SEV_6 = (-sinkGas.gasPortIn.Ndot) > 0.0; ($RES_EVT_1400)
(259|321)      [SCAL] (1) column1.balanceEquations.hu_stat[1] = 0.033 * $FUN_29; ($RES_SIM_100)
(260|322)      [SCAL] (1) column1.holdup.Vdot[1] = max(0.0, 50.0 * column1.eps_liq[1] * column1.holdup.geometry.eps * column1.holdup.geometry.A * (($FUN_27 / 0.555) ^ 3.0 * 9.80665 * column1.holdup.geometry.eps ^ 4.65 / column1.holdup.geometry.a) ^ 0.5 * CAST(Real, $FUN_28)); ($RES_SIM_101)
(261|323)      [SCAL] (1) column1.p_v_in - column1.p_hyd[1] = 0.5 * column1.pressureLoss.K * column1.Vdot_v_in; ($RES_SIM_102)
(262|324)      [SCAL] (1) max(0.0, CAST(Real, sign(column1.p_hyd[1] - column1.p_hyd[2])) * max(0.0, abs(column1.p_hyd[1] - column1.p_hyd[2]))) = (0.5 * column1.pressureLoss.K) * column1.Vdot_v[1]; ($RES_SIM_105)
(263|325)      [SCAL] (1) column1.Qdot_reac[1] = 0.0; ($RES_SIM_108)
(264|326)      [ARRY] (2) column1.c_v_feed_used[1, :] = column1.c_v_feed[1, :]; ($RES_SIM_61)
(265|328)      [SCAL] (1) column1.internalFeedPort[1].feedVapourInternal.h_outflow = column1.h_v[1]; ($RES_SIM_62)
(266|329)      [ARRY] (2) column1.internalFeedPort[1].feedVapourInternal.x_outflow = column1.x_v_star[1, :]; ($RES_SIM_63)
(267|331)      [SCAL] (1) column1.mediumVapourLink[1].mediumConIn.rho = 1.0; ($RES_SIM_65)
(268|332)      [SCAL] (1) column1.mediumVapourLink[1].mediumConIn.MM = 1.0; ($RES_SIM_66)
(269|333)      [SCAL] (1) column1.internalFeedPort[1].feedVapourInternal.p = column1.p_hyd[1]; ($RES_SIM_67)
(270|334)      [SCAL] (1) column1.internalFeedPort[1].feedLiquidInternal.h_outflow = column1.h_l[1]; ($RES_SIM_68)
(271|335)      [FOR-] (2) ($RES_SIM_950)
(271|335)      [----] for $f3 in 1:2 loop
(271|335)      [----]   [SCAL] (1) sinkGas.gasPortIn.x_outflow[$f3] = 1.0; ($RES_SIM_951)
(271|335)      [----] end for;
(272|337)      [ARRY] (2) column1.internalFeedPort[1].feedLiquidInternal.x_outflow = column1.x_l[1, :]; ($RES_SIM_69)
(273|339)      [FOR-] (2) ($RES_SIM_952)
(273|339)      [----] for $f4 in 1:2 loop
(273|339)      [----]   [SCAL] (1) column1.mediumLiquid[1].eta_comp[$f4] = column1.mediumLiquid[1].eta; ($RES_SIM_953)
(273|339)      [----] end for;
(274|341)      [FOR-] (2) ($RES_SIM_954)
(274|341)      [----] for $f5 in 1:2 loop
(274|341)      [----]   [SCAL] (1) column1.mediumLiquidIn.eta_comp[$f5] = column1.mediumLiquidIn.eta; ($RES_SIM_955)
(274|341)      [----] end for;
(275|343)      [RECD] (18) column1.reaction[1].propsLiq = column1.propsLiq[1]; ($RES_SIM_1188)
(276|361)      [ARRY] (2) column1.balanceEquations.c_l_star[1, :] = column1.c_l[1, :]; ($RES_SIM_110)
(277|363)      [ARRY] (2) column1.balanceEquations.c_v_star[1, :] = column1.c_v[1, :]; ($RES_SIM_111)
(278|365)      [ARRY] (2) column1.balanceEquations.p_sat[1, :] = column1.p_sat_bulk[1, :]; ($RES_SIM_112)
(279|367)      [SCAL] (1) column1.bool_eps[1] = $SEV_4; ($RES_SIM_113)
(280|368)      [SCAL] (1) 0.0 = ($FUN_25 * column1.Vdot_v_feed[1] * column1.h_v_feed[1] + column1.Qdot_reac[1] + $FUN_26 * column1.Vdot_l_feed[1] * column1.h_l_feed[1] + column1.Ndot_v_in * column1.h_upStreamIn_act + column1.Ndot_l_in * sourceLiquid.h + column1.Ndot_source_startUp[1] * column1.balanceEquations.h_v[1]) - ($FUN_21 * column1.Vdot_le[1] * column1.balanceEquations.h_l[1] + column1.Qdot_wall[1] + column1.Ndot_v[1] * column1.h_upStreamOut_act + column1.Ndot_l[1] * column1.h_downStreamOut_act); ($RES_SIM_115)
(281|369)      [ARRY] (2) {0.0 for $f12 in 1:2} = (column1.Ndot_reac[1, :] + column1.Ndot_source_startUp[1] * column1.x_v_star[1, :] + column1.Vdot_l_feed[1] * column1.c_l_feed[1, :] + column1.Vdot_v_feed[1] * column1.c_v_feed[1, :] + column1.Ndot_l_in * column1.x_downStreamIn_act + column1.Ndot_v_in * column1.x_upStreamIn_act) - (column1.Vdot_le[1] * column1.c_l[1, :] + column1.Ndot_l[1] * column1.x_downStreamOut_act + column1.Ndot_v[1] * column1.x_upStreamOut_act); ($RES_SIM_119)
(282|371)      [SCAL] (1) column1.mediumLink[1].mediumConIn.rho = 1.0; ($RES_SIM_71)
(283|372)      [SCAL] (1) column1.mediumLink[1].mediumConIn.MM = 1.0; ($RES_SIM_72)
(284|373)      [SCAL] (1) column1.internalFeedPort[1].feedLiquidInternal.p = column1.p_hyd[1]; ($RES_SIM_73)
(285|374)      [SCAL] (1) column1.Vdot_v_feed[1] = column1.internalFeedPort[1].feedVapourInternal.Ndot / $FUN_31; ($RES_SIM_74)
(286|375)      [SCAL] (1) column1.rho_v_feed[1] = column1.mediumVapourLink[1].mediumConIn.rho; ($RES_SIM_75)
(287|376)      [SCAL] (1) column1.MM_v_feed[1] = column1.mediumVapourLink[1].mediumConIn.MM; ($RES_SIM_76)
(288|377)      [SCAL] (1) column1.Vdot_l_feed[1] = column1.internalFeedPort[1].feedLiquidInternal.Ndot / $FUN_30; ($RES_SIM_77)
(289|378)      [SCAL] (1) column1.rho_l_feed[1] = column1.mediumLink[1].mediumConIn.rho; ($RES_SIM_78)
(290|379)      [SCAL] (1) column1.MM_l_feed[1] = column1.mediumLink[1].mediumConIn.MM; ($RES_SIM_79)
(291|380)      [FOR-] (2) ($RES_SIM_966)
(291|380)      [----] for $f14 in 1:2 loop
(291|380)      [----]   [SCAL] (1) column1.Ndot_reac[1, $f14] = 0.0; ($RES_SIM_967)
(291|380)      [----] end for;
(292|382)      [SCAL] (1) column1.heatTransferWall.Qdot_out[1] = 0.0; ($RES_SIM_968)
(293|383)      [SCAL] (1) column1.Vdot_le[1] = 0.0; ($RES_SIM_969)
(294|384)      [SCAL] (1) column1.balanceEquations.n_mol_V[1] = column1.eps_vap[1] * column1.balanceEquations.eps * column1.balanceEquations.H * column1.balanceEquations.A * $FUN_22; ($RES_SIM_120)
(295|385)      [SCAL] (1) column1.balanceEquations.n_mol_L[1] = column1.eps_liq[1] * column1.balanceEquations.eps * column1.balanceEquations.H * column1.balanceEquations.A * $FUN_21; ($RES_SIM_121)
(296|386)      [FOR-] (2) ($RES_SIM_122)
(296|386)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(296|386)      [----]   [SCAL] (1) column1.balanceEquations.n_mol[$i1, $i2] = column1.eps_vap[$i1] * column1.balanceEquations.eps * column1.balanceEquations.H * column1.balanceEquations.A * column1.c_v[$i1, ({{1, 1}, {2, 2}})[$i2, 1]] + column1.eps_liq[$i1] * column1.balanceEquations.eps * column1.balanceEquations.H * column1.balanceEquations.A * column1.c_l[$i1, ({{1, 1}, {2, 2}})[$i2, 1]]; ($RES_SIM_123)
(296|386)      [----] end for;
(297|388)      [SCAL] (1) sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat[2] = 309.2; ($RES_SIM_1175)
(298|389)      [SCAL] (1) sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat[1] = 231.1; ($RES_SIM_1174)
(299|390)      [FOR-] (2) ($RES_SIM_125)
(299|390)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(299|390)      [----]   [SCAL] (1) column1.x_v_star[$i1, ({{1, 1}, {2, 2}})[$i2, 1]] = column1.balanceEquations.K[$i1, $i2] * column1.x_l_star[$i1, ({{1, 1}, {2, 2}})[$i2, 2]]; ($RES_SIM_126)
(299|390)      [----] end for;
(300|392)      [SCAL] (1) column1.mediumVapour[1].p = column1.p_hyd[1]; ($RES_SIM_1171)
(301|393)      [FOR-] (2) ($RES_SIM_127)
(301|393)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(301|393)      [----]   [SCAL] (1) column1.x_l_star[$i1, $i2] = column1.balanceEquations.x_l[$i1, $i2]; ($RES_SIM_128)
(301|393)      [----] end for;
(302|395)      [ARRY] (2) column1.mediumVapour[1].x = column1.x_v_star[1, :]; ($RES_SIM_1170)
(303|397)      [SCAL] (1) column1.balanceEquations.K[j, i] = if column1.balanceEquations.homotopyMethod.bool_K then homotopy(column1.balanceEquations.thermoEquilibrium[j].K[i], column1.balanceEquations.homotopyMethod.K[i]) else column1.balanceEquations.thermoEquilibrium[j].K[i]; ($RES_SIM_129)
(304|398)      [SCAL] (1) column1.h_l_feed[1] = column1.internalFeedPort[1].feedLiquidInternal.h_outflow; ($RES_SIM_80)
(305|399)      [ARRY] (2) column1.c_l_feed[1, :] = (column1.rho_l_feed[1] .* {column1.internalFeedPort[1].feedLiquidInternal.x_outflow[1], column1.internalFeedPort[1].feedLiquidInternal.x_outflow[2]}) / column1.MM_l_feed[1]; ($RES_SIM_81)
(306|401)      [SCAL] (1) column1.h_v_feed[1] = column1.internalFeedPort[1].feedVapourInternal.h_outflow; ($RES_SIM_82)
(307|402)      [ARRY] (2) column1.c_v_feed[1, :] = (column1.rho_v_feed[1] .* {column1.internalFeedPort[1].feedVapourInternal.x_outflow[1], column1.internalFeedPort[1].feedVapourInternal.x_outflow[2]}) / column1.MM_v_feed[1]; ($RES_SIM_83)
(308|404)      [SCAL] (1) column1.Ndot_source_startUp[1] = 0.0; ($RES_SIM_972)
(309|405)      [FOR-] (2) ($RES_SIM_973)
(309|405)      [----] for i in 1:2 loop
(309|405)      [----]   [SCAL] (1) column1.p_sat_bulk[1, i] = column1.mediumLiquid[1].p_sat[i]; ($RES_SIM_974)
(309|405)      [----] end for;
(310|407)      [ARRY] (1) column1.propsLiq.d = column1.mediumLiquid.d; ($RES_SIM_300)
(311|408)      [ARRY] (1) column1.results.eps_liq = column1.eps_liq; ($RES_BND_600)
(312|409)      [ARRY] (2) column1.propsLiq[1].c[:] = column1.mediumLiquid[1].c; ($RES_SIM_301)
(313|411)      [ARRY] (2) column1.propsLiq[1].eta_comp[:] = column1.mediumLiquid[1].eta_comp; ($RES_SIM_302)
(314|413)      [ARRY] (1) column1.propsLiq.p = column1.mediumLiquid.p; ($RES_SIM_303)
(315|414)      [SCAL] (1) column1.mediumVapour[1].T = column1.T_v[1]; ($RES_SIM_1169)
(316|415)      [ARRY] (1) column1.propsLiq.u = column1.u_l; ($RES_SIM_304)
(317|416)      [ARRY] (2) column1.mediumVapour[1].c = column1.c_v[1, :]; ($RES_SIM_1168)
(318|418)      [ARRY] (1) column1.propsLiq.cp = column1.mediumLiquid.cp; ($RES_SIM_305)
(319|419)      [ARRY] (2) column1.mediumVapour[1].x_star = column1.x_v_star[1, :]; ($RES_SIM_1167)
(320|421)      [ARRY] (1) column1.propsLiq.lambda = column1.mediumLiquid.lambda; ($RES_SIM_306)
(321|422)      [SCAL] (1) column1.mediumVapour[1].mediumConOut.h = column1.mediumVapour[1].h; ($RES_SIM_1166)
(322|423)      [ARRY] (2) column1.propsLiq[1].x[:] = column1.mediumLiquid[1].x; ($RES_SIM_307)
(323|425)      [SCAL] (1) column1.mediumVapour[1].mediumConOut.rho = column1.mediumVapour[1].d; ($RES_SIM_1165)
(324|426)      [ARRY] (1) column1.propsLiq.v = column1.mediumLiquid.v; ($RES_SIM_308)
(325|427)      [SCAL] (1) column1.mediumVapour[1].mediumConOut.MM = column1.MM_v[1]; ($RES_SIM_1164)
(326|428)      [ARRY] (1) column1.propsLiq.MM = column1.MM_l; ($RES_SIM_309)
(327|429)      [SCAL] (1) column1.mediumVapour[1].calcSpecificEnthalpy.T = column1.mediumVapour[1].T; ($RES_SIM_1163)
(328|430)      [SCAL] (1) column1.propsVap[1].T = column1.T_star[1]; ($RES_SIM_135)
(329|431)      [ARRY] (2) column1.mediumVapour[1].calcSpecificEnthalpy.x = column1.mediumVapour[1].x; ($RES_SIM_1162)
(330|433)      [SCAL] (1) column1.propsLiq[1].T = column1.T_star[1]; ($RES_SIM_136)
(331|434)      [SCAL] (1) column1.mediumVapour[1].calcSpecificEnthalpy.p = column1.mediumVapour[1].p; ($RES_SIM_1161)
(332|435)      [SCAL] (1) column1.mediumVapour[1].calcSpecificEnthalpy.T_sat[2] = 309.2; ($RES_SIM_1160)
(333|436)      [SCAL] (1) column1.Vdot_l[1] = max(0.0, column1.holdup.Vdot[1]); ($RES_SIM_95)
(334|437)      [SCAL] (1) column1.Qdot_wall[1] = 6.283185307179586 * column1.heatTransferWall.geometry.r1 * column1.heatTransferWall.geometry.H * column1.heatTransferWall.alpha_inner * (column1.T[1] - column1.heatTransferWall.T_wall[1]); ($RES_SIM_96)
(335|438)      [SCAL] (1) 0.0 = column1.Qdot_wall[1] - column1.heatTransferWall.Qdot_out[1]; ($RES_SIM_97)
(336|439)      [SCAL] (1) column1.eps_liq[1] = column1.balanceEquations.hu_stat[1] + column1.pressureLoss.hu_dyn[1]; ($RES_SIM_99)
(337|440)      [ARRY] (1) column1.propsLiq.rho = column1.mediumLiquid.d; ($RES_SIM_310)
(338|441)      [ARRY] (1) column1.propsLiq.sigma = column1.mediumLiquid.sigma; ($RES_SIM_311)
(339|442)      [ARRY] (1) column1.propsLiq.eta = column1.mediumLiquid.eta; ($RES_SIM_312)
(340|443)      [ARRY] (1) column1.propsLiq.T = column1.mediumLiquid.T; ($RES_SIM_313)
(341|444)      [SCAL] (1) column1.mediumVapour[1].calcSpecificEnthalpy.T_sat[1] = 231.1; ($RES_SIM_1159)
(342|445)      [SCAL] (1) column1.balanceEquations.stateLiq[1].h = 55.55555555555556 * column1.mediumLiquid[1].h; ($RES_SIM_314)
(343|446)      [ARRY] (1) column1.balanceEquations.stateLiq.d = column1.mediumLiquid.d; ($RES_SIM_315)
(344|447)      [SCAL] (1) column1.balanceEquations.thermoEquilibrium[1].p_bubble = $FUN_20; ($RES_SIM_140)
(345|448)      [ARRY] (1) column1.balanceEquations.stateLiq.T = column1.mediumLiquid.T; ($RES_SIM_316)
(346|449)      [FOR-] (2) ($RES_SIM_141)
(346|449)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(346|449)      [----]   [SCAL] (1) column1.balanceEquations.thermoEquilibrium[$i1].K[$i2] = (column1.balanceEquations.thermoEquilibrium[$i1].p_sat[({{1, 1}, {2, 2}})[$i2, 2]] * column1.balanceEquations.thermoEquilibrium[$i1].factor_K[({{1, 1}, {2, 2}})[$i2, 2]]) / column1.balanceEquations.thermoEquilibrium[$i1].p; ($RES_SIM_142)
(346|449)      [----] end for;
(347|451)      [SCAL] (1) column1.mediumVapourIn.calcSpecificEnthalpy.T_sat[2] = 309.2; ($RES_SIM_1156)
(348|452)      [ARRY] (1) column1.balanceEquations.stateLiq.p = column1.mediumLiquid.p; ($RES_SIM_317)
(349|453)      [SCAL] (1) column1.mediumVapourIn.calcSpecificEnthalpy.T_sat[1] = 231.1; ($RES_SIM_1155)
(350|454)      [SCAL] (1) column1.mediumLiquid[1].cp = 2456.0 * column1.mediumLiquid[1].x[1] + 2206.0 * column1.mediumLiquid[1].x[2]; ($RES_SIM_318)
(351|455)      [FOR-] (2) ($RES_SIM_143)
(351|455)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(351|455)      [----]   [SCAL] (1) column1.balanceEquations.thermoEquilibrium[$i1].p_partial[$i2] = column1.balanceEquations.thermoEquilibrium[$i1].factor_K[({{1, 1}, {2, 2}})[$i2, 2]] * column1.balanceEquations.thermoEquilibrium[$i1].x_l[({{1, 1}, {2, 2}})[$i2, 2]] * column1.balanceEquations.thermoEquilibrium[$i1].p_sat[({{1, 1}, {2, 2}})[$i2, 2]]; ($RES_SIM_144)
(351|455)      [----] end for;
(352|457)      [SCAL] (1) column1.mediumLiquid[1].lambda = 0.108 * column1.mediumLiquid[1].x[1] + 0.122 * column1.mediumLiquid[1].x[2]; ($RES_SIM_319)
(353|458)      [SCAL] (1) column1.mediumLiquid[1].p = column1.p_hyd[1]; ($RES_SIM_1152)
(354|459)      [SCAL] (1) column1.mediumLiquid[1].T = column1.T[1]; ($RES_SIM_1151)
(355|460)      [FOR-] (2) ($RES_SIM_147)
(355|460)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(355|460)      [----]   [SCAL] (1) column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.He[$i2] = if column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.henry_temp then column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.henry_H[$i2] * $FUN_19[$i2, $i1] else column1.balanceEquations.thermoEquilibrium[$i1].henryCoefficient.henryCoeff.henry_H[$i2]; ($RES_SIM_148)
(355|460)      [----] end for;
(356|462)      [ARRY] (2) column1.mediumLiquid[1].x = column1.x_l[1, :]; ($RES_SIM_1150)
(357|464)      [SCAL] (1) sourceGas_Vdot.medium.properties.u = sourceGas_Vdot.medium.u; ($RES_SIM_400)
(358|465)      [SCAL] (1) sourceGas_Vdot.medium.properties.MM = sourceGas_Vdot.MM; ($RES_SIM_404)
(359|466)      [SCAL] (1) sourceGas_Vdot.medium.properties.rho = sourceGas_Vdot.rho; ($RES_SIM_405)
(360|467)      [ARRY] (1) column1.mdot_v = column1.Vdot_v .* column1.rho_v; ($RES_BND_706)
(361|468)      [ARRY] (1) column1.mdot_l = column1.Vdot_l .* column1.rho_l; ($RES_BND_707)
(362|469)      [ARRY] (2) sourceGas_Vdot.medium.state.X = sourceGas_Vdot.medium.X; ($RES_SIM_408)
(363|471)      [SCAL] (1) sourceGas_Vdot.medium.state.p = column1.p_v_in; ($RES_SIM_409)
(364|472)      [ARRY] (2) sourceLiquid.medium.c = sourceLiquid.rho * sourceLiquid.x / sourceLiquid.medium.MM; ($RES_BND_536)
(365|474)      [SCAL] (1) column1.mediumLiquid[1].v = column1.MM_l[1] / column1.mediumLiquid[1].d; ($RES_SIM_320)
(366|475)      [SCAL] (1) column1.u_l[1] = column1.mediumLiquid[1].h - (column1.mediumLiquid[1].p / column1.mediumLiquid[1].d) * column1.MM_l[1]; ($RES_SIM_322)
(367|476)      [SCAL] (1) $FUN_2 = sum({0.0441, 0.07215} .* sourceLiquid.x); ($RES_AUX_1021)
(368|477)      [ARRY] (2) column1.mediumLiquid[1].c = column1.mediumLiquid[1].x .* column1.mediumLiquid[1].d / column1.MM_l[1]; ($RES_SIM_1149)
(369|479)      [SCAL] (1) column1.mediumLiquid[1].R = 8.31446261815324 / column1.MM_l[1]; ($RES_SIM_324)
(370|480)      [SCAL] (1) column1.h_v_in = sum(sourceGas_Vdot.x .* sourceGas_Vdot.medium.calcSpecificEnthalpy.hX); ($RES_AUX_1020)
(371|481)      [SCAL] (1) column1.mediumLiquid[1].h = column1.h_l[1]; ($RES_SIM_1148)
(372|482)      [SCAL] (1) column1.mediumLiquid[1].eta = 1.1e-4 * column1.mediumLiquid[1].x[1] + 2.4e-4 * column1.mediumLiquid[1].x[2]; ($RES_SIM_325)
(373|483)      [SCAL] (1) column1.mediumLiquid[1].mediumConOut.h = column1.mediumLiquid[1].h; ($RES_SIM_1147)
(374|484)      [SCAL] (1) column1.mediumLiquid[1].d = 528.0 * column1.mediumLiquid[1].x[1] + 646.0 * column1.mediumLiquid[1].x[2]; ($RES_SIM_326)
(375|485)      [SCAL] (1) column1.mediumLiquid[1].mediumConOut.rho = column1.mediumLiquid[1].d; ($RES_SIM_1146)
(376|486)      [SCAL] (1) column1.MM_l[1] = max(1e-5, $FUN_10); ($RES_SIM_327)
(377|487)      [SCAL] (1) column1.mediumLiquid[1].mediumConOut.MM = column1.MM_l[1]; ($RES_SIM_1145)
(378|488)      [SCAL] (1) column1.mediumLiquid[1].sigma = 0.016; ($RES_SIM_328)
(379|489)      [SCAL] (1) column1.mediumLiquid[1].calcSpecificEnthalpy.p = column1.mediumLiquid[1].p; ($RES_SIM_1144)
(380|490)      [FOR-] (2) ($RES_SIM_329)
(380|490)      [----] for {$i1 in 1:1, $i2 in 1:2} loop
(380|490)      [----]   [SCAL] (1) $FUN_9[$i2, $i1] = (({369.8, 469.6})[$i2] / column1.mediumLiquid[$i1].T) * (column1.mediumLiquid[$i1].C[$i2] * (1.0 - column1.mediumLiquid[$i1].T / ({369.8, 469.6})[$i2]) ^ 3.0 + column1.mediumLiquid[$i1].B[$i2] * (1.0 - column1.mediumLiquid[$i1].T / ({369.8, 469.6})[$i2]) ^ 1.5 + column1.mediumLiquid[$i1].A[$i2] * (1.0 - column1.mediumLiquid[$i1].T / ({369.8, 469.6})[$i2]) + column1.mediumLiquid[$i1].D[$i2] * (1.0 - column1.mediumLiquid[$i1].T / ({369.8, 469.6})[$i2]) ^ 6.0); ($RES_SIM_330)
(380|490)      [----] end for;
(381|492)      [SCAL] (1) column1.mediumLiquid[1].calcSpecificEnthalpy.T = column1.mediumLiquid[1].T; ($RES_SIM_1143)
(382|493)      [ARRY] (2) column1.mediumLiquid[1].calcSpecificEnthalpy.x = column1.mediumLiquid[1].x; ($RES_SIM_1142)
(383|495)      [ARRY] (1) column1.n_total = column1.n_liq + column1.n_vap; ($RES_BND_711)
(384|496)      [SCAL] (1) sourceGas_Vdot.rho = column1.p_v_in / (sourceGas_Vdot.T * sourceGas_Vdot.medium.R); ($RES_SIM_414)
(385|497)      [SCAL] (1) column1.balanceEquations.thermoEquilibrium[1].T = column1.T_star[1]; ($RES_SIM_1058)
(386|498)      [SCAL] (1) sourceGas_Vdot.medium.u = column1.h_v_in - sourceGas_Vdot.T * sourceGas_Vdot.medium.R * sourceGas_Vdot.MM; ($RES_SIM_415)
(387|499)      [ARRY] (2) column1.balanceEquations.thermoEquilibrium[1].x_l = column1.x_l_star[1, :]; ($RES_SIM_1057)
(388|501)      [SCAL] (1) sourceGas_Vdot.medium.R = 8.31446261815324 / sourceGas_Vdot.MM; ($RES_SIM_416)
(389|502)      [ARRY] (2) column1.balanceEquations.thermoEquilibrium[1].x_v = column1.x_v_star[1, :]; ($RES_SIM_1056)
(390|504)      [ARRY] (2) column1.balanceEquations.thermoEquilibrium[1].x_vap_liq = column1.x_vap_liq[1, :]; ($RES_SIM_1055)
(391|506)      [SCAL] (1) sourceGas_Vdot.MM = max(1e-6, $FUN_4); ($RES_SIM_418)
(392|507)      [SCAL] (1) column1.balanceEquations.thermoEquilibrium[1].p = column1.p_hyd[1]; ($RES_SIM_1054)
(393|508)      [SCAL] (1) sourceGas_Vdot.medium.p_start = sourceGas_Vdot.medium.state.p; ($RES_SIM_419)
(394|509)      [ARRY] (2) column1.balanceEquations.thermoEquilibrium[1].p_sat = column1.balanceEquations.p_sat[1, :]; ($RES_SIM_1053)
(395|511)      [ARRY] (2) column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.x_l = column1.balanceEquations.thermoEquilibrium[1].x_l; ($RES_SIM_1052)
(396|513)      [SCAL] (1) column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.T = column1.balanceEquations.thermoEquilibrium[1].T; ($RES_SIM_1051)
(397|514)      [ARRY] (2) column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.henryCoeff.x_l = column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.x_l; ($RES_SIM_1050)
(398|516)      [SCAL] (1) $FUN_4 = sum({0.0441, 0.07215} .* sourceGas_Vdot.x); ($RES_AUX_1019)
(399|517)      [SCAL] (1) $FUN_5 = sum(column1.mediumVapour[1].calcSpecificEnthalpy.x .* column1.mediumVapour[1].calcSpecificEnthalpy.hX); ($RES_AUX_1018)
(400|518)      [SCAL] (1) $FUN_6 = sum({0.0441, 0.07215} .* column1.mediumVapour[1].x); ($RES_AUX_1017)
(401|519)      [SCAL] (1) column1.h_v_in = sum(column1.x_v_in .* column1.mediumVapourIn.calcSpecificEnthalpy.hX); ($RES_AUX_1016)
(402|520)      [SCAL] (1) $FUN_8 = sum({0.0441, 0.07215} .* column1.x_v_in); ($RES_AUX_1015)
(403|521)      [SCAL] (1) column1.mediumLiquid[1].h = (108.3096 * column1.mediumLiquid[1].calcSpecificEnthalpy.x[1] + 159.1629 * column1.mediumLiquid[1].calcSpecificEnthalpy.x[2]) * ((-1392.0) + column1.mediumLiquid[1].calcSpecificEnthalpy.T); ($RES_SIM_331)
(404|522)      [FOR-] (2) ($RES_AUX_1013)
(404|522)      [----] for {$i2 in 1:2, $i1 in 1:1} loop
(404|522)      [----]   [SCAL] (1) $FUN_9[$i2, $i1] = log(column1.mediumLiquid[$i1].p_sat[$i2] / ({4.24e6, 3.37e6})[$i2]); ($RES_AUX_1014)
(404|522)      [----] end for;
(405|524)      [SCAL] (1) column1.propsVapIn.h = column1.h_v_in; ($RES_SIM_332)
(406|525)      [SCAL] (1) $FUN_10 = sum({0.0441, 0.07215} .* column1.mediumLiquid[1].x); ($RES_AUX_1012)
(407|526)      [SCAL] (1) column1.propsVapIn.d = column1.rho_v_in; ($RES_SIM_333)
(408|527)      [ARRY] (2) column1.propsVapIn.c = column1.c_v_in; ($RES_SIM_334)
(409|529)      [FOR-] (2) ($RES_AUX_1010)
(409|529)      [----] for $i1 in 1:2 loop
(409|529)      [----]   [SCAL] (1) $FUN_11[$i1] = log(column1.mediumLiquidIn.p_sat[$i1] / ({4.24e6, 3.37e6})[$i1]); ($RES_AUX_1011)
(409|529)      [----] end for;
(410|531)      [SCAL] (1) column1.propsVapIn.u = column1.mediumVapourIn.u; ($RES_SIM_335)
(411|532)      [SCAL] (1) column1.mediumLiquid[1].Theta = min(1.0, 0.00154526068547764 * column1.mediumLiquid[1].T); ($RES_SIM_1137)
(412|533)      [SCAL] (1) column1.propsVapIn.cp = column1.mediumVapourIn.cp; ($RES_SIM_336)
(413|534)      [SCAL] (1) column1.mediumLiquid[1].A[2] = -7.27889; ($RES_SIM_1136)
(414|535)      [SCAL] (1) column1.propsVapIn.lambda = column1.mediumVapourIn.lambda; ($RES_SIM_337)
(415|536)      [SCAL] (1) column1.mediumLiquid[1].A[1] = -6.71791; ($RES_SIM_1135)
(416|537)      [ARRY] (2) column1.propsVapIn.x = column1.x_v_in; ($RES_SIM_338)
(417|539)      [SCAL] (1) column1.mediumLiquid[1].B[2] = 1.48754; ($RES_SIM_1134)
(418|540)      [SCAL] (1) column1.propsVapIn.MM = column1.MM_v_in; ($RES_SIM_339)
(419|541)      [SCAL] (1) column1.mediumLiquid[1].B[1] = 1.33932; ($RES_SIM_1133)
(420|542)      [SCAL] (1) column1.mediumLiquid[1].C[2] = -2.87593; ($RES_SIM_1132)
(421|543)      [SCAL] (1) column1.mediumLiquid[1].C[1] = -2.23017; ($RES_SIM_1131)
(422|544)      [SCAL] (1) column1.mediumLiquid[1].D[2] = -2.22146; ($RES_SIM_1130)
(423|545)      [FOR-] (2) ($RES_SIM_420)
(423|545)      [----] for $i1 in 1:2 loop
(423|545)      [----]   [SCAL] (1) sourceGas_Vdot.medium.X[$i1] = max(1e-5, min(1.0, (({0.0441, 0.07215})[$i1] * sourceGas_Vdot.x[$i1]) / sourceGas_Vdot.MM)); ($RES_SIM_421)
(423|545)      [----] end for;
(424|547)      [SCAL] (1) sourceGas_Vdot.medium.calcSpecificEnthalpy.hX[2] = 0.07215 * (358000.0 + 2206.0 * ((-1406.0) + sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat[2]) + 1599.0 * (sourceGas_Vdot.T - sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat[2])); ($RES_SIM_422)
(425|548)      [SCAL] (1) sourceGas_Vdot.medium.calcSpecificEnthalpy.hX[1] = 0.0441 * (426000.0 + 2456.0 * ((-1406.0) + sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat[1]) + 1549.0 * (sourceGas_Vdot.T - sourceGas_Vdot.medium.calcSpecificEnthalpy.T_sat[1])); ($RES_SIM_423)
(426|549)      [SCAL] (1) column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.henryCoeff.T = column1.balanceEquations.thermoEquilibrium[1].henryCoefficient.T; ($RES_SIM_1049)
(427|550)      [ARRY] (2) column1.reaction[1].c = column1.c_l[1, :]; ($RES_SIM_1048)
(428|552)      [ARRY] (2) sourceLiquid.X = {0.0441, 0.07215} .* sourceLiquid.x[:] / sourceLiquid.medium.MM; ($RES_BND_552)
(429|554)      [ARRY] (2) column1.reaction[1].x = column1.reaction[1].propsLiq.x; ($RES_SIM_1045)
(430|556)      [SCAL] (1) column1.p_bub[1] = $FUN_14; ($RES_SIM_253)
(431|557)      [SCAL] (1) column1.reaction[1].T = column1.reaction[1].propsLiq.T; ($RES_SIM_1044)


===================
  Scalar Matching
===================

variable to equation
**********************
	var 1 --> eqn 540
	var 2 --> eqn 76
	var 3 --> eqn 129
	var 4 --> eqn 185
	var 5 --> eqn 421
	var 6 --> eqn 533
	var 7 --> eqn 142
	var 8 --> eqn 367
	var 9 --> eqn 281
	var 10 --> eqn 202
	var 11 --> eqn 499
	var 12 --> eqn 500
	var 13 --> eqn 231
	var 14 --> eqn 154
	var 15 --> eqn 155
	var 16 --> eqn 12
	var 17 --> eqn 6
	var 18 --> eqn 393
	var 19 --> eqn 394
	var 20 --> eqn 178
	var 21 --> eqn 29
	var 22 --> eqn 112
	var 23 --> eqn 113
	var 24 --> eqn 40
	var 25 --> eqn 41
	var 26 --> eqn 419
	var 27 --> eqn 420
	var 28 --> eqn 368
	var 29 --> eqn 279
	var 30 --> eqn 332
	var 31 --> eqn 42
	var 32 --> eqn 39
	var 33 --> eqn 179
	var 34 --> eqn 262
	var 35 --> eqn 357
	var 36 --> eqn 168
	var 37 --> eqn 497
	var 38 --> eqn 97
	var 39 --> eqn 329
	var 40 --> eqn 330
	var 41 --> eqn 220
	var 42 --> eqn 120
	var 43 --> eqn 256
	var 44 --> eqn 92
	var 45 --> eqn 171
	var 46 --> eqn 80
	var 47 --> eqn 123
	var 48 --> eqn 508
	var 49 --> eqn 94
	var 50 --> eqn 95
	var 51 --> eqn 212
	var 52 --> eqn 213
	var 53 --> eqn 526
	var 54 --> eqn 132
	var 55 --> eqn 54
	var 56 --> eqn 253
	var 57 --> eqn 254
	var 58 --> eqn 186
	var 59 --> eqn 187
	var 60 --> eqn 485
	var 61 --> eqn 469
	var 62 --> eqn 470
	var 63 --> eqn 158
	var 64 --> eqn 239
	var 65 --> eqn 240
	var 66 --> eqn 200
	var 67 --> eqn 201
	var 68 --> eqn 7
	var 69 --> eqn 4
	var 70 --> eqn 78
	var 71 --> eqn 389
	var 72 --> eqn 388
	var 73 --> eqn 100
	var 74 --> eqn 31
	var 75 --> eqn 482
	var 76 --> eqn 326
	var 77 --> eqn 327
	var 78 --> eqn 199
	var 79 --> eqn 303
	var 80 --> eqn 304
	var 81 --> eqn 177
	var 82 --> eqn 130
	var 83 --> eqn 467
	var 84 --> eqn 134
	var 85 --> eqn 361
	var 86 --> eqn 362
	var 87 --> eqn 495
	var 88 --> eqn 261
	var 89 --> eqn 459
	var 90 --> eqn 496
	var 91 --> eqn 93
	var 92 --> eqn 537
	var 93 --> eqn 538
	var 94 --> eqn 264
	var 95 --> eqn 431
	var 96 --> eqn 432
	var 97 --> eqn 300
	var 98 --> eqn 301
	var 99 --> eqn 445
	var 100 --> eqn 471
	var 101 --> eqn 428
	var 102 --> eqn 46
	var 103 --> eqn 433
	var 104 --> eqn 379
	var 105 --> eqn 148
	var 106 --> eqn 430
	var 107 --> eqn 514
	var 108 --> eqn 515
	var 109 --> eqn 434
	var 110 --> eqn 38
	var 111 --> eqn 37
	var 112 --> eqn 88
	var 113 --> eqn 196
	var 114 --> eqn 197
	var 115 --> eqn 442
	var 116 --> eqn 147
	var 117 --> eqn 149
	var 118 --> eqn 22
	var 119 --> eqn 504
	var 120 --> eqn 505
	var 121 --> eqn 334
	var 122 --> eqn 535
	var 123 --> eqn 556
	var 124 --> eqn 278
	var 125 --> eqn 198
	var 126 --> eqn 36
	var 127 --> eqn 32
	var 128 --> eqn 444
	var 129 --> eqn 435
	var 130 --> eqn 89
	var 131 --> eqn 62
	var 132 --> eqn 49
	var 133 --> eqn 50
	var 134 --> eqn 66
	var 135 --> eqn 226
	var 136 --> eqn 227
	var 137 --> eqn 295
	var 138 --> eqn 296
	var 139 --> eqn 159
	var 140 --> eqn 141
	var 141 --> eqn 225
	var 142 --> eqn 549
	var 143 --> eqn 258
	var 144 --> eqn 440
	var 145 --> eqn 458
	var 146 --> eqn 131
	var 147 --> eqn 462
	var 148 --> eqn 463
	var 149 --> eqn 401
	var 150 --> eqn 377
	var 151 --> eqn 479
	var 152 --> eqn 172
	var 153 --> eqn 221
	var 154 --> eqn 492
	var 155 --> eqn 452
	var 156 --> eqn 370
	var 157 --> eqn 25
	var 158 --> eqn 23
	var 159 --> eqn 77
	var 160 --> eqn 243
	var 161 --> eqn 244
	var 162 --> eqn 153
	var 163 --> eqn 65
	var 164 --> eqn 61
	var 165 --> eqn 454
	var 166 --> eqn 272
	var 167 --> eqn 480
	var 168 --> eqn 176
	var 169 --> eqn 509
	var 170 --> eqn 510
	var 171 --> eqn 416
	var 172 --> eqn 417
	var 173 --> eqn 475
	var 174 --> eqn 215
	var 175 --> eqn 306
	var 176 --> eqn 307
	var 177 --> eqn 441
	var 178 --> eqn 465
	var 179 --> eqn 19
	var 180 --> eqn 5
	var 181 --> eqn 119
	var 182 --> eqn 60
	var 183 --> eqn 57
	var 184 --> eqn 139
	var 185 --> eqn 248
	var 186 --> eqn 408
	var 187 --> eqn 484
	var 188 --> eqn 490
	var 189 --> eqn 481
	var 190 --> eqn 276
	var 191 --> eqn 277
	var 192 --> eqn 174
	var 193 --> eqn 328
	var 194 --> eqn 405
	var 195 --> eqn 523
	var 196 --> eqn 266
	var 197 --> eqn 354
	var 198 --> eqn 218
	var 199 --> eqn 429
	var 200 --> eqn 138
	var 201 --> eqn 448
	var 202 --> eqn 108
	var 203 --> eqn 372
	var 204 --> eqn 53
	var 205 --> eqn 487
	var 206 --> eqn 255
	var 207 --> eqn 453
	var 208 --> eqn 451
	var 209 --> eqn 251
	var 210 --> eqn 259
	var 211 --> eqn 241
	var 212 --> eqn 242
	var 213 --> eqn 321
	var 214 --> eqn 81
	var 215 --> eqn 121
	var 216 --> eqn 122
	var 217 --> eqn 68
	var 218 --> eqn 474
	var 219 --> eqn 170
	var 220 --> eqn 82
	var 221 --> eqn 232
	var 222 --> eqn 233
	var 223 --> eqn 543
	var 224 --> eqn 542
	var 225 --> eqn 399
	var 226 --> eqn 400
	var 227 --> eqn 376
	var 228 --> eqn 222
	var 229 --> eqn 217
	var 230 --> eqn 437
	var 231 --> eqn 150
	var 232 --> eqn 346
	var 233 --> eqn 347
	var 234 --> eqn 536
	var 235 --> eqn 534
	var 236 --> eqn 333
	var 237 --> eqn 363
	var 238 --> eqn 364
	var 239 --> eqn 145
	var 240 --> eqn 282
	var 241 --> eqn 189
	var 242 --> eqn 464
	var 243 --> eqn 204
	var 244 --> eqn 205
	var 245 --> eqn 245
	var 246 --> eqn 424
	var 247 --> eqn 58
	var 248 --> eqn 320
	var 249 --> eqn 51
	var 250 --> eqn 552
	var 251 --> eqn 553
	var 252 --> eqn 52
	var 253 --> eqn 390
	var 254 --> eqn 397
	var 255 --> eqn 554
	var 256 --> eqn 555
	var 257 --> eqn 386
	var 258 --> eqn 387
	var 259 --> eqn 369
	var 260 --> eqn 27
	var 261 --> eqn 26
	var 262 --> eqn 33
	var 263 --> eqn 45
	var 264 --> eqn 43
	var 265 --> eqn 20
	var 266 --> eqn 183
	var 267 --> eqn 365
	var 268 --> eqn 366
	var 269 --> eqn 472
	var 270 --> eqn 473
	var 271 --> eqn 493
	var 272 --> eqn 521
	var 273 --> eqn 359
	var 274 --> eqn 360
	var 275 --> eqn 313
	var 276 --> eqn 314
	var 277 --> eqn 324
	var 278 --> eqn 263
	var 279 --> eqn 460
	var 280 --> eqn 461
	var 281 --> eqn 407
	var 282 --> eqn 486
	var 283 --> eqn 501
	var 284 --> eqn 162
	var 285 --> eqn 163
	var 286 --> eqn 164
	var 287 --> eqn 165
	var 288 --> eqn 166
	var 289 --> eqn 167
	var 290 --> eqn 144
	var 291 --> eqn 84
	var 292 --> eqn 446
	var 293 --> eqn 548
	var 294 --> eqn 547
	var 295 --> eqn 285
	var 296 --> eqn 353
	var 297 --> eqn 90
	var 298 --> eqn 182
	var 299 --> eqn 455
	var 300 --> eqn 456
	var 301 --> eqn 210
	var 302 --> eqn 211
	var 303 --> eqn 126
	var 304 --> eqn 206
	var 305 --> eqn 207
	var 306 --> eqn 59
	var 307 --> eqn 28
	var 308 --> eqn 236
	var 309 --> eqn 103
	var 310 --> eqn 104
	var 311 --> eqn 286
	var 312 --> eqn 229
	var 313 --> eqn 228
	var 314 --> eqn 270
	var 315 --> eqn 271
	var 316 --> eqn 331
	var 317 --> eqn 343
	var 318 --> eqn 322
	var 319 --> eqn 382
	var 320 --> eqn 175
	var 321 --> eqn 85
	var 322 --> eqn 79
	var 323 --> eqn 351
	var 324 --> eqn 352
	var 325 --> eqn 423
	var 326 --> eqn 494
	var 327 --> eqn 375
	var 328 --> eqn 269
	var 329 --> eqn 106
	var 330 --> eqn 107
	var 331 --> eqn 3
	var 332 --> eqn 418
	var 333 --> eqn 83
	var 334 --> eqn 349
	var 335 --> eqn 439
	var 336 --> eqn 273
	var 337 --> eqn 274
	var 338 --> eqn 348
	var 339 --> eqn 378
	var 340 --> eqn 151
	var 341 --> eqn 71
	var 342 --> eqn 519
	var 343 --> eqn 289
	var 344 --> eqn 290
	var 345 --> eqn 291
	var 346 --> eqn 257
	var 347 --> eqn 341
	var 348 --> eqn 342
	var 349 --> eqn 292
	var 350 --> eqn 234
	var 351 --> eqn 235
	var 352 --> eqn 293
	var 353 --> eqn 506
	var 354 --> eqn 63
	var 355 --> eqn 64
	var 356 --> eqn 398
	var 357 --> eqn 335
	var 358 --> eqn 336
	var 359 --> eqn 265
	var 360 --> eqn 124
	var 361 --> eqn 125
	var 362 --> eqn 73
	var 363 --> eqn 96
	var 364 --> eqn 193
	var 365 --> eqn 194
	var 366 --> eqn 527
	var 367 --> eqn 528
	var 368 --> eqn 356
	var 369 --> eqn 380
	var 370 --> eqn 381
	var 371 --> eqn 374
	var 372 --> eqn 325
	var 373 --> eqn 457
	var 374 --> eqn 427
	var 375 --> eqn 513
	var 376 --> eqn 13
	var 377 --> eqn 438
	var 378 --> eqn 21
	var 379 --> eqn 74
	var 380 --> eqn 75
	var 381 --> eqn 414
	var 382 --> eqn 384
	var 383 --> eqn 294
	var 384 --> eqn 114
	var 385 --> eqn 115
	var 386 --> eqn 383
	var 387 --> eqn 522
	var 388 --> eqn 491
	var 389 --> eqn 299
	var 390 --> eqn 520
	var 391 --> eqn 302
	var 392 --> eqn 305
	var 393 --> eqn 140
	var 394 --> eqn 518
	var 395 --> eqn 308
	var 396 --> eqn 350
	var 397 --> eqn 517
	var 398 --> eqn 184
	var 399 --> eqn 516
	var 400 --> eqn 315
	var 401 --> eqn 476
	var 402 --> eqn 316
	var 403 --> eqn 224
	var 404 --> eqn 252
	var 405 --> eqn 317
	var 406 --> eqn 402
	var 407 --> eqn 403
	var 408 --> eqn 507
	var 409 --> eqn 133
	var 410 --> eqn 477
	var 411 --> eqn 478
	var 412 --> eqn 55
	var 413 --> eqn 56
	var 414 --> eqn 436
	var 415 --> eqn 425
	var 416 --> eqn 422
	var 417 --> eqn 267
	var 418 --> eqn 268
	var 419 --> eqn 47
	var 420 --> eqn 48
	var 421 --> eqn 511
	var 422 --> eqn 512
	var 423 --> eqn 312
	var 424 --> eqn 391
	var 425 --> eqn 191
	var 426 --> eqn 192
	var 427 --> eqn 415
	var 428 --> eqn 466
	var 429 --> eqn 223
	var 430 --> eqn 230
	var 431 --> eqn 101
	var 432 --> eqn 181
	var 433 --> eqn 67
	var 434 --> eqn 98
	var 435 --> eqn 318
	var 436 --> eqn 319
	var 437 --> eqn 310
	var 438 --> eqn 311
	var 439 --> eqn 24
	var 440 --> eqn -1
	var 441 --> eqn 284
	var 442 --> eqn 10
	var 443 --> eqn 11
	var 444 --> eqn 102
	var 445 --> eqn 447
	var 446 --> eqn 9
	var 447 --> eqn 208
	var 448 --> eqn 209
	var 449 --> eqn 385
	var 450 --> eqn 525
	var 451 --> eqn 275
	var 452 --> eqn 14
	var 453 --> eqn 358
	var 454 --> eqn 287
	var 455 --> eqn 288
	var 456 --> eqn 283
	var 457 --> eqn 309
	var 458 --> eqn 246
	var 459 --> eqn 449
	var 460 --> eqn 450
	var 461 --> eqn 532
	var 462 --> eqn 111
	var 463 --> eqn 156
	var 464 --> eqn 137
	var 465 --> eqn 531
	var 466 --> eqn 143
	var 467 --> eqn 409
	var 468 --> eqn 410
	var 469 --> eqn 18
	var 470 --> eqn 127
	var 471 --> eqn 188
	var 472 --> eqn 157
	var 473 --> eqn 109
	var 474 --> eqn 110
	var 475 --> eqn 498
	var 476 --> eqn 529
	var 477 --> eqn 530
	var 478 --> eqn 118
	var 479 --> eqn 406
	var 480 --> eqn 550
	var 481 --> eqn 551
	var 482 --> eqn 355
	var 483 --> eqn 541
	var 484 --> eqn 539
	var 485 --> eqn 128
	var 486 --> eqn 260
	var 487 --> eqn 371
	var 488 --> eqn 99
	var 489 --> eqn 86
	var 490 --> eqn 87
	var 491 --> eqn 344
	var 492 --> eqn 135
	var 493 --> eqn 488
	var 494 --> eqn 502
	var 495 --> eqn 503
	var 496 --> eqn 30
	var 497 --> eqn 8
	var 498 --> eqn 136
	var 499 --> eqn 152
	var 500 --> eqn 426
	var 501 --> eqn 373
	var 502 --> eqn 91
	var 503 --> eqn 72
	var 504 --> eqn 247
	var 505 --> eqn 116
	var 506 --> eqn 17
	var 507 --> eqn 544
	var 508 --> eqn 34
	var 509 --> eqn 35
	var 510 --> eqn 339
	var 511 --> eqn 340
	var 512 --> eqn 169
	var 513 --> eqn 219
	var 514 --> eqn 323
	var 515 --> eqn 524
	var 516 --> eqn 216
	var 517 --> eqn 249
	var 518 --> eqn 250
	var 519 --> eqn 392
	var 520 --> eqn 180
	var 521 --> eqn 173
	var 522 --> eqn 411
	var 523 --> eqn 412
	var 524 --> eqn 345
	var 525 --> eqn 214
	var 526 --> eqn 160
	var 527 --> eqn 161
	var 528 --> eqn 489
	var 529 --> eqn 468
	var 530 --> eqn 337
	var 531 --> eqn 338
	var 532 --> eqn 15
	var 533 --> eqn 16
	var 534 --> eqn 146
	var 535 --> eqn 545
	var 536 --> eqn 546
	var 537 --> eqn 190
	var 538 --> eqn 280
	var 539 --> eqn 1
	var 540 --> eqn 2
	var 541 --> eqn 237
	var 542 --> eqn 238
	var 543 --> eqn 557
	var 544 --> eqn 297
	var 545 --> eqn 298
	var 546 --> eqn 413
	var 547 --> eqn 69
	var 548 --> eqn 70
	var 549 --> eqn 195
	var 550 --> eqn 117
	var 551 --> eqn 105
	var 552 --> eqn 483
	var 553 --> eqn 443
	var 554 --> eqn 203
	var 555 --> eqn 404
	var 556 --> eqn 44
	var 557 --> eqn 395
	var 558 --> eqn 396

equation to variable
**********************
	eqn 1 --> var 539
	eqn 2 --> var 540
	eqn 3 --> var 331
	eqn 4 --> var 69
	eqn 5 --> var 180
	eqn 6 --> var 17
	eqn 7 --> var 68
	eqn 8 --> var 497
	eqn 9 --> var 446
	eqn 10 --> var 442
	eqn 11 --> var 443
	eqn 12 --> var 16
	eqn 13 --> var 376
	eqn 14 --> var 452
	eqn 15 --> var 532
	eqn 16 --> var 533
	eqn 17 --> var 506
	eqn 18 --> var 469
	eqn 19 --> var 179
	eqn 20 --> var 265
	eqn 21 --> var 378
	eqn 22 --> var 118
	eqn 23 --> var 158
	eqn 24 --> var 439
	eqn 25 --> var 157
	eqn 26 --> var 261
	eqn 27 --> var 260
	eqn 28 --> var 307
	eqn 29 --> var 21
	eqn 30 --> var 496
	eqn 31 --> var 74
	eqn 32 --> var 127
	eqn 33 --> var 262
	eqn 34 --> var 508
	eqn 35 --> var 509
	eqn 36 --> var 126
	eqn 37 --> var 111
	eqn 38 --> var 110
	eqn 39 --> var 32
	eqn 40 --> var 24
	eqn 41 --> var 25
	eqn 42 --> var 31
	eqn 43 --> var 264
	eqn 44 --> var 556
	eqn 45 --> var 263
	eqn 46 --> var 102
	eqn 47 --> var 419
	eqn 48 --> var 420
	eqn 49 --> var 132
	eqn 50 --> var 133
	eqn 51 --> var 249
	eqn 52 --> var 252
	eqn 53 --> var 204
	eqn 54 --> var 55
	eqn 55 --> var 412
	eqn 56 --> var 413
	eqn 57 --> var 183
	eqn 58 --> var 247
	eqn 59 --> var 306
	eqn 60 --> var 182
	eqn 61 --> var 164
	eqn 62 --> var 131
	eqn 63 --> var 354
	eqn 64 --> var 355
	eqn 65 --> var 163
	eqn 66 --> var 134
	eqn 67 --> var 433
	eqn 68 --> var 217
	eqn 69 --> var 547
	eqn 70 --> var 548
	eqn 71 --> var 341
	eqn 72 --> var 503
	eqn 73 --> var 362
	eqn 74 --> var 379
	eqn 75 --> var 380
	eqn 76 --> var 2
	eqn 77 --> var 159
	eqn 78 --> var 70
	eqn 79 --> var 322
	eqn 80 --> var 46
	eqn 81 --> var 214
	eqn 82 --> var 220
	eqn 83 --> var 333
	eqn 84 --> var 291
	eqn 85 --> var 321
	eqn 86 --> var 489
	eqn 87 --> var 490
	eqn 88 --> var 112
	eqn 89 --> var 130
	eqn 90 --> var 297
	eqn 91 --> var 502
	eqn 92 --> var 44
	eqn 93 --> var 91
	eqn 94 --> var 49
	eqn 95 --> var 50
	eqn 96 --> var 363
	eqn 97 --> var 38
	eqn 98 --> var 434
	eqn 99 --> var 488
	eqn 100 --> var 73
	eqn 101 --> var 431
	eqn 102 --> var 444
	eqn 103 --> var 309
	eqn 104 --> var 310
	eqn 105 --> var 551
	eqn 106 --> var 329
	eqn 107 --> var 330
	eqn 108 --> var 202
	eqn 109 --> var 473
	eqn 110 --> var 474
	eqn 111 --> var 462
	eqn 112 --> var 22
	eqn 113 --> var 23
	eqn 114 --> var 384
	eqn 115 --> var 385
	eqn 116 --> var 505
	eqn 117 --> var 550
	eqn 118 --> var 478
	eqn 119 --> var 181
	eqn 120 --> var 42
	eqn 121 --> var 215
	eqn 122 --> var 216
	eqn 123 --> var 47
	eqn 124 --> var 360
	eqn 125 --> var 361
	eqn 126 --> var 303
	eqn 127 --> var 470
	eqn 128 --> var 485
	eqn 129 --> var 3
	eqn 130 --> var 82
	eqn 131 --> var 146
	eqn 132 --> var 54
	eqn 133 --> var 409
	eqn 134 --> var 84
	eqn 135 --> var 492
	eqn 136 --> var 498
	eqn 137 --> var 464
	eqn 138 --> var 200
	eqn 139 --> var 184
	eqn 140 --> var 393
	eqn 141 --> var 140
	eqn 142 --> var 7
	eqn 143 --> var 466
	eqn 144 --> var 290
	eqn 145 --> var 239
	eqn 146 --> var 534
	eqn 147 --> var 116
	eqn 148 --> var 105
	eqn 149 --> var 117
	eqn 150 --> var 231
	eqn 151 --> var 340
	eqn 152 --> var 499
	eqn 153 --> var 162
	eqn 154 --> var 14
	eqn 155 --> var 15
	eqn 156 --> var 463
	eqn 157 --> var 472
	eqn 158 --> var 63
	eqn 159 --> var 139
	eqn 160 --> var 526
	eqn 161 --> var 527
	eqn 162 --> var 284
	eqn 163 --> var 285
	eqn 164 --> var 286
	eqn 165 --> var 287
	eqn 166 --> var 288
	eqn 167 --> var 289
	eqn 168 --> var 36
	eqn 169 --> var 512
	eqn 170 --> var 219
	eqn 171 --> var 45
	eqn 172 --> var 152
	eqn 173 --> var 521
	eqn 174 --> var 192
	eqn 175 --> var 320
	eqn 176 --> var 168
	eqn 177 --> var 81
	eqn 178 --> var 20
	eqn 179 --> var 33
	eqn 180 --> var 520
	eqn 181 --> var 432
	eqn 182 --> var 298
	eqn 183 --> var 266
	eqn 184 --> var 398
	eqn 185 --> var 4
	eqn 186 --> var 58
	eqn 187 --> var 59
	eqn 188 --> var 471
	eqn 189 --> var 241
	eqn 190 --> var 537
	eqn 191 --> var 425
	eqn 192 --> var 426
	eqn 193 --> var 364
	eqn 194 --> var 365
	eqn 195 --> var 549
	eqn 196 --> var 113
	eqn 197 --> var 114
	eqn 198 --> var 125
	eqn 199 --> var 78
	eqn 200 --> var 66
	eqn 201 --> var 67
	eqn 202 --> var 10
	eqn 203 --> var 554
	eqn 204 --> var 243
	eqn 205 --> var 244
	eqn 206 --> var 304
	eqn 207 --> var 305
	eqn 208 --> var 447
	eqn 209 --> var 448
	eqn 210 --> var 301
	eqn 211 --> var 302
	eqn 212 --> var 51
	eqn 213 --> var 52
	eqn 214 --> var 525
	eqn 215 --> var 174
	eqn 216 --> var 516
	eqn 217 --> var 229
	eqn 218 --> var 198
	eqn 219 --> var 513
	eqn 220 --> var 41
	eqn 221 --> var 153
	eqn 222 --> var 228
	eqn 223 --> var 429
	eqn 224 --> var 403
	eqn 225 --> var 141
	eqn 226 --> var 135
	eqn 227 --> var 136
	eqn 228 --> var 313
	eqn 229 --> var 312
	eqn 230 --> var 430
	eqn 231 --> var 13
	eqn 232 --> var 221
	eqn 233 --> var 222
	eqn 234 --> var 350
	eqn 235 --> var 351
	eqn 236 --> var 308
	eqn 237 --> var 541
	eqn 238 --> var 542
	eqn 239 --> var 64
	eqn 240 --> var 65
	eqn 241 --> var 211
	eqn 242 --> var 212
	eqn 243 --> var 160
	eqn 244 --> var 161
	eqn 245 --> var 245
	eqn 246 --> var 458
	eqn 247 --> var 504
	eqn 248 --> var 185
	eqn 249 --> var 517
	eqn 250 --> var 518
	eqn 251 --> var 209
	eqn 252 --> var 404
	eqn 253 --> var 56
	eqn 254 --> var 57
	eqn 255 --> var 206
	eqn 256 --> var 43
	eqn 257 --> var 346
	eqn 258 --> var 143
	eqn 259 --> var 210
	eqn 260 --> var 486
	eqn 261 --> var 88
	eqn 262 --> var 34
	eqn 263 --> var 278
	eqn 264 --> var 94
	eqn 265 --> var 359
	eqn 266 --> var 196
	eqn 267 --> var 417
	eqn 268 --> var 418
	eqn 269 --> var 328
	eqn 270 --> var 314
	eqn 271 --> var 315
	eqn 272 --> var 166
	eqn 273 --> var 336
	eqn 274 --> var 337
	eqn 275 --> var 451
	eqn 276 --> var 190
	eqn 277 --> var 191
	eqn 278 --> var 124
	eqn 279 --> var 29
	eqn 280 --> var 538
	eqn 281 --> var 9
	eqn 282 --> var 240
	eqn 283 --> var 456
	eqn 284 --> var 441
	eqn 285 --> var 295
	eqn 286 --> var 311
	eqn 287 --> var 454
	eqn 288 --> var 455
	eqn 289 --> var 343
	eqn 290 --> var 344
	eqn 291 --> var 345
	eqn 292 --> var 349
	eqn 293 --> var 352
	eqn 294 --> var 383
	eqn 295 --> var 137
	eqn 296 --> var 138
	eqn 297 --> var 544
	eqn 298 --> var 545
	eqn 299 --> var 389
	eqn 300 --> var 97
	eqn 301 --> var 98
	eqn 302 --> var 391
	eqn 303 --> var 79
	eqn 304 --> var 80
	eqn 305 --> var 392
	eqn 306 --> var 175
	eqn 307 --> var 176
	eqn 308 --> var 395
	eqn 309 --> var 457
	eqn 310 --> var 437
	eqn 311 --> var 438
	eqn 312 --> var 423
	eqn 313 --> var 275
	eqn 314 --> var 276
	eqn 315 --> var 400
	eqn 316 --> var 402
	eqn 317 --> var 405
	eqn 318 --> var 435
	eqn 319 --> var 436
	eqn 320 --> var 248
	eqn 321 --> var 213
	eqn 322 --> var 318
	eqn 323 --> var 514
	eqn 324 --> var 277
	eqn 325 --> var 372
	eqn 326 --> var 76
	eqn 327 --> var 77
	eqn 328 --> var 193
	eqn 329 --> var 39
	eqn 330 --> var 40
	eqn 331 --> var 316
	eqn 332 --> var 30
	eqn 333 --> var 236
	eqn 334 --> var 121
	eqn 335 --> var 357
	eqn 336 --> var 358
	eqn 337 --> var 530
	eqn 338 --> var 531
	eqn 339 --> var 510
	eqn 340 --> var 511
	eqn 341 --> var 347
	eqn 342 --> var 348
	eqn 343 --> var 317
	eqn 344 --> var 491
	eqn 345 --> var 524
	eqn 346 --> var 232
	eqn 347 --> var 233
	eqn 348 --> var 338
	eqn 349 --> var 334
	eqn 350 --> var 396
	eqn 351 --> var 323
	eqn 352 --> var 324
	eqn 353 --> var 296
	eqn 354 --> var 197
	eqn 355 --> var 482
	eqn 356 --> var 368
	eqn 357 --> var 35
	eqn 358 --> var 453
	eqn 359 --> var 273
	eqn 360 --> var 274
	eqn 361 --> var 85
	eqn 362 --> var 86
	eqn 363 --> var 237
	eqn 364 --> var 238
	eqn 365 --> var 267
	eqn 366 --> var 268
	eqn 367 --> var 8
	eqn 368 --> var 28
	eqn 369 --> var 259
	eqn 370 --> var 156
	eqn 371 --> var 487
	eqn 372 --> var 203
	eqn 373 --> var 501
	eqn 374 --> var 371
	eqn 375 --> var 327
	eqn 376 --> var 227
	eqn 377 --> var 150
	eqn 378 --> var 339
	eqn 379 --> var 104
	eqn 380 --> var 369
	eqn 381 --> var 370
	eqn 382 --> var 319
	eqn 383 --> var 386
	eqn 384 --> var 382
	eqn 385 --> var 449
	eqn 386 --> var 257
	eqn 387 --> var 258
	eqn 388 --> var 72
	eqn 389 --> var 71
	eqn 390 --> var 253
	eqn 391 --> var 424
	eqn 392 --> var 519
	eqn 393 --> var 18
	eqn 394 --> var 19
	eqn 395 --> var 557
	eqn 396 --> var 558
	eqn 397 --> var 254
	eqn 398 --> var 356
	eqn 399 --> var 225
	eqn 400 --> var 226
	eqn 401 --> var 149
	eqn 402 --> var 406
	eqn 403 --> var 407
	eqn 404 --> var 555
	eqn 405 --> var 194
	eqn 406 --> var 479
	eqn 407 --> var 281
	eqn 408 --> var 186
	eqn 409 --> var 467
	eqn 410 --> var 468
	eqn 411 --> var 522
	eqn 412 --> var 523
	eqn 413 --> var 546
	eqn 414 --> var 381
	eqn 415 --> var 427
	eqn 416 --> var 171
	eqn 417 --> var 172
	eqn 418 --> var 332
	eqn 419 --> var 26
	eqn 420 --> var 27
	eqn 421 --> var 5
	eqn 422 --> var 416
	eqn 423 --> var 325
	eqn 424 --> var 246
	eqn 425 --> var 415
	eqn 426 --> var 500
	eqn 427 --> var 374
	eqn 428 --> var 101
	eqn 429 --> var 199
	eqn 430 --> var 106
	eqn 431 --> var 95
	eqn 432 --> var 96
	eqn 433 --> var 103
	eqn 434 --> var 109
	eqn 435 --> var 129
	eqn 436 --> var 414
	eqn 437 --> var 230
	eqn 438 --> var 377
	eqn 439 --> var 335
	eqn 440 --> var 144
	eqn 441 --> var 177
	eqn 442 --> var 115
	eqn 443 --> var 553
	eqn 444 --> var 128
	eqn 445 --> var 99
	eqn 446 --> var 292
	eqn 447 --> var 445
	eqn 448 --> var 201
	eqn 449 --> var 459
	eqn 450 --> var 460
	eqn 451 --> var 208
	eqn 452 --> var 155
	eqn 453 --> var 207
	eqn 454 --> var 165
	eqn 455 --> var 299
	eqn 456 --> var 300
	eqn 457 --> var 373
	eqn 458 --> var 145
	eqn 459 --> var 89
	eqn 460 --> var 279
	eqn 461 --> var 280
	eqn 462 --> var 147
	eqn 463 --> var 148
	eqn 464 --> var 242
	eqn 465 --> var 178
	eqn 466 --> var 428
	eqn 467 --> var 83
	eqn 468 --> var 529
	eqn 469 --> var 61
	eqn 470 --> var 62
	eqn 471 --> var 100
	eqn 472 --> var 269
	eqn 473 --> var 270
	eqn 474 --> var 218
	eqn 475 --> var 173
	eqn 476 --> var 401
	eqn 477 --> var 410
	eqn 478 --> var 411
	eqn 479 --> var 151
	eqn 480 --> var 167
	eqn 481 --> var 189
	eqn 482 --> var 75
	eqn 483 --> var 552
	eqn 484 --> var 187
	eqn 485 --> var 60
	eqn 486 --> var 282
	eqn 487 --> var 205
	eqn 488 --> var 493
	eqn 489 --> var 528
	eqn 490 --> var 188
	eqn 491 --> var 388
	eqn 492 --> var 154
	eqn 493 --> var 271
	eqn 494 --> var 326
	eqn 495 --> var 87
	eqn 496 --> var 90
	eqn 497 --> var 37
	eqn 498 --> var 475
	eqn 499 --> var 11
	eqn 500 --> var 12
	eqn 501 --> var 283
	eqn 502 --> var 494
	eqn 503 --> var 495
	eqn 504 --> var 119
	eqn 505 --> var 120
	eqn 506 --> var 353
	eqn 507 --> var 408
	eqn 508 --> var 48
	eqn 509 --> var 169
	eqn 510 --> var 170
	eqn 511 --> var 421
	eqn 512 --> var 422
	eqn 513 --> var 375
	eqn 514 --> var 107
	eqn 515 --> var 108
	eqn 516 --> var 399
	eqn 517 --> var 397
	eqn 518 --> var 394
	eqn 519 --> var 342
	eqn 520 --> var 390
	eqn 521 --> var 272
	eqn 522 --> var 387
	eqn 523 --> var 195
	eqn 524 --> var 515
	eqn 525 --> var 450
	eqn 526 --> var 53
	eqn 527 --> var 366
	eqn 528 --> var 367
	eqn 529 --> var 476
	eqn 530 --> var 477
	eqn 531 --> var 465
	eqn 532 --> var 461
	eqn 533 --> var 6
	eqn 534 --> var 235
	eqn 535 --> var 122
	eqn 536 --> var 234
	eqn 537 --> var 92
	eqn 538 --> var 93
	eqn 539 --> var 484
	eqn 540 --> var 1
	eqn 541 --> var 483
	eqn 542 --> var 224
	eqn 543 --> var 223
	eqn 544 --> var 507
	eqn 545 --> var 535
	eqn 546 --> var 536
	eqn 547 --> var 294
	eqn 548 --> var 293
	eqn 549 --> var 142
	eqn 550 --> var 480
	eqn 551 --> var 481
	eqn 552 --> var 250
	eqn 553 --> var 251
	eqn 554 --> var 255
	eqn 555 --> var 256
	eqn 556 --> var 123
	eqn 557 --> var 543
"
[Timeout remaining time 659]
[Calling sys.exit(0), Time elapsed: 3.19125073030591]