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