clc;clear all;close all

file='HadISSTyr_sst.nc'

sst=ncread(file,'sst');
lon=ncread(file,'longitude');
lat=ncread(file,'latitude');

lon_loc=lon>=-80 &lon<=-15;
lat_loc=lat>=45&lat<=65;

N_SST=squeeze(nanmean(nanmean(sst(lon_loc,lat_loc,:))));
time=linspace(1870,2021,151);



set(0,'defaultfigurecolor','w')
plot(time(110:end),N_SST(110:end),'Linewidth',3);
title('North Atlantic HadlSST')
xlabel('year');ylabel('degree')
set(gca,'XMinorTick','on','YMinorTick','on','TickDir','out');
ytickformat('%.2f');set(gca,'fontsize',16,'linewidth',2);grid on;











