4thVersion.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');

%Latitude
lat = ncread(FileNC,'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');

%Data having all the data in it
if type==1
Data = ncread(FileNC,'chimere_ozone') ;
elseif type==2
Data = ncread(FileNC,'emep_ozone') ;    
elseif type==3
Data = ncread(FileNC,'ensemble_ozone') ;    
elseif type==4
Data = ncread(FileNC,'eurad_ozone') ;    
elseif type==5
Data = ncread(FileNC,'lotoseuros_ozone') ;    
elseif type==6
Data = ncread(FileNC,'match_ozone') ;    
elseif type==7
Data = ncread(FileNC,'mocage_ozone') ;    
elseif type==8
Data = ncread(FileNC,'silam_ozone') ;   
    
end

% 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');