diff --git a/BayesValidRox/PostProcessing/PostProcessing.py b/BayesValidRox/PostProcessing/PostProcessing.py
index 7bfee9c3473a2e72e11c315c51ee33f0301cf576..5ef451a9c1586cbd22ddb20e5686cd5cb962f641 100644
--- a/BayesValidRox/PostProcessing/PostProcessing.py
+++ b/BayesValidRox/PostProcessing/PostProcessing.py
@@ -377,13 +377,16 @@ class PostProcessing:
         self.RMSE = {}
         # Loop over the keys and compute RMSE error.
         for key in list(ModelOutputs.keys())[1:]:
-            print("\n>>>>> Errors of %s <<<<<"%key)
-            self.RMSE[key] = mean_squared_error(ModelOutputs[key], PCEOutputs[key], squared=False)
-            print("\nRMSE of {0}: {1:.3e}".format(key,self.RMSE[key]))
+            print("\n>>>>> Errors of {} <<<<<".format(key))
+            self.RMSE[key] = mean_squared_error(ModelOutputs[key], PCEOutputs[key], squared=False, multioutput='raw_values')
+            print("\nIndex: RMSE")
+            print('\n'.join('{0}: {1:.3e}'.format(i+1,k) for i,k in enumerate(self.RMSE[key])))
+
+            validErr = mean_squared_error(ModelOutputs[key], PCEOutputs[key], multioutput='raw_values')/ \
+                np.var(ModelOutputs[key],ddof=1, axis=0)
+            print("\nIndex: Validation Error")
+            print('\n'.join('{0}: {1:.3e}'.format(i+1,k) for i,k in enumerate(validErr)))
 
-            validErr = mean_squared_error(ModelOutputs[key], PCEOutputs[key])/np.var(ModelOutputs[key],ddof=1)
-            print("\nvalidErr of {0}: {1:.3e}".format(key,validErr))
-            
     #--------------------------------------------------------------------------------------------------------
     def plotValidation(self, SaveFig=True):
         """