Model_Ensemble.m

ok=input("input 1 to continue and 0 to exit:");
while ok==1
    if ok==0 
       break
    end
%%
figure(1)
title('Meteorologic Model: Match')
    y = double(ncread('Individual_Models\W_fr-meteofrance,MODEL,ENSEMBLE+FORECAST+SURFACE+O3+0H24H_C_LFPW_20180701000000.nc','longitude'));
    zlabel('ozone')
    x = double(ncread('Individual_Models\W_fr-meteofrance,MODEL,ENSEMBLE+FORECAST+SURFACE+O3+0H24H_C_LFPW_20180701000000.nc','latitude'));
    time = ncread('Individual_Models\W_fr-meteofrance,MODEL,ENSEMBLE+FORECAST+SURFACE+O3+0H24H_C_LFPW_20180701000000.nc','time');
    ozone = ncread('Individual_Models\W_fr-meteofrance,MODEL,ENSEMBLE+FORECAST+SURFACE+O3+0H24H_C_LFPW_20180701000000.nc','unknown');
    h=input("input an integer between 0 and 26 representing 1 of the 25 hours:");
    clf
    z = ozone(:,:,h);
    mesh(x,y,z)
    
    
[x,y] = meshgrid(x, y);
%%
figure(2)
    clf
    % Create the map
    worldmap('Europe'); % set the part of the earth to show
    load coastlines
    plotm(coastlat,coastlon)

    %bordersm('countries','FaceColor', [0.5 0.7 0.5])
    colorbar
    land = shaperead('landareas', 'UseGeoCoords', true);
    geoshow(gca, land, 'FaceColor', [0.5 0.7 0.5])
    lakes = shaperead('worldlakes', 'UseGeoCoords', true);
    geoshow(lakes, 'FaceColor', 'blue', 'EdgeColor', 'cyan')
    rivers = shaperead('worldrivers', 'UseGeoCoords', true);
    geoshow(rivers, 'Color', 'blue')
    cities = shaperead('worldcities', 'UseGeoCoords', true);
    geoshow(cities, 'Marker', '.', 'Color', 'magenta')
    setm(gca,'ffacecolor','#0072BD')
    surfm(x,y,z, 'EdgeColor', 'none','FaceAlpha', 0.6) % edge colour outlines the edges, 'FaceAlpha', sets the transparency
%%    
figure(3)
    clf
    worldmap('Europe'); % set the part of the earth to show
    load coastlines
    plotm(coastlat,coastlon)

    %bordersm('countries','FaceColor', [0.5 0.7 0.5])
    setm(gca,'ffacecolor','#0072BD')
    land = shaperead('landareas', 'UseGeoCoords', true);
    geoshow(gca, land, 'FaceColor', [0.5 0.7 0.5])
    lakes = shaperead('worldlakes', 'UseGeoCoords', true);
    geoshow(lakes, 'FaceColor', 'blue', 'EdgeColor', 'cyan')
    rivers = shaperead('worldrivers', 'UseGeoCoords', true);
    geoshow(rivers, 'Color', 'blue')
    cities = shaperead('worldcities', 'UseGeoCoords', true);
    geoshow(cities, 'Marker', '.', 'Color', 'magenta')

    NumContours = 8;
    contourfm(x, y, z, NumContours)
    Plots = findobj(gca,'Type','Axes');
    Plots.SortMethod = 'depth';
    clear
    ok=input("input 1 to continue and 0 to exit:");
end