diff --git a/fsl_mrs_mce/mc_in.py b/fsl_mrs_mce/mc_in.py
index e5d3f7c..a40856c 100644
--- a/fsl_mrs_mce/mc_in.py
+++ b/fsl_mrs_mce/mc_in.py
@@ -30,7 +30,7 @@ def load_config(input_folder="input/", section="default"):
 
 
 def conc_df_from_file(input_path="input/", foundation="vivo_average",
-                      include_molecules=[]):
+                      include_molecules=[], save_df=True):
     sections = load_config(input_path, "")
     """Generates concentration dataframe for mc method.
     Supports different variation types (absolute/delta/relative)
@@ -118,4 +118,7 @@ def conc_df_from_file(input_path="input/", foundation="vivo_average",
         df_conc[parameter_name] = parameter_values
         i = i+1
 
+    if save_df:
+        df_conc.to_csv(str(input_path+"df_conc.csv"))
+
     return df_conc