CEMCECEM.m

ncfile = "o3_surface_20180701000000.nc";%ncfile chosen
ncinfo(ncfile);
ncdisp(ncfile);

lati = ncread(ncfile, 'lat');
longe = ncread(ncfile, 'lon');

fileDirectory = dir('24HR_CBE_*.csv');%shows the directory

for k =1 : length(fileDirectory) %loop has started
    file = fileDirectory(k).name;  
    matrix = readmatrix(file);
    matrix(400,:) = [0];
    matrix(:,700) = [0];
    showmap = pcolor(longe,lati,matrix); % creates the colour
    showmap.EdgeAlpha = 0;
    hold on;
    p = meshgrid(longe,lati) %this makes longe and lati equal to allow the program to run smoothly
    s = plot(p,'k'); % plots the values
    showmap;
    colorbar
    pause
end