1      SUBROUTINE VTDRAG(I)
2C
3C *** CALCULATES VERTICAL DRAG COEFFICIENT
4C
5      DIMENSION ROUTID(2)
6      DIMENSION Q15127(3),Q5128B(3),Q128BD(5)
7      COMMON /OPTION/ SREF,CBARR,ROUGFC,BLREF
8      COMMON /FLGTCD/ FLC(93)
9      COMMON /WHAERO/ C2(51),D2(55),CHT2(51),DHT2(55),DVT(55)
10      COMMON /VTDATA/ AVT(195)
11      COMMON /VTI/    VT(100)
12      DIMENSION D(55),A(195),B(93)
13      EQUIVALENCE (A(1),AVT(1)),(D(1),DVT(1)),(B(1),FLC(1))
14      DIMENSION WTYPE(4)
15      DIMENSION X228BD(9),X128BD(4),Y28BD(36)
16      DIMENSION X228B(11), X128B(4), Y28B(44)
17      DIMENSION X27M(4),X27I(4),C27(6)
18      DATA ROUTID/4HVTDR,4HAG  /
19      DATA Q15127/4H4.1.,4H5.1-,4H27  /,Q5128B/4H4.1.,4H5.1-,4H28-B/,
20     1     Q128BD/4H4.1.,4H5.1-,4H28-B,4H(DAS,4HHED)/
21      DATA X27M /0.0,1.0,2.0,3.0/
22      DATA X27I /1.57780,1.67221,1.98509,2.28874/
23      DATA I27  /0/
24C
25C---------          FIGURE 4.1.5.1-28B            ----
26C----     X228B=COS(SWEEP ANGLE(T/C)MAX)   X128B= MACH NUMBER
27C----     Y28B= (R)LS
28C
29      DATA WTYPE / 4HSTRA ,4HDOUB ,4HCRAN ,4HCURV /
30      DATA X228B /0.5,.55,.60,.65,.70,.75,.80,.85,.90,.95,1.0/
31      DATA X128B /0.9,.80,.60,.25/
32      DATA Y28B/ 1.1,1.13,1.17,1.20,1.24,1.27,1.3,1.33,1.34,1.35,1.36 ,
33     11.0,1.04,1.08,1.11,1.15,1.18,1.21,1.23,1.25,1.25,1.26 , 0.88, .92,
34     2.96,1.0,1.04,1.08,1.11,1.13,1.14,1.14,1.15 , 0.81,.85,.89,.925,
35     3.96,1.0,1.03,1.05,1.06,1.06,1.07/
36C
37C---------          FIGURE 4.1.5.1-28B-DASHED     ---------
38C----          X228BD= COS(T/C)MAX      X128BD= MACH NUMBER
39C----          Y28BD= (R)LS
40C
41      DATA X228BD /0.45,.65,.70,.75,.80,.85,.90,.95,1.0/
42      DATA X128BD /0.9,.80,.60,.25/
43      DATA Y28BD  /1.20,1.20,1.24,1.27,1.30,1.33,1.34,1.35,1.36 , 1.11,
44     11.11,1.15,1.18,1.21,1.23,1.25,1.25,1.26 , 1.0,1.0,1.04,1.08,1.11,
45     21.13,1.14,1.14,1.15 , 0.925,0.925,.96,1.0,1.03,1.05,2*1.06,1.07/
46C
47      D(2)=12.0*A(16)/ROUGFC
48C
49C---------     DETERMINE WING CONFIGURATION       ----------------------
50C              IF NOT STRAIGHT-TAPERED, GO TO PAGE 3-2
51C
52      IF(VT(15).NE.WTYPE(1)) GO TO 1010
53C
54C---------          PAGE 3-1 CALCULATIONS         ----------------------
55C*********     FOR STRAIGHT-TAPER WING CONFIGURATION        ************
56C---------     COMPUTE REYNOLDS NUMBER, (RN)
57C
58 1000 D(17)=A(16)*A(129)
59C
60C---------     COMPUTE RL FROM EQUATION FOR FIG.4.1.5.1-27
61C*********          FIGURE 4.1.5.1-27 EQUATION
62C
63      CALL TBFUNX(B(I+2),CEPT,DYDX,4,X27M,X27I,C27,I27,MI,NG,
64     1            0,0,Q15127,3,ROUTID)
65      D(14)=D(2)**1.0482*10.0**CEPT
66C
67C*********          FIGURE 4.1.5.1-26 EQUATIONS (SUBROUTINE FIG26 )
68C
69      IF( D(14).LT.D(17)) D(17)=D(14)
70      CALL FIG26  (D(17),B(I+2),D(10) )
71C
72C*********          FIGURE 4.1.5.1-28-B
73C
74      CALL TLINEX(X128B,X228B,Y28B,4,11,B(I+2),A(178),D(13),
75     1            0,2,0,2,Q5128B,3,ROUTID)
76      CAPL=2.0
77      IF(VT(18).GE.0.30) CAPL=1.20
78      D(20)=D(10)*(1.+CAPL*VT(16)+100.*VT(16)**4)*D(13)*(2.*A(3)/SREF)
79      GO TO 1030
80C
81C                  ***** INBOARD CALCUALTIONS *****
82C
83 1010 D(18)=A(15)*A(129)
84C
85C*********     FIGURE 4.1.5.1-27 EQUATION
86C
87      CALL TBFUNX(B(I+2),CEPT,DYDX,4,X27M,X27I,C27,I27,MI,NG,
88     1            0,0,Q15127,3,ROUTID)
89      D(15)=D(2)**1.0482*10.0**CEPT
90C
91C*********     IF RLI.LT.RNI  SET RNI=RLI
92C
93      IF( D(15).LT.D(18) ) D(18)= D(15)
94C
95C*********     FIGURE 4.1.5.1-26 SUBROUTINE
96C
97      CALL FIG26 ( D(18),B(I+2),D(11) )
98C
99C---------          FIGURE 4.1.5.1-28-B (DASHED)
100C
101      CALL TLINEX(X128BD,X228BD,Y28BD,4,9,B(I+2),A(190),D(23),
102     1            0,0,0,0,Q128BD,5,ROUTID)
103C
104C*********     IF XTI GREATER TAN 0.30, THEN CAPL= 2.0          *
105C
106 1020 CAPL= 1.20
107      IF(VT(18).LT.0.30) CAPL=2.00
108      D(21)=D(11)*(1.+CAPL*VT(16)+100.*VT(16)**4)*D(23)*(2.*A(1)/SREF)
109C
110C                   ***** OUTBOARD CALCULATIONS *****
111C
112      D(19)=A(17)*A(129)
113C
114C*********     FIGURE 4.1.5.1-27 EQUATION
115C
116      D(16)=D(15)
117C
118C*********   IF RLO .LT. RNO, SET RNO= RLO
119C
120      IF(D(16).LT.D(19)) D(19)=D(16)
121C
122C*********     FIGURE 4.1.5.1-26 SUBROUTINE
123C
124      CALL FIG26( D(19),B(I+2),D(12) )
125C
126C*********     FIGURE 4.1.5.1-28-B
127C
128      CALL TLINEX(X128B,X228B,Y28B,4,11,B(I+2),A(184),D(24),
129     1            0,2,0,2,Q5128B,3,ROUTID)
130C
131C*********     IF XTO GREATER THAN 0.30, THEN CAPL= 2.0
132C
133      CAPL=1.20
134      IF(VT(66).LT.0.30) CAPL=2.00
135      D(22)=D(12)*(1.+CAPL*VT(65)+100.*VT(65)**4)*D(24)*(2.*A(2)/SREF)
136C
137C*********     COMPUTE CDO FROM PAGE 3-2
138C
139      D(20)=D(21)+D(22)
140 1030 CONTINUE
141      RETURN
142      END
143