diff --git a/networkComplexityBigg/helper/applyMediaBoundsOnModel.R b/networkComplexityBigg/helper/applyMediaBoundsOnModel.R
new file mode 100644
index 0000000000000000000000000000000000000000..c51fb040385312f6ca7485e78dc995dc7af18a13
--- /dev/null
+++ b/networkComplexityBigg/helper/applyMediaBoundsOnModel.R
@@ -0,0 +1,17 @@
+applyMediaBoundsOnModel <- function(model, medium, uptakeLowBnd=-10){
+	
+	ex <- findExchReact(model)
+	
+	lowbnd(model)[react_pos(ex)] <- 0
+	
+	exOnly <- ex[grep("^EX_", react_id(ex))]
+	exMap <- react_pos(exOnly)
+	names(exMap) <- gsub("\\[\\w\\]$", "", met_id(exOnly))
+	
+	lb <- rep(0, length(exMap))
+	lb[names(exMap) %in% medium] <- uptakeLowBnd
+	
+	lowbnd(model)[exMap] <- lb
+	
+	return(model)
+}