subroutine propeller(altitude, Radius,Jar,CP,CT)


%	*****  Program  NDProp *****
%
%	Adaptation of the program
%	'Propeller 1.2.3' by Barry Young
%
%	***  Variable Usage   ***
%
%	A,Acousvel = speed of sound                                      （音速，計算馬赫數影響時用）
%	ACL = Cl correction constant                                     （升力系數修正因子）
%	Azero = lift curve slope                                         （升力線斜率）
%	AF = Airfoil designator (integer)                                （翼型  ？）
%	afname = airfoil designation                                     （翼型  ？）
%	AFds(i)= array of airfoil names                                  （所選翼型序號）
%	Ads = temporary register                                         （暫存）
%	AIRFOILNAME = name of airfoil used on the prop                   （所選翼型名號）
%	Alpha(i,j) = angle of attack at Xi & Jj                          （進距比Jj時Xi處的攻角--rad）
%	ALPHADEG(I,J) = angle of attack in degrees                       （進距比Jj時Xi處的攻角--degrees）
%	ANS = input integer 'decision' varibles
%	BCD = Cd correction constant                                     （阻力系數修正因子）
%	Beta(i) = blade angle at X(i) (DEG)                              （槳葉角--degrees）
%	C(i) = chord at X(i) (INCHES)                                    （弦長）
%	CF = sine of angle phi                                           （sin（φ））
%	CL = lift coefficient                                            （升力系數）
%	CP(j) = Cp at J(j)                                               （J(j)時功率系數）
%	CT(j) = Ct at J(j)                                               （J(j)時推力系數）
%	COLS = number of columns in a matrix                             （矩陣列數）
%	C0 = constant pi^3/8                                             （常數因子）
%	C1, C2 = constants                                               （常數）
%	DQ(i,j) = d(CQ)/dX at X(i) and J(j)                              （進距比Jj時Xi處的d(CQ)/dX）
%	DR = conversion, degrees to radians                              （度-- 變化-- 弧度）
%	DT(i,j) = d(CT)/dX at X(i) and J(j)                              （進距比Jj時Xi處的d(CT)/dX）
%	D1,D2 = temporary variables                                      （暫時變量）
%	ETA(j) = efficiency at j                                         （進距比J時槳的效率）
%	Fzero = thickness routine flag                                   （厚度……標志）
%	F1 = jmin flag                                                   （jmin標志）
%	F2 = jmax flag                                                   （jmax標志）
%	F3 = general flag                                                  ?
%	F4 = change routine flag                                           ?
% 	F5 = convergence flag                                             （收斂?標志）
%	F6 = airfoil selection input flag                                 （翼型輸入標志）
%	FX = input to subroutine                                          （子程序輸入變量）
%	FY = output from subroutine                                       （子程序輸出變量）
%	F(i,j) = matrix of spline fit coefficients                        （擬合系數）
%	GP(i) = geometric pitch at X(i)                                   （X(i)處的幾何螺距）
%	H = altitude in thousands of feet                                 （飛行高度--thousands of feet）
%	I = counter or dummy variable                                     （記數變量）
%	II= row number in matrix                                          （行號）
%	ikey= logical unit identifier for keyboard                        （鍵盤設備號）      
%	imon= logical unit identifier for monitor                         （顯示器設備號）        
%	INPUTDATAFILE = file name for input data                          （輸入文件名）
%	J, JJ = counter or dummy variable                                 （記數變量）
%	Jar(i) = advance ratio                                            （進距比）
%	Jmax = Advance ratio maximum                                      （最大進距比）
%	Jmin = Advance ratio minimum                                      （最小進距比）
%	Kgold = goldstein constant                                        （goldstein常數）
%	K(i,j) = array of goldstein constants                             （goldstein常數向量）
% 	KK=column number in matrix                                        （矩陣列號）
% 	LDA = indicator of lift/drag adjustments (integer) 
%		=1  no CL/CD corrections
%		=2  mach no. adjustments (glauert)                        （LDA=2馬赫數調整）
%		=3  Re adjustment                                         （LDA=3雷諾數調整）
%		=4  mach and Re adjustment                                （LDA=4 馬赫數，雷諾數調整）
%	MF = Conversion, MPH to FPS                                       （英里/小時--英尺/秒）
%	MPH = Miles per hour                                              （飛行速度--Miles per hour）
%	MT = mach number at tip (corrected for induced velocity)          （槳尖馬赫數，經誘導速度修正后）
%C	MTUC = mach number at tip uncorrected（槳尖馬赫數，經誘導速度修正前）
%	N0= mach/reynolds number correction flag
%	N = number of radial data points
%	NB = number of blades
%	NAIRFOILS = number of airfoils in the AIRFOIL data file
%	NU = Kinematic viscosity（動黏度）
%	PHI = angle of flow incident to blade（氣流角）
%	PI= 3.14159....
%	PROPNAME = string variable identifier for prop（螺旋槳名稱）
%	Pn = 10^n; weird little variable used in rounding numbers
%	R0-R8 = variables used in aerodynamic coefficient corrections
%	Radius = propeller radius (inches)（直徑--inches）
%	RA = indication of analysis refinement (integer)
%		=1  simple blade element theory（計算方法標志，RA=1，用simple blade element theory）
%		=2  include induced velocity（計算方法標志，RA=2，用include induced velocity）
%		=3  include induced velocity and tip losses（計算方法標志，RA=2，用 include induced velocity and tip losses）
%	RE(i,j) = Reynolds number（雷諾數矩陣）
%	REAFD = Reynolds number of the airfoil data（翼型臨界雷諾數）
%	RD = conversion, radians to degrees	（弧度-- 變化-- 度）
%	R(i) = radius inputs NONDIMENSIONAL(i.e.fraction of radius)（無量綱半徑輸入）
%	RHO = atmospheric density（空氣密度）
%	ROWS= number of rows in a matrix
%	RPM = revolutions per minute（轉速--  revolutions per minute）
%	S(i) = solidity at i       ？
%	SF = sine of angle phi（sin(φ)）
%	T(i) = thickness ratio at X(i) NONDIMENSIONAL（相對厚度）
%		(i.e. thickness/local chord)
%	Theta = induced angle of attack（誘導速度修正角）
%	Ufixed = RPM - MPH indicator (integer)
%           =1 MPH fixed
%           =2 RPM fixed
%	X(i) = % radial location
%	Z(ijk) = table of airfoil aerodynamic coefficients（翼型氣動系數表）