1      SUBROUTINE BODOWG(ALPHA,XCBO4,RCREO2,BWO2,TRAT,IVBW,GO2PAV,NA)
2C
3C***  COMPUTES BODY VORTEX EFFECTS ON THE WINGS
4C
5      COMMON /CONSNT/CONST(4)
6      COMMON /BDATA/ BD(762)
7      COMMON /BODYI/ XNX,X(20),S(20),P(20),R(20),ZU(20),ZL(20),
8     1               BNOSE,BTAIL,BLN,BLA,DS
9      EQUIVALENCE (RAD,CONST(4)),(PI,CONST(1))
10      DIMENSION ROUTID(2),Q1313A(3),Q1313B(3),Q31314(3),Q31315(3)
11C
12C     ----BODY EFFECT ON WING
13C
14      REAL   IVBW
15      DIMENSION ALPHA(1),IVBW(1),GO2PAV(1)
16      DIMENSION X1313B(11),Y1313B(11),C1313B(6)
17      DIMENSION X31314(15),Y31314(15),C31314(6)
18      DIMENSION X31315(15),Y31315(15),C31315(6)
19C
20C     ----FIG 4.3.1.3-13A
21C
22      DIMENSION X1313A(16),Y1313A(16),C1313A(6)
23      DATA Q1313A/4H4.3.,4H1.3-,4H13A /,Q1313B/4H4.3.,4H1.3-,4H13B /,
24     1     Q31314/4H4.3.,4H1.3-,4H14  /,Q31315/4H4.3.,4H1.3-,4H15  /
25     2     ,ROUTID/4HBODO,4HWG  /
26      DATA I1313A /0/
27      DATA X1313A /6.8,7.2,7.5,8.,8.4,9.,9.6,10.4,11.2,12.1,13.4,15.,16.
28     1             ,17.1,18.,20./
29      DATA Y1313A /20.,19.,18.,17.,16.,15.,14.,13.,12.,11.,10.,9.,8.5,8.
30     1             ,7.7,7./
31C
32C     ----FIG 4.3.13-13B
33C
34      DATA I1313B /0/
35      DATA X1313B /0.,.5,1.,2.,2.5,3.3,4.,4.7,5.5,6.,7./
36      DATA Y1313B /.86,1.02,1.2,1.52,1.65,1.8,1.92,2.,2.09,2.14,2.23/
37C
38C     ----FIG 4.3.1.3-14
39C
40      DATA I31314 /0/
41      DATA X31314 /0.,.5,1.,1.5,2.,2.5,3.,3.5,4.,4.5,5.,5.5,6.,6.5,7./
42      DATA Y31314 /.5,.57,.62,.66,.69,.71,.72,.735,.75,.755,.76,4*.765/
43C
44C     ----4.3.1.3-15
45C
46      DATA I31315 /0/
47      DATA X31315 /0.,.5,1.,1.5,2.,2.5,3.,3.5,4.,4.5,5.,5.5,6.,6.5,7./
48      DATA Y31315 /.4,.5,.6,.7,.81,.91,1.02,1.13,1.24,1.35,1.47,1.6,1.74
49     1             ,1.88,2.02/
50C
51C***  GETMAX CROSS SECTIONAL AREA
52C
53      NX=XNX+.5
54      CALL GETMAX(X,S,NX,BD(2),BD(3),IMAX)
55      RCREO2=SQRT(BD(3)/PI)
56      XRT=XCBO4/RCREO2
57      DO 1010 I=1,NA
58         ARG=ABS(ALPHA(I))
59         IF(ARG.GE.6.)GO TO 1000
60         IVBW(I)=0.
61         GO2PAV(I)=0.
62         GO TO 1010
63C
64C     ----FIG. 4.3.1.3-13A
65C
66 1000    CALL TBFUNX(ARG,XOR,DYDX,16,X1313A,Y1313A,C1313A,I1313A,MI,N,
67     1               0,0,Q1313A,3,ROUTID)
68         XD=ARG*(XRT-XOR)/RAD
69         IF(XD.LE.0.)GO TO 1005
70C
71C     ----FIG. 4.3.1.3-13B
72C
73         CALL TBFUNX(XD,ZOR,DYDX,11,X1313B,Y1313B,C1313B,I1313B,MI,N,
74     1               0,0,Q1313B,3,ROUTID)
75C
76C     ----FIG  4.3.1.3-14
77C
78         CALL TBFUNX(XD,YOR,DYDX,15,X31314,Y31314,C31314,I31314,MI,N,
79     1               0,1,Q31314,3,ROUTID)
80C
81C     ----FIG  4.3.1.3-15
82C
83         CALL TBFUNX(XD,GO2PAV(I),DYDX,15,X31315,Y31315,C31315,I31315,
84     1               MI,N,0,1,Q31315,3,ROUTID)
85         Z0=ZOR*RCREO2
86         Y0=YOR*RCREO2
87         IVBW(I)=ALI(Z0,Y0,BWO2,RCREO2,TRAT)
88 1005    IF(XD.LE.0.)IVBW(I)=0.0
89 1010 CONTINUE
90      RETURN
91      END
92