!   TITLE:  airfoil.mac
!
mod_assemb_prep
!
sf=5.0              ! airfoil scaling factor
ar=10.0             ! aspect ratio (wingspan/trunk wing chord)
tip2trnk=2.0        ! scaling factor from of wing tip chord to wing trunk
!
csys,0
!
create_cs,0,'aircs',0,sf*(-0.5),0,0,0,0,0
!
init_entity_num
lep=km+1
!
k,,sf*0.00000,sf*(-0.00006)
!
k,,sf*0.00240,sf*0.00917                ! Upper Camber
k,,sf*0.00909,sf*0.01947
k,,sf*0.02004,sf*0.03027
k,,sf*0.03527,sf*0.04120
k,,sf*0.05469,sf*0.05201
!
k,,sf*0.07816,sf*0.06250
k,,sf*0.10546,sf*0.07247
k,,sf*0.13635,sf*0.08175
k,,sf*0.17050,sf*0.09019
k,,sf*0.20758,sf*0.09761
!
k,,sf*0.24720,sf*0.10389
k,,sf*0.28894,sf*0.10887
k,,sf*0.33237,sf*0.11240
k,,sf*0.37702,sf*0.11428
k,,sf*0.42253,sf*0.11427
!
k,,sf*0.46864,sf*0.11219
k,,sf*0.51524,sf*0.10784
k,,sf*0.56247,sf*0.10147
k,,sf*0.61010,sf*0.09373
k,,sf*0.65752,sf*0.08513
!
k,,sf*0.70408,sf*0.07603
k,,sf*0.74914,sf*0.06673
k,,sf*0.79206,sf*0.05746
k,,sf*0.83222,sf*0.04844
k,,sf*0.86902,sf*0.03983
!
k,,sf*0.90193,sf*0.03175
k,,sf*0.93044,sf*0.02428
k,,sf*0.95409,sf*0.01737
k,,sf*0.97285,sf*0.01082
k,,sf*0.98710,sf*0.00507
!
k,,sf*0.99658,sf*0.00126
k,,sf*1.00000,sf*0.00000
!
!
init_entity_num
bsplin,all,,,,,,0,-1
init_entity_num
!
k,,sf*0.99630,sf*0.00086                ! Lower Camber
k,,sf*0.98504,sf*0.00323
k,,sf*0.96588,sf*0.00655
k,,sf*0.93862,sf*0.00990
k,,sf*0.90359,sf*0.01198
!
k,,sf*0.86198,sf*0.01155
k,,sf*0.81565,sf*0.00804
k,,sf*0.76645,sf*0.00167
k,,sf*0.71614,sf*(-0.00705)
k,,sf*0.66662,sf*(-0.01663)
!
k,,sf*0.61825,sf*(-0.02441)
k,,sf*0.56920,sf*(-0.02925)
k,,sf*0.51867,sf*(-0.03225)
k,,sf*0.46731,sf*(-0.03415)
k,,sf*0.41576,sf*(-0.03519)
!
k,,sf*0.36468,sf*(-0.03554)
k,,sf*0.31473,sf*(-0.03531)
k,,sf*0.26654,sf*(-0.03456)
k,,sf*0.22073,sf*(-0.03334)
k,,sf*0.17788,sf*(-0.03166)
!
k,,sf*0.13854,sf*(-0.02954)
k,,sf*0.10324,sf*(-0.02699)
k,,sf*0.07242,sf*(-0.02399)
k,,sf*0.04653,sf*(-0.02052)
k,,sf*0.02592,sf*(-0.01656)
!
k,,sf*0.01099,sf*(-0.01211)
k,,sf*0.00245,sf*(-0.00704)
!
init_entity_num
!
ksel,s,kp,,km_prev,km
ksel,a,kp,,lep
bsplin,all,,,,,,,,,0,1
!
create_area-lines
trnkarea=am
!
csys,0
vext,trnkarea,,,0,0,sf*ar/2,tip2trnk,tip2trnk
init_entity_num
!
lplot