%% 

function [] = plot_stn_daily(soil_stn, soil_wrf, stn_number)

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


plot(1:92, soil_stn(1,:, stn_number),'k-o');hold on
plot(1:92, soil_wrf(1,:,stn_number,1),'r-x'); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,2)); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,3)); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,4)); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,5)); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,6)); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,7)); hold on; 
% plot(1:92, soil_wrf(1,:,stn_number,8),'r-s'); hold on; 

% plot(1:92, soil_wrf(1,:,stn_number,9),'b-x'); hold on; 
plot(1:92, soil_wrf(1,:,stn_number,10),'b-s'); hold on; 



% legend('STAMP','Runoption1', 'Runoption2', 'Runoption3', 'Runoption4',...
%     'Runoption5','Runoption6','Runoption7','Runoption8')
legend('STAMP','Runoption1', 'Runoption06')

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

ylabel('soil moisture (m^3 m^{-3})')
xlabel('day')
xlim([1 92])  

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