function ep_chl_plot(NumDays,ext)
% Brian Schlining
% 27 Apr 1999
if nargin < 2
   ext = '';
end

PaperPosition = [.25 2.5 4.5 3.25];

DateStart = now-NumDays;
DateEnd   = now;
%===================================
% Create plots of standard CHL products
%===================================
YLim = [0 1];
ep_plot(1,62,DateStart,DateEnd); % Morel[0m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['morel1_0' ext]);
close all
ep_plot(1,63,DateStart,DateEnd); % OC2V2[0m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['oc2v21_0' ext]);
close all
ep_plot(1,64,DateStart,DateEnd); % OC2V2[20m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['oc2v21_20' ext]);
close all
ep_plot(1,65,DateStart,DateEnd); % lu443/lu555[0m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['rat1_0' ext]);
close all
ep_plot(1,66,DateStart,DateEnd); % lu443/lu555[20m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['rat1_20' ext]);
close all
ep_plot(1,69,DateStart,DateEnd); % CZCS
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['czcs1_0' ext]);
close all

ep_plot(2,62,DateStart,DateEnd); % Morel[0m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['morel2_0' ext]);
close all
ep_plot(2,63,DateStart,DateEnd); % OC2V2[0m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['oc2v22_0' ext]);
close all
ep_plot(2,64,DateStart,DateEnd); % OC2V2[20m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['oc2v22_20' ext]);
close all
ep_plot(2,65,DateStart,DateEnd); % lu443/lu555[0m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['rat2_0' ext]);
close all
ep_plot(2,66,DateStart,DateEnd); % lu443/lu555[20m]
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['rat2_20' ext]);
close all
ep_plot(2,69,DateStart,DateEnd); % CZCS
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['czcs2_0' ext]);
close all

% Put all chl on one graph
ep_plot(1,[62:66 69],DateStart,DateEnd,[1.0 0.4 0.4;1.0 0.8 0.0; 0.8 1.0 0.0;0.4 1.0 .4;0.0 0.8 1.0;0.4 0.4 0.8])
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['ep1chl' ext]);
close all

ep_plot(2,[62:66 69],DateStart,DateEnd,[1.0 0.4 0.4;1.0 0.8 0.0; 0.8 1.0 0.0;0.4 1.0 .4;0.0 0.8 1.0;0.4 0.4 0.8])
set(gcf,'PaperPosition',PaperPosition);
set(gca,'YLim',YLim);
makegif(['ep2chl' ext]);
close all

%=================================================
% Create a gif. Matlab no longer supports gifs so 
%  1. Create a temporary png image
%  2. Change it to a gif using a java program
%=================================================
function makegif(Filename)
java on;
at_setfontsize
png = ['d:\apps\matlab\work\' Filename '.png'];
gif = ['d:\www\bog\images\' Filename '.gif'];
print('-dpng',png);
java_method('convert', 'PngToGif', png, gif);
delete(png);