clc;close all;clear all

filename='MOCtot.957-1108.eof.atl.nodt.nc'
PC=ncread(filename,'eofts');
PC1=PC(:,1);


plot(PC1,'Linewidth',2);
set(gca,'Fontsize',16,'Linewidth',2);
xlabel('time(years)')
title('b40.1850.track1.1deg.006a')
set(gcf,'color','w');grid on
ylim([-3 3])

