1      SUBROUTINE TRANCD
2C
3C***  CALCULATES TRANSONIC WING AND WING-BODY CD
4C
5      COMMON /OVERLY/ NLOG,NMACH,I,NALPHA,IG,NF
6      COMMON /CONSNT/ PI,DEG,UNUSED,RAD
7      COMMON /FLGTCD/ FLC(95), TR
8      COMMON /BDATA/  BD(762)
9      COMMON /WINGD/  A(195)
10      COMMON /HTDATA/ AHT(195), BHT(49)
11      COMMON /HTI/    HTIN(154)
12      COMMON /IBH/    PBH, BH(380)
13      COMMON /WINGI/  WINGIN(100)
14      COMMON /SYNTSS/ XCG, XW, ZW, ALIW, ZCG, XH, ZH, ALIH, XV,
15     1                VERTUP, HINAX, XVF, SCALE, ZV, ZVF, YV, YF,
16     2                PHIV, PHIF
17      COMMON /IBW/    PBW, BW(380)
18      COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC,
19     1                HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN,SUPERS,SUBSON,
20     2                TRANSN,HYPERS,SYMFP,ASYFP,TRIMC,TRIM,DAMP,
21     3                HYPEF,TRAJET,BUILD,FIRST
22      LOGICAL  FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC,
23     1         HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN,SUPERS,SUBSON,
24     2         TRANSN,HYPERS,SYMFP,ASYFP,TRIMC,TRIM,DAMP,
25     3         HYPEF,TRAJET,BUILD,FIRST
26      REAL LN, LA
27      LOGICAL NA
28      DIMENSION CDL(20)
29      DATA STRA / 4HSTRA /
30C
31C***  WING-BODY DRAG
32C
33      IF(NF .LT. 0) GO TO 1040
34      IF(WINGIN(15) .NE. STRA .OR. .NOT. (BO .AND. WGPL)) GO TO 1020
35      CD0 = BW(1)
36      RN = FLC(I+42)*A(122)
37      DB = 2.0*(WINGIN(4)-WINGIN(3))
38      LN = (XW + 0.5*DB*A(38))/DB
39      LA = (XW + WINGIN(6) + 0.5*DB*A(56))/DB
40      IF(BD(2) .LT. (LN*DB)) LN = LN*DB/BD(85)
41      IF(BD(2) .GT. (LA*DB)) LA = LA*DB/BD(85)
42      LA = BD(1)/DB-LA
43      CALL WBCDL(A(120), A(38), WINGIN(16), LN, LA, A(118),
44     1          WINGIN(62), WINGIN(11), WINGIN(93), WINGIN(94), RN,
45     1          TR, FLC(I+2), FLC(23), CDL, NA)
46      IF(NA) GO TO 1020
47        DO 1010 J=1,NALPHA
48           BW(J) = CD0 + CDL(J)
49           IF(CDL(J) .NE. UNUSED) GO TO 1010
50           BW(J) = -UNUSED
51 1010 CONTINUE
52 1020 CONTINUE
53C
54C***  H.T.-BODY DRAG
55C
56      IF(HTIN(15) .NE. STRA .OR. .NOT. HTPL) GO TO 1040
57      CD0 = BH(1)
58      RN = FLC(I+42)*AHT(122)
59      DB = 2.0*(HTIN(4)-HTIN(3))
60      LN = (XH + 0.5*DB*AHT(38))/DB
61      LA = (XH + HTIN(6) + 0.5*DB*AHT(56))/DB
62      IF(BD(2) .LT. (LN*DB)) LN = LN*DB/BD(85)
63      IF(BD(2) .GT. (LA*DB)) LA = LA*DB/BD(85)
64      LA = BD(1)/DB-LA
65      CALL WBCDL(AHT(120), AHT(38), HTIN(16), LN, LA, AHT(118),
66     1          HTIN(62), HTIN(11), HTIN(93), HTIN(94), RN,
67     1          TR, FLC(I+2), FLC(23), CDL, NA)
68      IF(NA) GO TO 1040
69        DO 1030 J=1,NALPHA
70           BH(J) = CD0 + CDL(J)
71           IF(CDL(J) .NE. UNUSED) GO TO 1030
72           BH(J) = -UNUSED
73 1030 CONTINUE
74 1040 CONTINUE
75      RETURN
76      END
77