
 load /eddy/s0/nmiles/marcellus/output_files/processed_July2013/marcellus.mat

% load marcellus.mat

dates=146:7:213;
str_dates=datestr(dates+1,'dd-mmm');
sites=[1 2];
  symbol_dot.site01='.b';
    symbol_dot.site02='.g';

    for i=1:length(sites)  
tmp=num2str(sites(i));
    if length(tmp)==1
        site_short=strcat('site0',tmp);
    else
        site_short=strcat('site',tmp);
    end
    figure(100)
    subplot(2,1,1)
    plot(doy,co2_sites.(site_short)(1:length(co2_sites.site01)),symbol_dot.(site_short));
    set(gca,'FontSize',14)
    set(gca,'FontWeight','bold');
    legend('Site01 all','Site02 all');
    title('CO2 Dry Mole Fraction')
ylabel('[CO2] (ppm)')
set(gca,'Xtick',dates);
set(gca,'XTickLabel',str_dates);
hold on;grid on
    subplot(2,1,2)
    plot(doy,aft_co2_sites.(site_short)(1:length(aft_co2_sites.site01)),symbol_dot.(site_short));
     set(gca,'FontSize',14)
    set(gca,'FontWeight','bold');
    % ax=axis;axis([1 121 370 440]);
legend('Site01 aft','Site02aft');
title('CO2 Dry Mole Fraction')
ylabel('[CO2] (ppm)')
set(gca,'Xtick',dates);
set(gca,'XTickLabel',str_dates);
hold on
    legend('Site01 aft','Site02 aft');
    title('CO2 Dry Mole Fraction')
    grid on;
% print -djpeg co2_time_series
% print -depsc co2_time_series
    
    
    figure(101)
    subplot(2,1,1)
    plot(doy,ch4_sites.(site_short)(1:length(ch4_sites.site01)),symbol_dot.(site_short));
    set(gca,'FontSize',14)
    set(gca,'FontWeight','bold');
    legend('Site01 all','Site02 all');
title('CH4 Dry Mole Fraction')
ylabel('[CH4] (ppb)')
set(gca,'Xtick',dates);
set(gca,'XTickLabel',str_dates);
hold on;grid on;
    subplot(2,1,2)
    plot(doy,aft_ch4_sites.(site_short)(1:length(aft_ch4_sites.site01)),symbol_dot.(site_short));
    % ax=axis;axis([1 121 1800 2150]);
     set(gca,'FontSize',14)
    set(gca,'FontWeight','bold');
hold on
legend('Site01 aft','Site02 aft');
title('CH4 Dry Mole Fraction')
ylabel('[CH4] (ppb)')
set(gca,'Xtick',dates);
set(gca,'XTickLabel',str_dates);
grid on;
% print -djpeg ch4_time_series
% print -depsc ch4_time_series
    
end



figure(200)
plot(doy,aft_co2_sites.site02(1:length(aft_co2_sites.site01))-aft_co2_sites.site01,'.b');
title('Daytime CO2 Difference')
ylabel('CO2_{Site02} - CO2_{Site01} (ppm)')
set(gca,'Xtick',dates);
set(gca,'XTickLabel',str_dates);
grid on

figure(201)
plot(doy,aft_ch4_sites.site02(1:length(aft_ch4_sites.site01))-aft_ch4_sites.site01,'.g');
title('Daytime CH4 Difference')
ylabel('CH4_{Site02} - CH4_{Site01} (ppb)')
set(gca,'Xtick',dates);
set(gca,'XTickLabel',str_dates);
grid on


