From 63bdd86a41a1b27c3d1a3f6ebc9e29db26b2f90b Mon Sep 17 00:00:00 2001
From: Eduard Kerkhoven <eduardk@chalmers.se>
Date: Mon, 27 May 2024 00:50:50 +0200
Subject: [PATCH] feat: constrainEnzConcs keep prot pool draw

addresses #375
---
 src/geckomat/limit_proteins/constrainEnzConcs.m | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/geckomat/limit_proteins/constrainEnzConcs.m b/src/geckomat/limit_proteins/constrainEnzConcs.m
index 736f30c28..66094b53b 100644
--- a/src/geckomat/limit_proteins/constrainEnzConcs.m
+++ b/src/geckomat/limit_proteins/constrainEnzConcs.m
@@ -43,7 +43,7 @@
 
 %If non-NaN in model.ec.concs, then constrain by UB
 if any(protCons)
-    model.S(protPoolIdx, usageRxnsIdx(protCons)) = 0;
+%    model.S(protPoolIdx, usageRxnsIdx(protCons)) = 0; % Since GECKO 3.2.0
     model.lb(usageRxnsIdx(protCons)) = -model.ec.concs(protCons);
 end
 end