From d9f9d2fe8d2c01d4e591306c1cfe309f16f05764 Mon Sep 17 00:00:00 2001 From: "Faiza Bibi (bibif5)" Date: Sun, 2 Aug 2020 19:10:24 +0100 Subject: [PATCH] 4thVersion.m This has if-statements added for first user inputs --- 4thVersion.m | 45 +++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 45 insertions(+) create mode 100644 4thVersion.m 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');