ind = 0
labels = []
for s in stgvals:
CLaAlt = []
AVLCLa = []
for g in gapvals:
biwing.Stagger = s
biwing.Gap = g
# AVL CLalpha
CL = biwing.GetAVLCL(a, 'g' + str(g) + 's' + str(s), RunDir, Execute)
CLalpha = (CL[-1] - CL[0])
AVLCLa.append(CLalpha)
# Altman CLalpha
CLaAlt.append(biwing.CL2(1 * ARCDEG) - biwing.CL2(0 * ARCDEG))
# pyl.plot(gapvals,CLaAlt,clrstr[ind] + '--')
pyl.plot(gapvals, AVLCLa, clrstr[ind])
ind = ind + 1
# labels.append(r'Altman $CL_{\alpha}$ S = '+str(s))
labels.append(r'AVL $CL_{\alpha}$ S = ' + str(s))
#
# Find the BWCF to match the AVL output
#
biwing.BWCFInterp = [0.745, 0.845, 0.88, 0.915] #chord=12.5IN
#biwing.BWCFInterp = [0.69,0.785,0.835,0.87] #chord=15.0IN
biwing.Stagger = 0
gapvals = npy.linspace(0.1, 0.4, 4)
#a = alpha2d / (ARCDEG)
pyl.figure(1)
clrstr = 'bgrycmk'
ind = 0
labels=[]
for s in stgvals:
CL2_CL3 = []
CL2_CL1 = []
for g in gapvals:
BoxWing.Stagger = s
BoxWing.Gap = g
alpha2d=0*ARCDEG
cl1 = BoxWing.CL(alpha2d)
cl2 = BoxWing.CL2(alpha2d)
cl3 = (BoxWing.UpperWing.CL(alpha2d) + BoxWing.LowerWing.CL(alpha2d))/2
CL2_CL1.append(cl2/cl1)
CL2_CL3.append(cl2/cl3)
pyl.subplot(121)
pyl.plot(gapvals,CL2_CL3,clrstr[ind])
pyl.subplot(122)
pyl.plot(gapvals,CL2_CL1,clrstr[ind])
ind = ind + 1
labels.append('Stagger = '+str(s))
pyl.subplot(121)
pyl.xlabel('Gap/Span')