1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
|
C
C----- CIRCLE.INC include file for circle-plane operations
C
C n
C ICX number of circle-plane points for complex mapping ( 2 + 1 )
C IMX number of complex mapping coefficients Cn
PARAMETER (ICX=257)
PARAMETER (IMX=(ICX-1)/4)
C
COMPLEX ZCOLDW, DZTE, CHORDZ, ZLEOLD, ZC, ZC_CN, PIQ, CN, EIW
C
COMMON/CPI01/ NC,MC,MCT
C
COMMON/CPR01/ PI,AGTE,AG0,QIM0,QIMOLD,
& DWC,WC(ICX),SC(ICX),
& SCOLD(ICX),XCOLD(ICX),YCOLD(ICX)
C
COMMON/CPC01/ DZTE, CHORDZ, ZLEOLD, ZCOLDW(ICX),
& ZC(ICX), ZC_CN(ICX,IMX/4),
& PIQ(ICX), CN(0:IMX), EIW(ICX,0:IMX)
C
C NC number of circle plane points, must be 2**n + 1
C MC number of Fourier harmonics of P(w) + iQ(w)
C MCT number of Fourier harmonics for which dZC/dCN are calculated
C
C PI 3.1415926
C AGTE trailing edge angle/pi
C AG0 angle of airfoil surface at first point
C QIM0 Q(w) offset = Q(0)
C QIMOLD Q(w) offset for old airfoil
C DWC increment of circle-plane coordinate w, DWC = 2 pi/(NC-1)
C WC(.) circle plane coordinate w for Fourier operations
C SC(.) normalized arc length array s(w)
C SCOLD(.) normalized arc length s(w) of old airfoil
C XCOLD(.) x coordinate x(w) of old airfoil
C YCOLD(.) y coordinate y(w) of old airfoil
C
C DZTE trailing edge gap specified in the complex plane
C CHORDZ airfoil chord specified in the complex plane
C ZLEOLD leading edge of old airfoil
C ZCOLDW(.) d(x+iy)/dw of old airfoil
C ZC(.) complex airfoil coordinates derived from P(w) + iQ(w)
C ZC_CN(..) sensitivities dZC/dCN for driving geometry constraints
C PIQ(.) complex harmonic function P(w) + iQ(w)
C CN(.) Fourier coefficients of P(w) + iQ(w)
C EIW(..) complex number exp(inw) array on the unit circle
|