From dc31c808030a427fdd04ca761c5d5e2661492e3d Mon Sep 17 00:00:00 2001
From: "Jibraan Keyani (keyanij)" <keyanij@coventry.ac.uk>
Date: Tue, 28 Jul 2020 12:46:58 +0100
Subject: [PATCH] Add files via upload

---
 OzoneVisualisation5011CEM.m | 76 +++++++++++++++++++++++++++++++++++
 Untitled.asv                | 79 +++++++++++++++++++++++++++++++++++++
 2 files changed, 155 insertions(+)
 create mode 100644 OzoneVisualisation5011CEM.m
 create mode 100644 Untitled.asv

diff --git a/OzoneVisualisation5011CEM.m b/OzoneVisualisation5011CEM.m
new file mode 100644
index 0000000..5bedfb3
--- /dev/null
+++ b/OzoneVisualisation5011CEM.m
@@ -0,0 +1,76 @@
+%%
+ncfile = 'o3_surface_20180701000000.nc';
+
+ncinfo(ncfile) %- gives info on nc file
+ncdisp(ncfile) % displays nc file
+
+%%select an ozone map 
+disp('1: emep ozone map')
+disp('2: chimere ozone map')
+disp('3: ozone map')
+disp('4: lotoseuros ozone map')
+disp('5: eurad ozone map')
+disp('6: ensemble ozone map')
+disp('7: silam ozone map')
+disp('8: mocage ozone map')
+selected_map = input('Select ozone map')
+if selected_map == 1
+    select = 'emep_ozone'
+elseif selected_map == 2
+    select = 'chimere_ozone'
+elseif selected_map == 3
+    select = 'match_ozone'
+elseif selected_map == 4
+    select = 'lotoseuros_ozone'
+elseif selected_map == 5
+    select = 'eurad_ozone'
+elseif selected_map == 6
+    select = 'ensemble_ozone'
+elseif selected_map == 7
+    select = 'silam_ozone'
+elseif selected_map == 8
+    select = 'mocage_ozone'
+end
+X = ncread(ncfile, 'lat');
+Y = ncread(ncfile, 'lon');
+time = ncread(ncfile, 'hour'); 
+selection = ncread(ncfile, select);
+chimere_ozone = ncread(ncfile, "chimere_ozone");
+
+%colour blindness selection
+disp('Do you require colourblind assistance?')
+cBlind = input('input 1 for yes or 2 for no ')
+if cBlind == 1
+    colormap bone;
+elseif cBlind == 2
+    colormap default;
+end
+
+%% main loop
+for i = 1:length(time) %iterates through 1 hour intervals
+    Z = chimere_ozone(:,:,i);
+    % Create the map
+    worldmap('Europe'); % set the part of the earth to show
+    
+    load coastlines %plots coastlines on the map
+    plotm(coastlat,coastlon)
+
+    land = shaperead('landareas', 'UseGeoCoords', true);
+    geoshow(gca, land, 'FaceColor', [0.5 0.7 0.5])
+
+    lakes = shaperead('worldlakes', 'UseGeoCoords', true);
+    geoshow(lakes, 'FaceColor', 'blue')
+
+    rivers = shaperead('worldrivers', 'UseGeoCoords', true);
+    geoshow(rivers, 'Color', 'blue')
+
+    cities = shaperead('worldcities', 'UseGeoCoords', true);
+    geoshow(cities, 'Marker', '.', 'Color', 'red')
+
+    % Plot the data
+    surfm(double(X), double(Y), double(Z), 'EdgeColor', 'none',...
+        'FaceAlpha', 0.5) % edge colour outlines the edges, 'FaceAlpha', sets the transparency
+    title(sprintf('Current time:  %.0f', time(i)));
+    pause(0.05)
+    
+end
diff --git a/Untitled.asv b/Untitled.asv
new file mode 100644
index 0000000..9c29f5e
--- /dev/null
+++ b/Untitled.asv
@@ -0,0 +1,79 @@
+%%
+ncfile = 'o3_surface_20180701000000.nc';
+
+&ncinfo(fileName) %- gives info on nc file
+ncdisp(fileName) % displays nc file
+
+%%select a ozone map of your choice
+disp('1: emep ozone map')
+disp('2: chimere ozone map')
+disp('3: ozone map')
+disp('4: lotoseuros ozone map')
+disp('5: eurad ozone map')
+disp('6: ensemble ozone map')
+disp('7: silam ozone map')
+disp('8: mocage ozone map')
+selected_map = input('Select ozone map')
+if selected_map == 1
+    select = 'emep_ozone'
+elseif selected_map == 2
+    select = 'chimere_ozone'
+elseif selected_map == 3
+    select = 'match_ozone'
+elseif selected_map == 4
+    select = 'lotoseuros_ozone'
+elseif selected_map == 5
+    select = 'eurad_ozone'
+elseif selected_map == 6
+    select = 'ensemble_ozone'
+elseif selected_map == 7
+    select = 'silam_ozone'
+elseif selected_map == 8
+    select = 'mocage_ozone'
+end
+X = ncread(ncfile, 'lat');
+Y = ncread(ncfile, 'lon');
+time = ncread(ncfile, 'hour'); 
+selection = ncread(ncfile, select);
+chimere_ozone = ncread(ncfile, "chimere_ozone");
+
+%colour blindness selection
+disp('Do you require colourblind assistance?')
+cBlind = input('input 1 for yes or 2 for no')
+if cBlind == 1
+    colormap bone;
+elseif cBlind == 2
+    colormap default;
+end
+
+%% Create a display of the data from the NetCDF files like this
+figure(2);
+clf
+%
+for i = 1:length(time)
+    Z = chimere_ozone(:,:,i);
+    % Create the map
+    worldmap('Europe'); % set the part of the earth to show
+    
+    load coastlines
+    plotm(coastlat,coastlon)
+
+    land = shaperead('landareas', 'UseGeoCoords', true);
+    geoshow(gca, land, 'FaceColor', [0.5 0.7 0.5])
+
+    lakes = shaperead('worldlakes', 'UseGeoCoords', true);
+    geoshow(lakes, 'FaceColor', 'blue')
+
+    rivers = shaperead('worldrivers', 'UseGeoCoords', true);
+    geoshow(rivers, 'Color', 'blue')
+
+    cities = shaperead('worldcities', 'UseGeoCoords', true);
+    geoshow(cities, 'Marker', '.', 'Color', 'red')
+
+    % Plot the data
+    surfm(double(X), double(Y), double(Z), 'EdgeColor', 'none',...
+        'FaceAlpha', 0.5) % edge colour outlines the edges, 'FaceAlpha', sets the transparency
+    title(sprintf('Current time:  %.0f', time(i)));
+    pause(0.05)
+    
+end