%% 

function [] = plot_stn_daily(lh_stn, lh_wrf)

figure
% stn_list = {'E11','E21','E12', 'E13','E15','E31', 'E32', 'E33'...
%             'E34','E35','E36','E37','E38','E39','E40', 'E9'};

ndays = size(lh_stn,1);
lh_stn_mean = nanmean(lh_stn, 2);
lh_wrf_mean = squeeze(nanmean(lh_wrf, 2));

plot(1:ndays, lh_stn_mean,'k-o');hold on
% plot(1:ndays, lh_wrf_mean(1,:),'r-x'); hold on; 
% plot(1:ndays, lh_wrf_mean(2,:),'r-s'); hold on; 
plot(1:ndays, lh_wrf_mean(3,:),'b-x'); hold on; 
plot(1:ndays, lh_wrf_mean(4,:),'r-s'); hold on; 
plot(1:ndays, lh_wrf_mean(5,:),'g-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(6,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(7,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(8,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(9,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(10,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(11,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(12,:),'b-s'); hold on; 
% plot(1:ndays, lh_wrf_mean(13,:),'b-s'); hold on; 



% legend('STAMP','Runoption1', 'Runoption2', 'Runoption3', 'Runoption4',...
%     'Runoption5','Runoption6','Runoption7','Runoption8')
% legend('Best estimate','Runoption1', 'Runoption2', 'Runoption3', 'Runoption4',...
%        'Runoption5','Runoption6','Runoption7','Runoption8',...
%        'Runoption9','Runoption10','Runoption11','Runoption12','Runoption13');

legend('Best estimate','Runoption1', 'Runoption4',...
       'Runoption5');


%  title (['station ', char(stn_list(stn_number)) ])

ylabel('LH (W m^{-2})')
xlabel('day')
xlim([1 ndays])

% print('-depsc', ['station_number_1_8_9_10_',num2str(stn_number)])
end