% EP_FRAN1   - Plot of CHl vs seawifs for Francisco

% Brian Schlining
% 01 Sep 1999

% loads s0n150w and s2s170w. The first column of each variable is days since 1 Jan 1997.
% Note the dat is in 89 day bins with the date stored as the beginning of the bin. Add 4 to 
% get the center of the bin
% The second column is [chl] in ug/L
load d:\brian\matlab\working\seawifs_dat 
ds = jul2date(s0n150w(:,1)+4,ones(size(s0n150w(:,1)))*1997);
cs = s0n150w(:,2);

y97 = jul2date(1,1997);
y98 = jul2date(1,1998);
y99 = jul2date(1,1999);

load d:\brian\matlab\working\ep1hflf   % load d and dat
c      = morel_chl(dat(:,4),dat(:,32),20,490);
bad    = [173:175 586:676 887:890];		% Remove bad data
c(957) = (c(958)+c(956))./2;				% Patch a small gap
d(bad) = NaN;
c7 = specmean(c,7);
d7 = specmean(d,7);
c21 = specmean(c,21);
d21 = specmean(d,21);
plot(d,c,'g')
hold on
plot(d7,c7,'b');
plot(d21,c21,'r')
%plot(ds,cs,'r')
plot(ds,cs,'k.')
plot([y97 y97], [0 1.5],'k:')
plot([y98 y98], [0 1.5],'k:')
plot([y99 y99], [0 1.5],'k:')
hold off
datetick('x',4)
set(gca,'YLim',[0 1.5],'YTick',[0:0.1:1.5],'YTickLabel',strvcat('0',' ',' ',' ',' ','0.5',' ',' ',' ',' ','1.0',' ',' ',' ',' ','1.5'))
title('CHL at 0, 155W');
legend('Daily Avg','7day Avg','21day Avg','SeaWiFS',1)
tp = jul2date([182 182 182],[1997:1999]);
th = text(tp,ones(size(tp))*-0.12,['1997';'1998';'1999']);
set(th,'HorizontalAlignment','Center');

ylabel('CHL (\mug L^-^1');

set(gcf,'PaperPosition',[0 0 5 4])
writegif('d:\brian\matlab\working\ep1seawif_chl.gif')
xlim = get(gca,'XLim');

figure
ds = jul2date(s2s170w(:,1)+4,ones(size(s2s170w(:,1)))*1997);	% get date of seawifs values
cs = s2s170w(:,2);					
load d:\brian\matlab\working\ep2hflf   % load d and dat
i = find(isnan(dat(:,39)));				% If no depth is found set it to 21 m
dat(i,39) = 21;
bad    = [1:631];		% Remove bad data
d(bad) = NaN;
bad = find(ds < datenum('01 Nov 1998'));
ds(bad) = NaN;
%c      = morel_chl(dat(:,4),dat(:,32),dat(:,39),490);
c      = morel_chl(dat(:,4),dat(:,32),20,490);

%c      = morel_chl(dat(:,5),dat(:,33),20,510);

c7 = specmean(c,7);
d7 = specmean(d,7);
c21 = specmean(c,21);
d21 = specmean(d,21);
plot(d,c,'g')
hold on
plot(d7,c7,'b');
plot(d21,c21,'r')
plot(ds,cs,'k.')
%plot([y97 y97], [0 1.5],'k:')
%plot([y98 y98], [0 1.5],'k:')
plot([y99 y99], [0 1.5],'k:')

hold off
datetick('x',4)
set(gca,'YLim',[0 0.5],'YTick',[0:0.1:1.5],...
   'YTickLabel',strvcat('0',' ',' ',' ',' ','0.5',' ',' ',' ',' ','1.0',' ',' ',' ',' ','1.5'),...
   'XLim',[datenum('01 Nov 1998') xlim(2)])
tH = title('CHL at 2S, 170W');
xlabel('1999')
ylabel('CHL (\mug L^-^1');
legend('Daily Avg','7day Avg','21day Avg','SeaWiFS',1)
set(gcf,'PaperPosition',[0 0 5 4])
writegif('d:\brian\matlab\working\ep2seawif_chl.gif')