diff --git a/createTestData/fitting_all_test_data_profit.m b/createTestData/fitting_all_test_data_profit.m
index 7b7bf3386f31f4c3265826c0866d67a969cd8684..3564864e4ff18653ff3228a6a654f0b05b0185b6 100644
--- a/createTestData/fitting_all_test_data_profit.m
+++ b/createTestData/fitting_all_test_data_profit.m
@@ -96,6 +96,7 @@ for indexSubj = 1:numberOfSubjects
     %     load([filenameWater '.mat'], 'aWater', 'maxWaterPeak');
     
     concentrations = cell(numberOfTEs,length(namingSuffixesAve));
+    concentrations_H2O = cell(numberOfTEs,length(namingSuffixesAve));
     if ~fitLCModel
         FQNs = cell(numberOfTEs,length(namingSuffixesAve));
         FQNs2 = cell(numberOfTEs,length(namingSuffixesAve));
@@ -136,6 +137,8 @@ for indexSubj = 1:numberOfSubjects
                     end
                     [concentrations{indexTE,indexAverageDelete, indexCoilDelete}, activeMetabolites, FQNs{indexTE,indexAverageDelete, indexCoilDelete}, FQNs2{indexTE,indexAverageDelete, indexCoilDelete}] = ...
                         plot_ProFit_Metabolites_Fit([filenameData, '_profit.mat'],dataExportPathTE, 'Met', true);
+                    load([dataExportPathTE, filenameData, '_profit_H2O.mat'],'concH2O');
+                    concentrations_H2O{indexTE,indexAverageDelete, indexCoilDelete}=concH2O;
                 end
                 close all
             end
@@ -145,6 +148,7 @@ for indexSubj = 1:numberOfSubjects
         save([dataExportPath{indexSubj}, 'concentrations' truncSuffix '.mat'], 'concentrations');
     else
         save([dataExportPath{indexSubj}, 'concentrations', truncSuffix '_profit.mat'], 'concentrations', 'activeMetabolites', 'FQNs', 'FQNs2');
+        save([dataExportPath{indexSubj}, 'concentrations', truncSuffix '_H2O_profit.mat'], 'concentrations_H2O');
     end
 end