1      SUBROUTINE TRSONJ(NZ)
2C
3C*** COMPUTES TRANSONIC H.T. LIFT SLOPE,CLMAX,ALPHA CLMAX,AND CD0
4C***      BODY LIFT AND MOMENT SLOPES,DRAG AT ANGLE OF ATTACK
5C
6      DIMENSION ROUTID(2)
7      DIMENSION Q3243A(3),Q3243B(3),Q3244B(3),Q3244A(3),Q3244C(3),Q34171
8     1  (4),Q34172(4),Q3418A(3),Q34181(3),Q34182(3),Q3419B(3),Q3419C(3),
9     2  Q3419A(3),Q15127(4),Q5129R(3),Q5129L(3),Q23124(3),Q23126(3)
10      REAL K,MACH,MFB,MT(5)
11      REAL MFB0
12      DIMENSION T429L(15),D429L(56),T429R(14),D429R(35),T424(16),
13     1          D424(39)
14      DIMENSION WTYPE(4),C17(6) ,TR(11),DR(11)
15      DIMENSION T426(16),D426(55)
16      DIMENSION X(11),Y(11),C(6)
17      DIMENSION T418A(10),D418A(10), T18B1(3),D18B1(20),C18B1(60),
18     1          T18B2(10), D18B2(10) ,C18B2(100)
19      DIMENSION T419A(22),D419A(77) ,T419B(13),D419B(13),T419C(26),
20     1          D419C(130)
21      DIMENSION T43A(15),D43A(56),T43B(13),D43B(40),T44A(13),D44A(42),
22     1          T44B(13),D44B(42),T44C(6),D44C(6)
23      DIMENSION XMTD(15),CDW1(15),CDW2(15),VAR(4),LGH(4)
24      DIMENSION CLAMT(5),ALPHA(20) ,CDJB(20)
25      DIMENSION X27M(4),X27I(4),C27(6)
26      COMMON /CONSNT/ PI,DEG,UNUSED,RAD
27      COMMON /OVERLY/ NLOG,NMACH,M,NALPHA,IG,NF
28      COMMON /OPTION/ SR,CBARR,RUFF,BLREF
29      COMMON /FLGTCD/ FLC(93)
30      COMMON /BDATA/  BD(762)
31      COMMON /HTDATA/ A(195),B(49)
32      COMMON /HTI/    WINGIN(77)
33      COMMON /SUPBOD/ SBD(227)
34      COMMON /SBETA/  STB(243),TRA(108)
35      COMMON /IBODY/  PB, BODY(400)
36      COMMON /IHT/    PW, WING(380)
37      COMMON /IBH/    PBW, BW(380)
38      COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC,
39     1                HEAD,PRPOWR,JEQPOW,LOASRT,TVTPAN,
40     2                SUPERS,SUBSON,TRANSN,HYPERS
41      LOGICAL         FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC,
42     1                HEAD,PRPOWR,JEQPOW,LOASRT,TVTPAN,
43     2                SUPERS,SUBSON,TRANSN,HYPERS
44      EQUIVALENCE (TOC,WINGIN(16)),(ARSTAR,A(7)),(SWEPE2,A(71))
45     1,(COSLE,A(61)),(TAPR,A(27)),(CBAR,A(16)),(RNFS,A(129)),
46     2(XOC,WINGIN(18)),(SRSTAR,A(3)),(COSC4,A(67)),(SS,BD(93))
47     3,(DB,SBD(6)),(DMAX,SBD(120)),(XLB,SBD(2)),(SB,BD(57))
48     4,(TANC2,A(74))
49      EQUIVALENCE (CLA14,TRA(1)),(CLAMO,WINGIN(69)),(K,TRA(3))
50     1,(MFB0,TRA(5)),(MFB,TRA(6)),(AOC,TRA(7)),(CFBCT,TRA(8)),(BETAFB,
51     2TRA(9)),(CLAFBT,TRA(10)),(MACH,TRA(4)),(CLAFB,TRA(12)),(CLAA,TRA
52     3(13)),(BOC,TRA(14)),(CLAB,TRA(15)),(MT(1),TRA(16)),(CLAMT(1),TRA
53     4(21)),(CLA,WING(101)),(C1,TRA(27)),(ARATIO,TRA(28)),(BU4,TRA(29)),
54     5(CLMAX6,TRA(30)),(ACLBAS,TRA(31)),(DACMA6,TRA(32)),(C3,TRA(33))
55      EQUIVALENCE (CLAB14,SBD(18)),(CMAB14,SBD(110)),(CD014,SBD(124)),
56     1(CMAB6,BODY(121)),(CFB6,BD(92)),(CD0B6,BD(61)),(CDBB6,BD(60))
57      EQUIVALENCE (DALCM,TRA(34)),(DCLMAX,TRA(35)),(ALCLM6,TRA(36)),
58     1(ALCLMT,TRA(37)),(CLMAXT,TRA(38)),(RLCOFF,TRA(39)),(RNN,TRA(40)),
59     2(RL,TRA(41)),(CF,TRA(42)),(CDW2(1),TRA(43)),
60     3(CDW,TRA(67)),(CDF,TRA(68)),(CD0W,WING(1)),(CLAW6,TRA(70))
61     4,(CLABM6,BODY(101)),(CLABD,BODY(101)),(CMAB,BODY(121))
62     5,(CDPB,TRA(80)),(CDBFIG,TRA(81)),(CDBB,TRA(76)),(CDWB,TRA(77))
63     6,(CD0B,TRA(78)),(CDJB(1),BODY(1)),(CD0WB,BW(1))
64     7,(CDFB,TRA(79)),(DCLA14,TRA(82))
65      EQUIVALENCE (TANLE,A(62))
66      LOGICAL LOWAR
67      DATA ROUTID /4HTRSO,4HNI  /
68      DATA Q3243A/4H4.1.,4H3.2-,4H53A /,Q3243B/4H4.1.,4H3.2-,4H53B /,
69     1Q3244B/4H4.1.,4H3.2-,4H54B /,Q3244A/4H4.1.,4H3.2-,4H54A /,Q3244C/
70     24H4.1.,4H3.2-,4H54C /,QCLA/4HCLA /,Q34171/4H4.1.,4H3.4-,4H24B , 4H
71     3C1  /,Q34172/4H4.1.,4H3.4-,4H24B ,4HC2  /,Q3418A/4H4.1.,4H3.4-, 4H
72     425A /,Q34181/4H4.1.,4H3.4-,4H25B1/,Q34182/4H4.1.,4H3.4-,4H25B2/,
73     5Q3419B/4H4.1.,4H3.4-,4H26B /,Q3419C/4H4.1.,4H3.4-,4H26C /,Q3419A
74     6/4H4.1.,4H3.4-,4H26A /,Q15127/4H4.1.,4H5.1-,4H27  ,4HCEPT/,Q5129R
75     7/4H4.1.,4H5.1-,4H29R /,Q5129L/4H4.1.,4H5.1-,4H29L /,QCDW/4HCDW /,
76     8Q23124/4H4.2.,4H3.1-,4H24  /,Q23126/4H4.2.,4H3.1-,4H26  /
77C
78C              FIGURE 4.1.5.1-27
79C
80      DATA X27I  /1.57780,1.67221,1.98509,2.28874/
81      DATA X27M  /0.0,1.0,2.0,3.0/
82      DATA I27   /0/
83      DATA WTYPE /4HSTRA,4HDOUB,4HCRAN,4HCURV/
84      DATA I17,IN/0,0/
85      DATA NP/11/
86C
87C              FIGURE 4.1.3.4-24B
88C
89      DATA TR /0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1./
90      DATA DR /0.,.21,.5,.9,1.08,1.05,1.,.94,.9,.86,.85/
91      DATA Y  /0.,.225,.47,.496,.43,.32,.21,.125,.075,.0475,.0/
92      DATA X  /0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1./
93C
94C     ----4.1.3.4-25(A)
95C
96      DATA T418A /0.,.4,.8,1.2,1.6,2.0,2.4,2.8,3.2,3.6/
97      DATA D418A /3*35.,32.,28.,25.,23.2,22.,21.5,21./
98C
99C     ----4.1.3.4-25(B)
100C
101      DATA T18B1/.2,.4,.6/
102      DATA D18B1/4.5,5.,5.5,6.,6.5,7.,7.5,8.,8.5,9.,9.5,10.,10.5,11.,
103     111.5,12.,12.5,13.,13.5,14./
104      DATA C18B1/0.,.5,.9,1.4,1.9,2.5,3.3,4.,4.6,5.6,6.4,7.3,8.2,9.2,10.
105     1,11.,12.,13.,14.,15. , 0.,.2,.4,.7,1.2,1.7,2.4,3.,3.7,4.6,5.2,6.,
106     26.9,7.8,8.6,9.5,10.4,11.4,12.3,13.5 , 0.,0.,.1,.2,.5,.7,1.,1.3,1.6
107     3,2.,2.5,3.,3.6,4.3,4.9,5.5,6.2,7.,7.6,8.5/
108      DATA T18B2/0.,2.,3.,4.,5.,6.,7.,8.,9.,30./
109      DATA C18B2 /                                    10.,8.5,6.9,5.5,4.
110     1,2.6,1.5,.7,.1,0. , 8.7,7.3,5.3,4.2,2.6,1.4,0.5,-.2,-.5,0. , 7.5,
111     25.9,4.2,2.5,1.2,0.,-.7,-1.1,-.8,0. , 5.5,3.4,1.6,0.,-1.3,-2.1,-2.5
112     3,-2.,-.8,0. , 3.,.7,-1.4,-3.3,-4.3,-4.3,-3.1,-2.,-.8,0. , .3,-2.5,
113     4-4.7,-5.8,-5.3,-4.3,-3.1,-2.,-.8,0. ,-2.2,-5.,-6.7,-6.3,-5.3,-4.3,
114     5-3.1,-2.,-0.8,0. , -3.3,-6.6,-7.2,-6.3,-5.3,-4.3,-3.1,-2.,-.8,0.,
115     6  -4.2,-7.,-7.2,-6.3,-5.3,-4.3,-3.1,-2.,-.8,0. , -8.5,-7.9,-7.2,
116     7 -6.3,-5.3,-4.3,-3.1,-2.,-.8,0. /
117      DATA D18B2  /.0,.5,1.,1.5,2.,2.5,3.,3.5,4.,4.5/
118C
119C              FIGURE 4.1.3.4-26B
120C
121      DATA T419B
122     1 / 0., 4., 5., 6., 7., 8., 9., 10., 11., 12., 13., 14., 16. /
123      DATA D419B
124     1 / 1., 1., .93, .8, .55, .4, .3, .24, .2, .16, .14, .12, .12 /
125C
126C              FIGURE 4.1.3.2-53A
127C
128      DATA T43A/
129     1 0.,     2.,     4.,     6.,     8.,     10.,    12.,    14. ,
130     2 0., 1., 2., 3., 4., 6., 8./
131      DATA D43A
132     1/8*1.0,
133     2 6*1.0,.98,.9,
134     3 5*1.0,.983,.89,.855,
135     4 4*1.0,.939,.888,.857,.835,
136     5 2*1.0,.958,.904,.869,.842,.82,.802,
137     6 1.0,    .952,   .908,   .872,   .842,   .816,   .80,    .787,
138     7 1.0,    .952,   .908,   .872,   .842,   .816,   .80,    .787 /
139C
140C              FIGURE 4.1.3.2-53B
141C
142      DATA T43B
143     1  / 0.0, 0.17453, 0.34907 , 0.52360, 0.69813,
144     2         1.04720, 1.22173 , 1.57080    ,
145     3 .799,.85,.9,.95,1.000/
146      DATA D43B/
147     1 .8,     .807,   .82,    .841,   .867,   .918,   .945,   1.0,
148     2 .85,    .855,   .865,   .88,    .9,     .938,   .958,   1.0,
149     3 .90,    .903,   .91,    .92,    .932,   .958,   .974,   1.0,
150     4 .95,    .952,   .955,   .96,    .967,   .98,    .987,   1.0,
151     5    1.0,  1.0,  1.0,  1.0,  1.0,  1.0,1.0,  1.0  /
152C
153C              FIGURE 4.1.3.2-54A
154C
155      DATA T44A/
156     1 0.,     2.,     4.,     6.,     8.,     12.,    14. ,
157     2 1.,     2.,     3.,     4.,     6.,     8. /
158      DATA D44A/
159     1 1.09,   1.088,  1.056,  1.008,  .962,   .872,   .828,
160     2 .98,    1.088,  1.1,    1.05,   .992,   .882,   .83,
161     3 .85,    1.,     1.075,  1.06,   1.005,  .89,    .835,
162     4 .74,    .9,     1.015,  1.07,   1.025,  .895,   .837,
163     54*1.15,1.07,.915,.915,
164     65*1.12,.94,.94/
165C
166C              FIGURE 4.1.3.2-54B
167C
168      DATA T44B/
169     1 0.,     4.,     6.,     8.,     10.,    12.,    14. ,
170     2 1.,     2.,     3.,     4.,     6.,     8. /
171      DATA D44B/
172     1 -.14,   -.14,   -.14,   -.14,   -.14,   0.,     .44,
173     2 -.07,   -.07,   -.07,   -.07,   0.,     .31,    .64,
174     3 .05,    .05,    .05,    .06,    .23,    .49,    .75,
175     4 .06,    .06,    .06,    .15,    .35,    .58,    .79,
176     5 .08,    .08,    .14,    .29,    .48,    .66,    .84,
177     6 .09,    .15,    .24,    .39,    .555,   .72,    .89 /
178C
179C              FIGURE 4.1.3.2-54C
180C
181      DATA T44C
182     1/0.,     2.,     4.,     6.,     8.,     16. /
183      DATA D44C
184     1  / -.04, .01, .075, .13, .15, .15 /
185C
186C              FIGURE 4.1.3.4-26A
187C
188      DATA T419A
189     1 /          0.8, 1.0, 1.2, 1.6, 2.0, 2.4, 3.0, 4*0.,
190     2   .6, .7, .8, .9, .95, 1.0, 1.05, 1.1, 1.15, 1.2, 1.4 /
191      DATA D419A
192     1 /              0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 ,
193     2               -.050,-.050,-.050,-.030,-.015,-.010, 0.0 ,
194     3               -.115,-.100,-.083,-.047,-.018, 0.0 , .025,
195     4               -.145,-.110,-.080,-.020, .022, .050, .097,
196     5               -.120,-.080,-.030, .050, .100, .140, .204,
197     6               -.060, 0.0 , .060, .150, .210, .250, .325,
198     7                .090, .130, .180, .257, .310, .330, .360,
199     8                .113, .165, .210, .280, .320, .345, .360,
200     9                .108, .160, .205, .284, .320, .344, .360,
201     A                .096, .145, .190, .273, .320, .343, .355,
202     B                .020, .050, .100, .200, .300, .330, .340/
203C
204C              FIGURE 4.1.3.4-26C
205C
206      DATA T419C
207     1 / 0.0, 0.4, 0.8, 1.0, 1.2, 1.6, 2.0, 2.4, 3.0, 5.0  ,3*0.,
208     2   .6,.7,.75,.8,.85,.9,.95,1.0,1.05,1.1,1.15,1.2,1.4  /
209      DATA D419C
210     1 /   0.0, 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0 , 0.0,  0.0 , 0.0 ,
211     2   -1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,-1.00,
212     3   -1.25,-1.22,-1.20,-1.20,-1.20,-1.10,-1.10,-1.00,-1.00,-0.85,
213     4   -2.05,-1.75,-1.50,-1.40,-1.30,-1.20,-1.10,-1.00,-1.00,-0.95,
214     5   -2.20,-1.90,-1.60,-1.50,-1.40,-1.20,-1.10,-1.00,-0.90,-0.63,
215     6   -1.85,-1.60,-1.40,-1.30,-1.20,-1.00,-0.80,-0.70,-0.60,-0.40,
216     7   -1.80,-1.38,-1.00,-0.90,-0.70,-0.50,-0.30,-0.10, 0.0 , 0.10,
217     8   -0.60,-0.20, 0.20, 0.30, 0.50, 0.80, 1.00, 1.20, 1.40, 2.00,
218     9    1.90, 2.40, 3.00, 3.20, 3.50, 3.90, 4.30, 4.60, 5.00, 6.20,
219     A    2.40, 4.05, 5.50, 6.20, 7.00, 8.00, 8.90, 10.0, 10.9, 13.9,
220     B    1.60, 4.04, 6.90, 8.00, 9.10, 10.6, 11.9, 13.0, 14.3, 18.2,
221     C   -1.20, 3.90, 7.60, 8.90, 10.2, 12.1, 13.5, 14.7, 15.9, 18.8,
222     D   -11.0,-0.60, 8.10, 10.5, 12.0, 14.4, 16.2, 17.5, 18.8, 21.3 /
223C
224C              FIGURE 4.1.5.1-29 LEFT SIDE (SUBSONIC)
225C
226      DATA T429L/
227     10.0,0.4,0.6,0.8,1.0,1.2,1.4,1.475,
228     20.,.5,1.0,1.5,2.0,3.0,4.0/
229      DATA D429L/
230     A8*0.0,
231     10.90,0.71,0.55,0.39,0.21,0.10,0.00,0.0,
232     22.34,2.03,1.70,1.05,0.51,0.21,0.01,0.0,
233     33.01,2.92,2.68,1.99,1.00,0.32,0.01,0.0,
234     43.30,3.27,3.12,2.50,1.26,0.42,0.03,0.0,
235     53.49,3.43,3.38,2.80,1.52,0.51,0.06,0.0,
236     63.61,3.58,3.52,3.14,1.83,0.52,0.09,0.0/
237C
238C              FIGURE 4.1.5.1-29 RIGHT SIDE (SUPERSONIC)
239C
240      DATA T429R/
241     10.0,.4,.8,1.2,1.4,0.0,0.0,
242     20.,.5,1.0,1.5,2.0,3.0,4.0/
243      DATA D429R/
244     A5*0.0,
245     10.90,1.08,1.18,1.28,1.31,
246     22.34,2.53,2.68,2.75,2.78,
247     33.01,3.08,3.10,3.12,3.12,
248     43.30,3.32,3.32,3.30,3.28,
249     53.49,3.49,3.48,3.42,3.39,
250     63.61,3.61,3.60,3.56,3.51/
251C
252C              FIGURE 4.2.3.1-24
253C
254      DATA T424/
255     1.75,.85,.90,.95,1.0,1.015,1.030,1.050,1.075,1.125,1.15,1.2,1.3,
256     2.011,.067,.122/
257      DATA D424/
258     1.011,.014,.018,.031,.069,.040,.020,6*0.0,
259     2.067,.071,.078,.092,.140,.150,.120,.086, .067, .050,.047,.045,.05,
260     3.122,.130,.139,.155,.204,.220,.236,.210, .168, .110, .101,.100,.1/
261      DATA T426/
262     16.0,7.0,8.0,10.,12.,14.,16.,18.,20.,24.,26.,
263     21.0,1.025,1.05,1.1,1.2/
264      DATA D426/
265     1.118,.084,.065,.041,.029,.020,.016,.012,.010,.007,.006,
266     2.147,.110,.088,.061,.043,.033,.026,.022,.019,.015,.013,
267     3.166,.120,.097,.068,.050,.039,.031,.026,.022,.015,.013,
268     4.186,.138,.110,.076,.057,.044,.034,.028,.022,.015,.013,
269     5.200,.152,.120,.080,.059,.046,.036,.028,.022,.015,.013/
270C
271C  *** WING ALONE ***
272C *** WING LIFT CURVE SLOPE ***
273C
274      MACH=FLC(NZ+2)
275      IF(WINGIN(15).NE.WTYPE(1))GO TO 1100
276      CALL TRNHT(CLA14, DCLA14)
277      K=CLAMO*RAD/(2.0*PI)
278      BETA6=0.80
279C
280C              FIGURE 4.1.3.2-53A FORCE BREAK MACH NUMBER=MFB ( 0 SWEEP)
281C
282      VAR(1)=TOC*100.
283      VAR(2)=ARSTAR
284      LGH(1)=8
285      LGH(2)=7
286      CALL INTERX(2,T43A,VAR,LGH,D43A,MFB0,8,56,
287     1            2,2,0,0,2,-1,0,0,Q3243A,3,ROUTID)
288C
289C              FIGURE 4.1.3.2-53B MFB0 CORRECTED FOR SWEEP
290C
291      VAR(1)=SWEPE2
292      VAR(2)=MFB0
293      LGH(1)=8
294      LGH(2)=5
295      CALL INTERX(2,T43B,VAR,LGH,D43B,MFB,8,40,
296     1            0,2,0,0,2,0,0,0,Q3243B,3,ROUTID)
297      IF(NF .LT. 0) RETURN
298C
299C              FIGURE 4.1.3.2-54B (A/C)
300C
301      VAR(1)=TOC*100.
302      VAR(2)=ARSTAR
303      LGH(1)=7
304      LGH(2)=6
305      CALL INTERX(2,T44B,VAR,LGH,D44B,AOC,7,42,
306     1            0,2,0,0,1,2,0,0,Q3244B,3,ROUTID)
307C
308C              FIGURE 4.1.3.2-54A CORRECTION TO LIFT CURVE SLOPE
309C
310      CALL INTERX(2,T44A,VAR,LGH,D44A,CFBCT,7,42,
311     1            2,1,0,0,2,2,0,0,Q3244A,3,ROUTID)
312C
313C
314C
315      BETAFB = SQRT(1.0-MFB**2)
316      IF(MFB .GT. 0.98) BETAFB = 0.0
317      BB   = BETAFB
318      ARG1 = 2.0*PI*ARSTAR/RAD
319      ARG2=(ARSTAR/K)**2
320      ARG3=BB**2+TANC2**2
321      ARG4 = 2.0+SQRT(ARG2*ARG3+4.0)
322      CLAFBT = ARG1/ARG4
323C
324      BB   = 0.8
325      ARG2 = (BB*ARSTAR/K)**2
326      ARG3 = 0.0
327      IF(ARG2 .GT. 0.0) ARG3 = 1.0+(TANC2/BB)**2
328      ARG4 = 2.0+SQRT(ARG2*ARG3+4.0)
329      CLAW6 = ARG1*SRSTAR/(ARG4*SR)
330C
331      XM = 0.75
332      IF(XM .GT. (MFB-0.1)) XM = MFB-0.1
333      IF(XM .GT. MACH) XM=MACH
334      BB   = SQRT(1.0-XM**2)
335      ARG2=(ARSTAR/K)**2
336      ARG3=BB**2+TANC2**2
337      ARG4 = 2.0+SQRT(ARG2*ARG3+4.0)
338      CLAW7 = ARG1*SRSTAR/(ARG4*SR)
339      DCLA7 = XM*ARG1*ARG2*SRSTAR/SR
340      DCLA7 = DCLA7/(ARG4**2*SQRT(ARG2*ARG3+4.0))
341      CLAFB=CLAFBT*CFBCT *SRSTAR/SR
342      CLAA=(1.-AOC)*CLAFB
343C
344C              FIGURE 4.1.3.2-54C (B/C)
345C
346      VAR(1)=TOC*100.
347      LGH(1)=6
348      CALL INTERX(1,T44C,VAR,LGH,D44C,BOC,6,6,
349     1            2,0,0,0,1,0,0,0,Q3244C,3,ROUTID)
350      CLAB=(1.-BOC)*CLAFB
351      MT(1)=XM
352      MT(2)=MFB
353      MT(3)=MFB+.07
354      MT(4)=MFB+.14
355      MT(5)=1.40
356      CLAMT(1)=CLAW7
357      CLAMT(2)=CLAFB
358      CLAMT(3)=CLAA
359      CLAMT(4)=CLAB
360      CLAMT(5)=CLA14
361      CALL TRANF(5, MT, CLAMT, DCLA7, DCLA14, MACH, CLA)
362C
363C  *** CLASSIFICATION OF ASPECT RATIO ***
364C
365      CALL TBFUNX(A(27), C1   ,DYDX,NP,X,Y,C,IN,MI,NG,
366     1            0,0,Q34171,4,ROUTID)
367      ARG1=(C1+1.)*COSLE
368      ARATIO=A(128)/ARG1
369C
370C              FIGURE 4.1.3.4-24B (C1) AND (C2)
371C
372      VAR(1)=TAPR
373      LGH(1)=11
374      CALL INTERX(1,TR,VAR,LGH,DR,C2,11,11,
375     1           0,0,0,0,0,0,0,0,Q34172,4,ROUTID)
376      LOWAR=.FALSE.
377      IF(ARSTAR.LT.ARATIO)LOWAR=.TRUE.
378C
379C  *** IF LOW ASPECT RATIO CALCULATE CLMAX AND (ALPHA)CLMAX ***
380C
381      IF(.NOT.LOWAR)GO TO 1070
382      BU4=(C1+1.)*ARSTAR*COSLE/BETA6
383      ARG=(C2+1.)*ARSTAR*TANLE
384      A(160)=ARG
385      CALL CLMXB1(BU4,CLMAX6,A,WINGIN)
386C
387C  *** CALCULATE (ALPHA)CLMAX AT MACH=0.60 ***
388C              FIGURE 4.1.3.4-25A
389C
390      CALL TBFUNX(BU4,ACLBAS,DYDX,10,T418A,D418A,C17,I17,MI,NG,
391     1            0,0,Q3418A,3,ROUTID)
392      IF(ARG.LE.4.5)GO TO 1050
393C
394C              FIGURE 4.1.3.4-25B
395C
396      CALL TLINEX(T18B1,D18B1,C18B1,3,20,MT(1),ARG,DACMA6,
397     1            0,0,2,1,Q34181,3,ROUTID)
398      GO TO 1060
399 1050 ARG1=ARSTAR*COSLE*(1.+4.*TAPR**2)
400      CALL TLINEX(T18B2,D18B2,C18B2,10,10,ARG1,ARG,DACMA6,
401     1            0,2,2,0,Q34182,3,ROUTID)
402 1060 ALCLM6=ACLBAS+DACMA6
403C
404C  *** ADJUST PARAMETERS TO TRANSONIC MACH NUMBER ***
405C              FIGURE 4.1.3.4-26B
406C
407      VAR(1)=ARG
408      LGH(1)=13
409      CALL INTERX(1,T419B,VAR,LGH,D419B,C3,13,13,
410     1            0,0,0,0,0,0,0,0,Q3419B,3,ROUTID)
411C
412C              FIGURE 4.1.3.4-26C  DELTA (ALPHA)CLMAX
413C
414      VAR(1)=BU4*BETA6
415      VAR(2)=MACH
416      LGH(1)=10
417      LGH(2)=13
418      CALL INTERX(2,T419C,VAR,LGH,D419C,DALCM ,13,130,
419     1            2,0,0,0,2,2,0,0,Q3419C,3,ROUTID)
420C
421C              FIGURE 4.1.3.4-26A DELTA CLMAX
422C
423      VAR(1)=C3*VAR(1)
424      LGH(1)=7
425      LGH(2)=11
426      CALL INTERX(2,T419A,VAR,LGH,D419A,DCLMAX,11,77,
427     1            2,0,0,0,2,2,0,0,Q3419A,3,ROUTID)
428      ALCLMT=ALCLM6+DALCM
429      CLMAXT=CLMAX6+DCLMAX*(SRSTAR/SR)
430 1070 CONTINUE
431C
432C  *** ZERO ANGLE OF ATTACK DRAG,CD0 ***
433C              FIGURE 4.1.5.1-27 (EQUATION FOR RLCOFF)
434C
435      ARG=12.*CBAR/RUFF
436      CALL TBFUNX(0.6,CEPT,DYDX,4,X27M,X27I,C27,I27,MI,NG,
437     1            0,0,Q15127,4,ROUTID)
438      RLCOFF=ARG**1.0482*10.0**CEPT
439      RNN=CBAR*RNFS
440      IF(RLCOFF.LT.RNN)RNN=RLCOFF
441      CALL FIG26(RNN,0.60,CF)
442      RL=1.2
443      IF(XOC.LT.0.30)RL=2.0
444      CDF=CF*(1.+RL*TOC)*2.*SRSTAR/SR
445      ARG=SQRT(COSC4)
446      XMTD(1)=0.60/ARG
447      XMTD(2)=0.65/ARG
448      XMTD(3)=0.70/ARG
449      XMTD(4)=0.75/ARG
450      XMTD(5)=0.80/ARG
451      XMTD(6)=0.85/ARG
452      XMTD(7)=0.90/ARG
453      XMTD(8)=0.925/ARG
454      XMTD(9)=0.950/ARG
455      XMTD(10)=0.975/ARG
456      XMTD(11)=1.000/ARG
457      XMTD(12)=1.050/ARG
458      XMTD(13)=1.100/ARG
459      XMTD(14)=1.200/ARG
460      XMTD(15)=1.400/ARG
461      ARG2=TOC**.3333
462      VAR(2)=ARSTAR*ARG2
463      ARG3=COSC4**2.5
464      DO 1090 I=1,15
465         ARG4=ABS((XMTD(I)*ARG)**2-1.)
466         VAR(1)=SQRT(ARG4)/ARG2
467C
468C              FIGURE 4.1.5.1-29
469C
470         IF(XMTD(I)*ARG.LE.1.)GO TO 1080
471         LGH(1)=5
472         LGH(2)=7
473         CALL INTERX(2,T429R,VAR,LGH,D429R,CDW1(I),7,35,
474     1               0,0,0,0,0,2,0,0,Q5129R,3,ROUTID)
475         GO TO 1090
476 1080    LGH(1)=8
477         LGH(2)=7
478         CALL INTERX(2,T429L,VAR,LGH,D429L,CDW1(I),8,56,
479     1               0,0,0,0,-1,2,0,0,Q5129L,3,ROUTID)
480 1090 CDW2(I)=CDW1(I)*TOC**1.666*ARG3*SRSTAR/SR
481      CALL TRANF(15, XMTD, CDW2, 0.0, 0.0, MACH, CDW)
482      CD0W=CDW+CDF
483C
484C  *** BODY ALONE ***
485C
486 1100 IF(.NOT.BO)RETURN
487C
488C  ***BODY LIFT AND MOMENT SLOPES ***
489C
490      ARG1=(MACH-.60)/.80
491      CLABD=CLABM6+(CLAB14-CLABM6)*ARG1
492      CMAB=CMAB6+(CMAB14-CMAB6)*ARG1
493C
494C  *** BODY DRAG ***
495C
496      CDFB=CFB6*SS/SR
497      CDPB=CD0B6-CDBB6-CDFB
498      IF(MACH.GE.1.0.AND.MACH.LE.1.2)CDPB=CDPB*(1.-(MACH-1.0)/.2)
499      IF(MACH.GT.1.2)CDPB=0.0
500      IF(DB.LT.0.3*DMAX)DB=0.3*DMAX
501      ARG=(DB/DMAX)**2
502C
503C              FIGURE 4.2.3.1-24 (CDBT)
504C
505      VAR(1)=MACH
506      VAR(2)=CDBB6/ARG*4.0*SR/(PI*DMAX**2)
507      LGH(1)=13
508      LGH(2)=3
509      CALL INTERX(2,T424,VAR,LGH,D424,CDBFIG,13,39,
510     1            0,2,0,0,0,2,0,0,Q23124,3,ROUTID)
511      CDBB=CDBFIG*ARG*PI*DMAX**2/(4.0*SR)
512      IF(MACH.LT.1.0)GO TO 1110
513      VAR(1)=XLB/DMAX
514      VAR(2)=MACH
515      LGH(1)= 11
516      LGH(2)=5
517      CALL INTERX(2,T426,VAR,LGH,D426,CDWB,11,55,
518     1            0,0,0,0,2,2,0,0,Q23126,3,ROUTID)
519      CDWB=CDWB*PI*DMAX**2/(4.*SR)
520      GO TO 1120
521 1110 CDWB=0.0
522 1120 CD0B=CDFB+CDPB+CDBB+CDWB
523      IF(MACH.GT.1.2)GO TO 1130
524      GO TO 1140
525 1130 CD0B=CD0B+(CD014-CD0B)*(MACH-1.2)/0.2
526C
527C  *** BODY DRAG AT ANGLE OF ATTACK ***
528C
529 1140 DO 1150 J=1,NALPHA
530 1150 CDJB(J)=CD0B+(FLC(J+22)/RAD)**2*SB/SR
531      CD0WB=CD0B+CD0W
532      TRA(73)=CD0WB
533      RETURN
534      END
535