Edit two examples, remove some commented lines

sintefmath · Oct 9, 2022 · 17d9e21 · 17d9e21
1 parent bda74f3
commit 17d9e21
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 20 deletions.
diff --git a/darcy/co2_brine_2d_vertical.jl b/darcy/co2_brine_2d_vertical.jl
@@ -1,7 +1,3 @@
-#activate("../JutulExamples.jl/darcy/.")
-
-#using Pkg
-#Pkg.activate("./../JutulExamples.jl/darcy")
 
 using MultiComponentFlash
 h2o = MolecularProperty(0.018015268, 22.064e6, 647.096, 5.595e-05, 0.3442920843)
@@ -53,12 +49,7 @@ day = 24*3600.0
 dt0 = repeat([1]*day, 26)
 dt1 = repeat([10.0]*day, 180)
 dt = append!(dt0, dt1)
-#reservoir = reservoir_model(model)
-#pv = pore_volume(model)
-#inj_rate = sum(pv)/sum(dt)
-#rate_target = TotalRateTarget(inj_rate)
 rate_target = TotalRateTarget(9.5066e-06)
-#ϵ = 0
 I_ctrl = InjectorControl(rate_target, [0, 1], density = rhoVS)
 bhp_target = BottomHolePressureTarget(50*bar)
 P_ctrl = ProducerControl(bhp_target)
@@ -78,4 +69,4 @@ display(f)
 ## Plot the wells
 wd = full_well_outputs(sim.model, states, forces)
 time = report_times(reports)
-plot_well_results(wd, time)
+plot_well_results(wd, time)
diff --git a/darcy/compositional_5components.jl b/darcy/compositional_5components.jl
@@ -1,7 +1,3 @@
-#activate("../JutulExamples.jl/darcy/.")
-
-#using Pkg
-#Pkg.activate("./../JutulExamples.jl/darcy")
 
 using MultiComponentFlash
 n2_ch4 = MolecularProperty(0.0161594, 4.58e6, 189.515, 9.9701e-05, 0.00854)
@@ -55,12 +51,7 @@ state0 = setup_reservoir_state(model, Pressure = 225*bar, OverallMoleFractions =
 
 day = 24*3600.0
 dt = repeat([2.0]*day, 365)
-#reservoir = reservoir_model(model)
-#pv = pore_volume(model)
-#inj_rate = sum(pv)/sum(dt)
-#rate_target = TotalRateTarget(inj_rate)
 rate_target = TotalRateTarget(0.0015)
-#ϵ = 0
 I_ctrl = InjectorControl(rate_target, [0, 1, 0, 0, 0], density = rhoVS)
 bhp_target = BottomHolePressureTarget(100*bar)
 P_ctrl = ProducerControl(bhp_target)
@@ -80,4 +71,4 @@ display(f)
 ## Plot the wells
 wd = full_well_outputs(sim.model, states, forces)
 time = report_times(reports)
-plot_well_results(wd, time)
+plot_well_results(wd, time)