diff --git a/RBBGCMuso/R/musoMonte.R b/RBBGCMuso/R/musoMonte.R
index 5e36263..fe94833 100644
--- a/RBBGCMuso/R/musoMonte.R
+++ b/RBBGCMuso/R/musoMonte.R
@@ -33,6 +33,7 @@ musoMonte <- function(settings=NULL,
                      constrains = NULL,
                      skipZero = TRUE,
                      postProcString=NULL,
+                     modifyOut=TRUE,
                      ...){
 
 
@@ -149,6 +150,7 @@ musoMonte <- function(settings=NULL,
                              skipSpinup = skipSpinup,
                              keepEpc = keepEpc,
                              debugging = debugging,
+                             modifyOriginal = modifyOut,
                              outVars = outVars,postProcString=postProcString), error = function (e) NA)
             
             if(length(dim(tmp))>=1){
diff --git a/RBBGCMuso/R/musoSensi.R b/RBBGCMuso/R/musoSensi.R
index 674d96b..4f4130b 100644
--- a/RBBGCMuso/R/musoSensi.R
+++ b/RBBGCMuso/R/musoSensi.R
@@ -32,6 +32,7 @@ musoSensi <- function(monteCarloFile = NULL,
                      skipSpinup = TRUE,
                      skipZero = TRUE,
                      postProcString=NULL,
+                     modifyOut=TRUE,
                      dpi=300){
 
     if(is.null(parameters)){
@@ -96,7 +97,8 @@ musoSensi <- function(monteCarloFile = NULL,
                       varIndex = varIndex,
                       skipSpinup = skipSpinup,
                       skipZero=skipZero,
-                      postProcString=postProcString
+                      postProcString=postProcString,
+                      modifyOut=modifyOut
                       )
         M <- cbind(seq_along(M[,1]),M)
         yInd <-  grep("mod.", colnames(M))[varIndex]