function varargout = AFname(N)
% function [x xu xl] = AFname(N)
%
% Put comments related to this airfoil here.  This will be viewed when typing
% 'help <AFname>' at the MATLAB prompt.
%
% Replace the comments in the definitions of xu and xl with the airfoil
% coordinates of the UPPER and LOWER surfaces (from leading to trailing edges)
% and save the file as <AFname>.m.  The SOUSSA MATLAB geometry preprocessors 
% then recognize this airfoil as <AFname>.  Also, change the function name 
% above and make the appropriate adjustments to these help comments.  


TEfix = 0;
% Set TEfix to 1 if the trailing edge of the airfoil geometry is blunt (thick).  
% This will sharpen the edge (by altering the end and neighbouring points) to
% allow for the modelling of the wake generation from a sharp trailing edge
% point using the Kutta condition.  

xu = flipdim([
   % Put your x and y coordinates of the UPPER surface of the airfoil here
   % starting from the leading edge (0.0 0.0) and ending at the trailing edge
   % (1.0 0.0).  Then remove these comments and leave no blank lines.  e.g.,
   %
   %  0.000000  0.000000
   %  0.003990  0.012500
   %  0.006370  0.015280
   %  0.011200  0.019800
   %  0.023480  0.027920
   %  0.048290  0.039600
   %  0.073230  0.048470
   %  0.098230  0.055690
   %      
   %    etc.      etc.
   %      
   %  0.900300  0.013680
   %  0.950140  0.006160
   %  1.000000  0.000000
],1);

xl = [
   % Put your x and y coordinates of the LOWER surface of the airfoil here
   % starting from the leading edge (0.0 0.0) and ending at the trailing edge
   % (1.0 0.0).  Then remove these comments and leave no blank lines.
];

x = [xu; xl];

if TEfix 
   xu(1,2) = 0; 
   xl(end-1,2) = xl(end-1,2)+(xl(end-1,1)-xl(end-2,1))/(xl(end,1)-xl(end-2,1))*abs(xl(end,2)); 
   xl(end,2) = 0; 
end;

[xu, xl] = AFinterp(xu, xl, N);
x = [xu ; xl(2:end,:)];

varargout = {x xu xl};