overlay data to map

Visualisation.m

@@ -77,3 +77,20 @@ elseif colourblind == 3
 end
 
 clc
+
+disp('The map in the figure will load shortly ..')
+for k = 1 : length(h)
+    clf; 
+    matrix = pcolor(X,Y,ozoneChoice(:,:,k)');
+    matrix.EdgeAlpha = 0;
+    matrix.FaceAlpha = 0.8;
+    load coast 
+    cities = shaperead('worldcities', 'UseGeoCoords', true);
+    bar = colorbar; 
+    bar.Label.String = 'Ozone Level';
+    geoshow(cities, 'Marker', '+', 'Color', 'red') 
+    hold on;
+    plot(long,lat,'k')
+    title(sprintf('Ensemble Data on Europe'))
+    pause(0.3)
+end