close all;
clear all;
clc;
hold on;

% equi-dimensional MFD
mfd = csvread('mfd/mfd_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(mfd(:, 6), mfd(:, 1), 'Color', 'k', 'LineWidth', 1);

% Box-DFM
dfm = csvread('boxdfm/boxdfm_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(dfm(:, 14), dfm(:, 4), 'Color', 'g', 'LineWidth', 1);

% CC-DFM
ccdfm = csvread('ccdfm/ccdfm_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(ccdfm(:, 4), ccdfm(:, 1), 'Color', 'c', 'LineWidth', 1);

% MPFA
mpfa = csvread('ccdfm/mpfa_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(mpfa(:, 4), mpfa(:, 1), 'Color', [0 .5 0], 'LineWidth', 1);

% EDFM
edfm = csvread('edfm/edfm_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(edfm(:, 5), edfm(:, 1), 'Color', 'r', 'LineWidth', 1);

% Mortar-DFM
mortar = csvread('mortardfm/mortardfm_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(mortar(:, 4), mortar(:, 1), 'Color', [1 0.7 0], 'LineWidth', 1);

% P-XFEM
p_xfem = csvread('pxfem/pxfem_geiger_conductive_x_eq_0.5.csv', 1, 0);
plot(p_xfem(:, 3), p_xfem(:, 1), 'Color', 'b', 'LineWidth', 1);

% D-XFEM
d_xfem_fine = csvread('dxfem/dxfem_geiger_conductive_x_eq_0.5.csv', 1, 0);
[blabla,I]=sort(d_xfem_fine(:, 3));
plot(d_xfem_fine(I, 3), d_xfem_fine(I, 1), 'Color', 'm', 'LineWidth', 1);

% format the plot
xlabel('arc length', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
ylabel('pressure', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
set(gca, 'FontSize', 12);
legend('Reference', 'Box', 'TPFA', 'MPFA', 'EDFM', 'Flux-Mortar', 'P-XFEM', 'D-XFEM');