1 SUBROUTINE SUBRYW 2C 3C THIS ROUTINE YIELDS THE SUBSONIC WING,WING-BODY 4C ROLLING DERIVATIVES,CLP CNP CNR AND 5C YAWING DERIVATIVE , CYP 6C 7 COMMON /OVERLY/ NLOG,NMACH,IM,NALPHA 8 COMMON /OPTION/ SR,CBARR,RUFF,BLREF 9 COMMON /CONSNT/ PI,DR,UNUSED,RAD 10 COMMON /FLGTCD/ FLC(93) 11 COMMON /SYNTSS/ SYNA(19) 12 COMMON /WINGD/ A(195),B(49) 13 COMMON /WINGI/ WINGIN(77) 14 COMMON /IWING/ PWING, WING(400) 15 COMMON /POWR/ DYN(213) 16 COMMON /WHAERO/ C(51) 17 COMMON /SBETA/ STB(135) 18 COMMON /FLOLOG/ FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 19 1 HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN, 20 2 SUPERS,SUBSON,TRANSN,HYPERS, 21 3 SYMFP,ASYFP,TRIMC,TRIM 22 LOGICAL FLTC,OPTI,BO,WGPL,WGSC,SYNT,HTPL,HTSC,VTPL,VTSC, 23 1 HEAD,PRPOWR,JETPOW,LOASRT,TVTPAN, 24 2 SUPERS,SUBSON,TRANSN,HYPERS, 25 3 SYMFP,ASYFP,TRIMC,TRIM 26 INTEGER Y71336,D11336,D21336,D31336,Y71337,D11337,D21337,D31337, 27 1 Y2220A,D1220A,D2220A,Y2220B,Y2220C,Y2220D,D1220D,D2220D,Y12224, 28 2 Y12312,Y1219B,Y1219C 29 INTEGER Y7133A 30 REAL LAMCO4,MACOE 31 REAL MACH,LAMB,KAY(20),K,LAMDA 32 DIMENSION ROUTID(2),ALMALO(20) 33C AJT DIMENSION CL(20),ALSCHD(20),Z(6),CLSUB(20),ALSCH2(20),CD(20), 34C AJT1 DCLPD(20),CLPWBS(20),CLTAN(20),CDCDO(20),CL2PAW(20) 35 DIMENSION CL(20),ALSCHD(20),Z(6),CLSUB(21),ALSCH2(21),CD(20), 36 1 DCLPD(20),CLPWBS(20),CLTAN(20),CDCDO(20),CL2PAW(20) 37 DIMENSION Y(6),D(6),E(6),G(6),H(6),O(6),P(6) 38 DIMENSION CLA(20),DCDDA(20),DCADA(20) 39 DIMENSION CLPTAN(20),DCLDA(20),CDOO(20) 40 DIMENSION CYPWBS(20),I2A219(4),QEXPD(5), 41 1 QIMPR1(4),QIMPR2(4),QIMPR3(4),Q71336(3),Q71337(3) 42 DIMENSION X11336(4),X21336(15),X31336(3),Y71336( 90), 43 1 D11336(30),D21336(30),D31336(30) 44 DIMENSION X1133A(2),X2133A(2),Y7133A(2) 45 DIMENSION X11337(4),X21337(9),X31337(3),Y71337( 54), 46 1 D11337(18),D21337(18),D31337(18) 47 DIMENSION X1220A( 13),X2220A( 17),Y2220A( 111),Q2220A(3) 48 DIMENSION D1220A( 97),D2220A(14) 49 DIMENSION X1220B( 13),X2220B( 10),Y2220B( 65),Q2220B(3) 50 DIMENSION X1220C( 13),X2220C( 10),Y2220C( 65),Q2220C(3) 51 DIMENSION X1220D( 13),X2220D( 18),Y2220D( 117),Q2220D(3) 52 DIMENSION D1220D( 97),D2220D(20) 53 DIMENSION X12224( 8),X22224( 10),Y12224( 40),Q12224(3) 54 DIMENSION X12312( 6),X22312( 8),Y12312( 24),Q12312(3) 55 DIMENSION X1219A( 10),Y1219A( 10),Q1219A(3) 56 DIMENSION X1219B( 8),X2219B( 2),Y1219B( 8),Q1219B(3) 57 DIMENSION X1219C( 3),X2219C( 2),Y1219C( 3),Q1219C(3) 58 DIMENSION CNPWBS(20) 59 DIMENSION CNRWBS(20) 60 EQUIVALENCE(DCMDCL,DYN(21)),(MACOE,A(16)) 61 EQUIVALENCE(CL(1),WING(21)),(THETA,WINGIN(11)) 62 EQUIVALENCE(ZW,SYNA(3)),(ZCG,SYNA(5)),(GEG,STB(122)), 63 1 (ALSCHD(1),FLC(23)),(SSPN,WINGIN(4)),(SSPNE,WINGIN(3)), 64 2 (LAMCO4,A(64)),(SW,A(4)),(TANCO4,A(68)), 65 3 (LAMDA,A(27)),(COSLO4,A(67)),(CD(1),WING(1)),(ASTRW,A(120)) 66 EQUIVALENCE(ALSD,DYN(44)),(CLACLO,DYN(45)),(CDO,DYN(94)), 67 1 (CLA(1),DYN(47)),(ZEE,DYN(67)),(CLPCLP,DYN(68)), 68 2 (CLPCL2,DYN(69)),(BAOK,DYN(70)),(BCLPCL,DYN(71)), 69 3 (DCLPD(1),DYN(72)),(CLPWBS(1),WING(281)),(CNPCLO,DYN(92)), 70 4 (BEE,DYN(93)),(CNPCLM,DYN(46)),(CNPTHE,DYN(95)), 71 5 (DCLDA(1),DYN(96)),(DCDDA(1),DYN(116)),(DCADA(1),DYN(136)), 72 6 (KAY(1),DYN(156)),(CNPWBS(1),WING(321)) 73 EQUIVALENCE(CLPG,DYN(176)),(DCYPG,DYN(177)),(TRANS,DYN(178)), 74 1 (CHANGE,DYN(179)),(CYPCLM,DYN(180)),(CYPWBS(1),WING(301)), 75 2 (TRADE,DYN(181)),(CNRCL2,DYN(182)),(CNRCDO,DYN(183)), 76 3 (CDOO(1),DYN(184)),(CNRWBS(1),WING(341)) 77 EQUIVALENCE (D11336(1),Y71336(1)),(D21336(1),Y71336(31)), 78 1 (D31336(1),Y71336( 61)) 79 EQUIVALENCE (D11337(1),Y71337(1)),(D21337(1),Y71337(19)), 80 1 (D31337(1),Y71337(37)) 81 EQUIVALENCE(D1220A(1),Y2220A(1)),(D2220A(1),Y2220A( 98)) 82 EQUIVALENCE(D1220D(1),Y2220D(1)),(D2220D(1),Y2220D( 98)) 83 DATA ROUTID /4HSUBR,4HYW / 84 DATA I2A219 /10,3*0/ 85 DATA IN/0/, IX/0/, IO/0/, IP/0/, IQ/0/, IR/0/, IS/0/, IT/0/ 86 DATA QEXPD /4HEXPE,4HRIME,4HNTAL,4H DAT,4HA / 87 DATA Q12224 /4H7.1.,4H2.2-,4H24 / 88 DATA Q12312 /4H7.1.,4H2.3-,4H12 / 89 DATA QIMPR1 /4HCLTA,4HN VS,4H ALS,4HCHD / 90 DATA QIMPR2 /4HCDCD,4HO VS,4H ALS,4HCHD / 91 DATA QIMPR3 /4HCL2P,4HAW V,4HS AL,4HSCHD/ 92 DATA Q71336 /4H7.1.,4H3.3-,4H6 / 93 DATA Q71337 /4H7.1.,4H3.3-,4H7 / 94C 95C ----FIGURE 7.1.2.2-20A 96C 97 DATA Q2220A /4H7.1.,4H2.2-,4H20A / 98 DATA X1220A / .100E 02, .900E 01, .800E 01, .700E 01, 99 1 .600E 01, .500E 01, .450E 01, .400E 01, .350E 01, .300E 01, 100 2 .250E 01, .200E 01, .150E 01/ 101 DATA X2220A / -.193E 02, -.989E 01, -.437E 01, .134E-01, 102 1 .104E 02, .143E 02, .203E 02, .251E 02, .300E 02, .355E 02, 103 2 .402E 02, .445E 02, .505E 02, .547E 02, .598E 02, .651E 02, 104 3 .702E 02/ 105C 106C ----FIGURE 7.1.2.2-20B 107C 108 DATA Q2220B/4H7.1.,4H2.2-,4H20B / 109 DATA X1220B/ .100E 02, .900E 01, .800E 01, .700E 01, 110 1 .600E 01, .500E 01, .450E 01, .400E 01, .350E 01, .300E 01, 111 2 .250E 01, .200E 01, .150E 01/ 112 DATA X2220B/ -.201E 02, -.105E 02, .230E 00, .994E 01, 113 1 .197E 02, .300E 02, .403E 02, .499E 02, .604E 02, .703E 02/ 114C 115C ----FIGURE 7.1.2.2-20C 116C 117 DATA Q2220C/4H7.1.,4H2.2-,4H20C / 118 DATA X1220C/ .100E 02, .900E 01, .800E 01, .700E 01, 119 1 .600E 01, .500E 01, .450E 01, .400E 01, .350E 01, .300E 01, 120 2 .250E 01, .200E 01, .150E 01/ 121 DATA X2220C/ -.193E 02, -.937E 01, .893E-01, .105E 02, 122 1 .208E 02, .304E 02, .405E 02, .505E 02, .607E 02, .704E 02/ 123C 124C ----FIGURE 7.1.2.2-20D 125C 126 DATA Q2220D/4H7.1.,4H2.2-,4H20D / 127 DATA X1220D/ .100E 02, .900E 01, .800E 01, .700E 01, 128 1 .600E 01, .500E 01, .450E 01, .400E 01, .350E 01, .300E 01, 129 2 .250E 01, .200E 01, .150E 01/ 130 DATA X2220D/ -.200E 02, -.965E 01, -.458E 01, .293E-01, 131 1 .556E 01, .995E 01, .148E 02, .203E 02, .245E 02, .301E 02, 132 2 .354E 02, .400E 02, .447E 02, .502E 02, .553E 02, .602E 02, 133 3 .646E 02, .698E 02/ 134C 135C ----FIGURE 7.1.2.2-24 136C 137 DATA X12224/ .700E 02, .600E 02, .500E 02, .400E 02, 138 1 .300E 02, .200E 02, .100E 02, .000E 00/ 139 DATA X22224/ .196E 01, .252E 01, .302E 01, .350E 01, 140 1 .402E 01, .505E 01, .599E 01, .700E 01, .804E 01, .100E 02/ 141C 142C ----FIGURE 7.1.2.3-12 143C 144 DATA X12312/ .100E 01, .800E 00, .600E 00, .400E 00, 145 1 .200E 00, .000E 00/ 146 DATA X22312/ .202E 01, .307E 01, .404E 01, .502E 01, 147 1 .603E 01, .808E 01, .100E 02, .120E 02/ 148C 149C ----FIGURE 7.1.2.1-9A 150C 151 DATA Q1219A/4H7.1.,4H2.1-,4H9A / 152 DATA X1219A/ .822E 00, .109E 01, .144E 01, .207E 01, 153 1 .257E 01, .316E 01, .398E 01, .462E 01, .531E 01, .600E 01/ 154 DATA Y1219A/ .452E 01, .381E 01, .297E 01, .215E 01, 155 1 .147E 01, .966E 00, .486E 00, .279E 00, .117E 00, .154E-06/ 156C 157C ----FIGURE 7.1.2.1-9B 158C 159 DATA Q1219B/4H7.1.,4H2.1-,4H9B / 160 DATA X1219B/ .700E 02, .650E 02, .600E 02, .450E 02, 161 1 .300E 02, .000E 00, -.300E 02, -.450E 02/ 162 DATA X2219B/ 0.0,4.5/ 163C 164C ----FIGURE 7.1.2.1-9C 165C 166 DATA Q1219C/4H7.1.,4H2.1-,4H9C / 167 DATA X1219C/ .000E 00, .500E 00, .100E 01/ 168 DATA X2219C/ .453E 01, -.963E 00/ 169 DATA X11336/ .600E 02, .500E 02, .400E 02, .000E 00/ 170 DATA X21336/ .101E 01, .126E 01, .151E 01, .178E 01, 171 1 .201E 01, .244E 01, .299E 01, .352E 01, .400E 01, .501E 01, 172 2 .599E 01, .698E 01, .799E 01, .902E 01, .994E 01/ 173 DATA X31336/0.0,0.2,0.4/ 174C 175C ----FIGURE 7.1.3.3-6D 176C 177 DATA X1133A /0.0,1.0/ 178 DATA X2133A /0.0,4.0/ 179 DATA X11337 / .000E 00, .400E 02, .500E 02, .600E 02/ 180 DATA X21337 / .103E 01, .144E 01, .197E 01, .297E 01, 181 1 .397E 01, .503E 01, .606E 01, .706E 01, .803E 01/ 182 DATA X31337 /0.0,0.2,0.4/ 183C 184C IN THE FOLLOWING DATA STATEMENTS 2 FLOATING POINT NUMBERS,WITH 185C 3 SIGNIFICANT DIGITS EACH,ARE PACKED IN 1 FIXED POINT WORD. 186C 187C THE FORMAT IS AS FOLLOWS- 188C 189C SIGN1,NS2,NSPS,NP1,NP2,ND1,ND2 190C 191C WHERE- 192C 193C SIGN1 IS THE ACTUAL SIGN OF THE 1ST NUMBER 194C NS2 REPRESENTS THE SIGN OF THE 2ND NUMBER.0 IF + , 1 IF - 195C NSPS REPRESENTS THE SIGNS OF THE POWERS OF 10 FOR BOTH NUMBERS 196C 0 IF ++ , 1 IF -- , 2 IF +- , 3 IF -+ 197C NP1 IS THE POWER OF 10 FOR THE 1ST NUMBER. LIMITED TO 1 DIGIT 198C NP2 IS THE POWER OF 10 FOR THE 2ND NUMBER. LIMITED TO 1 DIGIT 199C ND1 3 DIGITS OF THE 1ST NUMBER (DECIMAL POINT AT END) 200C ND2 3 DIGITS OF THE 2ND NUMBER (DECIMAL POINT AT END) 201C 202 DATA D11336/ 122392378, 122366359, 122352345, 122335325, 203 1 122320318, 122314315, 122313316, 122317157, 122167177, 204 2 122181187, 122204212, 122218226, 122234241, 122241242, 205 3 122240239, 133758974, 122116132, 122142160, 122176185, 206 4 122200209, 122216221, 122224222, 124225341, 133325655, 207 5 132872102, 122127147, 122160173, 122189201, 122206206, 208 6 122207205/ 209 DATA D21336/ 122762654, 122584539, 122499463, 122436427, 210 1 122418408, 122404404, 122402403, 122403381, 122349331, 211 2 122322312, 122300299, 122297295, 122290289, 122282283, 212 3 122279279, 122253253, 122253253, 122253256, 122258261, 213 4 122265270, 122273273, 122272272, 122268175, 122180185, 214 5 122191196, 122209218, 122222229, 122245253, 122255260, 215 6 122261261/ 216 DATA D31336/ 111160130, 122950630, 122572509, 122462437, 217 1 122416388, 122379371, 122370371, 121373150, 122450447, 218 2 122369347, 112112630, 122291284, 122281278, 122280281, 219 3 122280283, 122351293, 122272247, 122242230, 122225225, 220 4 122224223, 122222228, 122236242, 123259951, 132994111, 221 5 122125130, 122141159, 122169178, 122196204, 122206206, 222 6 122206205/ 223 DATA D11337/-1133563462,-1133383332,-1133310299,-1133293293, 224 1-1133293569,-1133501468,-1133451456,-1133456456,-1133456451, 225 2-1133670644,-1133608608,-1133619631,-1133631637,-1133642975, 226 3-1133953952,-1132979101,-1122102103,-1122103103/ 227 DATA D21337/-1133743503,-1133440371,-1133331297,-1133291291, 228 1-1133291759,-1133560480,-1133457451,-1133457457,-1133463469, 229 2-1133806669,-1133623606,-1133611623,-1133634640,-1132646106, 230 3-1133985949,-1133971988,-1122101101,-1122102102/ 231 DATA D31337/-1133978712,-1133518417,-1133364346,-1133340334, 232 1-1132340102,-1133712538,-1133500472,-1133465465,-1133471471, 233 2-1123107796,-1133664614,-1133626638,-1133644650,-1132656121, 234 3-1123104979,-1132995101,-1122101103,-1122104105/ 235 DATA D1220A/-1133347366,-1133375377,-1133379378,-1133376372, 236 1-1133361354,-1133338321,-1133298280,-1133248215,-1133177339, 237 2-1133357361,-1133365369,-1133368363,-1133359353,-1133339325, 238 3-1133308288,-1133267237,-1133211175,-1133326340,-1133345353, 239 4-1133355355,-1133353347,-1133339329,-1133316301,-1133281258, 240 5-1133230206,-1133170310,-1133323329,-1133334338,-1133338336, 241 6-1133332326,-1133319305,-1133290268,-1133251227,-1133198169, 242 7-1133289303,-1133307311,-1133318318,-1133316310,-1133305298, 243 8-1133287276,-1133256241,-1133214192,-1133160267,-1133283286, 244 9-1133289289,-1133292289,-1133286284,-1133280268,-1133261247, 245 A-1133230208,-1133184157,-1133254266,-1133273277,-1133281281, 246 B-1133279277,-1133270267,-1133259250,-1133233222,-1133198179, 247 C-1133151242,-1133250255,-1133261263,-1133263263,-1133262256, 248 D-1133253247,-1133241228,-1133214193,-1133175150,-1133227232, 249 E-1133238242,-1133245246,-1133245243,-1133239233,-1133230224, 250 F-1133212203,-1133184169,-1133140205,-1133211214,-1133216220, 251 G-1133222222,-1133222219,-1133218213,-1133209198,-1133192176, 252 H-1133159139,-1133184190,-1133191193,-1133194194,-1133195195, 253 I-1133194193,-1133193191,-1133184176,-1133161151,-1133131159, 254 J-1133162164,-1133164168,-1133169169/ 255 DATA D2220A/-1133169169,-1133165164,-1133161157,-1133156150, 256 1-1133137120,-1133133132,-1133133134,-1133133134,-1133133133, 257 2-1133133134,-1133133132,-1133129124,-1133120114, -330102000/ 258 DATA Y2220B/-1133464486,-1133500494,-1133484457,-1133419364, 259 1-1133298218,-1133445465,-1133479473,-1133463441,-1133407354, 260 2-1133291215,-1133428444,-1133450449,-1133437419,-1133388344, 261 3-1133287209,-1133407414,-1133422422,-1133415398,-1133372331, 262 4-1133276202,-1133375386,-1133391392,-1133385377,-1133353320, 263 5-1133270199,-1133340353,-1133358358,-1133353342,-1133325296, 264 6-1133255189,-1133325330,-1133334334,-1133331324,-1133310299, 265 7-1133247190,-1133300304,-1133312307,-1133305301,-1133288270, 266 8-1133237183,-1133278280,-1133282282,-1133281276,-1133268251, 267 9-1133221178,-1133245248,-1133251251,-1133250248,-1133241232, 268 A-1133208167,-1133221221,-1133222221,-1133220217,-1133215206, 269 B-1133192161,-1133183182,-1133182181,-1133180180,-1133178174, 270 C-1133162140,-1133138139,-1133141141,-1133141141,-1133142141, 271 D-1133134124/ 272 DATA Y2220C/-1133523535,-1133539534,-1133519496,-1133456407, 273 1-1133337245,-1133503509,-1133513505,-1133499475,-1133441390, 274 2-1133317234,-1133470483,-1133486483,-1133470454,-1133421375, 275 3-1133312228,-1133445451,-1133453451,-1133445428,-1133400355, 276 4-1133298223,-1133410415,-1133420414,-1133409392,-1133370338, 277 5-1133287213,-1133371376,-1133377375,-1133371358,-1133338314, 278 6-1133269204,-1133343349,-1133354353,-1133349339,-1133321299, 279 7-1133260204,-1133320320,-1133322322,-1133319316,-1133300280, 280 8-1133245199,-1133293297,-1133298297,-1133295286,-1133278264, 281 9-1133238189,-1133264264,-1133263264,-1133262259,-1133248239, 282 A-1133220180,-1133226226,-1133225226,-1133224222,-1133219207, 283 B-1133191165,-1133188188,-1133189187,-1133186186,-1133180172, 284 C-1133163145,-1133149149,-1133149149,-1133148147,-1133146142, 285 D-1133134121/ 286 DATA D1220D/-1133560573,-1133575575,-1133573567,-1133560546, 287 1-1133539520,-1133502481,-1133461430,-1133398367,-1133335280, 288 2-1133538545,-1133546546,-1133544542,-1133537525,-1133518501, 289 3-1133484464,-1133446415,-1133384353,-1133319269,-1133508517, 290 4-1133519517,-1133515512,-1133507499,-1133487475,-1133462442, 291 5-1133425396,-1133367339,-1133305266,-1133469477,-1133480480, 292 6-1133477471,-1133469463,-1133454448,-1133434418,-1133401376, 293 7-1133350321,-1133290246,-1133431434,-1133435435,-1133433431, 294 8-1133428423,-1133416409,-1133397383,-1133371350,-1133322299, 295 9-1133273232,-1133385390,-1133391390,-1133389386,-1133384378, 296 A-1133375368,-1133360353,-1133340325,-1133304284,-1133260224, 297 B-1133359361,-1133361362,-1133362360,-1133360356,-1133354347, 298 C-1133339334,-1133323307,-1133288271,-1133249216,-1133332334, 299 D-1133334334,-1133333333,-1133331329,-1133325322,-1133317310, 300 E-1133302290,-1133276257,-1133237208,-1133265266,-1133266266, 301 F-1133266265,-1133265264,-1133262261,-1133258255,-1133250246, 302 G-1133234222,-1133209183,-1133298300,-1133302302,-1133302301, 303 H-1133300298,-1133298293,-1133289284,-1133280267,-1133255243, 304 I-1133225197,-1133230230,-1133230230,-1133230231,-1133230230, 305 J-1133229228,-1133224223,-1133219216/ 306 DATA D2220D/-1133208198,-1133188170,-1133189190,-1133190190, 307 1-1133190189,-1133188187,-1133187184,-1133183181,-1133179175, 308 2-1133171163,-1133156143,-1133148146,-1133146145,-1133145145, 309 3-1133146146,-1133145145,-1133144142,-1133142141,-1133141139, 310 4-1133135122/ 311 DATA Y12224/-1133407331,-1133278240,-1133216175,-1133146127, 312 1-1134112891,-1133171137,-1133116103,-1144901734,-1144620532, 313 2-1144456365,-1144918736,-1144606535,-1144467381,-1144344280, 314 3-1144239183,-1144567428,-1144382349,-1144297232,-1144183177, 315 4-1144160125,-1144400303,-1144267223,-1144193165,-1144149120, 316 5-1144114103,-1144302225,-1144199187,-1144159131,-1145103856, 317 6-1155683804,-1144230201,-1144165141,-1145136843,-1155686624, 318 7-1155338463,-1144221179,-1144153141,-1145113850,-1155687398, 319 8-1155341456/ 320 DATA Y12312/-1166968781,-1166626494,-1166388197, -176319113, 321 1-1166941724,-1166541382, -168238339, 166218422,-1166885655, 322 2-1166450256, -176949195, 166454705,-1166777520,-1166300113, 323 3 176533337, 166590825,-1166548303, -177979781, 166231499, 324 4 166724919, -176387127, 166264371, 166473661, 166827994/ 325 DATA Y1219B/ 122468642 , 122368554 , 122300494 , 122176413, 326 1 122111372 , 132160292 , -132680236 ,-122168172/ 327 DATA Y1219C/ 1123153677, 1123161573 , 1123182380 / 328 DATA Y7133A/-133110100,-134127830/ 329 BETA=B(2) 330 IF( WGPL.AND.(.NOT.BO))GO TO 1000 331 DB=(SSPN-SSPNE)/SSPN 332 IF(DB.GT.0.3)RETURN 333 1000 K=BETA*WINGIN(IM+20)*RAD/(2.*PI) 334 MACH=FLC(IM+2) 335 XYZ=0. 336C 337C WING-BODY ROLLING DERIVATIVE,CLP--- 338C 339 CALL TBFUNX(XYZ,ALSD,DYDX,NALPHA,CL,ALSCHD,Y,IN,MI,NG,1,1,QEXPD,5, 340 1 ROUTID) 341 CALL TBFUNX(XYZ,CDO,DYDX,NALPHA,CL,CD,D,IO,MI,NG,2,2, 342 1 QEXPD,5,ROUTID) 343 IF(CL(1).LE.0..AND.CL(NALPHA).GE.0.) GO TO 1050 344 NNALPH=NALPHA+1 345 IF(CL(1).LT.0.) GO TO 1020 346 CLSUB(1)=0. 347 ALSCH2(1)=ALSD 348 DO 1010 J=1,NALPHA 349 CLSUB(J+1)=CL(J) 350 1010 ALSCH2(J+1)=ALSCHD(J) 351 GO TO 1040 352 1020 CLSUB(NNALPH)=0. 353 ALSCH2(NNALPH)=ALSD 354 DO 1030 J=1,NALPHA 355 CLSUB(J)=CL(J) 356 1030 ALSCH2(J)=ALSCHD(J) 357 1040 CALL TBFUNX(ALSD,ABCD,CLACLO,NNALPH,ALSCH2,CLSUB,E,IX,MI, 358 1 NG,1,1,QEXPD,5,ROUTID) 359 GO TO 1060 360 1050 CALL TBFUNX(ALSD,ABCD,CLACLO,NALPHA,ALSCHD,CL,Z,IX,MI, 361 1 NG,1,1,QEXPD,5,ROUTID) 362 1060 CONTINUE 363 XBOCB=-DCMDCL*CBARR/MACOE 364 DO 1070 I=1,NALPHA 365 CALL TBFUNX(ALSCHD(I),CDEF,CLA(I),NALPHA,ALSCHD,CL,G,IQ,MI,NG, 366 1 0,0,QEXPD,5,ROUTID) 367 1070 CONTINUE 368 ZEE=ZW-ZCG 369 SINGEG=SIN(GEG/RAD) 370 CLPCLP=1.-2.*ZEE*SINGEG/SSPN+3.*(ZEE/SSPN)**2*SINGEG**2 371 LAMB=RAD*ATAN(TANCO4/BETA) 372C 373C -- FIGURE 7.1.2.2-24 -- 374C 375 CALL TLIP2X(X12224,X22224,Y12224,8,10,LAMCO4,ASTRW,CLPCL2, 376 1 1,1,1,1,Q12224,3,ROUTID) 377 SRSW=SR/SW 378 AAAA=4.*SSPN**2/(RAD*SRSW*BLREF**2) 379 BAOK =BETA*ASTRW/K 380C 381C FIGURE 7.1.2.2-20 A-D 382C 383 CALL INTEP3(BAOK,LAMB,LAMDA,X1220A,X2220A,Y2220A,13,17,Q2220A, 384 1 X1220B,X2220B,Y2220B,13,10,Q2220B, 385 2 X1220C,X2220C,Y2220C,13,10,Q2220C, 386 3 1,1,1,1,1,1, 387 4 X1220D,X2220D,Y2220D,13,18,Q2220D,BCLPCL) 388 BBBB=AAAA*BCLPCL*K*CLPCLP/(BETA*CLACLO) 389 DO 1080 I=1,NALPHA 390 CLT=CL(I)*SRSW 391 DCLPD(I)=CLPCL2*CLT**2-CDO/8.*SRSW 392 CLPWBS(I)=BBBB*CLA(I)+AAAA*DCLPD(I) 393 1080 CONTINUE 394C 395C ---WING-BODY ROLLING DERIVATIVE,CNP--- 396C 397 DO 1090 I=1,NALPHA 398 1090 ALMALO(I) = ALSCHD(I)-ALSD 399 BEE=SQRT(1.-MACH**2*COSLO4**2) 400 CCCC=(ASTRW+4.*COSLO4)/(ASTRW*BEE+4.*COSLO4) 401 DDDD=(ASTRW*BEE+.5*(ASTRW*BEE+COSLO4)*TANCO4**2) 402 EEEE=ASTRW+.5*(ASTRW+COSLO4)*TANCO4**2 403C AJT IF NALPHA=20 IT WILL NOW BE 21 AFTER DO LOOP AND CL(I) WILL BE 404C AJT OUTSIDE ARRAY BOUND 405 406C AJT IF(CL(I).EQ.0.0)GO TO 1100 407 408C AJT PAGE 7.1.2.3-1 OF DATCOM SAYS NO ESTIMATE OF CNP VARIATION WITH CL 409C AJT SHOULD BE MADE IF CL AND CD DATA NOT AVAILABLE 410 IF(CL(1).EQ.0.0)GO TO 1100 411 CNPCLO=-(ASTRW+6.*(ASTRW+COSLO4)*(XBOCB*TANCO4/ASTRW+TANCO4**2/12. 412 1 ))/((ASTRW+4.*COSLO4)*6.) 413 CNPCLM =CCCC*DDDD*CNPCLO/EEEE 414 GO TO 1110 415 1100 CNPCLM=0. 416 1110 CONTINUE 417 DO 1120 I=1,NALPHA 418 CLTAN(I)=CL(I)*SRSW*TAN(ALMALO(I)/RAD) 419 CDCDO(I)=(CD(I)-CDO)*SRSW 420 CL2PAW(I)=(CL(I)*SRSW)**2/(ASTRW*PI) 421 1120 CONTINUE 422C 423C -- FIGURE 7.1.2.3-12 -- 424C 425 CALL TLIP2X(X12312,X22312,Y12312,6, 8,LAMDA,ASTRW,CNPTHE,2,1,2,1, 426 1 Q12312,3,ROUTID) 427 CNP=CNPTHE*THETA 428 FFFF=AAAA*BLREF/(2.*SSPN) 429 DO 1150 I=1,NALPHA 430 CLPTAN(I)=CLPWBS(I)*TAN(ALMALO(I)/RAD) 431 IF(CL(I).EQ.0.0)GO TO 1130 432 CALL TBFUNX(ALSCHD(I),ABC,DCLDA(I),NALPHA,ALSCHD,CLTAN,H,IR, 433 1 MI,NG,1,1,QIMPR1,4,ROUTID) 434 CALL TBFUNX(ALSCHD(I),ABC,DCDDA(I),NALPHA,ALSCHD,CDCDO,O,IS, 435 1 MI,NG,1,1,QIMPR2,4,ROUTID) 436 CALL TBFUNX(ALSCHD(I),ABC,DCADA(I),NALPHA,ALSCHD,CL2PAW,P,IT, 437 1 MI,NG,1,1,QIMPR3,4,ROUTID) 438 KAY(I)=(DCLDA(I)-DCDDA(I))/(DCLDA(I)-DCADA(I)) 439 GO TO 1140 440 1130 CONTINUE 441 KAY(I)=1. 442 1140 CONTINUE 443 CNPWBS(I)=(KAY(I)-1.)*CLPTAN(I)+(KAY(I)*(CNPCLM*CL(I)*SRSW) 444 1 +CNP)*AAAA 445 1150 CONTINUE 446C 447C WING-BODY YAWING DERIVATIVE,CYP--- 448C 449 CLPG=BCLPCL*K/BETA 450 DCYPG=(3.*SINGEG*(1.-2.*ZEE*SINGEG/SSPN))*CLPG 451C 452C -- FIGURE 7.1.2.1-9 -- 453C 454 CALL INTERX(1,X1219A,ASTRW,I2A219,Y1219A,TRANS,10,10,1,0,0,0,0,0 455 1 ,0,0,Q1219A,3,ROUTID) 456 CALL TLIP2X(X1219B,X2219B,Y1219B,8,2,LAMCO4,TRANS,CHANGE,1 457 1 ,1,1,1,Q1219B,3,ROUTID) 458 CALL TLIP2X(X1219C,X2219C,Y1219C,3,2,LAMDA,CHANGE,CYPCLO, 459 1 0,0,1,1,Q1219C,3,ROUTID) 460 AABB=ASTRW+4.*COSLO4 461 BBCC=ASTRW*BEE+4.*COSLO4 462 CCDD=ASTRW*BEE+COSLO4 463 DDEE=ASTRW+COSLO4 464 CYPCLM=AABB*CCDD*CYPCLO/(BBCC*DDEE) 465 FFPG=FFFF*DCYPG 466 DO 1160 I=1,NALPHA 467 CYPWBS(I)=FFPG+FFFF*CYPCLM*SRSW*KAY(I)*CL(I) 468 1160 CONTINUE 469C 470C WING-BODY ROLLING DERIVATIVE,CNR--- 471C --- FIGURE 7.1.3.3-7 -- 472C -- FIGURE 7.1.3.3-6 -- 473C 474 CALL TLIP3X(X11336,X21336,X31336,Y71336,4,15,3,LAMCO4,ASTRW,XBOCB, 475 1 TRADE,0,1,0,2,1,1,Q71336,3,ROUTID) 476 CALL TLIP2X(X1133A,X2133A,Y7133A,2,2,LAMDA,TRADE,CNRCL2,0,1,0,1, 477 1 Q71336,3,ROUTID) 478 CALL TLIP3X(X11337,X21337,X31337,Y71337,4, 9,3,ABS(LAMCO4), 479 1 ASTRW,XBOCB,CNRCDO,0,1,0,2,1,1,Q71337,3,ROUTID) 480 DO 1170 I=1,NALPHA 481 CDOO(I)=CD(I)*SRSW-CL2PAW(I) 482 CNRWBS(I)=AAAA*((CNRCL2*(CL(I)*SRSW)**2)+CNRCDO*CDOO(I)) 483 1170 CONTINUE 484 RETURN 485 END 486