diff --git a/4thVersion.m b/4thVersion.m
new file mode 100644
index 0000000..7d64832
--- /dev/null
+++ b/4thVersion.m
@@ -0,0 +1,45 @@
+%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');