ncf = '/Users/iljagolikov/Documents/ssh/study_Ilja/baka/gor/joonis1.nc'

data = ncread(ncf,'siconc');
lon = ncread(ncf,'longitude');
lat = ncread(ncf,'latitude');
data = transpose(data);

%figure(1)
f = figure(1);
f.Position = [100 100 1200 500];

subplot(1,2,1)
p1 = pcolor(lon,lat,data)
text(0.01, 0.9, '(a)', 'Units', 'normalized', 'FontSize', 24)
shading flat
grid on
%title('Liivi lahe piirkond')
xlabel('Pikkuskraadid [\circE]','fontweight','bold','FontSize', 20) 
ylabel('Laiuskraadid [\circN]','fontweight','bold','FontSize', 20) 
axis([21.5 25 56.5 59.5])
ax = gca;
ax.XAxis.FontSize = 16;
ax.YAxis.FontSize = 16;
ax.LineWidth = 1;


% Joonis 2

ncf2 = '/Users/iljagolikov/Documents/ssh/study_Ilja/baka/gor/joonis2.nc'

data2 = ncread(ncf2,'siconc');
lon2 = ncread(ncf2,'longitude');
lat2 = ncread(ncf2,'latitude');
data2 = transpose(data2);


subplot(1,2,2)
p2 = pcolor(lon2,lat2,data2)
text(0.01, 0.9, '(b)', 'Units', 'normalized', 'FontSize', 24)
shading flat
grid on
%title('Liivi lahe piirkond')
xlabel('Pikkuskraadid [\circE]','fontweight','bold','FontSize', 20) 
ylabel('Laiuskraadid [\circN]','fontweight','bold','FontSize', 20) 
axis([21.5 25 56.5 59.5])
ax = gca;
ax.XAxis.FontSize = 16;
ax.YAxis.FontSize = 16;
ax.LineWidth = 1;
%ax.GridLineWidth = 1;

%% jaa tn

ncf3 = '/Users/iljagolikov/Documents/ssh/study_Ilja/baka/period/prob.nc'

ncdisp(ncf3)

data3 = ncread(ncf3,'siconc');
lon3 = ncread(ncf3,'longitude');
lat3 = ncread(ncf3,'latitude');
data3 = transpose(data3);

figure(2)
pcolor(lon3,lat3,data3)
cb = colorbar();
%cb.Label.FontSize = 12;
cb.Label.String = 'Tõenäosus'
cb.Label.FontWeight = 'bold'
cb.FontSize = 12
shading flat

grid on
%title('Jää tõenäosus')
xlabel('Pikkuskraadid [\circE]','fontweight','bold','FontSize', 12) 
ylabel('Laiuskraadid [\circN]','fontweight','bold','FontSize', 12) 
axis([21.5 25 56.5 59.5])
ax = gca;
ax.XAxis.FontSize = 10;
ax.YAxis.FontSize = 10;
ax.LineWidth = 1;


