Matlab plot -- annual lake wide average surface temperature

clc;clear;

yr=2019;
%% import observation data
filename='/Users/chuyanzhao/Desktop/gls_observation_data/obs_glsea_SWT/glsea-temps2019_1024.dat';
tmp=importdata(filename);
temp_obs=tmp.data;
temp_obs_superior=temp_obs(:,3);


%% import fvcom data
ncload(['/Users/chuyanzhao/Desktop/FVCOM_output/Lake_Superior/2019/raw_data/gls_lst_2019_LowRes.nc'],'temp','art1');

clear lst
temp_mod=squeeze(temp(:,1,:));
[nt,nn]=size(temp_mod);

for t=1:nt
    temp_mod_avg(t)=sum(temp_mod(t,:).*art1')./sum(art1);% art1 is the area of each element
end

temp_mod_superior=mean_by_segments(temp_mod_avg',[0:24:8760],1);


%% plot
figure(101)
ff1=plot(1:365,temp_obs_superior,'Color','r','LineWidth',1)
hold on;
ff2=plot(1:365,temp_mod_superior,'Color','b','LineWidth',1)
set(gca,'FontSize',18,'FontName','Times New Roman','tickdir','out');

xlabel(['Julian day (' num2str(yr) ')']);
ylabel('LST (\circC)');
axis([0 365 -5 20]);
legend([ff1,ff2],'obs','fvcom');

saveas(gcf,['lake-wide_avg_lst_2019.png'],'png');