% EP_FRAN2

load d:\brian\matlab\working\seawifs_dat 
ds = jul2date(s0n150w(:,1)+4,ones(size(s0n150w(:,1)))*1997);
cs = s0n150w(:,2);
fid = fopen('s0n150w.txt','wt');
fprintf(fid,'%% Col 1: Date\n')
fprintf(fid,'%% Col 2: Chl from SeaWifs at 2N 150W (ug/L)\n')
[Y M D] = datevec(ds);
for i = 1:length(ds)
   fprintf(fid,'%02i/%02i/%04i %8.5f\n',M(i),D(i),Y(i),cs(i));
end
fclose(fid);

ds = jul2date(s2s170w(:,1)+4,ones(size(s2s170w(:,1)))*1997);	% get date of seawifs values
cs = s2s170w(:,2);
fid = fopen('s2s170w.txt','wt');
fprintf(fid,'%% Col 1: Date\n')
fprintf(fid,'%% Col 2: Chl from SeaWifs at 2S 170W (ug/L)\n')

[Y M D] = datevec(ds);
for i = 1:length(ds)
   fprintf(fid,'%02i/%02i/%04i %8.5f\n',M(i),D(i),Y(i),cs(i));
end
fclose(fid);

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);
j = find(isnan(d));
d(j) = 0;
[Y M D] = datevec(d);
Y(j) = 0;M(j) = 0; D(j) = 0;

j = find(isnan(d7));
d7(j) = 0;
[Y7 M7 D7] = datevec(d7);
Y7(j) = 0;M7(j) = 0; D7(j) = 0;

j = find(isnan(d21));
d21(j) = 0;
[Y21 M21 D21] = datevec(d21);
Y21(j) = 0;M21(j) = 0; D21(j) = 0;

fid = fopen('ep1chl.txt','wt');
fprintf(fid,'%% Col 1: Date\n')
fprintf(fid,'%% Col 2: Morel Chl (490 nm) from ep1 (ug/L)\n')
fprintf(fid,'%% Col 3: Date of 7 day average data\n')
fprintf(fid,'%% Col 4: Morel Chl (490 nm) from ep1 (ug/L), 7 day average\n')
fprintf(fid,'%% Col 3: Date of 21 day average data\n')
fprintf(fid,'%% Col 4: Morel Chl (490 nm) from ep1 (ug/L), 21 day average\n')
for i = 1:length(d)
   fprintf(fid,'%02i/%02i/%04i %8.5f %02i/%02i/%04i %8.5f %02i/%02i/%04i %8.5f\n',...
      M(i),D(i),Y(i),c(i),M7(i),D7(i),Y7(i),c7(i),M21(i),D21(i),Y21(i),c21(i));
end
fclose(fid);

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);

c7 = specmean(c,7);
d7 = specmean(d,7);
c21 = specmean(c,21);
d21 = specmean(d,21);

j = find(isnan(d));
d(j) = 0;
[Y M D] = datevec(d);
Y(j) = 0;M(j) = 0; D(j) = 0;

j = find(isnan(d7));
d7(j) = 0;
[Y7 M7 D7] = datevec(d7);
Y7(j) = 0;M7(j) = 0; D7(j) = 0;

j = find(isnan(d21));
d21(j) = 0;
[Y21 M21 D21] = datevec(d21);
Y21(j) = 0;M21(j) = 0; D21(j) = 0;

fid = fopen('ep2chl.txt','wt');
fprintf(fid,'%% Col 1: Date\n')
fprintf(fid,'%% Col 2: Morel Chl (490 nm) from ep2 (ug/L)\n')
fprintf(fid,'%% Col 3: Date of 7 day average data\n')
fprintf(fid,'%% Col 4: Morel Chl (490 nm) from ep2 (ug/L), 7 day average\n')
fprintf(fid,'%% Col 3: Date of 21 day average data\n')
fprintf(fid,'%% Col 4: Morel Chl (490 nm) from ep2 (ug/L), 21 day average\n')
for i = 1:length(d)
   fprintf(fid,'%02i/%02i/%04i %8.5f %02i/%02i/%04i %8.5f %02i/%02i/%04i %8.5f\n',...
      M(i),D(i),Y(i),c(i),M7(i),D7(i),Y7(i),c7(i),M21(i),D21(i),Y21(i),c21(i));
end
fclose(fid);