Update plot_blocking_diag.m with mpfa

geiger/results/plot_blocking_diag.m

@@ -15,6 +15,11 @@ plot(dfm(:, 15), dfm(:, 4), 'Color', 'g', 'LineWidth', 1);
 ccdfm = csvread('ccdfm/ccdfm_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
 plot( ccdfm(:, 3), ccdfm(:, 1), 'Color', 'c', 'LineWidth', 1);
 
+% MPFA
+ccdfm = csvread('ccdfm/mpfa_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
+plot( ccdfm(:, 3), ccdfm(:, 1), 'Color', [0 .5 0], 'LineWidth', 1);
+
+
 % EDFM
 edfm = csvread('edfm/edfm_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
 plot( edfm(:, 3), edfm(:, 1), 'Color', 'r', 'LineWidth', 1);
@@ -36,4 +41,4 @@ xlabel('arc length', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
 ylabel('pressure', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
 xlim([0 1.2728]);
 set(gca, 'FontSize', 12);
-legend('reference', 'Box-DFM', 'CC-DFM', 'EDFM', 'mortar-DFM', 'P-XFEM', 'D-XFEM', 'location', 'northeastoutside');
+legend('reference', 'Box-DFM', 'TPFA', 'MPFA', 'EDFM', 'flux-mortar', 'P-XFEM', 'D-XFEM', 'location', 'northeastoutside');