diff --git a/src/optimization/__pycache__/model_builder.cpython-313.pyc b/src/optimization/__pycache__/model_builder.cpython-313.pyc
index b1842a4..e3cdeb7 100644
Binary files a/src/optimization/__pycache__/model_builder.cpython-313.pyc and b/src/optimization/__pycache__/model_builder.cpython-313.pyc differ
diff --git a/src/optimization/model_builder.py b/src/optimization/model_builder.py
index 7bf9e45..b6b560c 100644
--- a/src/optimization/model_builder.py
+++ b/src/optimization/model_builder.py
@@ -835,19 +835,6 @@ def build_model(
                         >= model.bunker_target[j] - model.bunker_total[j, w_prev, d_prev]
                     )
 
-    model.lambda_J = pyo.Param(initialize=100_000, mutable=True, within=pyo.NonNegativeReals)
-
-    def welzow_to_J(m):
-        return sum(
-            m.x["Welzow", "J", w, d, s]
-            for w in m.W
-            for d in m.D
-            for s in m.S
-            if ("Welzow", "J", w, d, s) in m.x
-        )
-
-    model.bonus_welzow_J = pyo.Expression(rule=welzow_to_J)
-
     def demand_rule(m, j, w, d):
         return m.dV_N[w, d] + m.dV_W[w, d] if j == "V" else m.d[j, w, d]