aboutsummaryrefslogtreecommitdiff
path: root/src/CIRCLE.INC
diff options
context:
space:
mode:
Diffstat (limited to 'src/CIRCLE.INC')
-rw-r--r--src/CIRCLE.INC48
1 files changed, 48 insertions, 0 deletions
diff --git a/src/CIRCLE.INC b/src/CIRCLE.INC
new file mode 100644
index 0000000..4a540a8
--- /dev/null
+++ b/src/CIRCLE.INC
@@ -0,0 +1,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