3rdVersion.m

%Reading the .nc file
FileNC = 'o3_surface_20180701000000.nc' ;

%Printing information of the file in command window
ncinfo(FileNC)
%Longitude
long = ncread(FileNC,'lon');
%Length of longitude
nx = length(long) ;
%Latitude
lat = ncread(FileNC,'lat'); 
%Length of Latitude
ny = length(lat) ;
%Total number of time hours in a day
time = ncread(FileNC,'hour');

%Longitude and latitude for plotting. Converting it to double values
lon=double(long);
lat=double(lat);

%Selecting the data type
disp('Please select which climate model you would like to view: ');
type=input('1- chimere_ozone \n2- emep_ozone \n3- ensemble_ozone \n4- eurad_ozone \n5- lotoseuros_ozone \n6- match_ozone \n7- mocage_ozone \n8- silam_ozone \n');



% user interface - allows the user to select an option to their needs
disp('Please select your desired format: ')
userInput = ('a - Standard Display (Trichromacy)\nb - Monochromacy Display \nc - Protanopia Display  \n');