diff --git a/tools/MATLAB scripts/outputfiles_merge.m b/tools/MATLAB scripts/outputfiles_merge.m
deleted file mode 100644
index 5523f65a..00000000
--- a/tools/MATLAB scripts/outputfiles_merge.m	
+++ /dev/null
@@ -1,56 +0,0 @@
-% outputfiles_merge.m
-% Script to merge gprMax output files of A-scans (traces) into a single
-% HDF5 file
-%
-% Craig Warren
-
-clear all, clc
-
-[filenames, pathname] = uigetfile('*.out', 'Select gprMax A-scan output files to merge', 'MultiSelect', 'on');
-
-% Combine A-scans (traces) into a single HDF5 output file
-if filenames{1} ~= 0
-    [pathstr, basefilename, ext] = fileparts(filenames{1});
-    outputfile = strcat(fullfile(pathname, basefilename(1:end-1)), '_all.out');
-    modelruns = length(filenames);
-    filename = fullfile(pathname, filenames{1});
-    iterations = double(h5readatt(filename, '/', 'Iterations'));
-    dt = h5readatt(filename, '/', 'dt');
-    h5create(outputfile, '/rxs/rx1/Ex', [modelruns iterations]);
-    h5create(outputfile, '/rxs/rx1/Ey', [modelruns iterations]);
-    h5create(outputfile, '/rxs/rx1/Ez', [modelruns iterations]);
-    h5create(outputfile, '/rxs/rx1/Hx', [modelruns iterations]);
-    h5create(outputfile, '/rxs/rx1/Hy', [modelruns iterations]);
-    h5create(outputfile, '/rxs/rx1/Hz', [modelruns iterations]);
-    h5writeatt(outputfile, '/', 'Iterations', iterations);
-    h5writeatt(outputfile, '/', 'dt', dt);
-    Ex = zeros(iterations, modelruns);
-    Ey = zeros(iterations, modelruns);
-    Ez = zeros(iterations, modelruns);
-    Hx = zeros(iterations, modelruns);
-    Hy = zeros(iterations, modelruns);
-    Hz = zeros(iterations, modelruns);
-    for rx=1:modelruns
-        filename = fullfile(pathname, filenames{rx});
-        Ex(:, rx) = h5read(filename, '/rxs/rx1/Ex');
-        Ey(:, rx) = h5read(filename, '/rxs/rx1/Ey');
-        Ez(:, rx) = h5read(filename, '/rxs/rx1/Ez');
-        Hx(:, rx) = h5read(filename, '/rxs/rx1/Hx');
-        Hy(:, rx) = h5read(filename, '/rxs/rx1/Hy');
-        Hz(:, rx) = h5read(filename, '/rxs/rx1/Hz');
-    end
-    h5write(outputfile, '/rxs/rx1/Ex', Ex');
-    h5write(outputfile, '/rxs/rx1/Ey', Ey');
-    h5write(outputfile, '/rxs/rx1/Ez', Ez');
-    h5write(outputfile, '/rxs/rx1/Hx', Hx');
-    h5write(outputfile, '/rxs/rx1/Hy', Hy');
-    h5write(outputfile, '/rxs/rx1/Hz', Hz');
-    prompt = 'Do you want to remove the multiple individual output files? [y] or n: ';
-    check = input(prompt,'s');
-    if isempty(check) || check == 'y'
-        for f=1:length(filenames)
-            filename = fullfile(pathname, filenames{f});
-            delete(filename);
-        end
-    end
-end
\ No newline at end of file
diff --git a/tools/MATLAB scripts/plot_Ascan.m b/tools/MATLAB scripts/plot_Ascan.m
deleted file mode 100644
index 9b22b5ed..00000000
--- a/tools/MATLAB scripts/plot_Ascan.m	
+++ /dev/null
@@ -1,89 +0,0 @@
-% plot_Ascan.m
-% Script to save and plot EM fields from a gprMax A-scan
-%
-% Craig Warren
-
-clear all, clc
-
-[filename, pathname] = uigetfile('*.out', 'Select gprMax A-scan output file to plot');
-fullfilename = strcat(pathname, filename);
-
-if filename ~= 0
-    header.title = h5readatt(fullfilename, '/', 'Title');
-    header.iterations = double(h5readatt(fullfilename,'/', 'Iterations'));
-    tmp = h5readatt(fullfilename, '/', 'dx, dy, dz');
-    header.dx = tmp(1);
-    header.dy = tmp(2);
-    header.dz = tmp(3);
-    header.dt = h5readatt(fullfilename, '/', 'dt');
-    header.nsrc = h5readatt(fullfilename, '/', 'nsrc');
-    header.nrx = h5readatt(fullfilename, '/', 'nrx');
-    
-    % Time vector for plotting
-    time = linspace(0, (header.iterations)*(header.dt)*1E9, header.iterations);
-    
-    % Initialise structure for field arrays
-    fields.ex = zeros(header.iterations, header.nrx);
-    fields.ey = zeros(header.iterations, header.nrx);
-    fields.ez = zeros(header.iterations, header.nrx);
-    fields.hx = zeros(header.iterations, header.nrx);
-    fields.hy = zeros(header.iterations, header.nrx);
-    fields.hz = zeros(header.iterations, header.nrx);
-    
-    prompt = 'Would you like to save a PDF of the plot(s)? y or [n]: ';
-    pdf = input(prompt,'s');
-    
-    % Save and plot fields from each receiver
-    for n=1:header.nrx
-        path = strcat('/rxs/rx', num2str(n));
-        tmp = h5readatt(fullfilename, path, 'Position');
-        header.rx(n) = tmp(1);
-        header.ry(n) = tmp(2);
-        header.rz(n) = tmp(3);
-        path = strcat(path, '/');
-        fields.ex(:,n) = h5read(fullfilename, strcat(path, 'Ex'));
-        fields.ey(:,n) = h5read(fullfilename, strcat(path, 'Ey'));
-        fields.ez(:,n) = h5read(fullfilename, strcat(path, 'Ez'));
-        fields.hx(:,n) = h5read(fullfilename, strcat(path, 'Hx'));
-        fields.hy(:,n) = h5read(fullfilename, strcat(path, 'Hy'));
-        fields.hz(:,n) = h5read(fullfilename, strcat(path, 'Hz'));
-        
-        fh1=figure('Name', strcat('rx', num2str(n)));
-        ax1 = subplot(2,3,1); plot(time, fields.ex(:,n), 'r', 'LineWidth', 2), grid on, xlabel('Time [ns]'), ylabel('Field strength [V/m]'), title('E_x')
-        ax2 = subplot(2,3,2); plot(time, fields.ey(:,n), 'r', 'LineWidth', 2), grid on, xlabel('Time [ns]'), ylabel('Field strength [V/m]'), title('E_y')
-        ax3 = subplot(2,3,3); plot(time, fields.ez(:,n), 'r', 'LineWidth', 2), grid on, xlabel('Time [ns]'), ylabel('Field strength [V/m]'), title('E_z')
-        ax4 = subplot(2,3,4); plot(time, fields.hx(:,n), 'b', 'LineWidth', 2), grid on, xlabel('Time [ns]'), ylabel('Field strength [A/m]'), title('H_x')
-        ax5 = subplot(2,3,5); plot(time, fields.hy(:,n), 'b', 'LineWidth', 2), grid on, xlabel('Time [ns]'), ylabel('Field strength [A/m]'), title('H_y')
-        ax6 = subplot(2,3,6); plot(time, fields.hz(:,n), 'b', 'LineWidth', 2), grid on, xlabel('Time [ns]'), ylabel('Field strength [A/m]'), title('H_z')
-        ax1.FontSize = 16;
-        ax2.FontSize = ax1.FontSize;
-        ax3.FontSize = ax1.FontSize;
-        ax4.FontSize = ax1.FontSize;
-        ax5.FontSize = ax1.FontSize;
-        ax6.FontSize = ax1.FontSize;
-        
-        % Options to create a nice looking figure for display and printing
-        set(fh1,'Color','white','Menubar','none');
-        X = 60;   % Paper size
-        Y = 30;   % Paper size
-        xMargin = 0; % Left/right margins from page borders
-        yMargin = 0;  % Bottom/top margins from page borders
-        xSize = X - 2*xMargin;    % Figure size on paper (width & height)
-        ySize = Y - 2*yMargin;    % Figure size on paper (width & height)
-        
-        % Figure size displayed on screen
-        set(fh1, 'Units','centimeters', 'Position', [0 0 xSize ySize])
-        movegui(fh1, 'center')
-        
-        % Figure size printed on paper
-        set(fh1,'PaperUnits', 'centimeters')
-        set(fh1,'PaperSize', [X Y])
-        set(fh1,'PaperPosition', [xMargin yMargin xSize ySize])
-        set(fh1,'PaperOrientation', 'portrait')
-        
-        % Export to PDF
-        if pdf == 'y'
-            print(fh1, '-dpdf', '-r0', strcat(fullfilename(1:end-4), '_rx', num2str(n), '.pdf'));
-        end
-    end
-end
\ No newline at end of file
diff --git a/tools/MATLAB scripts/plot_Bscan.m b/tools/MATLAB scripts/plot_Bscan.m
deleted file mode 100644
index f100a7a1..00000000
--- a/tools/MATLAB scripts/plot_Bscan.m	
+++ /dev/null
@@ -1,65 +0,0 @@
-% plot_Bscan.m
-% Script to plot EM fields from a gprMax B-scan
-%
-% Craig Warren
-
-clear all, clc
-
-[filename, pathname] = uigetfile('*.out', 'Select gprMax output file to plot B-scan', 'MultiSelect', 'on');
-filename = fullfile(pathname, filename);
-
-% Open file and read fields
-if filename ~= 0
-    iterations = double(h5readatt(filename, '/', 'Iterations'));
-    dt = h5readatt(filename, '/', 'dt');
-    Ex = h5read(filename, '/rxs/rx1/Ex')';
-    Ey = h5read(filename, '/rxs/rx1/Ey')';
-    Ez = h5read(filename, '/rxs/rx1/Ez')';
-    Hx = h5read(filename, '/rxs/rx1/Hx')';
-    Hy = h5read(filename, '/rxs/rx1/Hy')';
-    Hz = h5read(filename, '/rxs/rx1/Hz')';
-end
-
-prompt = 'Would you like to save a PDF of the plot? y or [n]: ';
-pdf = input(prompt,'s');
-
-prompt = 'Which field do you want to view? Ex, Ey, or Ez: ';
-field = input(prompt,'s');
-field = eval(field);
-time = 0:dt:iterations * dt;
-traces = 0:size(field,2);
-
-fh1=figure('Name', filename);
-clims = [-max(max(abs(field))) max(max(abs(field)))];
-im = imagesc(traces, time, field, clims);
-xlabel('Trace number');
-ylabel('Time [s]');
-c = colorbar;
-c.Label.String = 'Field strength [V/m]';
-ax = gca;
-ax.FontSize = 16;
-xlim([0 traces(end)]);
-
-% Options to create a nice looking figure for display and printing
-set(fh1,'Color','white','Menubar','none');
-X = 60;   % Paper size
-Y = 30;   % Paper size
-xMargin = 0; % Left/right margins from page borders
-yMargin = 0;  % Bottom/top margins from page borders
-xSize = X - 2*xMargin;    % Figure size on paper (width & height)
-ySize = Y - 2*yMargin;    % Figure size on paper (width & height)
-
-% Figure size displayed on screen
-set(fh1, 'Units','centimeters', 'Position', [0 0 xSize ySize])
-movegui(fh1, 'center')
-
-% Figure size printed on paper
-set(fh1,'PaperUnits', 'centimeters')
-set(fh1,'PaperSize', [X Y])
-set(fh1,'PaperPosition', [xMargin yMargin xSize ySize])
-set(fh1,'PaperOrientation', 'portrait')
-
-% Export to PDF
-if pdf == 'y'
-    print(fh1, '-dpdf', '-r0', strcat(filename(1:end-4), '.pdf'));
-end