
function Plot_acs(filenames)

awl = 10;
cwl = 30;
hf1 = figure
hf2 = figure
hf3 = figure
for kf = 1:length(filenames)
    
    if ~exist(filenames{kf})
        disp([filenames{kf} ' NOT FOUND! SKIPPING.'])
        continue
    end
    
    load(filenames{kf})
    
    figure(hf1)
    subplot(211), hold on
    plot(acs_a.time+datenum(0,0,0,0,0,30), acs_a.med(:,awl), 'bo-', 'LineWidth',3.0)
    plot([(acs_a.time+datenum(0,0,0,0,0,30)) acs_a.time+datenum(0,0,0,0,0,30)]', ...
        [acs_a.med(:,awl)+acs_a.std(:,awl) acs_a.med(:,awl)-acs_a.std(:,awl)]', ...
        'b.-','LineWidth',2.0)
    
    subplot(212), hold on
    plot(acs_c.time+datenum(0,0,0,0,0,30), acs_c.med(:,cwl), 'ro-', 'LineWidth',3.0)
    plot([(acs_c.time+datenum(0,0,0,0,0,30)) acs_c.time+datenum(0,0,0,0,0,30)]', ...
        [acs_c.med(:,cwl)+acs_c.std(:,cwl) acs_c.med(:,cwl)-acs_c.std(:,cwl)]', ...
        'r.-','LineWidth',2.0)
    
    
    figure(hf2)
    subplot(211), hold on
    plot(acs_a.time+datenum(0,0,0,0,0,30), acs_a.std(:,awl), 'b-', 'LineWidth',3.0)
    subplot(212), hold on
    plot(acs_c.time+datenum(0,0,0,0,0,30), acs_c.std(:,cwl), 'r-', 'LineWidth',3.0)
    
    figure(hf3)
    subplot(211), hold on
    plot(acs_a.time+datenum(0,0,0,0,0,30), acs_a.ns, 'k-', 'LineWidth',1.5)
    plot(acs_a.time+datenum(0,0,0,0,0,30), acs_a.nsclean, 'b--', 'LineWidth',1.5)
    subplot(212), hold on
    plot(acs_a.time+datenum(0,0,0,0,0,30), acs_c.ns, 'k-', 'LineWidth',1.5)
    plot(acs_a.time+datenum(0,0,0,0,0,30), acs_c.nsclean, 'r--', 'LineWidth',1.5)
    
end

figure(hf1)
subplot(211)
    grid on
    datetick('keeplimits')
    ylabel(['Measured Absorption(' num2str(acs_a.wavel(awl)) 'nm) [1/m]'])
    title([filenames{1} ' to ' filenames{end}])
subplot(212)
    grid on
    datetick('keeplimits')
    ylabel(['Measured Attenuation(' num2str(acs_c.wavel(cwl)) 'nm) [1/m]'])
    title([filenames{1} ' to ' filenames{end}])

    
figure(hf2)
subplot(211)
    grid on
    datetick('keeplimits')
    ylabel(['StdDev Absorption(' num2str(acs_a.wavel(awl)) 'nm) [1/m]'])
    title([filenames{1} ' to ' filenames{end}])
subplot(212)
    grid on
    datetick('keeplimits')
    ylabel(['StdDev Attenuation(' num2str(acs_c.wavel(cwl)) 'nm) [1/m]'])
    title([filenames{1} ' to ' filenames{end}])
figure(hf3)
subplot(211)
    grid on
    datetick('keeplimits')
    ylabel('Nspectra Absorption')
    title([filenames{1} ' to ' filenames{end}])
subplot(212)
    grid on
    datetick('keeplimits')
    ylabel('Nspectra Attenuation')
    title([filenames{1} ' to ' filenames{end}])
