1      SUBROUTINE SUBLAT(STB,WING,WINGIN,BW,A,B,AHT,AVT,VT,VTIN,
2     1                  VTPL,TVTPAN,HTPL,WGPL,ITY)
3C
4C  THIS SUBROUTINE COMPUTES THE SUBSONIC LATERAL STABILITY DERIVATIVES
5C
6      COMMON /OVERLY/ NLOG,NMACH,I,NALPHA,IG
7      COMMON /BDATA/  BD(762)
8      COMMON /IBODY/  PBODY,BODY(400)
9      COMMON /SYNTSS/ SYNA(19)
10      COMMON /VTI/    XXXX(154), TVTIN(8)
11      COMMON /HTI/    HTIN(154)
12      COMMON /FLOLOG/ FLTC,OPTI,BO,XX(14),TRANSN
13      COMMON /FLGTCD/ FLC(160)
14      COMMON /CONSNT/ PI,DEG,UNUSED,RAD
15      COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF
16C
17      LOGICAL FLTC,OPTI,BO,WGPL,TVTPAN,TRANSN,HTPL,VTPL
18C
19      DIMENSION STB(135),WING(380),WINGIN(100),BW(380),A(195)
20      DIMENSION B(49), AHT(195), AVT(195), VT(380), VTIN(100)
21      DIMENSION ROUTID(2),Q12127(3),Q2128A(3),Q2128B(3),Q2130B(3),
22     1          Q2130A(3),Q12129(3),Q12131(3),Q52318(3),Q22126(3),
23     2          Q1122A(3),Q1122B(3),Q1122C(3),Q1122D(3),Q1124A(3),
24     3          Q1124B(3),Q1124C(3)
25      DIMENSION WTYPE(4)
26      DIMENSION C4(6),C5(6),C6(6),C7(6),C8(6),C9(6),C1(6),C2(6)
27      DIMENSION X327(3),X227(8),X127(5),Y27(120)
28      DIMENSION X128A(7),X228A(10),Y28A(70)
29      DIMENSION X128B(3),X228B(9),Y28B(27)
30      DIMENSION X129(3),X229(9),X329(3),Y29(81)
31      DIMENSION X130A(5),X230A(9),Y30A(45)
32      DIMENSION X130B(4),X230B(9),Y30B(36)
33      DIMENSION X1526(7),X2526(9),Y526(63)
34      DIMENSION X131(4),X231(11),X331(3),X431(3),Y31(396),T31A(132),
35     1          T31B(132),T31C(132)
36      DIMENSION X158A(10),X258A(2),Y58A(20),X158B(5),X258B(3),Y58B(15),
37     1          X158C(5),X258C(2),Y58C(10)
38      DIMENSION X5322D(13),Y5322D(13),X122A(2),X222A(19),Y22A(38),X122B(
39     14),X222B(13),Y22B(52),X5322C(11),Y5322C(11),X5324A(11),Y5324A(11),
40     2X124B(2),X224B(10),Y24B(20),X124C(5),X224C(11),Y24C(55)
41C
42      EQUIVALENCE(Y31(1),T31A(1)),(Y31(133),T31B(1)),(Y31(265),T31C(1))
43C
44      DATA ROUTID/4HSUBL,4HAT  /,Q12127/4H5.1.,4H2.1-,4H27  /,
45     1 Q2128A/4H5.1.,4H2.1-,4H28-A/,Q2130B/4H5.1.,4H2.1-,4H30-B/,
46     2 Q2128B/4H5.1.,4H2.1-,4H28-B/,Q2130A/4H5.1.,4H2.1-,4H30-A/,
47     3 Q12129/4H5.1.,4H2.1-,4H29  /,Q12131/4H5.1.,4H2.1-,4H31  /,
48     4 Q52318/4H5.2.,4H3.1-,4H8   /,Q22126/4H5.2.,4H2.1-,4H26  /,
49     5 Q1122A/4H5.3.,4H1.1-,4H22A /,Q1124A/4H5.3.,4H1.1-,4H24A /,
50     6 Q1122B/4H5.3.,4H1.1-,4H22B /,Q1124B/4H5.3.,4H1.1-,4H24B /,
51     7 Q1122C/4H5.3.,4H1.1-,4H22C /,Q1124C/4H5.3.,4H1.1-,4H24C /,
52     8 Q1122D/4H5.3.,4H1.1-,4H22D /
53      DATA WTYPE/4HSTRA,4HDOUB,4HCRAN,4HCURV/
54C
55C---------          FIGURE 5.1.2.1-27
56C
57      DATA X327/0.0,0.5,1.0/
58      DATA X127/1.0,2.0,4.0,6.0,8.0/
59      DATA X227/-20.,0.0,20.,30.,40.,50.,55.,60./
60      DATA Y27/.0014,0.0,-.00125,-.002,-.0027,-.0036,-.004,-.0044 ,
61     1.0015,0.0,-.00145,-.0022,-.003,-.0041,-.005,-.00595 , .0016,0.0,
62     2-.0016,-.0024,-.0033,-.0047,-.0057,-.0071 , .0016,0.0,-.0016,
63     3-.0024,-.0035,-.0049,-.006,-.0074 , .0016,0.0,-.0016,-.0027,-.0035
64     4,-.0049,-.006,-.0074 , .0012,0.0,-.0012,-.0019,-.0026,-.0034,
65     5-.0039,-.0044 , .0013,0.0,-.0013,-.0021,-.003,-.0043,-.00515,
66     6-.0064 , .0015,0.0,-.0014,-.0024,-.0036,-.005,-.00605,-.0075 ,
67     7.00165,0.0,-.0016,-.0025,-.0038,-.0054,-.0066,-.0082 , .0018,0.0,
68     8-.00175,-.0027,-.004,-.0058,-.007,-.0089 , .00105,0.0,-.001,-.0016
69     9,-.0023,-.003,-.0035,-.0038 , .0012,0.0,-.0013,-.0021,-.0031,
70     A-.00435,-.00505,-.0062 , .0014,0.0,-.00165,-.00245,-.0036,-.0052,
71     B-.0061,-.0078 , .00167,0.0,-.0017,-.0028,-.004,-.00595,-.00715,
72     C-.009 , .0018,0.0,-.0018,-.00295,-.0042,-.0062,-.0078,-.010/
73C
74C---------          FIGURE 5.1.2.1-28-A                     ------------
75C
76      DATA X128A/2.,3.,4.,5.,6.,8.,10./
77      DATA X228A/0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,0.95/
78      DATA Y28A/8*1.0,0.995,0.990 , 6*1.0,1.01,2*1.03,1.02 , 2*1.0,1.01,
79     11.015,1.025,1.05,1.08,1.09,1.10,1.10 , 1.0,1.01,1.015,1.02,1.05,
80     21.09,1.115,1.16,1.20,1.21 , 1.0,1.01,1.02,1.04,1.07,1.12,1.17,1.24
81     3,1.32,1.36 , 1.0,1.01,1.05,1.07,1.12,1.18,1.27,1.40,1.58,1.70 ,
82     41.0,1.02,1.05,1.10,1.15,1.23,1.37,1.54,1.84,2.08/
83C
84C---------          FIGURE 5.1.2.1-28-B                     ------------
85C
86      DATA X128B/0.0,0.50,1.0/
87      DATA X228B/1.0,1.5,2.0,2.5,3.0,4.0,5.0,6.0,8.0/
88      DATA Y28B/-.0058,-.00345,-.00235,-.00145,-.001,-.00045,-.00025,
89     1.00005,.0004 , -.008,-.00555,-.004,-.003,-.00235,-.0014,-.001,
90     2-.00065,-.0002 , -.0113,-.008,-.00595,-.00465,-.0037,-.00255,
91     3-.00182,-.00147,-.00097/
92C
93C---------          FIGURE 5.1.2.1-29                       ------------
94C---------     NOTE%  INPUT ABSOLUTE VALUE OF THE           ------------
95C---------            SWEEP ANGLE FOR FIG. 5.1.2.1-29       ------------
96C---------                 ( THE X129 ARGUMENT )            ------------
97C
98      DATA X329/0.0,0.5,1.0/
99      DATA X129/0.0,40.0,60.0/
100      DATA X229/0.0,1.0,2.0,3.0,4.0,5.0,6.0,8.0,10./
101      DATA Y29/0.0,-.000052,-.000088,-.00011,-.000134,-.000153,-.000168,
102     1-.00019,-.0002 , 0.0,-.000048,-.000085,-.000108,-.000128,-.000141,
103     2-.000153,-.000173,-.000178 , 0.0,-.00004,-.000073,-.000095,
104     3-.000108,-.000119,-.000127,-.000135,-.000138 , 0.0,-.000052,
105     4-.000098,-.000132,-.000162,-.000186,-.000208,-.00024,-.00026 ,
106     50.0,-.00005,-.000096,-.000124,-.00015,-.00017,-.000188,-.000217,
107     6-.00023 , 0.0,-.00005,-.000087,-.000111,-.000129,-.000142,
108     7-.000153,-.000166,-.000170 , 0.0,-.00005,-.000096,-.000133,
109     8-.000167,-.000193,-.000216,-.000252,-.00028 , 0.0,-.00005,-.000095
110     9,-.000129,-.000155,-.000178,-.000197,-.000225,-.000245 , 0.0,
111     A-.00005,-.000088,-.000113,-.000132,-.000147,-0.000159,-.000172,
112     B-.00018/
113C
114C---------          FIGURE 5.1.2.1-30-A
115C
116      DATA X130A/2.,4.,6.,8.,10./
117      DATA X230A/0.0,.2,.4,.5,.6,.7,.8,.9,.95/
118      DATA Y30A/1.0,1.01,1.018,1.02,1.023,1.03,1.04,1.05,1.057 , 1.0,
119     11.012,1.03,1.045,1.06,1.085,1.118,1.16,1.19 , 1.0,1.015,1.045,
120     21.07,1.1,1.14,1.197,1.27,1.33 , 1.0,1.018,1.05,1.085,1.125,1.19,
121     31.26,1.39,1.485 , 1.0,1.02,1.058,1.097,1.148,1.215,1.325,1.495,
122     41.635/
123C
124C---------          FIGURE 5.1.2.1-30-B                     ------------
125C
126      DATA X130B/0.0,0.4,0.6,1.0/
127      DATA X230B/3.,4.,5.,6.,7.,8.,9.,10.,11./
128      DATA Y30B/-.0000192,-.0000222,-.0000238,-.0000231,-.0000230,
129     1-.0000241,-.000026,-.0000284,-.0000328 , -.000022,-.0000287,
130     2-.0000323,-.0000335,-.0000339,-.0000342,-.000035,-.000037,-.000042
131     3 , -.0000233,-.00003,-.0000335,-.000035,-.0000366,-.000037,
132     4-.0000375,-.00004,-.000047 , -.0000233,-.00003,-.0000335,-.000035,
133     5-.0000366,-.000037,-.0000375,-.00004,-.000047/
134C
135C---------          FIGURE 5.1.2.1-31                       ------------
136C
137      DATA X431 /0.0,0.5,1.0/
138      DATA X331 /2.,4.,8./
139      DATA X131 /-40.,0.0,40.,60./
140      DATA X231 /0.,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/
141      DATA T31A/0.,-.0001,-.00026,-.00056,-.001,-.0015,-.00217,-.0028,
142     1-.00336,-.0039,-.00435 , 0.,-.0001,-.00032,-.0006,-.0011,-.00176,
143     2-.0025,-.0032,-.00385,-.00438,-.0048 , 0.,-.0001,-.0003,-.00063,
144     3-.00117,-.0018,-.00256,-.0033,-.00395,-.00448,-.00496 , 0.0,
145     4-.0001,-.0003,-.0006,-.00115,-.0017,-.00236,-.0031,-.0037,-.00417,
146     5-.00450 , .0,-.0001,-.00025,-.00078,-.00153,-.0024,-.00335,-.0043,
147     6-.0051,-.00586,-.0066 , 0.,-.0001,-.00038,-.001,-.00183,-.0029,
148     7-.00405,-.0052,-.00633,-.00726,-.0079 , 0.,-.0001,-.0004,-.00102,
149     8-.00186,-.00303,-.00425,-.0054,-.00655,-.00736,-.00775 , 0.,-.0001
150     9,-.00040,-.001,-.00177,-.00275,-.0038,-.00476,-.00565,-.00635,
151     A-.00653 , 0.,-.00008,-.00057,-.00145,-.00255,-.00374,-.00484,
152     B-.00585,-.00682,-.0077,-.00845 , 0.,-.0001,-.00063,-.00165,-.003,
153     C-.00445,-.00625,-.00786,-.0093,-.01040,-.01110 , 0.,-.00013,-.0009
154     D,-.00195,-.00325,-.00483,-.00655,-.00815,-.0095,-.01035,-.0106 ,
155     E0.,-.00013,-.0009,-.00195,-.00325,-.0046,-.00585,-.00697,-.00787,
156     F-.0084,-.00842/
157      DATA T31B/0.,-.00015,-.0003,-.00065,-.0011,-.0017,-.0025,-.00325,
158     1-.004,-.00453,-.005 , 0.0,-.00015,-.00037,-.00076,
159     2-.00128,-.0019,-.00276,-.00355,-.00426,-.00493,-.00536 , 0.,
160     3-.00015,-.00037,-.00076,-.00128,-.0019,-.00276,-.00355,-.00426,
161     4-.0048,-.0052 , 0.0,-.00014,-.00036,-.0007,-.00125,-.0019,-.0027,
162     5-.0034,-.004,-.00452,-.00508 , 0.,-.0001,-.0004,-.00096,-.00163,
163     6-.00255,-.0038,-.0051,-.00637,-.00743,-.00825 , 0.,-.0002,-.00065,
164     7-.00125,-.00215,-.0032,-.0045,-.00595,-.0072,-.00845,-.00935 , 0.,
165     8-.0002,-.00065,-.00125,-.00215,-.0032,-.0045,-.00595,-.0072,
166     9-.00823,-.009 , 0.,-.0001,-.0005,-.00113,-.00185,-.00285,-.004,
167     A-.0051,-.00605,-.0069,-.00745 , 0.,-.00015,-.0006,-.00125,-.00215,
168     B-.00327,-.0047,-.00645,-.0084,-.0100,-.01140 , 0.,-.0003,-.00094,
169     C-.00192,-.00315,-.0046,-.0068,-.0089,-.0109,-.01255,-.0141 , 0.,
170     D-.0003,-.0009,-.00195,-.00325,-.00482,-.00675,-.00855,-.01015,
171     E-.01155,-.01265 , 0.,-.00015,-.0009,-.00176,-.00285,-.00405,-.0053
172     F,-.00665,-.0078,-.00885,-.0097/
173      DATA T31C/0.,-.00015,-.00038,-.0007,-.00126,-.0019,-.00268,-.00344
174     1,-.00418,-.00484,-.0053 , 0.,-.00016,-.0004,-.00076,-.0013,-.002,
175     2-.00288,-.0037,-.00445,-.0051,-.0056 , 0.,-.00015,-.00035,-.0007,
176     3-.00125,-.0019,-.00278,-.00356,-.0043,-.005,-.0054 , 0.,-.00012,
177     4-.0003,-.00067,-.0012,-.0018,-.00257,-.0033,-.00396,-.00458,
178     5-.00505 , 0.,-.00015,-.0005,-.00115,-.0018,-.00273,-.00385,-.00515
179     6,-.0066,-.0079,-.0088 , 0.,-.00015,-.00052,-.0012,-.0021,-.0033,
180     7-.00486,-.00627,-.00756,-.0087,-.00972 , 0.,-.00015,-.00052,-.0012
181     8,-.0021,-.0033,-.00465,-.006,-.00728,-.00836,-.00925 , 0.,-.00015,
182     9-.0005,-.00116,-.0019,-.00285,-.004,-.0051,-.0061,-.00705,-.0079 ,
183     A-0.,-.00018,-.0006,-.00125,-.0021,-.00335,-.00505,-.0071,-.00915,
184     B-.01100,-.01255 , 0.,-.00025,-.0009,-.0018,-.0031,-.0048,-.007,
185     C-.0094,-.01145,-.0133,-.01475 , 0.,-.0003,-.0009,-.00195,-.00327,
186     D-.0048,-.00697,-.0089,-.01057,-.01205,-.0132 , 0.,-.00012,-.00073,
187     E-.00155,-.00265,-.004,-.00536,-.0067,-.00795,-.009,-.010/
188C
189C                        FIGURE 5.2.2.1-26
190C
191      DATA X1526/4.0,4.5,5.0,5.5,6.0,7.0,8.0/
192      DATA X2526/0.0,.2,.4,.6,.8,1.0,1.2,1.4,1.6/
193      DATA Y526/1.0,1.0,1.0,1.0,1.0,1.0,1.0,.99,.97,
194     11.0,1.0,1.0,1.0,1.0,1.0,.98,.948,.911,
195     21.0,1.0,1.0,1.0,.997,.971,.933,.883,.827,
196     31.0,1.0,1.0,.991,.963,.922,.870,.811,.746,
197     41.0,1.0,.995,.970,.932,.884,.829,.764,.695,
198     51.0,1.0,.977,.944,.899,.845,.780,.715,.641,
199     61.0,.985,.960,.921,.870,.812,.745,.670,.592/
200C
201C                         FIGURE 5.2.3.1-8A
202C
203      DATA X158A/20.,14.,10.,8.,7.,6.,5.,4.,3.,2.5/
204      DATA X258A/.2,.8/
205      DATA Y58A/.1,1.88,
206     1.40,2.21,
207     2.74,2.60,
208     3.98,2.80,
209     41.30,3.13,
210     51.61,3.50,
211     62.00,3.88,
212     72.50,4.40,
213     82.99,5.00,
214     93.45,5.40/
215C
216C                        FIGURE 5.2.3.1-8B
217C
218      DATA X158B/.8,1.0,1.2,1.4,1.6/
219      DATA X258B/0.0,3.0,6.0/
220      DATA Y58B/0.0,2.35,4.68,
221     10.0,3.00,6.00,
222     20.0,3.60,7.25,
223     30.0,4.18,8.50,
224     40.0,4.79,9.50/
225C
226C                        FIGURE 5.2.3.1-8C
227C
228      DATA X158C/.5,.6,.8,1.,2./
229      DATA X258C/0.0,6.0/
230      DATA Y58C/-.00048,.00251,
231     1-.00048,.0035,
232     2-.00048,.00477,
233     3-.00048,.00559,
234     4-.00048,.00641/
235C
236C                        FIGURE 5.3.1.1-22A
237C
238      DATA X122A/1.0,.6/
239      DATA X222A/0.,.125,.25,.5,.75,1.,1.25,1.5,1.75,2.,2.25,2.5,
240     13.00,3.25,3.50,3.75,4.00,5.00,7.00/
241      DATA Y22A/0.,.40,.720,.990,1.19,1.32,1.40,1.46,1.50,1.51,1.48,1.42
242     1,1.27,1.21,1.17,1.13,1.10,1.04,1.02,
243     20.,.70,.940,1.18,1.35,1.46,1.54,1.60,1.63,1.64,1.60,1.53,1.36,
244     31.28,1.21,1.16,1.13,1.06,1.02/
245C
246C                        FIGURE 5.3.1.1-22B
247C
248      DATA X122B/.5,.6,.7,.8/
249      DATA X222B/0.0,-.2,-.3,-.4,-.5,-.6,-.65,-.7,-.75,-.8,-.85,-.9,-1./
250      DATA Y22B/1.05,.94,.90,.87,.86,.87,.90,.93,.98,1.06,1.16,1.29,1.7,
251     11.15,1.00,.95,.90,.89,.90,.92,.96,1.01,1.08,1.18,1.31,1.70,
252     21.22,1.05,.99,.94,.92,.92,.95,.98,1.03,1.10,1.20,1.33,1.70,
253     31.29,1.09,1.02,.97,.94,.94,.96,1.0,1.06,1.12,1.22,1.36,1.70/
254C
255C                        FIGURE 5.3.1.1-22C
256C
257      DATA X5322C/0.0,.2,.4,.6,.7,.8,.9,1.2,1.4,1.6,2.0/
258      DATA Y5322C/0.0,.29,.52,.70,.77,.83,.87,.98,1.04,1.07,1.13/
259C
260C                        FIGURE 5.3.1.1-22D
261C
262      DATA X5322D/0.0,.50,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0/
263      DATA Y5322D/5*.75,.835,.92,6*1.0/
264C
265C                        FIGURE 5.3.1.1-24A
266C
267      DATA X5324A/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/
268      DATA Y5324A/1.5,1.33,1.197,1.098,1.03,1.0,1.03,1.098,1.197,1.33,
269     11.5/
270C
271C                        FIGURE 5.3.1.1-24B
272C
273      DATA X124B/20.0,0.0/
274      DATA X224B/.5,1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0/
275      DATA Y24B/.90,1.64,2.19,2.63,2.95,3.22,3.44,3.61,3.75,3.86,
276     11.0,1.72,2.30,2.75,3.11,3.39,3.62,3.82,3.98,4.11/
277C
278C                        FIGURE 5.3.1.1-24C
279C
280      DATA X124C/.2,.4,.6,.8,1.0/
281      DATA X224C/0.0,.1,.2,.3,.4,.5,.6,.7,.8,.9,1.0/
282      DATA Y24C/1.0,.90,.818,.741,.68,.631,.593,.563,.542,.527,.510,
283     11.0,.93,.87,.825,.784,.752,.730,.710,.696,.688,.680,
284     21.0,.95,.911,.882,.859,.844,.830,.820,.812,.805,.800,
285     31.0,.968,.941,.923,.908,.898,.890,.883,.881,.880,.879,
286     41.0,.979,.962,.954,.949,.943,.938,.935,.935,.931,.931/
287      DATA I4,I5,I6/3*0/
288C
289      IF(.NOT. WGPL) GO TO 1255
290      YA27 =UNUSED
291      YA28A=UNUSED
292      YA28B=UNUSED
293      YA29 =UNUSED
294      YA30A=UNUSED
295      YA30B=UNUSED
296      YA31I=UNUSED
297      YA31O=UNUSED
298C
299C                   LATERAL STABILITY DERIVATIVES
300C
301      IF(TRANSN)GO TO 1190
302C
303C     ***** WING ALONE CALCULATIONS *****
304C
305C     WING ALONE SIDEFORCE DERIVATIVE
306C
307      BDAT=2.*WINGIN(4)
308      ANUM1=6.0*A(44)*A(42)
309      DEN1=PI*A(120)*(A(120)+4.*A(43))
310      BB=SQRT(1.-(FLC(I+2)*A(43))**2)
311      ANUM2=A(120)+4.0*A(43)
312      DEN2=A(120)*BB+4.0*A(43)
313      DO 1050 J=1,NALPHA
314         STB(J+15)=(1./RAD)*(B(J+2)**2)*(ANUM1/DEN1)-.0001*ABS(STB(122))
315         IF(B(J+2).EQ.0.0)GO TO 1000
316         STB(J+35)=(ANUM2/DEN2)*(STB(J+15)/B(J+2))
317         WING(J+140)=STB(J+35)*WING(J+20)
318         GO TO 1010
319 1000    STB(J+35)=0.0
320         WING(J+140)=0.0
321 1010    CONTINUE
322C
323C     WING ALONE YAWING MOMENT DERIVATIVE
324C
325         IF(J.GT.1)GO TO 1020
326         APART=1.0/(4.0*PI*A(120))
327         BPART=ABS(A(44))/(PI*A(120)*(A(120)+4.0*A(43)))
328         CPART=A(43)-0.5*A(120)-A(120)**2/(8.*A(43))
329         EPART=(A(120)+4.*A(43))/(A(120)*BB+4.*A(43))
330         ANUM3=(A(120)*BB)**2+4.*A(120)*BB*A(43)-8.0*A(43)*A(43)
331         DEN3=A(120)*A(120)+4.0*A(120)*   A(43)-8.0*A(43)*A(43)
332 1020    CONTINUE
333         IF(WING(J+60).EQ.0.0)GO TO 1030
334         DPART=6.0*(WING(J+40)/WING(J+60)/CBARR)*(ABS(A(42))/A(120))
335         GO TO 1040
336 1030    DPART=0.0
337 1040    CONTINUE
338         STB(J+75)=(APART-BPART*(CPART+DPART))/RAD
339         WING(J+160)=EPART*(ANUM3/DEN3)*STB(J+75)*WING(J+20)**2
340         WING(J+160)=WING(J+160)*BDAT*SREF/(BLREF*A(3))
341 1050 CONTINUE
342C
343C---------     TEST FOR CONFIGURATION AND UNIFORM DIHEDRAL ANGLE -------
344C
345C     WING ALONE ROLLING MOMENT DERIVATIVE
346C     (FOR STRAIGHT WING CONFIGURATION WITH UNIFORM DIHEDRAL)
347C
348      IF(WINGIN(15).NE.WTYPE(1)) GO TO 1120
349      IF(A(120).GE.1.0) GO TO 1070
350      DO 1060 J=1,NALPHA
351         WING(J+180)=-2./(RAD*3.*A(120))*WING(J+20)-STB(122)*A(120)/
352     1           (6.*RAD**2)*(A(4)/SREF)
353 1060 WING(J+180)=WING(J+180)*BDAT/ BLREF
354      GO TO 1180
355 1070 CONTINUE
356      XA2= FLC(I+2)*A(73)
357      X1ARG=ABS(A(70))
358      XA1=A(120)/A(73)
359C
360C     ---------     FIGURE 5.1.2.1-27        ----------
361C
362      CALL TLIN3X(X127,X227,X327,Y27,5,8,3,A(120),X1ARG,A(118),YA27,
363     1            0,1,0,2,2,0,Q12127,3,ROUTID)
364      IF(A(70).LT.0.0) YA27=ABS(YA27)
365C
366C     ---------     FIGURE 5.1.2.1-28-A      ----------
367C
368      CALL TLINEX(X128A,X228A,Y28A,7,10,XA1,XA2,YA28A,
369     1            0,2,1,2,Q2128A,3,ROUTID)
370C
371C     ---------     FIGURE 5.1.2.1-28-B      ----------
372C
373      CALL TLINEX(X128B,X228B,Y28B,3,9,A(118),A(120),YA28B,
374     1         0,0,0,0,Q2128B,3,ROUTID)
375C
376C     ---------     FIGURE 5.1.2.1-30-B      ----------
377C
378      CALL TLINEX(X130B,X230B,Y30B,4,9,A(118),A(120),YA30B,
379     1            0,2,0,2,Q2130B,3,ROUTID)
380      DO 1080 J=1,NALPHA
381 1080 STB(J+95)=(YA27*YA28A+YA28B)*WING(J+20)*SREF/A(4)+WINGIN(11)*
382     1          (A(68))*YA30B
383      IF(WINGIN(12).NE.UNUSED.AND.WINGIN(12).NE.0.0)GO TO 1100
384C
385C     ---------     FIGURE 5.1.2.1-30-A      ----------
386C
387      CALL TLINEX(X130A,X230A,Y30A,5,9,XA1,XA2,YA30A,
388     1            2,2,2,2,Q2130A,3,ROUTID)
389C
390C     ---------     FIGURE 5.1.2.1-29        ----------
391C
392      CALL TLIN3X(X129,X229,X329,Y29,3,9,3,X1ARG,A(120),A(118),YA29,
393     1            0,0,0,2,2,0,Q12129,3,ROUTID)
394      DO 1090 J=1,NALPHA
395         WING(J+180)=(STB(J+95)+WINGIN(13)*YA29*YA30A)*A(4)/SREF
396         WING(J+180)=WING(J+180)*BDAT/BLREF
397 1090 CONTINUE
398      GO TO 1180
399C
400C     WING ALONE ROLLING MOMENT DERIVATIVE
401C     (FOR STRAIGHT WING CONFIGURATION WITH NON-UNIFORM DIHEDRAL)
402C
403 1100 CONTINUE
404      STB(7)=A(131)/(2.0*PI)*RAD
405      X1ARG=ATAN(A(68)/B(2))*RAD
406      X3ARG=B(2)*A(120)/STB(7)
407      AKOVB=STB(7)/B(2)
408C
409C     ---------     FIGURE 5.1.2.1-31        ----------  (INBOARD)
410C
411      CALL TLIN4X(X131,X231,X331,X431,Y31,4,11,3,3,X1ARG,STB(2),X3ARG,
412     1            A(118),YA31I,0,0,1,0,2,2,1,0,Q12131,3,ROUTID)
413C
414C     ---------     FIGURE 5.1.2.1-31        ----------  (OUTBOARD)
415C
416      CALL TLIN4X(X131,X231,X331,X431,Y31,4,11,3,3,X1ARG,STB(3),X3ARG,
417     1            A(118),YA31O,0,0,1,0,2,2,1,0,Q12131,3,ROUTID)
418      DO 1110 J=1,NALPHA
419 1110 WING(J+180)=(STB(J+95)+(YA31O-YA31I)*AKOVB*WINGIN(14)/RAD+YA31I/
420     1            RAD*AKOVB*WINGIN(13))*A(4)/SREF*BDAT/BLREF
421      GO TO 1180
422C
423C     WING ALONE ROLLING MOMENT DERIVATIVE
424C              DOUBLE DELTA OR CRANKED WING CONFIGURATION
425C     (INBOARD PANEL ROLLING MOMENT)
426C
427 1120 CONTINUE
428      X1ARG=ABS(A(70))
429      STB(135)=-2.0/(3.0*RAD*A(163))
430C
431C     ---------     FIGURE 5.1.2.1-27        -----------
432C
433      CALL TLIN3X(X127,X227,X327,Y27,5,8,3,A(163),X1ARG,A(25),STB(132),
434     1            0,1,0,2,2,0,Q12127,3,ROUTID)
435      IF(A(70).LT.0.0) YA27=ABS(YA27)
436C
437C     ---------     FIGURE 5.1.2.1-28-A      ----------
438C
439      XA1I=A(163)/A(73)
440      XA2I=A(73)*FLC(I+2)
441      CALL TLINEX(X128A,X228A,Y28A,7,10,XA1I,XA2I,STB(134),
442     1            0,2,1,2,Q2128A,3,ROUTID)
443C
444C     ---------     FIGURE 5.1.2.1-28-B      ----------
445C
446      CALL TLINEX(X128B,X228B,Y28B,3,9,A(25),A(163),STB(133),
447     1         0,0,0,0,Q2128B,3,ROUTID)
448C
449C     (OUTBOARD PANEL ROLLING MOMENT)
450C
451      STB(131)=-2.0/(3.*RAD*A(168))
452      X1ARG=ABS(A(94))
453C
454C     ---------     FIGURE 5.1.2.1-27        ----------
455C
456      CALL TLIN3X(X127,X227,X327,Y27,5,8,3,A(168),X1ARG,A(169),STB(128),
457     1            0,1,0,2,2,0,Q12127,3,ROUTID)
458      IF(A(94).LT.0.0) YA27=ABS(YA27)
459C
460C     ---------     FIGURE 5.1.2.1-28-A      ----------
461C
462      XA1O=A(168)/A(97)
463      XA2O=A(97)*FLC(I+2)
464      CALL TLINEX(X128A,X228A,Y28A,7,10,XA1O,XA2O,STB(130),
465     1            0,2,1,2,Q2128A,3,ROUTID)
466C
467C     ---------     FIGURE 5.1.2.1-28-B      ----------
468C
469      CALL TLINEX(X128B,X228B,Y28B,3,9,A(169),A(168),STB(129),
470     1         0,0,0,0,Q2128B,3,ROUTID)
471      STB(58)=A(171)*(A(1)/SREF)+A(172)*(A(167)/SREF)
472      IF(A(163).GE.1.0) GO TO 1130
473      STB(127)=A(171)*(A(1)/SREF)*STB(135)*(WINGIN(3)-WINGIN(2))/
474     1  WINGIN(4)
475      GO TO 1140
476 1130 STB(127)=A(171)*(A(1)/SREF)*(STB(132)*STB(134)+STB(133))*
477     1(WINGIN(3)-WINGIN(2))/WINGIN(4)
478 1140 IF(A(168).GE.1.0)GO TO 1150
479      STB(126)=A(172)*(A(167)/SREF)*STB(131)
480      GO TO 1160
481 1150 STB(126)=A(172)*(A(167)/SREF)*(STB(128)*STB(130)+STB(129))
482 1160 CONTINUE
483      DO 1170 J=1,NALPHA
484 1170 WING(J+180)=(WING(J+20)/STB(58))*(STB(127)+STB(126))*BDAT/BLREF
485 1180 CONTINUE
486      STB(68)=YA27
487      STB(71)=YA28A
488      STB(70)=YA28B
489      STB(67)=YA29
490      STB(66)=YA30A
491      STB(69)=YA30B
492      STB(64)=YA31I
493      STB(65)=YA31O
494      IF(BO)GO TO 1190
495      GO TO 1330
496C
497C     ***** WING BODY CALCULATIONS *****
498C
499C     WING-BODY SIDEFORCE DERIVATIVE(ALL PLANFORMS)
500C
501 1190 CONTINUE
502      ARG=2.0*STB(1)/STB(72)
503      IF(ARG.LT.0.0)GO TO 1200
504C
505C                        FIGURE 5.2.1.1-7
506C
507      STB(57)=1.0+.49*ARG
508      GO TO 1210
509 1200 STB(57)=1.0-.85*ARG
510 1210 BW(141)=-STB(57)*BODY(101)-0.0001*ABS(STB(122))
511C
512C     WING-BODY YAWING MOMENT DERIVATIVE(ALL PLANFORMS)
513C
514      ARG1=FLC(I+42)*BD(1)
515      STB(56)=1.+ALOG(1.E-6*ARG1)/4.86
516      ARG11=STB(59)/STB(123)
517      ARG12=SQRT(STB(61)/STB(60))
518      ARG13=BD(1)**2/STB(62)
519      ARG14=SYNA(1)/BD(1)
520C                        FIGURE 5.2.3.1-8
521C
522      CALL TLINEX(X158A,X258A,Y58A,10,2,ARG13,ARG14,YDUMY,
523     1            2,1,2,1,Q52318,3,ROUTID)
524      CALL TLINEX(X158B,X258B,Y58B,5,3,ARG12,YDUMY,YDUMY2,
525     1            2,0,2,1,Q52318,3,ROUTID)
526      CALL TLINEX(X158C,X258C,Y58C,5,2,ARG11,YDUMY2,STB(15),
527     1            2,0,2,1,Q52318,3,ROUTID)
528      BW(161)=-STB(15)*STB(56)*(STB(62)*BD(1))/(SREF*BLREF)
529      IF(.NOT.TRANSN)GO TO 1220
530      RETURN
531 1220 CONTINUE
532C
533C     WING-BODY ROLLING MOMENT DERIVATIVE(ALL PLANFORMS,UNIFORM AND NON-
534C     UNIFORM DIHEDRAL)
535C
536      ANUM=1.2*STB(4)*STB(72)*SQRT(A(120))
537      ADEN=RAD*2.0*  WINGIN(4)**2
538      STB(14) = ANUM/ADEN
539      STB(13)=-.0005*SQRT(A(120))*STB(122)*(STB(72)/(2.0*WINGIN(4)))**2
540C
541C                        FIGURE 5.2.2.1-26
542C
543      ARG1=A(120)/A(49)
544      ARG2=STB(63)/(2.0*  WINGIN(4))
545      CALL TLINEX(X1526,X2526,Y526,7,9,ARG1,ARG2,STB(8),
546     1            0,0,0,1,Q22126,3,ROUTID)
547C
548      IF(WINGIN(15).EQ.WTYPE(1))GO TO 1230
549      APART= STB(8)*(STB(132)*STB(134)+STB(128)*STB(130))+STB(133)
550     1      +STB(129)
551      FPART=0.0
552      EPART=STB(13)
553      GO TO 1240
554 1230 CONTINUE
555      APART=STB(68)*STB(71)*STB(8)+STB(70)
556      FPART=WINGIN(11)*A(68)*STB(69)
557      EPART=WING(181)*BLREF/BDAT-(STB(96)-STB(13))*A(4)/SREF
558 1240 CONTINUE
559      DO 1250 J=1,NALPHA
560 1250 BW(J+180)=(BW(J+20)*APART+EPART+(STB(14)+FPART)*A(4)/SREF)*
561     1          BDAT/BLREF
562 1255 CONTINUE
563C
564C   *****WING-BODY-TAIL CALCULATIONS*****
565C
566      IF(VTPL.OR.TVTPAN)GO TO 1260
567      GO TO 1330
568 1260 CONTINUE
569      VT(141) = 0.0
570      VT(161) = 0.0
571      DO 1270 J=1,NALPHA
572 1270 VT(J+180) = 0.0
573      IF(.NOT.VTPL) GO TO 1310
574C
575C              SINGLE VERTICAL PANEL ON PLANE OF SYMMETRY
576C                        FIGURE 5.3.1.1-22A
577C
578      ARG2=VTIN(4)/(2.0*(VTIN(4)-VTIN(3)))
579      CALL TLINEX(X122A,X222A,Y22A,2,19,AVT(118),ARG2,STB(120),
580     1            0,0,0,0,Q1122A,3,ROUTID)
581C
582C                        FIGURE 5.3.1.1-22D
583C
584      CALL TBFUNX(ARG2,STB(118),DYDX,13,X5322D,Y5322D,C6,I6,MI,NG,
585     1            0,0,Q1122D,3,ROUTID)
586      ZH = SYNA(7)-((HTIN(4)-HTIN(3))*AHT(62)+AHT(30)-AHT(16)/4.0)*
587     1     SIN(SYNA(8)/RAD)
588C
589      IF(.NOT.HTPL.OR.ZH.LT.0.0) GO TO 1280
590C
591C                        FIGURE 5.3.1.1-22B
592C
593      ARG1=STB(9)/STB(10)
594      ARG2=-ZH/ VTIN(4)
595      CALL TLINEX(X122B,X222B,Y22B,4,13,ARG1,ARG2,STB(121),
596     1            0,0,2,0,Q1122B,3,ROUTID)
597C
598C                        FIGURE 5.3.1.1-22C
599C
600      ARG=(AHT(2)+AHT(119))/AVT(4)
601      CALL TBFUNX(ARG,STB(119),DYDX,11,X5322C,Y5322C,C5,I5,MI,NG,
602     1            0,0,Q1122C,3,ROUTID)
603C
604      GO TO 1290
605 1280 STB(121)=0.0
606      STB(119)=0.0
607 1290 STB(116)=STB(120)*AVT(120)*(1.+STB(119)*(STB(121)-1.))
608      BPART=3.06*((AVT(4)/A(4))/(1.+A(43)))
609      CPART=.4* STB(4)/STB(59)+ .009*A(120)
610      STB(117)=0.724+BPART+CPART
611      ANUM=2.*PI*STB(116)
612      APART=STB(116)**2/(VTIN(I+20)*RAD/(2.*PI))**2
613      BPART= 1.+AVT(50)**2/B(2)**2
614      STB(5)=ANUM/(2.+SQRT(APART*BPART+4.))
615      VT(141)=-STB(118)*STB(5)*STB(117)*AVT(4)/(RAD*SREF)
616      VT(161)=-VT(141)*STB(11)/BLREF
617      DO 1300 J=1,NALPHA
618 1300 VT(J+180)=(VT(141)*(STB(12)*COS(FLC(J+22)/RAD)-STB(11)*
619     1          SIN(FLC(J+22)/RAD)))/BLREF
620      IF(SYNA(18).EQ.UNUSED.AND.SYNA(19).EQ.UNUSED) GO TO 1310
621C
622C     TWIN VERTICAL TAIL CORRECTION
623C
624      IF(ITY.EQ.0) TVCORR=2*(COS(SYNA(18)/RAD))**2
625      IF(ITY.EQ.1) TVCORR=2*(COS(SYNA(19)/RAD))**2
626      VT(141)=TVCORR*VT(141)
627      VT(161)=TVCORR*VT(161)
628      DO 1305 J=1,NALPHA
629         VT(J+180)=TVCORR*VT(J+180)
630 1305 CONTINUE
631 1310 CONTINUE
632      IF(.NOT.TVTPAN) GO TO 1330
633C
634C                        FIGURE 5.3.1.1-24A
635C
636      ARG =TVTIN(1)/TVTIN(2)
637      CALL TBFUNX(ARG ,STB(75),DYDX,11,X5324A,Y5324A,C4,I4,MI,NG,
638     1            0,0,Q1124A,3,ROUTID)
639      STB(6)=TVTIN(2)**2/TVTIN(5)
640C
641C                        FIGURE 5.3.1.1-24B
642C
643      ARG2=STB(75)*STB(6)
644      STB(116) = ARG2
645      ARG1=TVTIN(6)
646      CALL TLINEX(X124B,X224B,Y24B,2,10,ARG1,ARG2,STB(73),
647     1            0,2,0,2,Q1124B,3,ROUTID)
648C
649C                        FIGURE 5.3.1.1-24C
650C
651      ARG1=TVTIN(4)/BD(1)
652      ARG2=TVTIN(3)/TVTIN(2)
653      CALL TLINEX(X124C,X224C,Y24C,5,11,ARG1,ARG2,STB(74),
654     1            2,0,0,0,Q1124C,3,ROUTID)
655C
656      TVT=-STB(74)*STB(73)*2.*TVTIN(5)/(SREF*RAD)
657      VT(141)=VT(141)+TVT
658      VT(161)=VT(161)-TVT*TVTIN(7)/BLREF
659C
660      DO 1320 J=1,NALPHA
661 1320 VT(J+180)=VT(J+180)+TVT*(TVTIN(8)*COS(FLC(J+22)/RAD)-TVTIN(7)*
662     1  SIN(FLC(J+22)/RAD))/BLREF
663 1330 CONTINUE
664      RETURN
665      END
666