clear
clc
format short g
file_names = dir('300k');

num_foil = length(file_names);
clcdmax = ones(num_foil,3)*NaN;
disp('Reading Airfoils at Re = 300k')
prct0 = 0;
for k=3:num_foil
    prct = round(k/num_foil*10)*10;
    if prct ~= prct0
        prct0 = prct;
        disp(sprintf('%g%% Completed - %g of %g remaining',prct,num_foil-k,num_foil))
    end
    fname = file_names(k).name;
    foil300(k).filename = fname;
    [foil300(k).aoa,foil300(k).c_l,foil300(k).c_d,foil300(k).c_m] = parse_polar(['300k\',fname]);
    [c_lmax cl_loc] = max(foil300(k).c_l);
    foil300(k).clmax = foil300(k).c_l(cl_loc);
    foil300(k).cdmax = foil300(k).c_d(cl_loc);
    foil300(k).clcdmax = foil300(k).c_l(cl_loc)/foil300(k).c_d(cl_loc);
    if isempty(foil300(k).clcdmax)==0
        clcdmax300(k,:) = [foil300(k).clcdmax foil300(k).c_l(cl_loc) foil300(k).c_d(cl_loc)];
    end
    clear fname alfa c_l c_d c_m
end

clear file_names
file_names = dir('450k');

num_foil = length(file_names);
clcdmax = ones(num_foil,3)*NaN;
disp('Reading Airfoils at Re = 450k')
prct0 = 0;
for k=3:num_foil
    prct = round(k/num_foil*10)*10;
    if prct ~= prct0
        prct0 = prct;
        disp(sprintf('%g%% Completed - %g of %g remaining',prct,num_foil-k,num_foil))
    end
    fname = file_names(k).name;
    foil450(k).filename = fname;
    [foil450(k).aoa,foil450(k).c_l,foil450(k).c_d,foil450(k).c_m] = parse_polar(['450k\',fname]);
    [c_lmax cl_loc] = max(foil450(k).c_l);
    foil450(k).clmax = foil450(k).c_l(cl_loc);
    foil450(k).cdmax = foil450(k).c_d(cl_loc);
    foil450(k).clcdmax = foil450(k).c_l(cl_loc)/foil450(k).c_d(cl_loc);
    if isempty(foil450(k).clcdmax)==0
        clcdmax450(k,:) = [foil450(k).clcdmax foil450(k).c_l(cl_loc) foil450(k).c_d(cl_loc)];
    end
    clear fname alfa c_l c_d c_m
end

topfoils