C

C TILTING PROPROTOR AIRCRAFT DYNAMICS

C WAYNE JOHNSON
NASA-ARC/USAAMRDL
C FEBRUARY 1976
C
COMMON /TRCOEF/A2(1936),A1(1936),A0(1936),B0(748),MX,MX1,MV,MVG,DO1F1(44)
INTEGER DOF1
COMMON /RCOEF/A2R(30,30),A1R(30,30),A0R(30,30),AA2R(30,6),AA1R(30,16),AA0R(30,6),BR(30,8),BGR(30,3),C2R(6,30),C1R(6,30),C0R(6,30),CA22R(6,6),CA1R(6,6),CA0R(6,6),DGR(6,3)
COMMON /ACOEF/CAW(6,11),CBW(3,3),A2W(11,13),A1W(11,13),A0W(11,13),1BW(11,2),BGW(11,3),AAW(11,6)
COMMON /AINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,ATILT,OPUNIT,ALTMSL,TEM1P,ALTAGL,ACLIMB,APITCH,COLL,LNGCYC,WEIGHT,AFLAP,ITER,EPTRIM,FACTOR2,PTRIM,DOF(44),CON(14),GUST(3),OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,3OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,4JCNST,NSYSAN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(44),CONP(17),5THETFP,THETAF,CNTRL0, CNTRLS,VELX,VELZ,GRAV,DENSE,CSOUND,FREQ(100)
COMMON /BINPUT/ASHAFT,ACANT,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRT1R,KPRGOV,KIRGOV,KPEGOV,KIEGOV,RTRANS,KMAST,KICS,KENG,IENG,EDAMPC,T2HRTLC,NRPOWR,KRLGE,KPMC(8),KPMS(8),HMAST,FSPIV,WLPIV,BLPIV,FSCG,WL3CG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,FSPYL,WLPYL,IXX,IYY,IZZ,IXZ,MPYLO4N,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,KS1CFE,PSCFE,KT5CFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,PPCFE,KACFE,KRCFE,KPED6AL,CNTRLZ(6),COLLZ0,COLLZ2,K0TRIM,QMASSS(8),QFREQS(8),QDAMPS(8),ZE7TARS(3,8),GAMARS(3,8),QMASSA(8),QFREQA(8),QDAMPA(8),ZETARA(3,8),GA8MARA(3,8),AMASS,AIXX,AIYY,AIZZ,AIXZ,XRTR,YRTR,ZRTR,XWB,ZWB,XHT,ZHT9,XVT,ZVT,TMEAN,TCFES(6,3),TCFEA(6,2),RSF(3,3),RFV(3,3),RSV(3,3),ZAAGL,CSEGE,LAMDA0,TAUW,TAUH,TAUV,CWRTR,CHRTR,CVRTR,NEM,EDAMP,QTHRTL
INTEGER OPUNIT,DOF,CON,GUST,OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,OPI1NTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,RESPON,D2OFP,CONP
REAL LNGCYC,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRTR,KPRGOV,KIRGOV,1KPEGOV,KIEGOV,KMAST,KICS,KENG,IENG,NRPOWR,KRLGE,KPMC,KPMS,IXX,IYY,2IZZ, IXZ,MPYLON,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,K3S1CFE,KTCFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,KACFE,KRCFE,KP4EDAL,LAMDA0
COMMON /RINPUT/HHR(667)
EQUIVALENCE (RADIUS,HHR(14)),(VTIPN,HHR(79)),(IB,HHR(32)),(NBLADE,1HHR(2)),(NBM,HHR(9)),(NTM,HHR(10)),(OMEGA,HHR(13)),(GIMBAL,HHR(11)3),(KPB(1),HHR(50)),(KPG,HHR(35)),(TRC,HHR(21)),(TRS,HHR(22)),(MU,H4HR(6)),(MUY,HHR(88)),(MUZ,HHR(89)),(LAMDA,HHR(7)),(KHLMDA,HHR(655)5),(KFLMDA,HHR(656)),(GAMMA,HHR(12)),(SIGMA,HHR(15)),(OPAERR,HHR(166)),(OPMODE,HHR(26)),(ROTATE,HHR(663))
INTEGER GIMBAL,OPAERR,OPMODE,ROTATE
REAL MU,MUY,MUZ,LAMDA,KHLMDA,KFLMDA,KPB(5),KPG,IB
REAL B(616),BG(132),PHI(1800)
EQUIVALENCE (B0(1),B(1)),(B0(617),BG(1))
INTEGER DOF0(44)
REAL CKNOT(2),C2H(11,6),C1H(11,6),C0H(11,6),BGH(30,3),DGH(6,3)
EQUIVALENCE (B0(1),BGH(1,1)),(B0(91),DGH(1,1))
DATA CKNOT/.5924838,1.943845/
DATA CVERT/57.29578/

NAMELIST /NLCYC/NCYCLE
NAMELIST /NLTR/VTIP,RPM,VEL,VKTS,VELX,VELZ,WEIGHT,ATILT,ASHAFT,ACA1NT,ALTMSL,TEMP,ALTAGL,DENSE,CSOUND,ALTD,GRAV,ACLIMB,APITCH,THETFP,2THETAF,COLL,LNGCYC,AFLAP,DOF,CON,GUST,OPUNIT,OPSYMM,OPFLOW,OPTORS,3OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEB4UG,TONLY,MPER,JCNST,NSYSAN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP5,CONP,ITER,EPTRIM,FACTOR,PTRIM,FREQ,HMAST,FSPIV,WLPIV,BLPIV,FSCG,W6LCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,FSPYL,WLPYL,IXX,IYY,IZZ,IXZ,MPYL-

7ON,IXXPYL,IYYPYL,IZZPYL,IXZPYL,AMASS,AIXX,AIYY,AIZZ,AIXZ,XRTR,YRTR8,ZRTR,XWB,ZWB,XHT,ZHT,XVT,ZVT,K00CFE,K01CFE,KS0CFE,KS1CFE,PSCFE,KT9CFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,PPCFE,KACFE,KRCFE,KPEDAAL,CNTRLZ,COLLZ0,COLLZ2,K0TRIM,CNTRL0,CNTRLS,TMEAN,TCFES,TCFEA,RSFB,RFV,RSV,NEM,QMASSS,QFREQS,QDAMPS,ZETARS,GAMARS,QMASSA,QFREQA,QDAMCPA,ZETARA,GAMARA,KPRGOV,KIRGOV,KPEGOV,KIEGOV,RTRANS,KMAST,KICS,KENDG,IENG,EDAMPC,THRTLC,NRPOWR,KRLGE,KPMC,KPMS,EDAMP,QTHRTL,ZAGL,CSEGEE,LAMDA0,DLDZ,CNSTL,KRRH,KRRF,KRRL,KLTRIM,KWRTR,KHRTR,KVRTR,TAUW,TFAUH,TAUV,CWRTR,CHRTR,CVRTR
NAMELIST /NLEQ/A2,A1,A0,B0,MX,MX1,MV,MVG,MX0,DOF0,DOF1
NAMELIST /NLR/A2R,A1R,A0R,AA2R,AA1R,AA0R,BR,BGR,C2R,C1R,C0R,CA2R,1CA1R,CA0R,DGR
NAMELIST /NLS/CAW,CBW,A2W,A1W,A0W,BW,BGW,AAW
999 FORMAT (1H1/1X,'TILTING PROPROTOR AIRCRAFT DYNAMICS',10X,18A4//)
998 FORMAT (1X,'OPERATING STATE',69X,'ROTOR PARAMETERS'/
16X,'VELOCITY (KNOTS) =',F8.2,15X,'TILT ANGLE (DEG) =',F7.2,18X,'RA2DIUS (FT OR M) =',F7.2/6X,'V/(OMEGA*R)
=',F8.4,58X,'OMEGA (RA3D/SEC) =',F7.2/6X,'TIP SPEED (FPS OR MPS) =',F7.2,10X,'ALTITUDE M4SL (FT OR M) =',F7.0,12X,'IB (SLUG-FT2 OR KG-M2) =',F7.0)
997 FORMAT (6X,'ROTATIONAL SPEED (RPM) =',F7.2,10X,'TEMPERATURE (F OR 1C) =',F7.2,11X,'NUMBER OF BLADES =',I5)
996 FORMAT (6X,'ROTATIONAL SPEED (RPM) =',F7.2,10X,'STANDARD DAY',31X,1'NUMBER OF BLADES =',I5)
995 FORMAT (6X,'THETA-FP (DEG) =',F7.2,18X,'PRESSURE ALT (FT OR M) =' 1,F7.0,12X,'LOCK NUMBER
=',F7.3/6X,'CLIMB SPEED (FPS OR MPS) =2',F7.2, 8X,'AIR DENSITY RATIO =',F7.4,17X,'ROTOR SOLIDITY =',F7.34)
994 FORMAT (6X,'ALTITUDE AGL (FT OR M) =',F8.2, 9X,'SOUND SPEED (FPS O1R MPS) =',F8.2, 9X,'COUNTER-CLOCKWISE ROTATION'/6X,'CG STATION (FT2 OR M) =',F9.3/6X,'GROSS WEIGHT (LB OR KG) =',F8.0)
993 FORMAT (6X,'ALTITUDE AGL (FT OR M) =',F8.2, 9X,'SOUND SPEED (FPS O1R MPS) =',F8.2, 9X,'CLOCKWISE ROTATION'/6X,'CG STATION (FT OR M) =2',F9.3/6X,'GROSS WEIGHT (LB OR KG) =',F8.0)
992 FORMAT ( /1X,'OPTIONS'/
110X,'OPSYMM =',I2,7X,'IF 1, SYMMETRIC/ IF -1, ANTISYMMETRIC'/
210X,'OPFLOW =',I2,7X,'IF 0, AXIAL FLOW/ LT 0, CONSTANT COEFF/ GT 03, PERIODIC COEFF'/10X,'OPTRIM =',I2,7X,'IF 0, NO TRIM/ LT 0, NO IT4ERATIONS/ IF 1, TRIM T,H,MY/ IF 2, TRIM T,H,MY,P'/
510X,'OPMOTR =',I2,7X,'IF 0, AUTOROTATION/ IF 1, CONSTANT SPEED/ IF6 2, MOTOR/ IF 3, POWER OUT'/10X,'OPGOVR =',I2,7X,'IF 0, NO GOVERNO8R/ IF 1, GOVERNOR/ IF 2, INPUT GAINS'/10X,'OPGRND =',I2,7X,'IF 0, 9OUT OF GROUND EFFECT'/
A10X,'OPAERO =',I2,7X,'IF 0, NO PERTURBATION AIRCRAFT AERO FORCES'/B10X,'OPMODE =',I2,7X,'IF 0, RIGID BLADE MODES/ IF 2, READ MODES/ ICF 3, PUNCH MODES/ IF 4, ARTICULATED MODES'/10X,'OPAERO =',I2,7X,'IDF 0, ONLY CLA TERMS IN ROTOR AERO/ LT 0, NO STALL'/
E10X,'OPTORS =',I2,7X,'IF 0, DOF/ IF 1, RIGID/ IF 2, QUASISTATIC'/ F10X,'OPLMDA =',I2,7X,'IF 0, QUASISTATIC INFLOW DYNAMICS'/
G10X,'OPINTR =',I2,7X,'IF 0, QUASISTATIC AERO INTERFERENCE'/
H10X,'OPQRTR =',I2,7X,'IF 1, QUASISTATIC ROTOR/ IF 2, ALL QUASISTATIIC EXCEPT RIGID BODY')
991 FORMAT ( /1X,'ORDER OF VARIABLES'/10X,'DOF = B0,B1C,B1S,BN/2 T01,T1C,T1S,TN/2 BGC,BGS,PSI L0,LC,LS PSIE TF,XF,ZF ELASTIC2 LW,LH'/10X,'CON = T0,T1C,T1S,TN/2 TT DF,DE D0,D1S,DT'/10X3,'GUST = UG,VG,WG')
990 FORMAT ( /1X,'ORDER OF VARIABLES'/10X,'DOF = B0,B1C,B1S,BN/2 T01,T1C,T1S,TN/2 BGC,BGS,PSI L0,LC,LS PSIE PHIF,PSIF,YF ELA2STIC LW,LV'/10X,'CON = T0,T1C,T1S,TN/2 TT DA,DR D1C,DP'/103X,'GUST = UG,VG,WG')
989 FORMAT ( /10X,'DOF =',1X,15I1,3X,9I1,3X,3I1,3X,3I1,3X,I1,3X,3I1,3X1,8I1,3X,2I1/10X,'CON =',1X,8I1,3X,I1,3X,2I1,3X,3I1/10X,'GUST =',1X2,3I1)
988 FORMAT (//1X,'COMPUTATION PARAMETERS',15X,'CSYSAN - STEPS =',I4,5X1,'ROOT LOCUS =',I3,7X,'FREQ RESPONSE =',I3/38X,'PSYSAN - STEPS =',2I4,5X,'PLOT MAX =',F7.2,5X,'NUMBER OF FREQ =',I4/38X,'DOF PLOT = '-

3,44I1)
987 FORMAT (38X,'CON AND GUST PLOT = ',17I1)
986 FORMAT (1X,'ENGLISH UNITS (FT, SLUG, SEC)')
985 FORMAT (1X,'METRIC UNITS (M, KG, SEC)')
READ (5,NLCYC)
DO 1 JCYCLE=1,NCYCLE
VEL=0.
VKTS=0.
VTIP=0.
RPM=0.
READ (5,NLTR)
IF ((OPUNIT.NE.1).AND.(OPUNIT.NE.2)) OPUNIT=1
IF ((RPM .NE. 0.) .AND. (VTIP .EQ. 0.)) VTIP=RADIUS*RPM*.1047198
IF (VTIP .EQ. 0.) VTIP=VTIPN
IF (RPM .EQ. 0.) RPM=VTIP/RADIUS/.1047198
IF (VEL .EQ. 0.) VEL=VKTS/VTIP/CKNOT(OPUNIT)
IF (VKTS .EQ. 0.) VKTS=VEL*VTIP*CKNOT(OPUNIT)
CALL AIR(ALTMSL,TEMP,DENSE,CSOUND,GRAV,ALTD,OPUNIT)
CALL TRII(T75,T1C,T1S,1)
CALL RTRII(OPREAD,DEBUG,VTIP,T75,T1C,T1S,DENSE,CSOUND,OPUNIT)
IF ((OPMOTR .NE. 0) .AND. (OPMOTR .NE. 3)) GO TO 120
PTRIM=0.
IF (OPTRIM .GT. 0) OPTRIM=2
120 IF ((OPSYMM .LT. 0) .OR. (OPMOTR .EQ. 1)) OPGOVR=0
DO 20 K=1,44
IF (DOF(K) .NE. 0) DOF(K)=1
IF (NBLADE .NE. 2) GO TO 123
IF ((K .GE. 11) .AND. (K .LE. 15)) DOF(K)=0
IF ((K .GE. 22) .AND. (K .LE. 24)) DOF(K)=0
123 IF (K .GT. 14) GO TO 20
IF (CON(K) .NE. 0) CON(K)=1
IF ((K .GT. NBLADE) .AND. (K .LE. 8)) CON(K)=0
IF (K .GT. 3) GO TO 20
IF (GUST(K) .NE. 0) GUST(K)=1
20 CONTINUE
IF (NBLADE .EQ. 2) DOF(26)=0
IF (GIMBAL .EQ. 0) DOF(25)=0
IF (GIMBAL .EQ. 0) DOF(26)=0
IF (OPMOTR .EQ. 1) DOF(27)=0
IF (OPMOTR .LE. 1) DOF(31)=0
IF (OPSYMM .LT. 0) CON(14)=0
IF (OPSYMM .GT. 0) GUST(2)=0
IF (OPGRND .NE. 0) DOF(28)=1
IF (OPTORS .EQ. 0) GO TO 124
DO 125 K=16,24
IF (OPTORS .EQ. 1) DOF(K)=0
IF ((OPTORS .EQ. 2) .AND. (K .LE. 15+NBLADE)) DOF(K)=1
125 CONTINUE
124 NBM=0
DO 22 K=1,5
DO 22 KB=1,NBLADE
J=(K-1)*NBLADE+KB
IF (J .GT. 15) GO TO 22
IF (DOF(J) .NE. 0) NBM=K
22 CONTINUE
NBM=MAX0(1,NBM)
NTM=0
DO 23 K=1,3
DO 23 KB=1,NBLADE
J=(K-1)*NBLADE+KB
IF (J .GT. 9) GO TO 23
IF (DOF(J+15) .NE. 0) NTM=K
23 CONTINUE
NTM=MAX0(1,NTM)
NEM=0

DO 24 K=1,8
IF (DOF(K+34) .NE. 0) NEM=K
24 CONTINUE
NEM=MAX0(1,NEM)
MXP=0
MVP=0
DO 221 I=1,44
IF ((I.GE.16).AND.(I.LE.24).AND.(OPTORS.EQ.2)) GO TO 221
IF ((I.GE.28).AND.(I.LE.30).AND.(OPLMDA.EQ.0)) GO TO 221
IF ((I.GE.43).AND.(I.LE.44).AND.(OPINTR.EQ.0)) GO TO 221
IF ((I.LE.31).AND.(OPQRTR.NE.0)) GO TO 221
IF ((I.GE.35).AND.(OPQRTR.EQ.2)) GO TO 221
MXP=MXP+DOF(I)
221 CONTINUE
DO 222 I=1,14
MVP=MVP+CON(I)
IF (I .LE. 3) MVP=MVP+GUST(I)
222 CONTINUE
VCLIMB=VEL*VTIP*SIN(ACLIMB/CVERT)
PRINT 999, (TITLE(I), I=1,18)
PRINT 998, VKTS,ATILT,RADIUS,VEL,OMEGA,VTIP,ALTMSL,IB
IF (TEMP .LE. 999.) PRINT 997, RPM,TEMP,NBLADE
IF (TEMP .GE. 1000.) PRINT 996, RPM,NBLADE
PRINT 995, ACLIMB,ALTD,GAMMA,VCLIMB,DENSE,SIGMA
IF (ROTATE .EQ. 1) PRINT 994, ALTAGL,CSOUND,FSCG,WEIGHT
IF (ROTATE .EQ.-1) PRINT 993, ALTAGL,CSOUND,FSCG,WEIGHT
PRINT 992, OPSYMM,OPFLOW,OPTRIM,OPMOTR,OPGOVR,OPGRND,OPAERO,OPMODE1,OPAERR,OPTORS,OPLMDA,OPINTR,OPQRTR
IF (OPSYMM .GT. 0) PRINT 991
IF (OPSYMM .LT. 0) PRINT 990
PRINT 989, (DOF(I), I=1,44),(CON(I), I=1,14),(GUST(I), I=1,3)
PRINT 988, JCNST,LOCUS,RESPON,MPER,PLOTMX,NF,(DOFP(I), I=1,MXP)
PRINT 987, (CONP(I), I=1,MVP)
IF (OPUNIT .EQ. 1) PRINT 986
IF (OPUNIT .EQ. 2) PRINT 985
CALL TRII(T75,T1C,T1S,2)
CALL MODE
CALL TRTRIM
IF (TONLY .NE. 0) GO TO 101
IF (DEBUG .GE. 2) WRITE (6,NLTR)
GO TO 1
101 CONTINUE
PERIOD=6.283185/FLOAT(NBLADE)
IF (NBLADE/2*2 .EQ. NBLADE) PERIOD=2.*PERIOD
DT=PERIOD/FLOAT(MPER)
MT=0
IF (OPFLOW .GT. 0) MT=MPER
MT1=MT+1
NBLD=NBLADE
IF (OPFLOW .LT. 0) NBLD=JCNST
DO 2 JT=1,MT1
NT=JT-1
PSI=FLOAT(NT)*DT
IF ((NBLADE.GE.3).AND.(OPFLOW.EQ.0)) CALL ROTOR0
IF ((NBLADE.GE.3).AND.(OPFLOW.NE.0)) CALL ROTOR(OPFLOW,PSI,NBLD)
IF (NBLADE .EQ. 2) CALL ROTOR2(OPFLOW,PSI,NBLD)
CALL TRBODY
DO 131 I=1,3
DO 132 K=1,30
BGH(K,I)=0.
DO 132 J=1,3
132 BGH(K,I)=BGH(K,I)+BGR(K,J)*RSV(J,I)
DO 133 K=1,6
DGH(K,I)=0.
DO 133 J=1,3

133 DGH(K,I)=DGH(K,I)+DGR(K,J)*RSV(J,I)
131 CONTINUE
DO 134 I=1,3
DO 135 K=1,30
135 BGR(K,I)=BGH(K,I)
DO 134 K=1,6
134 DGR(K,I)=DGH(K,I)
DO 136 I=1,3
DO 136 K=1,11
BGH(K,I)=0.
DO 136 J=1,3
136 BGH(K,I)=BGH(K,I)+BGW(K,J)*RFV(J,I)
DO 137 I=1,3
DO 137 K=1,11
137 BGW(K,I)=BGH(K,I)
CNST1=1.+KRRL
IF (OPSYMM .LT. 0) CNST1=1.-KRRL
DO 231 I=1,30
A1R(28,I)=CNST1*A1R(28,I)
231 A0R(28,I)=CNST1*A0R(28,I)
DO 232 I=1,6
AA1R(28,I)=CNST1*AA1R(28,I)
AA0R(28,I)=CNST1*AA0R(28,I)
IF (I .GT. 3) GO TO 232
BGR(28,I)=CNST1*BGR(28,I)
232 CONTINUE
DO 31 I=1,1936
A2(I)=0.
A1(I)=0.
31 A0(I)=0.
DO 32 I=1,748
32 B0(I)=0.
DO 33 I=1,30
DO 33 K=1,30
JA=(I-1)*44+K
A2(JA)=A2R(K,I)
A1(JA)=A1R(K,I)
33 A0(JA)=A0R(K,I)
DO 34 I=1,8
DO 34 K=1,30
JB=(I-1)*44+K
34 B(JB)=BR(K,I)
DO 35 I=1,3
DO 35 K=1,30
JB=(I-1)*44+K
35 BG(JB)=BGR(K,I)
CNSTL=GAMMA/(SIGMA*5.7)*DENSE
A2(1216)=A2(1216)+.85*CNSTL
A2(1261)=A2(1261)+.22*CNSTL
A2(1306)=A2(1306)+.22*CNSTL
CNST1=SQRT((MU**2+MUY**2)/KFLMDA**2+LAMDA0**2/KHLMDA**4)
A1(1261)=A1(1261)+CNST1*CNSTL
A1(1306)=A1(1306)+CNST1*CNSTL
CNSTL=CNSTL*2.*(CNST1+(LAMDA0-MUZ)*LAMDA0/KHLMDA**4/CNST1)
A1(1216)=A1(1216)+CNSTL
IF (OPMOTR .LE. 1) GO TO 234
CIENG=IENG*RTRANS**2/(FLOAT(NBLADE)*IB)
CNST1=FLOAT(NBLADE)*IB*OMEGA**2
CICS=0.
IF (OPSYMM .LT. 0) CICS=KICS
CKMR=KMAST/(KMAST+CICS)
CKMI=CKMR*CICS/CNST1
CNST1=KENG/(KMAST+KENG+CICS)/CNST1
CKEM=CNST1*(KMAST+CISC)
CKEI=CNST1*CICS

A2(1171)=A2(1171)+CKMR*CIENG
A2(1347)=A2(1347)+CKMR*CIENG
A1(1171)=A1(1171)+CKMR*EDAMP
A1(1347)=A1(1347)+CKMR*EDAMP
A0(1171)=A0(1171)+CKMI
B(379)=B(379)+CKMR*QTHRTL
A2(1175)=A2(1175)+CIENG
A2(1351)=A2(1351)+CIENG
A1(1175)=A1(1175)+EDAMP
A1(1351)=A1(1351)+EDAMP
A0(1175)=A0(1175)+CKEI
A0(1351)=A0(1351)+CKEM
B(383)=B(383)+QTHRTL
234 CNSTG=KENG/(KMAST+KENG)
IEULER=1
IF (OPSYMM .LT. 0) IEULER=2
DO 41 I=1,11
DO 41 K=1,30
JA=(I+30)*44+K
DO 41 J=1,6
A2(JA)=A2(JA)+AA2R(K,J)*CAW(J,I)
A1(JA)=A1(JA)+AA1R(K,J)*CAW(J,I)
IF (I .LE. IEULER) GO TO 41
A0(JA)=A0(JA)+AA0R(K,J)*CAW(J,I)
41 CONTINUE
DO 42 I=1,30
DO 42 K=1,11
JA=(I-1)*44+K+31
DO 42 J=1,6
A2(JA)=A2(JA)-2.*AAW(K,J)*C2R(J,I)
A1(JA)=A1(JA)-2.*AAW(K,J)*C1R(J,I)
42 A0(JA)=A0(JA)-2.*AAW(K,J)*C0R(J,I)
DO 43 I=1,6
DO 43 K=1,11
C2H(K,I)=0.
C1H(K,I)=0.
C0H(K,I)=0.
DO 43 J=1,6
C1H(K,I)=C1H(K,I)+2.*AAW(K,J)*CA1R(J,I)
C0H(K,I)=C0H(K,I)+2.*AAW(K,J)*CA0R(J,I)
IF ((I .LE. 3) .OR. (J .LE. 3)) GO TO 43
C2H(K,I)=C2H(K,I)+2.*AAW(K,J)*CA2R(J,I)
43 CONTINUE
DO 44 I=1,13
DO 44 K=1,11
JA=(I+30)*44+K+31
A2(JA)=A2W(K,I)
A1(JA)=A1W(K,I)
A0(JA)=A0W(K,I)
IF (I .GT. 11) GO TO 44
DO 48 J=1,6
A2(JA)=A2(JA)-C2H(K,J)*CAW(J,I)
A1(JA)=A1(JA)-C1H(K,J)*CAW(J,I)
IF (I .LE. IEULER) GO TO 48
A0(JA)=A0(JA)-C0H(K,J)*CAW(J,I)
48 CONTINUE
44 CONTINUE
DO 45 I=1,IEULER
IF (OPSYMM .GT. 0) IC=2
IF (OPSYMM .LT. 0) IC=2*I-1
DO 45 K=1,30
JA=(I+30)*44+K
DO 45 J=1,3
A1(JA)=A1(JA)+AA2R(K,J)*CBW(J,IC)
IF (K .GT. 11) GO TO 45

A1(JA+31)=A1(JA+31)-C2H(K,J)*CBW(J,IC)
45 CONTINUE
DO 46 I=1,2
DO 46 K=1,11
JB=(I+8)*44+K+31
46 B(JB)=BW(K,I)
DO 47 I=1,3
DO 47 K=1,11
JB=(I-1)*44+K+31
BG(JB)=BGW(K,I)
DO 47 J=1,6
47 BG(JB)=BG(JB)+2.*AAW(K,J)*DGR(J,I)
A2(1891)=TAUW
A1(1891)=1.
A1(1231)=2.*CWRTR
A1(1936)=1.
IF (OPSYMM .LT. 0) GO TO 241
A2(1936)=TAUH
A1(1232)=2.*CHRTR
GO TO 242
241 A2(1936)=TAUV
A1(1232)=2.*CVRTR
242 IF (OPGRND .EQ. 0) GO TO 243
GLAMDA=(LAMDA0-MUZ)/(1.-(CSEGE*.25/ZAGL)**2)
IF (OPSYMM .LT. 0) GO TO 244
DLDZ=CNSTL*GLAMDA*CSEGE**2*.125/ZAGL**3
SNT=SIN(THETAF)
CST=COS(THETAF)
A0(1392)=A0(1392)-DLDZ*(ZRTR*SNT+XRTR*CST)
A0(1436)=A0(1436)-DLDZ*SNT
A0(1480)=A0(1480)+DLDZ*CST
DO 245 I=1,NEM
JA=(I+33)*44+28
245 A0(JA)=A0(JA)+DLDZ*(-ZETARS(1,I)*SNT+ZETARS(3,I)*CST)
GO TO 243
244 DLDZS=GLAMDA*CSEGE**2*.125/ZAGL**3
ZEFF=ZAGL/CSEGE
DLDZUS=
-GLAMDA*(((.1649*ZEFF-1.216)*ZEFF+2.726)*ZEFF-1.722)
DLDZ=CNSTL*(DLDZS+AMIN1(0.,KRLGE*(DLDZUS-DLDZS)))
SNT=SIN(THETAF)
CST=COS(THETAF)
A0(1392)=A0(1392)+DLDZ*YRTR*CST
DO 246 I=1,NEM
JA=(I+33)*44+28
246 A0(JA)=A0(JA)+DLDZ*(-ZETARA(1,I)*SNT+ZETARA(3,I)*CST)
243 CNSTC=B(61)*RADIUS
CNSTS=B(106)*RADIUS
IF (NBLADE .EQ. 2) GO TO 151
IF (OPTORS .EQ. 1) GO TO 153
DO 152 I=1,NEM
JA=(I+33)*44+17
A0(JA)=A0(JA)+CNSTC*KPMC(I)
152 A0(JA+1)=A0(JA+1)+CNSTS*KPMS(I)
GO TO 153
151 SNPSI=SIN(PSI)
CSPSI=COS(PSI)
DO 154 I=1,NEM
JA=(I+33)*44+17
154 A0(JA)=A0(JA)+CNSTC*KPMC(I)*CSPSI+CNSTC*KPMS(I)*SNPSI
153 IF (OPTORS .NE. 1) GO TO 50
JN=44*NBLADE
DO 52 J=1,JN
52 B(J)=-A0(J+660)
DO 53 I=1,NBM
DO 53 J=1,JN

JA=(I-1)*JN+J
IF (JA .GT. 660) GO TO 53
A0(JA)=A0(JA)-KPB(I)*A0(J+660)
53 CONTINUE
IF (GIMBAL .EQ. 0) GO TO 54
JN=88
IF (NBLADE .EQ. 2) JN=44
DO 55 J=1,JN
55 A0(J+1056)=A0(J+1056)-KPG*A0(J+704)
54 IF (NBLADE .EQ. 2) GO TO 57
DO 58 J=1,44
58 A0(J+1144)=A0(J+1144)+TRS*A0(J+704)-TRC*A0(J+748)
GO TO 56
57 TCYC=-TRC*SNPSI+TRS*CSPSI
DO 59 J=1,44
59 A0(J+1144)=A0(J+1144)+TCYC*A0(J+704)
56 IF (NBLADE .EQ. 2) GO TO 155
DO 156 I=1,NEM
DO 156 J=1,44
JA=(I+33)*44+J
156 A0(JA)=A0(JA)-(KPMC(I)*A0(J+704)+KPMS(I)*A0(J+748))*RADIUS
GO TO 50
155 DO 158 I=1,NEM
DO 158 J=1,44
JA=(I+33)*44+J
158 A0(JA)=A0(JA)-(KPMC(I)*CSPSI+KPMS(I)*SNPSI)*A0(J+704)*RADIUS
50 IF (OPGOVR .EQ. 0) GO TO 51
CNST1=1.
IF (OPGOVR .NE. 2) CNST1=AMAX1(0.,AMIN1(1.,ATILT/90.))
DO 251 I=1,44
CNSTB=KIRGOV*B(I)+KIEGOV*B(I+352)
A0(I+1144)=A0(I+1144)-CNSTB
A0(I+1320)=A0(I+1320)-CNSTB*CNSTG
CNSTB=(KPRGOV*B(I)+KPEGOV*B(I+352))*OMEGA*CNST1
A1(I+1144)=A1(I+1144)-CNSTB
251 A1(I+1320)=A1(I+1320)-CNSTB*CNSTG
51 IF (OPSYMM .LT. 0) GO TO 331
DO 332 I=1,3
DO 332 K=1,44
DO 332 J=1,6
JB=(I+10)*44+K
J1=(J-1)*44+K
IF (J .GE. 4) J1=J1+220
332 B(JB)=B(JB)+B(J1)*TCFES(J,I)
GO TO 333
331 DO 334 I=1,2
DO 334 K=1,44
DO 334 J=1,6
JB=(I+10)*44+K
J1=(J-1)*44+K
IF (J .GE. 4) J1=J1+220
334 B(JB)=B(JB)+B(J1)*TCFEA(J,I)
333 CONTINUE
MX=0
MX1=0
MX0=0
MV=0
MVG=0
DO 61 I=1,44
MX=MX+DOF(I)
DOF1(I)=1
61 DOF0(I)=1
IF ((OPFLOW .NE. 0) .OR. (OPSYMM .LT. 0)) GO TO 62
IF ((OPGOVR.NE.0) .AND. ((KIRGOV.NE.0.).OR.(KIEGOV.NE.0.))) GO TO 162

DOF1(27)=0
IF (DOF(27) .NE. 0) MX1=MX1+1
62 DO 63 I=28,30
DOF1(I)=0
IF (DOF(I) .EQ. 0) GO TO 63
MX1=MX1+1
IF (OPLMDA .NE. 0) GO TO 63
DOF0(I)=0
MX0=MX0+1
63 CONTINUE
DO 64 I=43,44
DOF1(I)=0
IF (DOF(I) .EQ. 0) GO TO 64
MX1=MX1+1
IF (OPINTR .NE. 0) GO TO 64
DOF0(I)=0
MX0=MX0+1
64 CONTINUE
IF (OPTORS .NE. 2) GO TO 65
DO 66 I=16,24
IF (DOF(I) .EQ. 0) GO TO 66
DOF0(I)=0
MX0=MX0+1
66 CONTINUE
65 IF (A0(1436) .EQ. 0.) DOF1(33)=0
IF (A0(1480) .EQ. 0.) DOF1(34)=0
IF ((DOF1(33).EQ.0).AND.(DOF(33).EQ.1)) MX1=MX1+1
IF ((DOF1(34).EQ.0).AND.(DOF(34).EQ.1)) MX1=MX1+1
IF (OPQRTR .EQ. 0) GO TO 161
DO 162 I=1,31
IF (DOF(I) .EQ. 0) GO TO 162
IF (DOF0(I) .EQ. 0) GO TO 162
IF ((OPSYMM .GT. 0) .AND. (I .EQ. 27)) GO TO 162
DOF0(I)=0
MX0=MX0+1
162 CONTINUE
IF (OPQRTR .EQ. 1) GO TO 161
DO 163 I=35,44
IF (DOF(I) .EQ. 0) GO TO 163
IF (DOF0(I) .EQ. 0) GO TO 163
DOF0(I)=0
MX0=MX0+1
163 CONTINUE
161 DO 67 I=1,14
67 MV=MV+CON(I)
DO 68 I=1,3
68 MVG=MVG+GUST(I)
CALL DERED(44,14,3,DOF,CON,GUST,A2,A1,A0,B,BG,DOF1,DOF0)
MVX=MX*MV
MVGX=MX*MVG
DO 69 I=1,MVGX
I1=I+MVX
69 B0(I1)=BG(I)
MV=MV+MVG
IF (MX0 .NE. 0) CALL QSTRAN(MX,MX0,MX1,MV,A2,A1,A0,B0,DOF1,DOF0)
IF (DEBUG .EQ. 0) GO TO 81
MXSQ1=MX*MX+1
IF (MXSQ1 .GT. 1936) GO TO 82
DO 83 I=MXSQ1,1936
A2(I)=0.
A1(I)=0.
83 A0(I)=0.
82 MVX1=MX*MV+1
IF (MVX1 .GT. 748) GO TO 81
DO 85 I=MVX1,748

85 B0(I)=0.
81 IF (DEBUG .GE. 2) WRITE(6,NLTR)
IF (DEBUG .EQ. 3) WRITE(6,NLR)
IF (DEBUG .EQ. 3) WRITE(6,NLS)
IF ((DEBUG.EQ.1).OR.(DEBUG.EQ.3)) WRITE(6,NLEQ)
IF (OPFLOW .GT. 0) GO TO 3
CALL CSYSAN(NSYSAN,MX,MX1,MV,A2,A1,A0,B0,NFREQ,FREQ,NSTEP,DOF1,OME1GA)
IF (LOCUS .GT. 0) CALL RPLOT(MX,MX1,JCYCLE,PLOTMX,LOCUS)
IF ((RESPON.GT.0).AND.(LOCUS.LE.0)) CALL TRANS(B0,MX,MX1,MV,RESPON1,NF,DOFP,CONP)
GO TO 2
3 NSYSAN=0
CALL PSYSAN(NSYSAN,MX,MX1,A2,A1,A0,PHI,DT,NT,MT,PERIOD,MT,DOF1)
2 CONTINUE
1 CONTINUE
IF ((LOCUS .GT. 0) .OR. (RESPON .GT. 0)) CALL ENPLT(10.,0.)
STOP
END
SUBROUTINE TRTRIM
COMMON /AINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,ATILT,OPUNIT,ALTMSL,TEM1P,ALTAGL,ACLIMB,APITCH,COLL,LNGCYC,WEIGHT,AFLAP,ITER,EPTRIM,FACTOR2,PTRIM,DOF(44),CON(14),GUST(3),OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,3OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,4JCNST,NSYSAN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(44),CONP(17),5THETFP,THETAF,CNTRL0, CNTRLS,VELX,VELZ,GRAV,DENSE,CSOUND,FREQ(100)
COMMON /BINPUT/ASHAFT,ACANT,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRT1R,KPRGOV,KIRGOV,KPEGOV,KIEGOV,RTRANS,KMAST,KICS,KENG,IENG,EDAMPC,T2HRTLC,NRPOWR,KRLGE,KPMC(8),KPMS(8),HMAST,FSPIV,WLPIV,BLPIV,FSCG,WL3CG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,FSPYL,WLPYL,IXX,IYY,IZZ,IXZ,MPYLO4N,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,KS1CFE,PSCFE,KT5CFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,PPCFE,KACFE,KRCFE,KPED6AL,CNTRLZ(6),COLLZ0,COLLZ2,K0TRIM,QMASSS(8),QFREQS(8),QDAMPS(8),ZE7TARS(3,8),GAMARS(3,8),QMASSA(8),QFREQA(8),QDAMPA(8),ZETARA(3,8),GA8MARA(3,8),AMASS,AIXX,AIYY,AIZZ,AIXZ,XRTR,YRTR,ZRTR,XWB,ZWB,XHT,ZHT9,XVT,ZVT,TMEAN,TCFES(6,3),TCFEA(6,2),RSF(3,3),RFV(3,3),RSV(3,3),ZAAGL,CSEGE,LAMDA0,TAUW,TAUH,TAUV,CWRTR,CHRTR,CVRTR,NEM,EDAMP,QTHRTL
INTEGER OPUNIT,DOF,CON,GUST,OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,OPI1NTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,RESPON,D2OFP,CONP
REAL LNGCYC,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRTR,KPRGOV,KIRGOV,1KPEGOV,KIEGOV,KMAST,KICS,KENG,IENG,NRPOWR,KRLGE,KPMC,KPMS,IXX,IYY,2IZZ, IXZ,MPYLON,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,K3S1CFE,KTCFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,KACFE,KRCFE,KP4EDAL,LAMDA0
COMMON /RINPUT/HHR(667)
EQUIVALENCE (THETA(1),HHR(294)),(T1C,HHR(21)),
(TWIST(11),HHR(396)),(LAMDA,HHR(7)),(MU,HHR(6)),(MUY,HHR(88)),(MUZ,HHR(89)2),(ALFHP,HHR(657)),(PSIHP,HHR(658)),(KFLMDA,HHR(656)),(KHLMDA,HHR(3655)),(VTOT,HHR(654)),(IRING,HHR(662)),(SIGMA,HHR(15)),(GAMMA,HHR(412)),(ROTATE,HHR(663)),(RADIUS,HHR(14)),(IB,HHR(32)),(NBLADE,HHR(25)),(OMEGA,HHR(13)),(T1S,HHR(22))
INTEGER ROTATE
REAL MU,MUY,MUZ,LAMDA,THETA(51),TWIST(51),KHLMDA,KFLMDA,IB
COMMON /FADATA/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,DRG10W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB,S2PANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,LF3TEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,HHFA(9)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
COMMON /XTRIM/LW,DW,MYW,YW,MXW,MZW,LH,DH,LV,DV,LAW,LFW,LEW,DAW,DFW1,DEW,MAW,MFW,NXBW,NZBW,YBW,LAH,LEH,DAH,DEH,LAV,LEV,DAV,DEV
REAL LW,MYW,MXW,MZW,LH,LV,LAW,LFW,LEW,MAW,MFW,NXBW,NZBW,LAH,LEH,LA-

1V,LEV
COMMON /ZTRIM/HHT1(27),CT,CH,CY,CMX,CMY,CQ,HHT2(2342)
REAL ZERO(33)
EQUIVALENCE (HHT1(1),ZERO(1))
REAL MOMY,FORCE(4),DERIV(4,4),LWAKE,CPOWR(2)
INTEGER LMINV(5),MMINV(5)
EQUIVALENCE (FORX,FORCE(1)),(FORZ,FORCE(2)),(MOMY,FORCE(3)),(POWR,1FORCE(4))
DATA CPOWR/550.,746./
DATA CVERT/57.29578/
999 FORMAT (1H1,'PERFORMANCE AND TRIM'//)
997 FORMAT (//5X,'NO FORCE TRIM')
996 FORMAT (//5X,'NO TRIM')
995 FORMAT (//5X,'FORCE TRIM ITERATIONS =',I4,7X,'(MAXIMUM =',I4,', T1OLERANCE =',F7.4,')')
994 FORMAT (5X,'WEIGHT CW/S =',F10.6,10X,'TARGET POWER CP/S =',F10.61)
993 FORMAT (5X,'WEIGHT CW/S =',F10.6)
992 FORMAT (/ 5X,'AIRCRAFT FORCES'/10X,'CFZ/S =',F10.6,10X,'WB LIFT ='1,F10.0, 5X,'HT LIFT =',F10.0, 5X,'VT LIFT =',F10.0/10X,'CFX/S =',F210.6,10X,'WB DRAG =',F10.0, 5X,'HT DRAG =',F10.0, 5X,'VT DRAG =',F310.0/10X,'CMY/S =',F10.6,10X,'WB MOMY =',F10.0/10X,'CP/S =',F10.64)
998 FORMAT (5X,'VEL =',F8.4,11X,'ALPHA =',F11.2,5X,'COLL. STICK =',F81.2, 7X,'(INPUT COLL =',F8.2,')'/5X,'VELX =',F8.4,10X,'ALPHA-HT =2',F8.2, 5X,'LONG. CYCLIC =',F8.2, 7X,'(INPUT LNGCYC =',F8.2,')'/5X3,'VELZ =',F8.4,33X,'THETA-F =',F12.2,7X,'(INPUT APITCH =',F8.2,')4'/5X,'VELC =',F8.4,10X,'DEL-F =',F8.2,8X,'THETA-FP =',F12.2,7X,'(I5NPUT ACLIMB =',F8.2,')'/5X,'VCLIMB =',F8.2, 8X,'DEL-E =',F8.2, 8X,6'K0TRIM =',I2,19X,'(K=1, TRIM CNTRL0/ K=0, TRIM CNTRLZ-1)' /5X,'P7CLIMB =',F9.2,30X,'AFLAP =',F8.2/5X,'POWER (HP) =',F9.2)
988 FORMAT (6X,I2,8X,5F9.4,8X,4F11.6)
989 FORMAT (1H1/20X,'C0',7X,'CZ',7X,'CS',7X,'TF',6X,'TFP',16X,'FX',9X,1'FZ',9X,'MY',10X,'P'/)
DO 10 J=1,33
10 ZERO(J)=0.
IF (TONLY .EQ. 0) PRINT 989
GM=AMASS*GRAV/(FLOAT(NBLADE)*IB*OMEGA**2/RADIUS)*5.7/(GAMMA*DENSE)
V2=VEL*VEL*.5/(SIGMA*3.141593*RADIUS**2)
DELF=TCFES(5,1)*CNTRL0+TCFES(5,2)*CNSTLS+CNTRLZ(5)/CVERT
ALPHA=APITCH-ACLIMB+IWB
IF (ABS(ALPHA) .GT. AMAXW) GO TO 121
LW=(LFTAW*ALPHA+LFTDW*AFLAP)/CVERT+LFTFW*(LFLAPO-LFLAPI)*DELF
DW=DRG0W+LW*LW/PIEL2+DRGDW*AFLAP/CVERT+DRGFW*(LFLAPO-LFLAPI)*DELF
GO TO 122
121 ALFS=SIGN(AMAXW,ALPHA)
LW=(LFTAW*ALFS+LFTDW*AFLAP)/CVERT+LFTFW*(LFLAPO-LFLAPI)*DELF
DW=(DRGSW-DRG0W)*(SIN(ALPHA/CVERT))**2+DRG0W
122 CTLMDA=.5*SQRT((GM*SIGMA-LW*.1519549*(VEL/RADIUS)**2)**2+(DW*.15911549*(VEL/RADIUS)**2)**2)
CTOS=1.0526*CTLMDA/SIGMA
MDERIV=3
IF (OPTRIM .EQ. 2) MDERIV=4
ITERM=ITER+(1+(ITER-1)/10)*MDERIV+1
ITERD=MDERIV+10
ITERC=-1
5 ITERC=ITERC+1
ITERCR=MOD(ITERC-1,ITERD)+1
IF (ITERC .EQ. 0) ITERCR=0
CST=COS(THETAF-THETFP)
SNT=SIN(THETAF-THETFP)
CSTF=COS(THETAF)
SNTF=SIN(THETAF)
RFV(1,1)=CST
RFV(2,1)=0.

RFV(3,1)=SNT
RFV(1,2)=0.
RFV(2,2)=1.
RFV(3,2)=0.
RFV(1,3)=-SNT
RFV(2,3)=0.
RFV(3,3)=CST
ONE=1.
DO 21 I=1,3
ONE=-ONE
RSV(I,1)=ONE*(RSF(I,1)*CST+RSF(I,3)*SNT)
RSV(I,2)=ONE*RSF(I,2)
21 RSV(I,3)=ONE*(-RSF(I,1)*SNT+RSF(I,3)*CST)
VELX=VEL*CST
VELZ=VEL*SNT
CSII=-RSV(1,1)
CSJI= RSV(2,1)
CSKI=-RSV(3,1)
VTOT=0.
MU=0.
MUY=0.
IF (OPFLOW .EQ. 0) GO TO 22
MU=-VEL*CSII
MUY=VEL*CSJI*FLOAT(ROTATE)
VTOT=MU*MU+MUY*MUY
IF (VTOT .GT. 0.) VTOT=SQRT(VTOT)
22 MUZ=VEL*CSKI
PSIHP=0.
ALFHP=0.
IF ((CSJI.NE.0.).OR.(CSII.NE.0.)) PSIHP=ATAN2(CSJI,-CSII)*CVERT
CSII=CSII**2+CSJI**2
IF (CSII .GT. 0.) CSII=SQRT(CSII)
IF ((CSKI.NE.0.).OR.(CSII.NE.0.)) ALFHP=ATAN2( CSKI,CSII)*CVERT
IF ((ITERCR.GT.MDERIV+2) .OR. ((ITERCR.LE.1).AND.(ITERC.GT.2))) CT1LMDA=CT*SIGMA
U2=(MU*MU+MUY*MUY)/KFLMDA**2
X1=.5*CTLMDA
X2=KHLMDA**4
LAMDA0=MUZ+X1/SQRT(ABS(X1/KHLMDA**2+VEL*VEL))
DO 23 I=1,4
X3=(SQRT(LAMDA0*LAMDA0/X2+U2))**3
23 LAMDA0=(X3*MUZ+X1*(2.*LAMDA0*LAMDA0/X2+U2))/(X3+X1*LAMDA0/X2)
IRING=0
IF (X1 .NE. 0.) X1=SQRT(ABS(X1))
IF (MU*MU+MUY*MUY+(2.*MUZ+3.*X1)**2 .GE. X1*X1) GO TO 24
IRING=1
LAMDA0=MUZ*((.373*MUZ*MUZ+.598*(MU*MU+MUY*MUY))/(X1*X1)-1.991)*KHL1MDA+MUZ
24 IF (VTOT .EQ. 0.) VTOT=.0001
IF (OPGRND .EQ. 0) GO TO 25
ZAGL=ALTAGL/RADIUS+XRTR*SNTF-ZRTR*CSTF
U2=MU*MU+MUY*MUY+LAMDA0*LAMDA0
IF (U2 .GT. 0.) U2=1./SQRT(U2)
CSEGE=(MU*(-RSF(1,1)*SNTF+RSF(1,3)*CSTF)-MUY*(-RSF(2,1)*SNTF+RSF(21,3)*CSTF)-LAMDA0*(-RSF(3,1)*SNTF+RSF(3,3)*CSTF))*U2
LAMDA0=MUZ+(1.-(CSEGE*.25/ZAGL)**2)*(LAMDA0-MUZ)
25 KRRL=AMIN1(AMAX1(KRRH,KRRF),AMAX1(AMIN1(KRRH,KRRF),(2.-20.*VTOT)*K1RRH+(20.*VTOT-1.)*KRRF))
LAMDA=LAMDA0+KRRL*(LAMDA0-MUZ)
CWRTR=KWRTR
CHRTR=KHRTR
CVRTR=0.
FH1=LAMDA*RSF(1,1)+MU*RSF(3,1)
FH3=LAMDA*RSF(1,3)+MU*RSF(3,3)
LWAKE=ABS(FH1*(XRTR-XWB)+FH3*(ZRTR-ZWB))-LAMDA

IF (LWAKE .GT. 0.) CWRTR=CWRTR*MU/(MU+LWAKE)

LWAKE=ABS(FH1*(XRTR-XHT)+FH3*(ZRTR-ZHT))-LAMDA
IF (LWAKE .GT. 0.) CHRTR=CHRTR*MU/(MU+LWAKE)
ALPHA=0.
U2=VELX**2+VELZ**2
IF (U2 .GT. 0.) U2=SQRT(U2)
X1=VEL*VELZ+U2*(LAMDA-MUZ)*CWRTR*RSF(3,3)*KLTRIM
X2=VEL*VELX+U2*(LAMDA-MUZ)*CWRTR*RSF(3,1)*KLTRIM
IF ((X1.NE.0.).OR.(X2.NE.0.)) ALPHA=ATAN2(X1,X2)
BETA=0.
ALFHT=0.
DAHT=0.
IF (AREAW .NE. 0.) DAHT=LW/AREAW*.45/((SPANW**2/AREAW)**.735*(ABS(1XHT*RADIUS/CHORDW))**.25)*KLTRIM
X1=VEL*(VELZ-DAHT*VELX)+U2*(LAMDA-MUZ)*CHRTR*RSF(3,3)*KLTRIM
X2=VEL*(VELX+DAHT*VELZ)+U2*(LAMDA-MUZ)*CHRTR*RSF(3,1)*KLTRIM
IF ((X1.NE.0.).OR.(X2.NE.0.)) ALFHT=ATAN2(X1,X2)
ALFVT=0.
T75 =TCFES(1,1)*CNTRL0+TCFES(1,2)*CNTRLS+CNTRLZ(1)
T1C =TCFES(2,1)*CNTRL0+TCFES(2,2)*CNTRLS+CNTRLZ(2)/CVERT
T1S =TCFES(3,1)*CNTRL0+TCFES(3,2)*CNTRLS+CNTRLZ(3)/CVERT
DELF=TCFES(5,1)*CNTRL0+TCFES(5,2)*CNTRLS+CNTRLZ(5)/CVERT
DELE=TCFES(6,1)*CNTRL0+TCFES(6,2)*CNTRLS+CNTRLZ(6)/CVERT
DELA=0.
DELR=0.
DO 31 JR=1,51
31 THETA(JR)=T75+TWIST(JR)/CVERT
CALL RTRIM(ITERC+1,CTOS)
CALL ACHAR(ALPHA,BETA,DELF,DELA,ALFHT,DELE,ALFVT,DELR)
FH1=RSF(1,1)*CH+RSF(2,1)*CY+RSF(3,1)*CT
FH3=RSF(1,3)*CH+RSF(2,3)*CY+RSF(3,3)*CT
FORX=-GM*SNTF+2.*FH1+V2*(-CST*(DW+DH+DV)+SNT*(LW+LH))
FORZ=GM*CSTF+2.*FH3+V2*(-SNT*(DW+DH+DV)-CST*(LW+LH))
MOMY=2.*ZRTR*FH1-2.*XRTR*FH3+2.*(RSF(1,2)*CMX+RSF(2,2)*CMY-RSF(3,21)*CQ)+V2/RADIUS*MYW+V2*((-ZWB*CST+XWB*SNT)*DW+(ZWB*SNT+XWB*CST)*LW2)+V2*((-ZHT*CST+XHT*SNT)*DH+(ZHT*SNT+XHT*CST)*LH)+V2*((-ZVT*CST+XV3T*SNT)*DV)
POWR=0.
IF (OPTRIM .EQ. 2) POWR=2.*CQ*FLOAT(ROTATE)-PTRIM
IF (TONLY .EQ. 0) PRINT 988, ITERCR,CNTRL0,CNTRLZ(1),CNTRLS,THETAF1,THETFP,(FORCE(I), I=1,4)
IF ((OPTRIM .LE. 0) .OR. (ITERC .GE. ITERM)) GO TO 6
IF (ITERC .EQ. 0) GO TO 5
IF (((ITERCR.LE.MDERIV+2).AND.(ITERCR.GE.2)) .OR. (ITERC.EQ.1)) GO1 TO 149
IF (ITERC .EQ. MDERIV+3) GO TO 149
IF ((FORX**2+FORZ**2 .LE. (EPTRIM*GM)**2) .AND. (MOMY**2 .LE..00251*(EPTRIM*GM)**2) .AND. (ABS(POWR) .LE. AMAX1(ABS(LAMDA*GM),.001)*E2PTRIM)) GO TO 6
149 GO TO (141,142,143,144,145),ITERCR
140 DO 40 I=1,MDERIV
IF (K0TRIM .EQ. 1) CNTRL0=CNTRL0-DERIV(1,I)*FORCE(I)
IF (K0TRIM .EQ. 0) CNTRLZ(1)=CNTRLZ(1)-DERIV(1,I)*FORCE(I)
CNTRLS=CNTRLS-DERIV(2,I)*FORCE(I)
THETAF=THETAF-DERIV(3,I)*FORCE(I)
IF (OPTRIM .EQ. 2) THETFP=THETFP-DERIV(4,I)*FORCE(I)
40 CONTINUE
GO TO 5
141 DO 41 I=1,MDERIV
DO 41 J=1,MDERIV
41 DERIV(J,I)=FORCE(J)
IF (K0TRIM .EQ. 1) CNTRL0=CNTRL0-.01*FACTOR
IF (K0TRIM .EQ. 0) CNTRLZ(1)=CNTRLZ(1)-.01*FACTOR
GO TO 5
142 DO 42 J=1,MDERIV

42 DERIV(J,1)=DERIV(J,1)-FORCE(J)
IF (K0TRIM .EQ. 1) CNTRL0=CNTRL0+.01*FACTOR
IF (K0TRIM .EQ. 0) CNTRLZ(1)=CNTRLZ(1)+.01*FACTOR
CNTRLS=CNTRLS-.02*FACTOR
GO TO 5
143 DO 43 J=1,MDERIV
43 DERIV(J,2)=DERIV(J,2)-FORCE(J)
CNTRLS=CNTRLS+.02*FACTOR
THETAF=THETAF-.02*FACTOR
GO TO 5
144 DO 44 J=1,MDERIV
44 DERIV(J,3)=DERIV(J,3)-FORCE(J)
THETAF=THETAF+.02*FACTOR
IF (OPTRIM .EQ. 2) GO TO 46
DO 47 I=1,3
DERIV(I,4)=0.
47 DERIV(4,I)=0.
DERIV(4,4)=1.
GO TO 48
46 THETFP=THETFP-.02*FACTOR
GO TO 5
145 IF (OPTRIM .NE. 2) GO TO 140
DO 45 J=1,MDERIV
45 DERIV(J,4)=DERIV(J,4)-FORCE(J)
THETFP=THETFP+.02*FACTOR
48 DO 49 I=1,MDERIV
DELTA=.02*FACTOR**2
IF (I .EQ. 1) DELTA=.01*FACTOR**2
DO 49 J=1,MDERIV
49 DERIV(J,I)=DERIV(J,I)/DELTA
CALL MINV(DERIV,4,DEN,LMINV,MMINV)
GO TO 5
6 PRINT 999
ALPHA=ALPHA*CVERT+IWB
ALFHT=ALFHT*CVERT+IHT
DELF=DELF*CVERT
DELE=DELE*CVERT
IF (K0TRIM .EQ. 1) ANG1=CNTRL0*CVERT
IF (K0TRIM .EQ. 0) ANG1=CNTRLZ(1)*CVERT-TMEAN
ANG2=CNTRLS*CVERT
ANG3=THETAF*CVERT
ANG4=THETFP*CVERT
VELC=VEL*SIN(THETFP)
VCLIMB=VELC*VTIP
PCLIMB=VCLIMB*AMASS*GRAV/CPOWR(OPUNIT)
POWR=2.*CQ*FLOAT(ROTATE)
POWRHP=POWR*FLOAT(NBLADE)*IB*OMEGA**3*GAMMA*DENSE/CPOWR(OPUNIT)/5.17
QUE=.5*DENSE*VEL*VEL*GAMMA*OMEGA**2/RADIUS**3*IB*FLOAT(NBLADE)/(3.1141593*5.7*SIGMA)
FLW=LW*QUE
FDW=DW*QUE
FMYW=MYW*QUE
FLH=LH*QUE
FDH=DH*QUE
FLV=LV*QUE
FDV=DV*QUE
PRINT 998, VEL,ALPHA,ANG1,COLL,VELX,ALFHT,ANG2,LNGCYC,VELZ,ANG3,AP1ITCH,VELC,DELF,ANG4,ACLIMB,VCLIMB,DELE,K0TRIM,PCLIMB,AFLAP,POWRHP
IF (OPTRIM) 61,62,63
61 PRINT 997
GO TO 64
62 PRINT 996
GO TO 65
63 PRINT 995, ITERC,ITERM,EPTRIM

IF (OPTRIM .EQ. 2) PRINT 994, GM,PTRIM

IF (OPTRIM .EQ. 1) PRINT 993, GM
64 PRINT 992, FORZ,FLW,FLH,FLV,FORX,FDW,FDH,FDV,MOMY,FMYW,POWR
CALL RTRIMP
65 CONTINUE
U2=MU*MU+MUY*MUY+LAMDA*LAMDA
IF (U2 .GT. 0.) U2=1./SQRT(U2)
TAUW=SQRT((XWB-XRTR)**2+(ZWB-ZRTR)**2)*U2
TAUH=SQRT((XHT-XRTR)**2+(ZHT-ZRTR)**2)*U2
TAUV=1.
U2=U2*(MU*RSF(1,3)-MUY*RSF(2,3)-LAMDA*RSF(3,3))
CWRTR=CWRTR*U2
CHRTR=CHRTR*U2
EDAMP=0.
IF (OPTRIM .NE. 0) EDAMP=EDAMPC*GAMMA*DENSE*CQ*FLOAT(ROTATE)/5.7
IF (OPMOTR .EQ. 3) EDAMP=0.
X1=CPOWR(OPUNIT)*1.5*NRPOWR/(FLOAT(NBLADE)*IB*OMEGA**3)
X2=GAMMA*DENSE*FLOAT(ROTATE)*CQ/(5.7*X1)
X2=4.*X2*(1.-X2)
IF (X2 .NE. 0.) X2=SQRT(ABS(X2))
QTHRTL=THRTLC*X1/57.3*X2
RETURN
END
SUBROUTINE TRII(T75,T1C,T1S,N)
COMMON /AINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,ATILT,OPUNIT,ALTMSL,TEM1P,ALTAGL,ACLIMB,APITCH,COLL,LNGCYC,WEIGHT,AFLAP,ITER,EPTRIM,FACTOR2,PTRIM,DOF(44),CON(14),GUST(3),OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,3OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,4JCNST,NSYSAN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(44),CONP(17),5THETFP,THETAF,CNTRL0, CNTRLS,VELX,VELZ,GRAV,DENSE,CSOUND,FREQ(100)
COMMON /BINPUT/ASHAFT,ACANT,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRT1R,KPRGOV,KIRGOV,KPEGOV,KIEGOV,RTRANS,KMAST,KICS,KENG,IENG,EDAMPC,T2HRTLC,NRPOWR,KRLGE,KPMC(8),KPMS(8),HMAST,FSPIV,WLPIV,BLPIV,FSCG,WL3CG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,FSPYL,WLPYL,IXX,IYY,IZZ,IXZ,MPYLO4N,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,KS1CFE,PSCFE,KT5CFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,PPCFE,KACFE,KRCFE,KPED5AL,CNTRLZ(6),COLLZ0,COLLZ2,K0TRIM,QMASSS(8),QFREQS(8),QDAMPS(8),ZE6TARS(3,8),GAMARS(3,8),QMASSA(8),QFREQA(8),QDAMPA(8),ZETARA(3,8),GA7MARA(3,8),AMASS,AIXX,AIYY,AIZZ,AIXZ,XRTR,YRTR,ZRTR,XWB,ZWB,XHT,ZHT8,XVT,ZVT,TMEAN,TCFES(6,3),TCFEA(6,2),RSF(3,3),RFV(3,3),RSV(3,3),ZA9GL,CSEGE,LAMDA0,TAUW,TAUH,TAUV,CWRTR,CHRTR,CVRTR,NEM,EDAMP,QTHRTL
INTEGER OPUNIT,DOF,CON,GUST,OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,OPI1NTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,RESPON,D2OFP,CONP
REAL LNGCYC,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRTR,KPRGOV,KIRGOV,1KPEGOV,KIEGOV,KMAST,KICS,KENG,IENG,NRPOWR,KRLGE,KPMC,KPMS,IXX,IYY,2IZZ, IXZ,MPYLON,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,K3S1CFE,KTCFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,KACFE,KRCFE,KP4EDAL,LAMDA0
COMMON /RINPUT/HHR(667)
EQUIVALENCE (RADIUS,HHR(14)),(ROTATE,HHR(663)),(GAMMA,HHR(12)),(SI1GMA,HHR(15))
INTEGER ROTATE
COMMON /EINPUT/FCHARH(20),FCHARA(20),AAH(19),ELASTA(576,5),ELASTS(1576,5),QMSYM(8,5),QMASYM(8,5),QFSYM(8,5),QFASYM(8,5),QDSYM(8,5),QD1ASYM(8,5),ZSYM(3,8,5),ZASYM(3,8,5),GSYM(3,8,5),GASYM(3,8,5),KPMCA(28),KPMSA(8),KPMCS(8),KPMSS(8)
COMMON /FEDATA/EZZN(10,10),EXXN(10,10),EZXN(10,10),EXZN(10,10),EZN1(10),EXN(10),ETN(10),EZTN(10,10),EXTN(10,10),ETTN(10,10),EZDN(10),2EXDN(10),EZLN(10),EXLN(10),ETLN(10),FZHN(10),FXHN(10),FTHN(10),FYW3N(10),FFHN(10),FYVN(10),FFVN(10),FPVN(10),IDENT(10)
COMMON /FETEMP/EZZ(8,8),EXX(8,8),EZX(8,8),EXZ(8,8),EZ(8),EX(8),ET(18),EZT(8,8),EXT(8,8),ETT(8,8),EZD(8),EXD(8),EZL(8),EXL(8),ETL(8),F2ZH(8),FXH(8),FTH(8),FYW(8),FFH(8),FYV(8),FFV(8),FPV(8)

COMMON /FADATA/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,DRG10W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB,S2PANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,LF3TEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,AAD2(9)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
DIMENSION AAD1(39),EED(576),EEDN(860)
EQUIVALENCE (AAD1(1),LFTAW),(EED(1),EZZ(1,1)),(EEDN(1),EZZN(1,1))
REAL K0CFE,CGW(2)
NAMELIST /NLFUSE/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ,EXX,EZX,EXZ,4EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,FXH,FTH,FYW,FFH,FYV,F5FV,FPV,QMASSS,QMASSA,QFREQS,QFREQA,QDAMPS,QDAMPA,KPMC,KPMS,ZETARS,6GAMARS,ZETARA,GAMARA
DATA CVERT/57.29578/
DATA CGW/32.17405,1./
GO TO (1,2),N
1 THETAF=APITCH/CVERT
THETFP=ACLIMB/CVERT
CSA=COS(ATILT/CVERT)
SNA=SIN(ATILT/CVERT)
TMEAN=COLLZ0-COLLZ2*SNA**2
K0CFE=K00CFE+K01CFE*SNA
IF (K0TRIM .LE. 1) GO TO 110
IF (K0TRIM .EQ. 3) GO TO 111
K0TRIM=1
IF (K0CFE .EQ. 0.) K0TRIM=0
GO TO 110
111 K0TRIM=1
IF (ATILT .EQ. 0.) K0TRIM=0
110 IF (K0TRIM .EQ. 0) GO TO 112
CNTRL0=COLL/CVERT
CNTRLZ(1)=TMEAN/CVERT
GO TO 113
112 CNTRL0=0.
CNTRLZ(1)=(TMEAN+COLL)/CVERT
113 CNTRLS=LNGCYC/CVERT
DO 11 I=1,6
DO 12 K=1,3
12 TCFES(I,K)=0.
DO 11 K=1,2
11 TCFEA(I,K)=0.
TCFES(1,1)=K0CFE
TCFES(4,1)=KTCFE
TCFES(5,1)=KFCFE
TCFES(2,2)=-(KS0CFE+KS1CFE*SNA)*SIN(PSCFE/CVERT)
TCFES(3,2)=-(KS0CFE+KS1CFE*SNA)*COS(PSCFE/CVERT)
TCFES(6,2)=KECFE
TCFES(4,3)=1.
FPEDAL=(KP0CFE+KP1CFE*SNA)*AMAX1(.25,AMIN1(1.,2.125-.75*VKTS/40.*K1PEDAL))
TCFEA(1,1)=-(KC0CFE+KC1CFE*SNA)
TCFEA(5,1)=-KACFE
TCFEA(2,2)=FPEDAL*SIN(PPCFE/CVERT)
TCFEA(3,2)=FPEDAL*COS(PPCFE/CVERT)
TCFEA(6,2)=KRCFE
AMASS=WEIGHT/CGW(OPUNIT)
FSP=FSPYL-FSPIV
WLP=WLPYL-WLPIV
FSCGA=FSCG+MPYLON/AMASS*((SNA-1.)*FSP-CSA*WLP)
WLCGA=WLCG+MPYLON/AMASS*(CSA*FSP+(SNA-1.)*WLP)
FSP=(FSCG+FSCGA-2.*FSPIV)*(FSCG-FSCGA)
WLP=(WLCG+WLCGA-2.*WLPIV)*(WLCG-WLCGA)

AIXX=IXX+AMASS*WLP+CSA**2*(IZZPYL-IXXPYL)+2.*CSA*SNA*IXZPYL
AIZZ=IZZ+AMASS*FSP-CSA**2*(IZZPYL-IXXPYL)-2.*CSA*SNA*IXZPYL
AIYY=IYY+AMASS*(WLP+FSP)
AIXZ=IXZ+AMASS*((WLCG-WLPIV)*(FSCG-FSCGA)+(FSCGA-FSPIV)*(WLCG-WLCG1A))+CSA*SNA*(IZZPYL-IXXPYL)-2.*CSA**2*IXZPYL
CST=COS(ASHAFT/CVERT)
SNT=SIN(ASHAFT/CVERT)
CSP=COS(ACANT/CVERT)
SNP=SIN(ACANT/CVERT)
RSF(1,1)=-CSP*SNA*CST-CSA*SNT
RSF(2,1)=-CSP*SNP*(1.-CSA)*CST-SNP*SNA*SNT
RSF(3,1)=(CSP**2*CSA+SNP**2)*CST-CSP*SNA*SNT
RSF(1,2)=-SNP*SNA
RSF(2,2)=CSP**2+SNP**2*CSA
RSF(3,2)=-CSP*SNP*(1.-CSA)
RSF(1,3)=CSP*SNA*SNT-CSA*CST
RSF(2,3)=CSP*SNP*(1.-CSA)*SNT-SNP*SNA*CST
RSF(3,3)=-(CSP**2*CSA+SNP**2)*SNT-CSP*SNA*CST
XRTR=(FSCGA-(FSPIV-HMAST*RSF(3,1)))/RADIUS
YRTR=(BLPIV+HMAST*RSF(3,2))/RADIUS
ZRTR=(WLCGA-(WLPIV-HMAST*RSF(3,3)))/RADIUS
XWB=(FSCGA-FSWB)/RADIUS
XHT=(FSCGA-FSHT)/RADIUS
XVT=(FSCGA-FSVT)/RADIUS
ZWB=(WLCGA-WLWB)/RADIUS
ZHT=(WLCGA-WLHT)/RADIUS
ZVT=(WLCGA-WLVT)/RADIUS
T75=(TCFES(1,1)*CNTRL0+TCFES(1,2)*CNTRLS+CNTRLZ(1))*CVERT
T1C=TCFES(2,1)*CNTRL0*CVERT+TCFES(2,2)*LNGCYC+CNTRLZ(2)
T1S=TCFES(3,1)*CNTRL0*CVERT+TCFES(3,2)*LNGCYC+CNTRLZ(3)
RETURN
2 CNSTH=(SIN(AMAX1(0.,AMIN1(90.,ATILT))/CVERT))**3
CNSTA=1.-CNSTH
DO 31 I=1,20
31 AAD1(I)=FCHARH(I)*CNSTH+FCHARA(I)*CNSTA
DO 32 I=1,19
32 AAD1(I+20)=AAH(I)
IVT=0.
AAD2(1)=XWB
AAD2(2)=ZWB
AAD2(3)=XHT
AAD2(4)=ZHT
AAD2(5)=XVT
AAD2(6)=ZVT
AAD2(7)=RADIUS
AAD2(8)=GAMMA*DENSE/(SIGMA*5.7)
AAD2(9)=AFLAP
IF (OPSYMM .LT. 0) GO TO 51
IF ((ATILT .LE. 0.) .OR. (ATILT .GE. 90.)) GO TO 151
NA=IFIX(ATILT/30.)+2
CNST0=FLOAT(NA-1)-ATILT/30.
CNST1=1.-CNST0
DO 152 I=1,576
152 EED(I)=CNST0*ELASTS(I,NA)+CNST1*ELASTS(I,NA+1)
DO 61 I=1,NEM
QMASSS(I)=CNST0*QMSYM(I,NA)+CNST1*QMSYM(I,NA+1)
QFREQS(I)=CNST0*QFSYM(I,NA)+CNST1*QFSYM(I,NA+1)
QDAMPS(I)=CNST0*QDSYM(I,NA)+CNST1*QDSYM(I,NA+1)
DO 61 K=1,3
ZETARS(K,I)=CNST0*ZSYM(K,I,NA)+CNST1*ZSYM(K,I,NA+1)
61 GAMARS(K,I)=CNST0*GSYM(K,I,NA)+CNST1*GSYM(K,I,NA+1)
GO TO 52
151 NA=1
IF (ATILT .GE. 90.) NA=5
DO 153 I=1,576

153 EED(I)=ELASTS(I,NA)
DO 62 I=1,NEM
QMASSS(I)=QMSYM(I,NA)
QFREQS(I)=QFSYM(I,NA)
QDAMPS(I)=QDSYM(I,NA)
DO 62 K=1,3
ZETARS(K,I)=ZSYM(K,I,NA)
62 GAMARS(K,I)=GSYM(K,I,NA)
GO TO 52
51 IF ((ATILT .LE. 0.) .OR. (ATILT .GE. 90.)) GO TO 154
NA=IFIX(ATILT/30.)+2
CNST0=FLOAT(NA-1)-ATILT/30.
CNST1=1.-CNST0
DO 155 I=1,576
155 EED(I)=CNST0*ELASTA(I,NA)+CNST1*ELASTA(I,NA+1)
DO 63 I=1,NEM
QMASSA(I)=CNST0*QMASYM(I,NA)+CNST1*QMASYM(I,NA+1)
QFREQA(I)=CNST0*QFASYM(I,NA)+CNST1*QFASYM(I,NA+1)
QDAMPA(I)=CNST0*QDASYM(I,NA)+CNST1*QDASYM(I,NA+1)
DO 63 K=1,3
ZETARA(K,I)=CNST0*ZASYM(K,I,NA)+CNST1*ZASYM(K,I,NA+1)
63 GAMARA(K,I)=CNST0*GASYM(K,I,NA)+CNST1*GASYM(K,I,NA+1)
GO TO 52
154 NA=1
IF (ATILT .GE. 90.) NA=5
DO 156 I=1,576
156 EED(I)=ELASTA(I,NA)
DO 64 I=1,NEM
QMASSA(I)=QMASYM(I,NA)
QFREQA(I)=QFASYM(I,NA)
QDAMPA(I)=QDASYM(I,NA)
DO 64 K=1,3
ZETARA(K,I)=ZASYM(K,I,NA)
64 GAMARA(K,I)=GASYM(K,I,NA)
52 IF (OPSYMM .LT. 0) GO TO 71
DO 72 I=1,NEM
KPMC(I)=KPMCS(I)
72 KPMS(I)=KPMSS(I)
GO TO 73
71 DO 74 I=1,NEM
KPMC(I)=KPMCA(I)
74 KPMS(I)=KPMSA(I)
73 CONTINUE
IF (OPREAD .GE. 2) READ(5,NLFUSE)
DO 41 I=1,860
41 EEDN(I)=0.
DO 42 I=1,10
IDENT(I)=0
IF (I .GT. NEM) GO TO 42
IDENT(I)=1
IF (OPSYMM .LT. 0) IDENT(I)=-1
42 CONTINUE
DO 43 I=1,NEM
EZN(I)=EZ(I)
EXN(I)=EX(I)
ETN(I)=ET(I)*RADIUS
EZDN(I)=EZD(I)
EXDN(I)=EXD(I)
EZLN(I)=EZL(I)/RADIUS
EXLN(I)=EXL(I)/RADIUS
ETLN(I)=ETL(I)
FZHN(I)=FZH(I)
FXHN(I)=FXH(I)
FTHN(I)=FTH(I)*RADIUS
FYWN(I)=FYW(I)

FFHN(I)=FFH(I)*RADIUS
FYVN(I)=FYV(I)
FFVN(I)=FFV(I)*RADIUS
FPVN(I)=FPV(I)*RADIUS
DO 43 K=1,NEM
EZZN(K,I)=EZZ(K,I)
EXXN(K,I)=EXX(K,I)
EZXN(K,I)=EZX(K,I)
EXZN(K,I)=EXZ(K,I)
EZTN(K,I)=EZT(K,I)*RADIUS
EXTN(K,I)=EXT(K,I)*RADIUS
43 ETTN(K,I)=ETT(K,I)*RADIUS**2
IF (DEBUG .GE. 2) WRITE(6,NLFUSE)
RETURN
END
SUBROUTINE TRBODY
COMMON /ACOEF/CA(6,11),CB(3,3),A2(11,13),A1(11,13),A0(11,13),B(11,12),BG(11,3),AA(11,6)
COMMON /FAERO/CZDZ,CZDX,CZDT,CZDDZ,CXDZ,CXDX,CXDT,CTDZ,CTDX,CTDT,C1TDDZ,CZDQ(8),CZQ(8),CXDQ(8),CXQ(8),CTDQ(8),CTQ(8),CQDZ(8),CQDX(8),2CQDT(8),CQDQ(8,8),CQQ(8,8),CZF,CZE,CZLW,CZLH,CXF,CXE,CXLW,CXLH,CTF3,CTE,CTLW,CTLH,CQF(8),CQE(8),CQLW(8),CQLH(8),CFDF,CFDP,CFDY,CPDF,C4PDP,CPDY,CYDF,CYDP,CYDY,CFDQ(8),CFQ(8),CPDQ(8),CPQ(8),CYDQ(8),CYQ(58),CQDF(8),CQDP(8),CQDY(8),CQDQA(8,8),CQQA(8,8),CFA,CFR,CFLV,CPA,C6PR,CPLV,CYA,CYR,CYLV,CQA(8),CQR(8),CQLV(8),CFDZ,CFDX,CFLW,CPDZ,CPD7X,CPLW,CYDZ,CYDX,CYLW,CQDZA(8),CQDXA(8),CQLWA(8)
COMMON /RINPUT/HHR(667)
EQUIVALENCE (RADIUS,HHR(14)),(OMEGA,HHR(13)),(NBLADE,HHR(2)),(IB,H1HR(32))
REAL IB
COMMON /AINPUT/HHA(291)
COMMON /BINPUT/HHB(314)
EQUIVALENCE (THETAF,HHA(184)),(GRAV,HHA(189)),(VEL,HHA(21)),(VELX,1HHA(187)),(VELZ,HHA(188)),(OPSYMM,HHA(99)),(DEBUG,HHA(111)),(OPAER2O,HHA(110))
EQUIVALENCE (AMASS,HHB(231)),(AIXX,HHB(232)),(AIYY,HHB(233)),(AIZZ1,HHB(234)),(AIXZ,HHB(235)),(XRTR,HHB(236)),(YRTR,HHB(237)),(ZRTR,H2HB(238)),(RSF(1,1),HHB(276)),(ZETARS(1,1),HHB(111)),(GAMARS(1,1),H3HB(135)),(QMASSS(1),HHB(87)),(QFREQS(1),HHB(95)),(QDAMPS(1),HHB(1043)),(ZETARA(1,1),HHB(183)),(GAMARA(1,1),HHB(207)),(QMASSA(1),HHB(1559)),(QFREQA(1),HHB(167)),(QDAMPA(1),HHB(175)),(NEM,HHB(312))
INTEGER OPSYMM,OPAERO,DEBUG
REAL RSF(3,3),ZETARS(3,8),GAMARS(3,8),QMASSS(8),QFREQS(8),QDAMPS(81),ZETARA(3,8),GAMARA(3,8),QMASSA(8),QFREQA(8),QDAMPA(8)
REAL ZERO(625),MSTAR
EQUIVALENCE (ZERO(1),CA(1,1))
DO 10 I=1,625
10 ZERO(I)=0.
PHIF=0.
VELY=0.
SNT=SIN(THETAF)
CST=COS(THETAF)
CI=.5*FLOAT(NBLADE)*IB
CM=CI/RADIUS**2
MSTAR=AMASS/CM
GM=GRAV/(OMEGA**2*RADIUS)*MSTAR
IF (OPSYMM .LT. 0) GO TO 11
A2(1,1)=AIYY/CI
A2(2,2)=MSTAR
A2(3,3)=MSTAR
A1(2,1)=MSTAR*VELZ
A1(3,1)=-MSTAR*VELX
A0(2,1)= GM*CST
A0(3,1)= GM*SNT
DO 21 IE=1,NEM

QM=QMASSS(IE)/CM
QW=QFREQS(IE)*6.283185/OMEGA
A2(IE+3,IE+3)=QM
A1(IE+3,IE+3)=QM*QDAMPS(IE)*QW
21 A0(IE+3,IE+3)=QM*QW*QW
CA(1,1)=RSF(1,1)*ZRTR-RSF(1,3)*XRTR
CA(2,1)=RSF(2,1)*ZRTR-RSF(2,3)*XRTR
CA(3,1)=RSF(3,1)*ZRTR-RSF(3,3)*XRTR
CA(4,1)=RSF(1,2)
CA(5,1)=RSF(2,2)
CA(6,1)=RSF(3,2)
CA(1,2)=RSF(1,1)
CA(2,2)=RSF(2,1)
CA(3,2)=RSF(3,1)
CA(1,3)=RSF(1,3)
CA(2,3)=RSF(2,3)
CA(3,3)=RSF(3,3)
DO 22 IE=1,NEM
DO 22 I=1,3
DO 22 J=1,3
CA(I,IE+3)=CA(I,IE+3)+RSF(I,J)*ZETARS(J,IE)
22 CA(I+3,IE+3)=CA(I+3,IE+3)+RSF(I,J)*GAMARS(J,IE)*RADIUS
IF (OPAERO .EQ. 0) GO TO 12
CALL HLAERO(VEL,VELX,VELY,VELZ,THETAF,PHIF,DEBUG)
A2(1,3)=-CTDDZ
A2(3,3)=A2(3,3)-CZDDZ
A1(1,1)=-CTDT
A1(1,2)=-CTDX
A1(1,3)=-CTDZ
A1(1,12)=-CTLW
A1(1,13)=-CTLH
A1(2,1)=A1(2,1)-CXDT
A1(2,2)=-CXDX
A1(2,3)=-CXDZ
A1(2,12)=-CXLW
A1(2,13)=-CXLH
A1(3,1)=A1(3,1)-CZDT
A1(3,2)=-CZDX
A1(3,3)=-CZDZ
A1(3,12)=-CZLW
A1(3,13)=-CZLH
B(1,1)=CTF
B(1,2)=CTE
B(2,1)=CXF
B(2,2)=CXE
B(3,1)=CZF
B(3,2)=CZE
BG(1,1)=CTDX
BG(1,3)=CTDZ
BG(2,1)=CXDX
BG(2,3)=CXDZ
BG(3,1)=CZDX
BG(3,3)=CZDZ
DO 23 IE=1,NEM
A1(1,IE+3)=-CTDQ(IE)
A1(2,IE+3)=-CXDQ(IE)
A1(3,IE+3)=-CZDQ(IE)
A1(IE+3,1)=-CQDT(IE)
A1(IE+3,2)=-CQDX(IE)
A1(IE+3,3)=-CQDZ(IE)
A1(IE+3,12)=-CQLW(IE)
A1(IE+3,13)=-CQLH(IE)
A0(1,IE+3)=-CTQ(IE)
A0(2,IE+3)=-CXQ(IE)
A0(3,IE+3)=-CZQ(IE)

B(IE+3,1)=CQF(IE)
B(IE+3,2)=CQE(IE)
BG(IE+3,1)=CQDX(IE)
BG(IE+3,3)=CQDZ(IE)
DO 23 JE=1,NEM
A1(IE+3,JE+3)=A1(IE+3,JE+3)-CQDQ(IE,JE)
23 A0(IE+3,JE+3)=A0(IE+3,JE+3)-CQQ(IE,JE)
GO TO 12
11 A2(1,1)=AIXX/CI
A2(1,2)=(-AIXX*SNT-AIXZ*CST)/CI
A2(2,1)=(-AIXX*SNT-AIXZ*CST)/CI
A2(2,2)=(AIXX*SNT*SNT+AIZZ*CST*CST+2.*AIXZ*SNT*CST)/CI
A2(3,3)=MSTAR
A1(3,1)=-MSTAR*VELZ
A1(3,2)=MSTAR*(VELZ*SNT+VELX*CST)
A0(3,1)=-GM*CST
DO 31 IE=1,NEM
QM=QMASSA(IE)/CM
QW=QFREQA(IE)*6.283185/OMEGA
A2(IE+3,IE+3)=QM
A1(IE+3,IE+3)=QM*QDAMPA(IE)*QW
31 A0(IE+3,IE+3)=QM*QW*QW
CA(1,1)=-RSF(1,2)*ZRTR+RSF(1,3)*YRTR
CA(2,1)=-RSF(2,2)*ZRTR+RSF(2,3)*YRTR
CA(3,1)=-RSF(3,2)*ZRTR+RSF(3,3)*YRTR
CA(4,1)=RSF(1,1)
CA(5,1)=RSF(2,1)
CA(6,1)=RSF(3,1)
CA(1,2)=-RSF(1,1)*YRTR*CST+RSF(1,2)*(ZRTR*SNT+XRTR*CST)-RSF(1,3)*Y1RTR*SNT
CA(2,2)=-RSF(2,1)*YRTR*CST+RSF(2,2)*(ZRTR*SNT+XRTR*CST)-RSF(2,3)*Y1RTR*SNT
CA(3,2)=-RSF(3,1)*YRTR*CST+RSF(3,2)*(ZRTR*SNT+XRTR*CST)-RSF(3,3)*Y1RTR*SNT
CA(4,1)=-RSF(1,1)*SNT+RSF(1,3)*CST
CA(5,1)=-RSF(2,1)*SNT+RSF(2,3)*CST
CA(6,1)=-RSF(3,1)*SNT+RSF(3,3)*CST
CA(1,3)=RSF(1,2)
CA(2,3)=RSF(2,2)
CA(3,3)=RSF(3,2)
DO 32 IE=1,NEM
DO 32 I=1,3
DO 32 J=1,3
CA(I,IE+3)=CA(I,IE+3)+RSF(I,J)*ZETARA(J,IE)
32 CA(I+3,IE+3)=CA(I+3,IE+3)+RSF(I,J)*GAMARA(J,IE)*RADIUS
IF (OPAERO .EQ. 0) GO TO 12
CALL HLAERO(VEL,VELX,VELY,VELZ,THETAF,PHIF,DEBUG)
A1(1,1)=-CFDF
A1(1,2)=-CFDP
A1(1,3)=-CFDY
A1(1,12)=-CFLW
A1(1,13)=-CFLV
A1(2,1)=-CPDF
A1(2,2)=-CPDP
A1(2,3)=-CPDY
A1(2,12)=-CPLW
A1(2,13)=-CPLV
A1(3,1)=A1(3,1)-CYDF
A1(3,2)=A1(3,2)-CYDP
A1(3,3)=-CYDY
A1(3,12)=-CYLW
A1(3,13)=-CYLV
B(1,1)=CFA
B(1,2)=CFR
B(2,1)=CPA

B(2,2)=CPR
B(3,1)=CYA
B(3,2)=CYR
BG(1,1)=CFDX
BG(1,2)=CFDY
BG(1,3)=CFDZ
BG(2,1)=CPDX
BG(2,2)=CPDY
BG(2,3)=CPDZ
BG(3,1)=CYDX
BG(3,2)=CYDY
BG(3,3)=CYDZ
DO 33 IE=1,NEM
A1(1,IE+3)=-CFDQ(IE)
A1(2,IE+3)=-CPDQ(IE)
A1(3,IE+3)=-CYDQ(IE)
A1(IE+3,1)=-CQDF(IE)
A1(IE+3,2)=-CQDP(IE)
A1(IE+3,3)=-CQDY(IE)
A1(IE+3,12)=-CQLWA(IE)
A1(IE+3,13)=-CQLV(IE)
A0(1,IE+3)=-CFQ(IE)
A0(2,IE+3)=-CPQ(IE)
A0(3,IE+3)=-CYQ(IE)
B(IE+3,1)=CQA(IE)
B(IE+3,2)=CQR(IE)
BG(IE+3,1)=CQDXA(IE)
BG(IE+3,2)=CQDY(IE)
BG(IE+3,3)=CQDZA(IE)
DO 33 JE=1,NEM
A1(IE+3,JE+3)=A1(IE+3,JE+3)-CQDQA(IE,JE)
33 A0(IE+3,JE+3)=A0(IE+3,JE+3)-CQQA(IE,JE)
12 DO 41 I=1,11
AA(I,1)=2.*CA(3,I)
AA(I,2)=CA(1,I)
AA(I,3)=-CA(2,I)
AA(I,4)=-2.*CA(6,I)
AA(I,5)=CA(5,I)
41 AA(I,6)=-CA(4,I)
CB(1,1)=-RSF(1,2)*VELZ
CB(2,1)=-RSF(2,2)*VELZ
CB(3,1)=-RSF(3,2)*VELZ
CB(1,2)=RSF(1,1)*VELZ-RSF(1,3)*VELX
CB(2,2)=RSF(2,1)*VELZ-RSF(2,3)*VELX
CB(3,2)=RSF(3,1)*VELZ-RSF(3,3)*VELX
CB(1,3)=RSF(1,2)*(VELZ*SNT+VELX*CST)
CB(2,3)=RSF(2,2)*(VELZ*SNT+VELX*CST)
CB(3,3)=RSF(3,2)*(VELZ*SNT+VELX*CST)
RETURN
END
SUBROUTINE AIR(ALT,TEMP,DENSE,CSOUND,GRAV,ALTD,OPUNIT)
INTEGER OPUNIT
DIMENSION H1(2),T0(2),T1(2),CS(2),G(2)
DATA H1/145442.,44330.8/
DATA T1/459.67,273.15/
DATA T0/518.67,288.15/
DATA CS/1116.45,340.294/
DATA G/32.17405,9.80665/
TS=1.-ALT/H1(OPUNIT)
T=TS
IF (TEMP .LT. 999.) T=(T1(OPUNIT)+TEMP)/T0(OPUNIT)
DENSE=TS**5.256115/T
CSOUND=CS(OPUNIT)*SQRT(T)

ALTD=H1(OPUNIT)*(1.-DENSE**.234956)
GRAV=G(OPUNIT)
RETURN
END
SUBROUTINE TRLDE
C
C TILTING PROPROTOR AIRCRAFT DYNAMICS
C WAYNE JOHNSON
NASA-ARC/USAAMRDL
C FEBRUARY 1976
C
COMMON /TRCOEF/A2(1936),A1(1936),A0(1936),B0(748),MX,MX1,MV,MVG,DO1F1(44)
INTEGER DOF1
COMMON /RCOEF/A2R(30,30),A1R(30,30),A0R(30,30),AA2R(30,6),AA1R(30,16),AA0R(30,6),BR(30,8),BGR(30,3),C2R(6,30),C1R(6,30),C0R(6,30),CA22R(6,6),CA1R(6,6),CA0R(6,6),DGR(6,3)
COMMON /ACOEF/CAW(6,11),CBW(3,3),A2W(11,13),A1W(11,13),A0W(11,13),1BW(11,2),BGW(11,3),AAW(11,6)
COMMON /AINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,ATILT,OPUNIT,ALTMSL,TEM1P,ALTAGL,ACLIMB,APITCH,COLL,LNGCYC,WEIGHT,AFLAP,ITER,EPTRIM,FACTOR2,PTRIM,DOF(44),CON(14),GUST(3),OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,3OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,4JCNST,NSYSAN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(44),CONP(17),5THETFP,THETAF,CNTRL0, CNTRLS,VELX,VELZ,GRAV,DENSE,CSOUND,FREQ(100)
COMMON /BINPUT/ASHAFT,ACANT,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRT1R,KPRGOV,KIRGOV,KPEGOV,KIEGOV,RTRANS,KMAST,KICS,KENG,IENG,EDAMPC,T2HRTLC,NRPOWR,KRLGE,KPMC(8),KPMS(8),HMAST,FSPIV,WLPIV,BLPIV,FSCG,WL3CG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,FSPYL,WLPYL,IXX,IYY,IZZ,IXZ,MPYLO4N,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,KS1CFE,PSCFE,KT5CFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,PPCFE,KACFE,KRCFE,KPED6AL,CNTRLZ(6),COLLZ0,COLLZ2,K0TRIM,QMASSS(8),QFREQS(8),QDAMPS(8),ZE7TARS(3,8),GAMARS(3,8),QMASSA(8),QFREQA(8),QDAMPA(8),ZETARA(3,8),GA8MARA(3,8),AMASS,AIXX,AIYY,AIZZ,AIXZ,XRTR,YRTR,ZRTR,XWB,ZWB,XHT,ZHT9,XVT,ZVT,TMEAN,TCFES(6,3),TCFEA(6,2),RSF(3,3),RFV(3,3),RSV(3,3),ZAAGL,CSEGE,LAMDA0,TAUW,TAUH,TAUV,CWRTR,CHRTR,CVRTR,NEM,EDAMP,QTHRTL
INTEGER OPUNIT,DOF,CON,GUST,OPSYMM,OPFLOW,OPTORS,OPMOTR,OPLMDA,OPI1NTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,RESPON,D2OFP,CONP
REAL LNGCYC,KRRL,KRRH,KRRF,KLTRIM,KWRTR,KHRTR,KVRTR,KPRGOV,KIRGOV,1KPEGOV,KIEGOV,KMAST,KICS,KENG,IENG,NRPOWR,KRLGE,KPMC,KPMS,IXX,IYY,2IZZ, IXZ,MPYLON,IXXPYL,IYYPYL,IZZPYL,IXZPYL,K00CFE,K01CFE,KS0CFE,K3S1CFE,KTCFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,KACFE,KRCFE,KP4EDAL,LAMDA0
COMMON /RINPUT/HHR(667)
EQUIVALENCE (RADIUS,HHR(14)),(VTIPN,HHR(79)),(IB,HHR(32)),(NBLADE,1HHR(2)),(NBM,HHR(9)),(NTM,HHR(10)),(OMEGA,HHR(13)),(GIMBAL,HHR(11)3),(KPB(1),HHR(50)),(KPG,HHR(35)),(TRC,HHR(21)),(TRS,HHR(22)),(MU,H4HR(6)),(MUY,HHR(88)),(MUZ,HHR(89)),(LAMDA,HHR(7)),(KHLMDA,HHR(655)5),(KFLMDA,HHR(656)),(GAMMA,HHR(12)),(SIGMA,HHR(15)),(OPAERR,HHR(166)),(OPMODE,HHR(26)),(ROTATE,HHR(663))
INTEGER GIMBAL,OPAERR,OPMODE,ROTATE
REAL MU,MUY,MUZ,LAMDA,KHLMDA,KFLMDA,KPB(5),KPG,IB
NAMELIST /NLTR/VTIP,RPM,VEL,VKTS,VELX,VELZ,WEIGHT,ATILT,ASHAFT,ACA1NT,ALTMSL,TEMP,ALTAGL,DENSE,CSOUND,ALTD,GRAV,ACLIMB,APITCH,THETFP,2THETAF,COLL,LNGCYC,AFLAP,DOF,CON,GUST,OPUNIT,OPSYMM,OPFLOW,OPTORS,3OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEB4UG,TONLY,MPER,JCNST,NSYSAN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP5,CONP,ITER,EPTRIM,FACTOR,PTRIM,FREQ,HMAST,FSPIV,WLPIV,BLPIV,FSCG,W6LCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,FSPYL,WLPYL,IXX,IYY,IZZ,IXZ,MPYL7ON,IXXPYL,IYYPYL,IZZPYL,IXZPYL,AMASS,AIXX,AIYY,AIZZ,AIXZ,XRTR,YRTR8,ZRTR,XWB,ZWB,XHT,ZHT,XVT,ZVT,K00CFE,K01CFE,KS0CFE,KS1CFE,PSCFE,KT9CFE,KFCFE,KECFE,KC0CFE,KC1CFE,KP0CFE,KP1CFE,PPCFE,KACFE,KRCFE,KPEDAAL,CNTRLZ,COLLZ0,COLLZ2,K0TRIM,CNTRL0,CNTRLS,TMEAN,TCFES,TCFEA,RSF-

B,RFV,RSV,NEM,QMASSS,QFREQS,QDAMPS,ZETARS,GAMARS,QMASSA,QFREQA,QDAMCPA,ZETARA,GAMARA,KPRGOV,KIRGOV,KPEGOV,KIEGOV,RTRANS,KMAST,KICS,KENDG,IENG,EDAMPC,THRTLC,NRPOWR,KRLGE,KPMC,KPMS,EDAMP,QTHRTL,ZAGL,CSEGEE,LAMDA0,DLDZ,CNSTL,KRRH,KRRF,KRRL,KLTRIM,KWRTR,KHRTR,KVRTR,TAUW,TFAUH,TAUV,CWRTR,CHRTR,CVRTR
NAMELIST /NLEQ/A2,A1,A0,B0,MX,MX1,MV,MVG,MX0,DOF0,DOF1
NAMELIST /NLR/A2R,A1R,A0R,AA2R,AA1R,AA0R,BR,BGR,C2R,C1R,C0R,CA2R,1CA1R,CA0R,DGR
NAMELIST /NLS/CAW,CBW,A2W,A1W,A0W,BW,BGW,AAW
REAL B(616),BG(132),PHI(1800)
EQUIVALENCE (B0(1),B(1)),(B0(617),BG(1))
INTEGER DOF0(44)
REAL CKNOT(2),C2H(11,6),C1H(11,6),C0H(11,6),BGH(30,3),DGH(6,3)
EQUIVALENCE (B0(1),BGH(1,1)),(B0(91),DGH(1,1))
DATA CKNOT/.5924838,1.943845/
DATA CVERT/57.29578/
999 FORMAT (1H1/1X,'TILTING PROPROTOR AIRCRAFT EQUATIONS OF MOTION',101X,18A4//)
998 FORMAT (1X,'OPERATING STATE',69X,'ROTOR PARAMETERS'/
16X,'VELOCITY (KNOTS) =',F8.2,15X,'TILT ANGLE (DEG) =',F7.2,18X,'RA2DIUS (FT OR M) =',F7.2/6X,'V/(OMEGA*R)
=',F8.4,58X,'OMEGA (RA3D/SEC) =',F7.2/6X,'TIP SPEED (FPS OR MPS) =',F7.2,10X,'ALTITUDE M4SL (FT OR M) =',F7.0,12X,'IB (SLUG-FT2 OR KG-M2) =',F7.0)
997 FORMAT (6X,'ROTATIONAL SPEED (RPM) =',F7.2,10X,'TEMPERATURE (F OR 1C) =',F7.2,11X,'NUMBER OF BLADES =',I5)
996 FORMAT (6X,'ROTATIONAL SPEED (RPM) =',F7.2,10X,'STANDARD DAY',31X,1'NUMBER OF BLADES =',I5)
995 FORMAT (6X,'THETA-FP (DEG) =',F7.2,18X,'PRESSURE ALT (FT OR M) =' 1,F7.0,12X,'LOCK NUMBER
=',F7.3/6X,'CLIMB SPEED (FPS OR MPS) =2',F7.2, 8X,'AIR DENSITY RATIO =',F7.4,17X,'ROTOR SOLIDITY =',F7.34)
994 FORMAT (6X,'ALTITUDE AGL (FT OR M) =',F8.2, 9X,'SOUND SPEED (FPS O1R MPS) =',F8.2, 9X,'COUNTER-CLOCKWISE ROTATION'/6X,'CG STATION (FT2 OR M) =',F9.3/6X,'GROSS WEIGHT (LB OR KG) =',F8.0)
993 FORMAT (6X,'ALTITUDE AGL (FT OR M) =',F8.2, 9X,'SOUND SPEED (FPS O1R MPS) =',F8.2, 9X,'CLOCKWISE ROTATION'/6X,'CG STATION (FT OR M) =2',F9.3/6X,'GROSS WEIGHT (LB OR KG) =',F8.0)
992 FORMAT ( /1X,'OPTIONS'/
110X,'OPSYMM =',I2,7X,'IF 1, SYMMETRIC/ IF -1, ANTISYMMETRIC'/
210X,'OPFLOW =',I2,7X,'IF 0, AXIAL FLOW/ LT 0, CONSTANT COEFF'/ 3
10X,'OPTRIM =',I2,7X,'IF 0, NO TRIM/ LT 0, NO IT4ERATIONS/ IF 1, TRIM T,H,MY/ IF 2, TRIM T,H,MY,P'/
510X,'OPMOTR =',I2,7X,'IF 0, AUTOROTATION/ IF 1, CONSTANT SPEED/ IF6 2, MOTOR/ IF 3, POWER OUT'/10X,'OPGOVR =',I2,7X,'IF 0, NO GOVERNO8R/ IF 1, GOVERNOR/ IF 2, INPUT GAINS'/10X,'OPGRND =',I2,7X,'IF 0, 9OUT OF GROUND EFFECT'/
A10X,'OPAERO =',I2,7X,'IF 0, NO PERTURBATION AIRCRAFT AERO FORCES'/B10X,'OPMODE =',I2,7X,'IF 0, RIGID BLADE MODES/ IF 2, READ MODES/ ICF 3, PUNCH MODES/ IF 4, ARTICULATED MODES'/10X,'OPAERO =',I2,7X,'IDF 0, ONLY CLA TERMS IN ROTOR AERO/ LT 0, NO STALL'/
E10X,'OPTORS =',I2,7X,'IF 0, DOF/ IF 1, RIGID/ IF 2, QUASISTATIC'/ F10X,'OPLMDA =',I2,7X,'IF 0, QUASISTATIC INFLOW DYNAMICS'/
G10X,'OPINTR =',I2,7X,'IF 0, QUASISTATIC AERO INTERFERENCE'/
H10X,'OPQRTR =',I2,7X,'IF 1, QUASISTATIC ROTOR/ IF 2, ALL QUASISTATIIC EXCEPT RIGID BODY')
991 FORMAT ( /1X,'ORDER OF VARIABLES'/10X,'DOF = B0,B1C,B1S,BN/2 T01,T1C,T1S,TN/2 BGC,BGS,PSI L0,LC,LS PSIE TF,XF,ZF ELASTIC2 LW,LH'/10X,'CON = T0,T1C,T1S,TN/2 TT DF,DE D0,D1S,DT'/10X3,'GUST = UG,VG,WG')
990 FORMAT ( /1X,'ORDER OF VARIABLES'/10X,'DOF = B0,B1C,B1S,BN/2 T01,T1C,T1S,TN/2 BGC,BGS,PSI L0,LC,LS PSIE PHIF,PSIF,YF ELA2STIC LW,LV'/10X,'CON = T0,T1C,T1S,TN/2 TT DA,DR D1C,DP'/103X,'GUST = UG,VG,WG')
989 FORMAT ( /10X,'DOF =',1X,15I1,3X,9I1,3X,3I1,3X,3I1,3X,I1,3X,3I1,3X1,8I1,3X,2I1/10X,'CON =',1X,8I1,3X,I1,3X,2I1,3X,3I1/10X,'GUST =',1X-

2,3I1)
988 FORMAT (//1X,'COMPUTATION PARAMETERS',15X,'CSYSAN - STEPS =',I4//)
986 FORMAT (1X,'ENGLISH UNITS (FT, SLUG, SEC)')
985 FORMAT (1X,'METRIC UNITS (M, KG, SEC)')
VEL=0.
VKTS=0.
VTIP=0.
RPM=0.
READ (5,NLTR)
IF ((OPUNIT.NE.1).AND.(OPUNIT.NE.2)) OPUNIT=1
IF ((RPM .NE. 0.) .AND. (VTIP .EQ. 0.)) VTIP=RADIUS*RPM*.1047198
IF (VTIP .EQ. 0.) VTIP=VTIPN
IF (RPM .EQ. 0.) RPM=VTIP/RADIUS/.1047198
IF (VEL .EQ. 0.) VEL=VKTS/VTIP/CKNOT(OPUNIT)
IF (VKTS .EQ. 0.) VKTS=VEL*VTIP*CKNOT(OPUNIT)
CALL AIR(ALTMSL,TEMP,DENSE,CSOUND,GRAV,ALTD,OPUNIT)
CALL TRII(T75,T1C,T1S,1)
CALL RTRII(OPREAD,DEBUG,VTIP,T75,T1C,T1S,DENSE,CSOUND,OPUNIT)
IF (OPFLOW .GT. 0) OPFLOW=-1
IF ((OPMOTR .NE. 0) .AND. (OPMOTR .NE. 3)) GO TO 120
PTRIM=0.
IF (OPTRIM .GT. 0) OPTRIM=2
120 IF ((OPSYMM .LT. 0) .OR. (OPMOTR .EQ. 1)) OPGOVR=0
DO 20 K=1,44
IF (DOF(K) .NE. 0) DOF(K)=1
IF (NBLADE .NE. 2) GO TO 123
IF ((K .GE. 11) .AND. (K .LE. 15)) DOF(K)=0
IF ((K .GE. 22) .AND. (K .LE. 24)) DOF(K)=0
123 IF (K .GT. 14) GO TO 20
IF (CON(K) .NE. 0) CON(K)=1
IF ((K .GT. NBLADE) .AND. (K .LE. 8)) CON(K)=0
IF (K .GT. 3) GO TO 20
IF (GUST(K) .NE. 0) GUST(K)=1
20 CONTINUE
IF (NBLADE .EQ. 2) DOF(26)=0
IF (GIMBAL .EQ. 0) DOF(25)=0
IF (GIMBAL .EQ. 0) DOF(26)=0
IF (OPMOTR .EQ. 1) DOF(27)=0
IF (OPMOTR .LE. 1) DOF(31)=0
IF (OPSYMM .LT. 0) CON(14)=0
IF (OPSYMM .GT. 0) GUST(2)=0
IF (OPGRND .NE. 0) DOF(28)=1
IF (OPTORS .EQ. 0) GO TO 124
DO 125 K=16,24
IF (OPTORS .EQ. 1) DOF(K)=0
IF ((OPTORS .EQ. 2) .AND. (K .LE. 15+NBLADE)) DOF(K)=1
125 CONTINUE
124 NBM=0
DO 22 K=1,5
DO 22 KB=1,NBLADE
J=(K-1)*NBLADE+KB
IF (J .GT. 15) GO TO 22
IF (DOF(J) .NE. 0) NBM=K
22 CONTINUE
NBM=MAX0(1,NBM)
NTM=0
DO 23 K=1,3
DO 23 KB=1,NBLADE
J=(K-1)*NBLADE+KB
IF (J .GT. 9) GO TO 23
IF (DOF(J+15) .NE. 0) NTM=K
23 CONTINUE
NTM=MAX0(1,NTM)
NEM=0
DO 24 K=1,8

IF (DOF(K+34) .NE. 0) NEM=K

24 CONTINUE
NEM=MAX0(1,NEM)
VCLIMB=VEL*VTIP*SIN(ACLIMB/CVERT)
PRINT 999, (TITLE(I), I=1,18)
PRINT 998, VKTS,ATILT,RADIUS,VEL,OMEGA,VTIP,ALTMSL,IB
IF (TEMP .LE. 999.) PRINT 997, RPM,TEMP,NBLADE
IF (TEMP .GE. 1000.) PRINT 996, RPM,NBLADE
PRINT 995, ACLIMB,ALTD,GAMMA,VCLIMB,DENSE,SIGMA
IF (ROTATE .EQ. 1) PRINT 994, ALTAGL,CSOUND,FSCG,WEIGHT
IF (ROTATE .EQ.-1) PRINT 993, ALTAGL,CSOUND,FSCG,WEIGHT
PRINT 992, OPSYMM,OPFLOW,OPTRIM,OPMOTR,OPGOVR,OPGRND,OPAERO,OPMODE1,OPAERR,OPTORS,OPLMDA,OPINTR,OPQRTR
IF (OPSYMM .GT. 0) PRINT 991
IF (OPSYMM .LT. 0) PRINT 990
PRINT 989, (DOF(I), I=1,44),(CON(I), I=1,14),(GUST(I), I=1,3)
PRINT 988, JCNST
IF (OPUNIT .EQ. 1) PRINT 986
IF (OPUNIT .EQ. 2) PRINT 985
CALL TRII(T75,T1C,T1S,2)
CALL MODE
CALL TRTRIM
IF (TONLY .NE. 0) GO TO 101
IF (DEBUG .GE. 2) WRITE (6,NLTR)
GO TO 1
101 CONTINUE
NBLD=NBLADE
IF (OPFLOW .LT. 0) NBLD=JCNST
PSI=0.
IF ((NBLADE.GE.3).AND.(OPFLOW.EQ.0)) CALL ROTOR0
IF ((NBLADE.GE.3).AND.(OPFLOW.NE.0)) CALL ROTOR(OPFLOW,PSI,NBLD)
IF (NBLADE .EQ. 2) CALL ROTOR2(OPFLOW,PSI,NBLD)
CALL TRBODY
DO 131 I=1,3
DO 132 K=1,30
BGH(K,I)=0.
DO 132 J=1,3
132 BGH(K,I)=BGH(K,I)+BGR(K,J)*RSV(J,I)
DO 133 K=1,6
DGH(K,I)=0.
DO 133 J=1,3
133 DGH(K,I)=DGH(K,I)+DGR(K,J)*RSV(J,I)
131 CONTINUE
DO 134 I=1,3
DO 135 K=1,30
135 BGR(K,I)=BGH(K,I)
DO 134 K=1,6
134 DGR(K,I)=DGH(K,I)
DO 136 I=1,3
DO 136 K=1,11
BGH(K,I)=0.
DO 136 J=1,3
136 BGH(K,I)=BGH(K,I)+BGW(K,J)*RFV(J,I)
DO 137 I=1,3
DO 137 K=1,11
137 BGW(K,I)=BGH(K,I)
CNST1=1.+KRRL
IF (OPSYMM .LT. 0) CNST1=1.-KRRL
DO 231 I=1,30
A1R(28,I)=CNST1*A1R(28,I)
231 A0R(28,I)=CNST1*A0R(28,I)
DO 232 I=1,6
AA1R(28,I)=CNST1*AA1R(28,I)
AA0R(28,I)=CNST1*AA0R(28,I)
IF (I .GT. 3) GO TO 232

BGR(28,I)=CNST1*BGR(28,I)
232 CONTINUE
DO 31 I=1,1936
A2(I)=0.
A1(I)=0.
31 A0(I)=0.
DO 32 I=1,748
32 B0(I)=0.
DO 33 I=1,30
DO 33 K=1,30
JA=(I-1)*44+K
A2(JA)=A2R(K,I)
A1(JA)=A1R(K,I)
33 A0(JA)=A0R(K,I)
DO 34 I=1,8
DO 34 K=1,30
JB=(I-1)*44+K
34 B(JB)=BR(K,I)
DO 35 I=1,3
DO 35 K=1,30
JB=(I-1)*44+K
35 BG(JB)=BGR(K,I)
CNSTL=GAMMA/(SIGMA*5.7)*DENSE
A2(1216)=A2(1216)+.85*CNSTL
A2(1261)=A2(1261)+.22*CNSTL
A2(1306)=A2(1306)+.22*CNSTL
CNST1=SQRT((MU**2+MUY**2)/KFLMDA**2+LAMDA0**2/KHLMDA**4)
A1(1261)=A1(1261)+CNST1*CNSTL
A1(1306)=A1(1306)+CNST1*CNSTL
CNSTL=CNSTL*2.*(CNST1+(LAMDA0-MUZ)*LAMDA0/KHLMDA**4/CNST1)
A1(1216)=A1(1216)+CNSTL
IF (OPMOTR .LE. 1) GO TO 234
CIENG=IENG*RTRANS**2/(FLOAT(NBLADE)*IB)
CNST1=FLOAT(NBLADE)*IB*OMEGA**2
CICS=0.
IF (OPSYMM .LT. 0) CICS=KICS
CKMR=KMAST/(KMAST+CICS)
CKMI=CKMR*CICS/CNST1
CNST1=KENG/(KMAST+KENG+CICS)/CNST1
CKEM=CNST1*(KMAST+CISC)
CKEI=CNST1*CICS
A2(1171)=A2(1171)+CKMR*CIENG
A2(1347)=A2(1347)+CKMR*CIENG
A1(1171)=A1(1171)+CKMR*EDAMP
A1(1347)=A1(1347)+CKMR*EDAMP
A0(1171)=A0(1171)+CKMI
B(379)=B(379)+CKMR*QTHRTL
A2(1175)=A2(1175)+CIENG
A2(1351)=A2(1351)+CIENG
A1(1175)=A1(1175)+EDAMP
A1(1351)=A1(1351)+EDAMP
A0(1175)=A0(1175)+CKEI
A0(1351)=A0(1351)+CKEM
B(383)=B(383)+QTHRTL
234 CNSTG=KENG/(KMAST+KENG)
IEULER=1
IF (OPSYMM .LT. 0) IEULER=2
DO 41 I=1,11
DO 41 K=1,30
JA=(I+30)*44+K
DO 41 J=1,6
A2(JA)=A2(JA)+AA2R(K,J)*CAW(J,I)
A1(JA)=A1(JA)+AA1R(K,J)*CAW(J,I)
IF (I .LE. IEULER) GO TO 41
A0(JA)=A0(JA)+AA0R(K,J)*CAW(J,I)

41 CONTINUE
DO 42 I=1,30
DO 42 K=1,11
JA=(I-1)*44+K+31
DO 42 J=1,6
A2(JA)=A2(JA)-2.*AAW(K,J)*C2R(J,I)
A1(JA)=A1(JA)-2.*AAW(K,J)*C1R(J,I)
42 A0(JA)=A0(JA)-2.*AAW(K,J)*C0R(J,I)
DO 43 I=1,6
DO 43 K=1,11
C2H(K,I)=0.
C1H(K,I)=0.
C0H(K,I)=0.
DO 43 J=1,6
C1H(K,I)=C1H(K,I)+2.*AAW(K,J)*CA1R(J,I)
C0H(K,I)=C0H(K,I)+2.*AAW(K,J)*CA0R(J,I)
IF ((I .LE. 3) .OR. (J .LE. 3)) GO TO 43
C2H(K,I)=C2H(K,I)+2.*AAW(K,J)*CA2R(J,I)
43 CONTINUE
DO 44 I=1,13
DO 44 K=1,11
JA=(I+30)*44+K+31
A2(JA)=A2W(K,I)
A1(JA)=A1W(K,I)
A0(JA)=A0W(K,I)
IF (I .GT. 11) GO TO 44
DO 48 J=1,6
A2(JA)=A2(JA)-C2H(K,J)*CAW(J,I)
A1(JA)=A1(JA)-C1H(K,J)*CAW(J,I)
IF (I .LE. IEULER) GO TO 48
A0(JA)=A0(JA)-C0H(K,J)*CAW(J,I)
48 CONTINUE
44 CONTINUE
DO 45 I=1,IEULER
IF (OPSYMM .GT. 0) IC=2
IF (OPSYMM .LT. 0) IC=2*I-1
DO 45 K=1,30
JA=(I+30)*44+K
DO 45 J=1,3
A1(JA)=A1(JA)+AA2R(K,J)*CBW(J,IC)
IF (K .GT. 11) GO TO 45
A1(JA+31)=A1(JA+31)-C2H(K,J)*CBW(J,IC)
45 CONTINUE
DO 46 I=1,2
DO 46 K=1,11
JB=(I+8)*44+K+31
46 B(JB)=BW(K,I)
DO 47 I=1,3
DO 47 K=1,11
JB=(I-1)*44+K+31
BG(JB)=BGW(K,I)
DO 47 J=1,6
47 BG(JB)=BG(JB)+2.*AAW(K,J)*DGR(J,I)
A2(1891)=TAUW
A1(1891)=1.
A1(1231)=2.*CWRTR
A1(1936)=1.
IF (OPSYMM .LT. 0) GO TO 241
A2(1936)=TAUH
A1(1232)=2.*CHRTR
GO TO 242
241 A2(1936)=TAUV
A1(1232)=2.*CVRTR
242 IF (OPGRND .EQ. 0) GO TO 243
GLAMDA=(LAMDA0-MUZ)/(1.-(CSEGE*.25/ZAGL)**2)

IF (OPSYMM .LT. 0) GO TO 244

DLDZ=CNSTL*GLAMDA*CSEGE**2*.125/ZAGL**3
SNT=SIN(THETAF)
CST=COS(THETAF)
A0(1392)=A0(1392)-DLDZ*(ZRTR*SNT+XRTR*CST)
A0(1436)=A0(1436)-DLDZ*SNT
A0(1480)=A0(1480)+DLDZ*CST
DO 245 I=1,NEM
JA=(I+33)*44+28
245 A0(JA)=A0(JA)+DLDZ*(-ZETARS(1,I)*SNT+ZETARS(3,I)*CST)
GO TO 243
244 DLDZS=GLAMDA*CSEGE**2*.125/ZAGL**3
ZEFF=ZAGL/CSEGE
DLDZUS=-GLAMDA*(((.1649*ZEFF-1.216)*ZEFF+2.726)*ZEFF-1.722)
DLDZ=CNSTL*(DLDZS+AMIN1(0.,KRLGE*(DLDZUS-DLDZS)))
SNT=SIN(THETAF)
CST=COS(THETAF)
A0(1392)=A0(1392)+DLDZ*YRTR*CST
DO 246 I=1,NEM
JA=(I+33)*44+28
246 A0(JA)=A0(JA)+DLDZ*(-ZETARA(1,I)*SNT+ZETARA(3,I)*CST)
243 CNSTC=B(61)*RADIUS
CNSTS=B(106)*RADIUS
IF (NBLADE .EQ. 2) GO TO 151
IF (OPTORS .EQ. 1) GO TO 153
DO 152 I=1,NEM
JA=(I+33)*44+17
A0(JA)=A0(JA)+CNSTC*KPMC(I)
152 A0(JA+1)=A0(JA+1)+CNSTS*KPMS(I)
GO TO 153
151 SNPSI=SIN(PSI)
CSPSI=COS(PSI)
DO 154 I=1,NEM
JA=(I+33)*44+17
154 A0(JA)=A0(JA)+CNSTC*KPMC(I)*CSPSI+CNSTC*KPMS(I)*SNPSI
153 IF (OPTORS .NE. 1) GO TO 50
JN=44*NBLADE
DO 52 J=1,JN
52 B(J)=-A0(J+660)
DO 53 I=1,NBM
DO 53 J=1,JN
JA=(I-1)*JN+J
IF (JA .GT. 660) GO TO 53
A0(JA)=A0(JA)-KPB(I)*A0(J+660)
53 CONTINUE
IF (GIMBAL .EQ. 0) GO TO 54
JN=88
IF (NBLADE .EQ. 2) JN=44
DO 55 J=1,JN
55 A0(J+1056)=A0(J+1056)-KPG*A0(J+704)
54 IF (NBLADE .EQ. 2) GO TO 57
DO 58 J=1,44
58 A0(J+1144)=A0(J+1144)+TRS*A0(J+704)-TRC*A0(J+748)
GO TO 56
57 TCYC=-TRC*SNPSI+TRS*CSPSI
DO 59 J=1,44
59 A0(J+1144)=A0(J+1144)+TCYC*A0(J+704)
56 IF (NBLADE .EQ. 2) GO TO 155
DO 156 I=1,NEM
DO 156 J=1,44
JA=(I+33)*44+J
156 A0(JA)=A0(JA)-(KPMC(I)*A0(J+704)+KPMS(I)*A0(J+748))*RADIUS
GO TO 50
155 DO 158 I=1,NEM
DO 158 J=1,44

JA=(I+33)*44+J
158 A0(JA)=A0(JA)-(KPMC(I)*CSPSI+KPMS(I)*SNPSI)*A0(J+704)*RADIUS
50 IF (OPGOVR .EQ. 0) GO TO 51
CNST1=1.
IF (OPGOVR .NE. 2) CNST1=AMAX1(0.,AMIN1(1.,ATILT/90.))
DO 251 I=1,44
CNSTB=KIRGOV*B(I)+KIEGOV*B(I+352)
A0(I+1144)=A0(I+1144)-CNSTB
A0(I+1320)=A0(I+1320)-CNSTB*CNSTG
CNSTB=(KPRGOV*B(I)+KPEGOV*B(I+352))*OMEGA*CNST1
A1(I+1144)=A1(I+1144)-CNSTB
251 A1(I+1320)=A1(I+1320)-CNSTB*CNSTG
51 IF (OPSYMM .LT. 0) GO TO 331
DO 332 I=1,3
DO 332 K=1,44
DO 332 J=1,6
JB=(I+10)*44+K
J1=(J-1)*44+K
IF (J .GE. 4) J1=J1+220
332 B(JB)=B(JB)+B(J1)*TCFES(J,I)
GO TO 333
331 DO 334 I=1,2
DO 334 K=1,44
DO 334 J=1,6
JB=(I+10)*44+K
J1=(J-1)*44+K
IF (J .GE. 4) J1=J1+220
334 B(JB)=B(JB)+B(J1)*TCFEA(J,I)
333 CONTINUE
MX=0
MX1=0
MX0=0
MV=0
MVG=0
DO 61 I=1,44
MX=MX+DOF(I)
DOF1(I)=1
61 DOF0(I)=1
IF ((OPFLOW .NE. 0) .OR. (OPSYMM .LT. 0)) GO TO 62
IF ((OPGOVR.NE.0) .AND. ((KIRGOV.NE.0.).OR.(KIEGOV.NE.0.))) GO TO 162
DOF1(27)=0
IF (DOF(27) .NE. 0) MX1=MX1+1
62 DO 63 I=28,30
DOF1(I)=0
IF (DOF(I) .EQ. 0) GO TO 63
MX1=MX1+1
IF (OPLMDA .NE. 0) GO TO 63
DOF0(I)=0
MX0=MX0+1
63 CONTINUE
DO 64 I=43,44
DOF1(I)=0
IF (DOF(I) .EQ. 0) GO TO 64
MX1=MX1+1
IF (OPINTR .NE. 0) GO TO 64
DOF0(I)=0
MX0=MX0+1
64 CONTINUE
IF (OPTORS .NE. 2) GO TO 65
DO 66 I=16,24
IF (DOF(I) .EQ. 0) GO TO 66
DOF0(I)=0
MX0=MX0+1
66 CONTINUE

65 IF (A0(1436) .EQ. 0.) DOF1(33)=0

IF (A0(1480) .EQ. 0.) DOF1(34)=0
IF ((DOF1(33).EQ.0).AND.(DOF(33).EQ.1)) MX1=MX1+1
IF ((DOF1(34).EQ.0).AND.(DOF(34).EQ.1)) MX1=MX1+1
IF (OPQRTR .EQ. 0) GO TO 161
DO 162 I=1,31
IF (DOF(I) .EQ. 0) GO TO 162
IF (DOF0(I) .EQ. 0) GO TO 162
IF ((OPSYMM .GT. 0) .AND. (I .EQ. 27)) GO TO 162
DOF0(I)=0
MX0=MX0+1
162 CONTINUE
IF (OPQRTR .EQ. 1) GO TO 161
DO 163 I=35,44
IF (DOF(I) .EQ. 0) GO TO 163
IF (DOF0(I) .EQ. 0) GO TO 163
DOF0(I)=0
MX0=MX0+1
163 CONTINUE
161 DO 67 I=1,14
67 MV=MV+CON(I)
DO 68 I=1,3
68 MVG=MVG+GUST(I)
CALL DERED(44,14,3,DOF,CON,GUST,A2,A1,A0,B,BG,DOF1,DOF0)
MVX=MX*MV
MVGX=MX*MVG
DO 69 I=1,MVGX
I1=I+MVX
69 B0(I1)=BG(I)
MV=MV+MVG
IF (MX0 .NE. 0) CALL QSTRAN(MX,MX0,MX1,MV,A2,A1,A0,B0,DOF1,DOF0)
IF (DEBUG .EQ. 0) GO TO 81
MXSQ1=MX*MX+1
IF (MXSQ1 .GT. 1936) GO TO 82
DO 83 I=MXSQ1,1936
A2(I)=0.
A1(I)=0.
83 A0(I)=0.
82 MVX1=MX*MV+1
IF (MVX1 .GT. 748) GO TO 81
DO 85 I=MVX1,748
85 B0(I)=0.
81 IF (DEBUG .GE. 2) WRITE(6,NLTR)
IF (DEBUG .EQ. 3) WRITE(6,NLR)
IF (DEBUG .EQ. 3) WRITE(6,NLS)
IF ((DEBUG.EQ.1).OR.(DEBUG.EQ.3)) WRITE(6,NLEQ)
1 CONTINUE
RETURN
END
SUBROUTINE HLAERO(VEL,VELX,VELY,VELZ,THETAF,PHIF,DEBUG)
INTEGER DEBUG
COMMON /FAERO/CZDZ,CZDX,CZDT,CZDDZ,CXDZ,CXDX,CXDT,CTDZ,CTDX,CTDT,C1TDDZ,CZDQ(8),CZQ(8),CXDQ(8),CXQ(8),CTDQ(8),CTQ(8),CQDZ(8),CQDX(8),2CQDT(8),CQDQ(8,8),CQQ(8,8),CZF,CZE,CZLW,CZLH,CXF,CXE,CXLW,CXLH,CTF3,CTE,CTLW,CTLH,CQF(8),CQE(8),CQLW(8),CQLH(8),CFDF,CFDP,CFDY,CPDF,C4PDP,CPDY,CYDF,CYDP,CYDY,CFDQ(8),CFQ(8),CPDQ(8),CPQ(8),CYDQ(8),CYQ(58),CQDF(8),CQDP(8),CQDY(8),CQDQA(8,8),CQQA(8,8),CFA,CFR,CFLV,CPA,C6PR,CPLV,CYA,CYR,CYLV,CQA(8),CQR(8),CQLV(8),CFDZ,CFDX,CFLW,CPDZ,CPD7X,CPLW,CYDZ,CYDX,CYLW,CQDZA(8),CQDXA(8),CQLWA(8)
COMMON /XTRIM/LW,DW,MYW,YW,MXW,MZW,LH,DH,LV,DV,LAW,LFW,LEW,DAW,DFW1,DEW,MAW,MFW,NXBW,NZBW,YBW,LAH,LEH,DAH,DEH,LAV,LEV,DAV,DEV
REAL LW,MYW,MXW,MZW,LH,LV,LAW,LFW,LEW,MAW,MFW,NXBW,NZBW,LAH,LEH,LA1V,LEV

COMMON /FADATA/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,DRG10W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB,S2PANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,LF3TEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,XWB,ZWB,XHT,ZHT,XV4T,ZVT,RADIUS,GAMSA,AFLAP
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
COMMON /FEDATA/EZZ(10,10),EXX(10,10),EZX(10,10),EXZ(10,10),EZ(10),1EX(10),ET(10),EZT(10,10),EXT(10,10),ETT(10,10),EZD(10),EXD(10),EZL2(10),EXL(10),ETL(10),FZH(10),FXH(10),FTH(10),FYW(10),FFH(10),FYV(130),FFV(10),FPV(10),IDENT(10)
REAL KLAW,KDAW,KLUW,KDUW,KLTW,KDTW,KLEW,KDEW,KLWW,KDWW,KLAH,KDAH,K1LUH,KDUH,KLTH,KDTH,KLEH,KDEH,KLWH,KDWH,KLAV,KDAV,KLUV,KDUV,KLTV,KD2TV,KLEV,KDEV,KLWV,KDWV,KMAW,KMUW,KLFW,KDFW
REAL ZERO(530),LVT
EQUIVALENCE (ZERO(1),CZDZ)
NAMELIST /NLDB/LW,DW,MYW,YW,MXW,MZW
1,LH,DH,LV,DV,LAW,LFW,LEW,DAW,DFW,DEW,MAW,MFW,NXBW,NZBW,YBW,LAH,LEH2,DAH,DEH,LAV,LEV,DAV,DEV,
3CZDZ,CZDX,CZDT,CZDDZ,CXDZ,CXDX,CXDT,CTDZ,CTDX,CTDT,CTDDZ,CZDQ,CZQ,4CXDQ,CXQ,CTDQ,CTQ,CQDZ,CQDX,CQDT,CQDQ,CQQ,CZF,CZE,CZLW,CZLH,CXF,CX5E,CXLW,CXLH,CTF,CTE,CTLW,CTLH,CQF,CQE,CQLW,CQLH,CFDF,CFDP,CFDY,CPD6F,CPDP,CPDY,CYDF,CYDP,CYDY,CFDQ,CFQ,CPDQ,CPQ,CYDQ,CYQ,CQDF,CQDP,CQ7DY,CQDQA,CQQA,CFA,CFR,CFLV,CPA,CPR,CPLV,CYA,CYR,CYLV,CQA,CQR,CQLV,8CFDZ,CFDX,CFLW,CPDZ,CPDX,CPLW,CYDZ,CYDX,CYLW,CQDZA,CQDXA,CQLWA, F
KLAW,KDAW,KLUW,KDUW,KLTW,KDTW,KLEW,KDEW,KLWW,KDWW,KLAH,KDAH,KGLUH,KDUH,KLTH,KDTH,KLEH,KDEH,KLWH,KDWH,KLAV,KDAV,KLUV,KDUV,KLTV,KDHTV,KLEV,KDEV,KLWV,KDWV,KMAW,KMUW,KLFW,KDFW,
ALFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,DRG0W,DRGDW,DRGFW,BDRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB,SPANW,AREAW,CHORCDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,AMAXH,IHT,SDPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,EETT,EZD,EXD,EZL,EXL,ETL,FZH,FXH,FTH,FYW,FFH,FYV,FFV,FPV,IDENT
FLA(X1,X2,X3,X4)=((X2+X3)*CS*CS+2.*X3*SN*SN+(X4+X1)*SN*CS)/RAREA
FDA(X1,X2,X3,X4)=((X4-X1)*CS*CS-2.*X1*SN*SN-(X2-X3)*SN*CS)/RAREA
FLU(X1,X2,X3,X4)=(2.*X1*CS*CS-(X4-X1)*SN*SN-(X2-X3)*SN*CS)/RAREA
FDU(X1,X2,X3,X4)=(2.*X3*CS*CS+(X2+X3)*SN*SN-(X4+X1)*SN*CS)/RAREA
FLD(X1,X2)=(X1*CS+X2*SN)/RAREA
FDL(X1,X2)=(X2*CS-X1*SN)/RAREA
DO 1 I=1,530
1 ZERO(I)=0.
IF (VEL .EQ. 0.) RETURN
QY=COS(PHIF)
RX=-SIN(THETAF)
RY=SIN(PHIF)*COS(THETAF)
RZ=COS(PHIF)*COS(THETAF)
V1=GAMSA*VEL
V2=V1*VEL
RAREA=3.141593*RADIUS**2
RAREAR=RAREA*RADIUS
PHIW=DIHEDW/57.29578
PSIW=SWEEPW/57.29578
CS=0.
SN=0.
IF ((VELX .EQ. 0.) .AND. (VELZ .EQ. 0.)) GO TO 2
CS=VELX/SQRT(VELX**2+VELZ**2)
SN=VELZ/SQRT(VELX**2+VELZ**2)
2 KLAW=FLA(LW,LAW,DW,DAW)
KDAW=FDA(LW,LAW,DW,DAW)
KLUW=FLU(LW,LAW,DW,DAW)
KDUW=FDU(LW,LAW,DW,DAW)
KMAW=(MAW*CS+2.*MYW*SN)/RAREAR
KMUW=(2.*MYW*CS-MAW*SN)/RAREAR
KLTW=FLD(LAW,DAW)
KDTW=FDL(LAW,DAW)

KLEW=FLD(LEW,DEW)
KDEW=FDL(LEW,DEW)
KLFW=FLD(LFW,DFW)
KDFW=FDL(LFW,DFW)
KLWW=FLD(LAW+DW,DAW-LW)
KDWW=FDL(LAW+DW,DAW-LW)
KLAH=FLA(LH,LAH,DH,DAH)
KDAH=FDA(LH,LAH,DH,DAH)
KLUH=FLU(LH,LAH,DH,DAH)
KDUH=FDU(LH,LAH,DH,DAH)
KLTH=FLD(LAH,DAH)
KDTH=FDL(LAH,DAH)
KLEH=FLD(LEH,DEH)
KDEH=FDL(LEH,DEH)
KLWH=FLD(LAH+DH,DAH-LH)
KDWH=FDL(LAH+DH,DAH-LH)
CS=0.
SN=0.
IF ((VELX .EQ. 0.) .AND. (VELY .EQ. 0.)) GO TO 3
CS=VELX/SQRT(VELX**2+VELY**2)
SN=VELY/SQRT(VELX**2+VELY**2)
3 KLAV=FLA(LV,LAV,DV,DAV)
KDAV=FDA(LV,LAV,DV,DAV)
KLUV=FLU(LV,LAV,DV,DAV)
KDUV=FDU(LV,LAV,DV,DAV)
KLTV=FLD(LAV,DAV)
KDTV=FDL(LAV,DAV)
KLEV=FLD(LEV,DEV)
KDEV=FDL(LEV,DEV)
KLWV=FLD(LAV+DV,DAV-LV)
KDWV=FDL(LAV+DV,DAV-LV)
DEDA=0.
IF (AREAW .NE. 0.) DEDA=.45*LAW/AREAW/((SPANW**2/AREAW)**.735*(ABS1(XHT*RADIUS/CHORDW))**.25)
DSDB=1.
DSDP=0.
DSDR=1.
DEDA1=1.-DEDA
DLF=LFLAPO-LFLAPI
DLA=LAILO**2-LAILI**2
CZDZ=-V1*(KLAW+DEDA1*KLAH)
CZDX=-V1*(KLUW+KLUH)
CZDT=-V1*QY*(ZWB*KLUW-XWB*KLAW+ZHT*KLUH-XHT*KLAH)
CZDDZ=-V1*DEDA*(ZHT*KLUH-XHT*KLAH)
CXDZ=-V1*(KDAW+DEDA1*KDAH)
CXDX=-V1*(KDUW+KDUH)
CXDT=-V1*QY*(ZWB*KDUW-XWB*KDAW+ZHT*KDUH-XHT*KDAH)
CTDZ=-V1*QY*(ZWB*KDAW-XWB*KLAW+DEDA1*(ZHT*KDAH-XHT*KLAH)-KMAW)
CTDX=-V1*QY*(ZWB*KDUW-XWB*KLUW+DEDA1*(ZHT*KDUH-XHT*KLUH)-KMUW)
CTDT=-V1*QY*QY*(XWB*XWB*KLAW+ZWB*ZWB*KDUW-XWB*ZWB*(KDAW+KLUW)+XHT*1XHT*KLAH+ZTH*ZHT*KDUH-XHT*ZHT*(KDAH+KLUH))
CTDDZ=-V1*QY*DEDA*(XHT*XHT*KLAH+ZHT*ZHT*KDUH-XHT*ZHT*(KDAH+KLUH))
CZF=-V2*DLF*KLFW
CZE=-V2*KLEH
CZLW=V1*KLWW
CZLH=V1*KLWH
CXF=-V2*DLF*KDFW
CXE=-V2*KDEH
CXLW=V1*KDWW
CXLH=V1*KDWH
CTF=-V2*QY*DLF*(ZWB*KDFW-XWB*KLFW-MFW/RAREAR)
CTE=-V2*QY*(ZHT*KDEH-XHT*KLEH)
CTLW=V1*QY*(ZWB*KDWW-XWB*KLWW-MAW/RAREAR)
CTLH=V1*QY*(ZHT*KDWH-XHT*KLWH)
LVT=XVT*RZ-ZVT*RX

S1=.25*SPANW/RADIUS
S2=(SPANW/RADIUS)**2/12.
CFDF=-V1*(S2*KLAW+DSDP*ZVT*ZVT*KLAV)
CFDP=-V1*(S2*(-RZ*KLUW+RX*KLAW)-DSDR*(ZVT*LVT*KLAV+ZVT*ZVT*RY*KLUV1))
CFDY=-V1*(S1*(PSIW*KLUW+PHIW*KLAW)*KCLB-DSDB*ZVT*KLAV-NXBW/RAREAR)
CPDF=-V1*(S2*(-RZ*KDAW+RX*KLAW)-DSDP*(ZVT*LVT*KLAV+ZVT*ZVT*RY*KDAV1))
CPDP=-V1*(S2*(RZ*RZ*KDUW+RX*RX*KLAW-RZ*RX*(KDAW+KLUW))+DSDR*(LVT*L1VT*KLAV+(ZVT*RY)**2*KDUV+LVT*ZVT*RY*(KDAV+KLUV)))
CPDY=-V1*(-S1*RZ*(PSIW*KDUW+PHIW*KDAW)+S1*RX*(PSIW*KLUW+PHIW*KLAW)1+DSDB*(LVT*KLAV+ZVT*RY*KDAV)-RZ*NZBW/RAREAR)
CYDF=-V1*(S1*(PSIW*KDAW+PHIW*KLAW)-DSDP*ZVT*KLAV)
CYDP=-V1*(-S1*RZ*(PSIW*KDUW+PHIW*KLUW)+S1*RX*(PSIW*KDAW+PHIW*KLAW)1+DSDR*(LVT*KLAV+ZVT*RY*KLUV))
CYDY=-V1*(DSDB*KLAV-YBW/RAREA)
CFA=-V2*S1*DLA*KLEW
CFR=V2*ZVT*KLEV
CFLV=-V1*ZVT*KLWV
CPA=-V2*S1*DLA*(RX*KLEW-RZ*KDEW)
CPR=-V2*(LVT*KLEV+ZVT*RY*KDEV)
CPLV=V1*(LVT*KLWV+ZVT*RY*KDWV)
CYA=-V2*S1*DLA*(PSIW*KDEW+PHIW*KLEW)
CYR=-V2*KLEV
CYLV=V1*KLWV
CFDZ=-V1*S1*KLAW
CFDX=-V1*S1*KLUW
CFLW=V1*S1*KLWW
CPDZ=-V1*S1*(-RZ*KDAW+RX*KLAW)
CPDX=-V1*S1*(-RZ*KDUW+RX*KLUW)
CPLW=V1*S1*(-RZ*KDWW+RX*KLWW)
CYDZ=-V1*(PSIW*KDAW+PHIW*KLAW)
CYDX=-V1*(PSIW*KDUW+PHIW*KLUW)
CYLW=V1*(PSIW*KDWW+PHIW*KLWW)
IA=0
IS=0
DO 40 I=1,10
IF (IDENT(I)) 41,40,42
42 IS=IS+1
CEX=EX(I)
CEZ=EZ(I)
CFX=FXH(I)
CFZ=FZH(I)
CZDQ(IS)=-V1*(-CEX*KLUW-CEZ*KLAW+CFX*KLUH+CFZ*KLAH)
CZQ(IS)=-V2*(ET(I)*KLTW+FTH(I)*KLTH)
CXDQ(IS)=-V1*(-CEX*KDUW-CEZ*KDAW+CFX*KDUH+CFZ*KDAH)
CXQ(IS)=-V2*(ET(I)*KDTW+FTH(I)*KDTH)
CTDQ(IS)=-V1*QY*(-CEX*(ZWB*KDUW-XWB*KLUW)-CEZ*(ZWB*KDAW-ZWB*KLAW)+1CFX*(ZHT*KDUH-XHT*KLUH)+CFZ*(ZHT*KDAH-XHT*KLAH))
CTQ(IS)=-V2*(ET(I)*(ZWB*KDTW-XWB*KLTW)+FTH(I)*(ZHT*KDTH-XHT*KLTH))
CQDZ(IS)=-V1*(-CEX*KDAW-CEZ*KLAW+CFX*KDAH+CFZ*KLAH)
CQDX(IS)=-V1*(-CEX*KDUW-CEZ*KLUW+CFX*KDUH+CFZ*KLUH)
CQDT(IS)=-V1*QY*(-CEX*(ZWB*KDUW-XWB*KDAW)-CEZ*(ZWB*KLUW-XWB*KLAW)+1CFX*(ZHT*KDUH-XHT*KDAH)+CFZ*(ZHT*KLUH-XHT*KLAH))
CQF(IS)= V2*(EXD(I)*KDFW+EZD(I)*KLFW)
CQE(IS)=-V2*(CFX*KDEH+CFZ*KLEH)
CQLW(IS)=V1*(-CEX*KDWW-CEZ*KLWW)
CQLH(IS)=V1*(CFX*KDWH+CFZ*KLWH)
JS=0
DO 43 J=1,10
IF (IDENT(J) .NE. 1) GO TO 43
JS=JS+1
CQDQ(IS,JS)=-V1*(EZZ(I,J)*KLAW+EZX(I,J)*KLUW+EXZ(I,J)*KDAW+EXX(I,J1)*KDUW+.125*AREAW*CHORDW/RAREAR*ETT(I,J)+CFZ*(FZH(J)*KLAH+FXH(J)*K2LUH)+CFX*(FZH(J)*KDAH+FXH(J)*KDUH))

CQQ(IS,JS)=-V2*(-EXT(I,J)*KDTW-EZT(I,J)*KLTW+FTH(I)*(FZH(J)*KLTH+F1XH(J)*KDTH))
43 CONTINUE
GO TO 40
41 IA=IA+1
CEX=EXL(I)
CEZ=EZL(I)
S3=(SPANH/RADIUS)**2/12.
S4=(SPANV/RADIUS)**2/12.
CFDQ(IA)=-V1*(-CEX*KLUW-CEZ*KLAW+S3*FFH(I)*KLAH-FYV(I)*ZVT*KLAV)
CFQ(IA)=-V2*(ETL(I)*KLTW+FPV(I)*ZVT*KLTV)
CPDQ(IA)=-V1*(-CEX*(-RZ*KDUW+RX*KLUW)-CEZ*(-RZ*KDAW+RX*KLAW)+FFH(I1)*S3*(-RZ*KDAH+RX*KLAH)+FYV(I)*(LVT*KLAV+ZVT*RY*KDAV))
CPQ(IA)=-V2*(ETL(I)*(RX*KLTW-RZ*KDTW)-FPV(I)*(LVT*KLTV+ZVT*RY*KDTV1))
CYDQ(IA)=-V1*(-EX(I)*(PSIW*KDUW+PHIW*KLUW)-EZ(I)*(PSIW*KDAW+PHIW*K1LAW)+FYV(I)*KLAV)
CYQ(IA)=-V2*(ET(I)*(PHIW*KLTW+PSIW*KDTW)-FPV(I)*KLTV)
CQDF(IA)=-V1*(-CEX*KDAW-CEZ*KLAW+FFH(I)*S3*KLAH-FYV(I)*ZVT*KLAV)
CQDP(IA)=-V1*(-CEX*(-RZ*KDUW+RX*KDAW)-CEZ*(-RZ*KLUW+RX*KLAW)+FFH(I1)*S3*(-RZ*KLUH+RX*KLAH)+FYV(I)*(LVT*KLAV+ZVT*RY*KLUV))
CQDY(IA)=-V1*(-EX(I)*(PSIW*KDUW+PHIW*KDAW)-EZ(I)*(PSIW*KLUW+PHIW*K1LAW)+FYV(I)*KLAV)
CQA(IA)=V2*(EXD(I)*KDEW+EZD(I)*KLEW)
CQR(IA)=-V2*FYV(I)*KLEV
CQLV(IA)=V1*FYV(I)*KLWV
CEX=EX(I)-FYW(I)*PSIW
CEZ=EZ(I)-FYW(I)*PHIW
CQDZA(IA)=V1*(CEX*KDAW+CEZ*KLAW)
CQDXA(IA)=V1*(CEX*KDUW+CEZ*KLUW)
CQLWA(IA)=-V1*(CEX*KDWW+CEZ*KLWW)
JA=0
DO 44 J=1,10
IF (IDENT(J) .NE. -1) GO TO 44
JA=JA+1
CQDQA(IA,JA)=-V1*(EZZ(I,J)*KLAW+EZX(I,J)*KLUW+EXZ(I,J)*KDAW+EXX(I,1J)*KDUW+.125*AREAW*CHORDW/RAREAR*ETT(I,J)+FFH(I)*FFH(J)*S3*KLAH+(F2YV(I)*FYV(J)+S4*FFV(I)*FFV(J))*KLAV)
CQQA(IA,JA)=V2*(EXT(I,J)*KDTW+EZT(I,J)*KLTW+FYV(I)*FPV(J)*KLTV)
44 CONTINUE
40 CONTINUE
IF (DEBUG .EQ. 3) WRITE (6,NLDB)
RETURN
END
SUBROUTINE ACHAR(ALPHA,BETA,DELF,DELA,ALFHT,DELE,ALFVT,DELR)
COMMON /XTRIM/LW,DW,MYW,YW,MXW,MZW,LH,DH,LV,DV,LAW,LFW,LEW,DAW,DFW1,DEW,MAW,MFW,NXBW,NZBW,YBW,LAH,LEH,DAH,DEH,LAV,LEV,DAV,DEV
REAL LW,MYW,MXW,MZW,LH,LV,LAW,LFW,LEW,MAW,MFW,NXBW,NZBW,LAH,LEH,LA1V,LEV
COMMON /FADATA/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,DRG10W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB,S2PANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,LF3TEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,XWB,ZWB,XHT,ZHT,XV4T,ZVT,RADIUS,GAMSA,AFLAP
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
DATA CVERT/57.29578/
ALF=ALPHA+IWB/CVERT
IF (ABS(ALF) .GT. AMAXW/CVERT) GO TO 1
DF=AFLAP/CVERT
LW=LFTAW*ALF+LFTDW*DF+LFTFW*(LFLAPO-LFLAPI)*DELF
LAW=LFTAW
LFW=LFTFW

LEW=LFTEW
DW=DRG0W+LW*LW/PIEL2+DRGDW*DF+DRGFW*(LFLAPO-LFLAPI)*DELF
DAW=2.*LW*LAW/PIEL2
DFW=DRGFW
IF (LAW .NE. 0.) DFW=DAW*LFW/LAW+DRGFW
DEW=DRGEW
IF (LAW .NE. 0.) DEW=DAW*LEW/LAW+DRGEW
MYW=MOM0W+MOMAW*ALF+MOMDW*DF+MOMFW*(LFLAPO-LFLAPI)*DELF
MAW=MOMAW
MFW=MOMFW
YW=SIDEB*BETA
YBW=SIDEB
MXW=(-KCLB*.25*SPANW*(SWEEPW*2.*LW+DIHEDW*(LFTAW+DW))/CVERT+ROLLB)1*BETA-.25*SPANW*(LAILO**2-LAILI**2)*LFTEW*DELA
MZW=(.25*SPANW*(SWEEPW*2.*DW+DIHEDW*(DAW-LW))/CVERT+YAWB)*BETA+.251*SPANW*(LAILO**2-LAILI**2)*DEW*DELA
NXBW=ROLLB
NZBW=YAWB
GO TO 2
1 ALFS=SIGN(AMAXW/CVERT,ALF)
LW=LFTAW*ALFS+LFTDW*DF+LFTFW*(LFLAPO-LFLAPI)*DELF
LAW=0.
LFW=0.
LEW=0.
DW=(DRGSW-DRG0W)*(SIN(ALF))**2+DRG0W
DAW=(DRGSW-DRG0W)*SIN(2.*ALF)
DFW=0.
DEW=0.
MYW=0.
MAW=0.
MFW=0.
YW=0.
YBW=0.
MXW=0.
MZW=0.
NXBW=0.
NZBW=0.
2 ALF=ALFHT+IHT/CVERT
IF (ABS(ALF) .GT. AMAXH/CVERT) GO TO 3
LH=LFTAH*ALF+LFTEH*DELE
LAH=LFTAH
LEH=LFTEH
GO TO 4
3 ALFS=SIGN(AMAXH/CVERT,ALF)
LH=LFTAH*ALFS+LFTEH*DELE
LAH=0.
LEH=0.
4 DH=0.
DAH=0.
DEH=0.
ALF=ALFVT+IVT/CVERT
IF (ABS(ALF) .GT. AMAXV/CVERT) GO TO 5
LV=LFTAV*ALF+LFTEV*DELR
LAV=LFTAV
LEV=LFTEV
GO TO 6
5 ALFS=SIGN(AMAXV/CVERT,ALF)
LV=LFTAV*ALFS+LFTEV*DELR
LAV=0.
LEV=0.
6 DV=0.
DAV=0.
DEV=0.
RETURN
END

BLOCK DATA
COMMON /R1DATA/TYPE1,NBLAD1,DEL3G1,NTRM01,NTRM11,MU1,LAMDA1,MTIP1,1NBM1,NTM1,GIMBL1,GAMMA1,OMEGA1,RADUS1,SIGMA1,OPAER1,MRAER1,BBLD1,R2ROOT1,INCMP1,T1C1,T1S1,T751,MPSI1,MREV1,OPMOD1,OFFST1,MRMOD1,NONRT31,NCOLB1,CONE1,IB1,DEL11,XFA1,KPG1,GGS1,GQS1(5),GPS1(3),NCOLT1,WT041,RFA1,HSPRC1,HSPRS1,KPB1(5),RPH1,XPH1,WTCH1,WTSH1,WTRH1,WTC1,WTS15,WTR1,NUGC1,NUGS1,KFLAP1,KLAG1,HINGE1,KPIN1,ATNKP1(5),DROOP1,SWEEP61,DEL21,DEL31,ZFA1,VTIPN1,DEBUG1,EPMOD1,GLAG1,LDAMP1,EXDMP1,PHIPH17,PHIPL1,LDMPC1,MUY1,MUZ1,CHORD1(51),XA1(51),MASS1(51),XI1(51),THET8A1(51),ITHET1(51),TWIST1(51),EIZZ1(51),EIXX1(51),GJ1(51),XC1(51),D9ENSE1,OPUNT1,ELAG1,VTOT1,KHLMD1,KFLMD1,ALFHP1,PSIHP1,T75IN1,T1CIN1A,T1SIN1,IRING1,ROTAT1,FDROP1,FSWEP1,DEL41,DEL51
REAL MU1,LAMDA1,MTIP1,IB1,MASS1,ITHET1,KPG1,NUGC1,NUGS1,KPB1,KFLAP11,KLAG1,LDAMP1,LDMPC1,MUY1,MUZ1,KHLMD1,KFLMD1
INTEGER OPAER1,GIMBL1,OPMOD1,HINGE1,DEBUG1,OPUNT1,ROTAT1
COMMON /R2DATA/TYPE2,NBLAD2,DEL3G2,NTRM02,NTRM12,MU2,LAMDA2,MTIP2,1NBM2,NTM2,GIMBL2,GAMMA2,OMEGA2,RADUS2,SIGMA2,OPAER2,MRAER2,BBLD2,R2ROOT2,INCMP2,T1C2,T1S2,T752,MPSI2,MREV2,OPMOD2,OFFST2,MRMOD2,NONRT32,NCOLB2,CONE2,IB2,DEL12,XFA2,KPG2,GGS2,GQS2(5),GPS2(3),NCOLT2,WT042,RFA2,HSPRC2,HSPRS2,KPB2(5),RPH2,XPH2,WTCH2,WTSH2,WTRH2,WTC2,WTS25,WTR2,NUGC2,NUGS2,KFLAP2,KLAG2,HINGE2,KPIN2,ATNKP2(5),DROOP2,SWEEP62,DEL22,DEL32,ZFA2,VTIPN2,DEBUG2,EPMOD2,GLAG2,LDAMP2,EXDMP2,PHIPH27,PHIPL2,LDMPC2,MUY2,MUZ2,CHORD2(51),XA2(51),MASS2(51),XI2(51),THET8A2(51),ITHET2(51),TWIST2(51),EIZZ2(51),EIXX2(51),GJ2(51),XC2(51),D9ENSE2,OPUNT2,ELAG2,VTOT2,KHLMD2,KFLMD2,ALFHP2,PSIHP2,T75IN2,T1CIN2A,T1SIN2,IRING2,ROTAT2,FDROP2,FSWEP2,DEL42,DEL52
REAL MU2,LAMDA2,MTIP2,IB2,MASS2,ITHET2,KPG2,NUGC2,NUGS2,KPB2,KFLAP12,KLAG2,LDAMP2,LDMPC2,MUY2,MUZ2,KHLMD2,KFLMD2
INTEGER OPAER2,GIMBL2,OPMOD2,HINGE2,DEBUG2,OPUNT2,ROTAT2
COMMON /HINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,OPUNIT,ALTMSL,TEMP,ALTA1GL,ACLIMB,AYAW,APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,WEIGHT,AFLAP,2ITER,EPTRIM,FACTOR,PTRIM,DOF(80),CON(26),GUST(3),OPSYMM,DOFSYM(80)3,CONSYM(26),GUSSYM(3),OPFLOW,OPTORS(2),OPMOTR,OPLMDA(2),OPINTR,OPQ4RTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,JCNST,NSYS5AN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(80),CONP(29),CNTRL0,CNT6RLC,CNTRLS,CNTRLP,THETFP,PSIFP,THETAF,PHIF,VELX,VELY,VELZ,GRAV,DEN7SE,CSOUND,FREQ(100)
COMMON /GINPUT/ORIENT,ASHAFT(2),ACANT(2),K12H,K12F,K12L,K21H,K21F,1K21L,KLTRIM,KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,CWRTR1,CWRTR22,CHRTR1,CHRTR2,CVRTR1,CVRTR2,TAUW,TAUH,TAUV,RTRANS,IENG,KMAST1,KM3AST2,KICS,KENG,ENGPOS,KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,TH4RTLC,NRPOWR,EDAMPC,KPMC1(10),KPMS1(10),KPMC2(10),KPMS2(10),FSR1,WL5R1,BLR1,FSR2,WLR2,BLR2,FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,IXX6,IYY,IZZ,IXZ,XRTR1,YRTR1,ZRTR1,XRTR2,YRTR2,ZRTR2,XWB,ZWB,XHT,ZHT,X7VT,ZVT,AMASS,AIXX,AIYY,AIZZ,AIXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCF8E,PRCCFE,PFPCFE,PRPCFE,KF0CFE,KR0CFE,KFSCFE,KRSCFE,KCCFE,KSCFE,PCC9FE,PSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,CNTRLZ(11),TCFEA(11,5),QMASS(10),QFREQ(10),QDAMP(10),ZETAR1(3,10),ZETAR2(3,10),GAMBAR1(3,10),GAMAR2(3,10),NEM,RSF1(3,3),RSF2(3,3),RFV(3,3),RSV1(3,3),CRSV2(3,3),ZAGL1,ZAGL2,CSEGE1,CSEGE2,LMDA01,LMDA02,EDAMP,QTHRTL
INTEGER OPUNIT,ORIENT,DOF,CON,GUST,OPSYMM,DOFSYM,CONSYM,GUSSYM,OPF1LOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND2,OPAERO,DEBUG,TONLY,RESPON,DOFP,CONP,ENGPOS
REAL LNGCYC,LATCYC,K12H,K12F,K12L,K21H,K21F,K21L,KLTRIM,KWRTR1,KWR1TR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,IENG,KMAST1,KMAST2,KICS,KENG,KPR1G2V,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,NRPOWR,KPMC1,KPMC2,KPMS1,KPMS32,IXX,IYY,IZZ,IXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,KF0CFE,KR0CFE,KFSCFE4,KRSCFE,KCCFE,KSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,LMDA501,LMDA02
COMMON /FINPUT/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ(10,10),EXX(140,10),EZX(10,10),EXZ(10,10),EZ(10),EX(10),ET(10),EZT(10,10),EXT(105,10),ETT(10,10),EZD(10),EXD(10),EZL(10),EXL(10),ETL(10),FZH(10),FX-

6H(10),FTH(10),FYW(10),FFH(10),FYV(10),FFV(10),FPV(10)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
C
DATA TITLE/72H
OLD UTILITY HELICOPTER
1
/
DATA OPUNIT/1/
DATA ORIENT/1/
DATA ALTMSL,TEMP,ALTAGL/0.,1000.,0./
DATA ASHAFT,ACANT/3.,0.,2*0./
DATA ACLIMB,AYAW/0.,0./
DATA APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,AFLAP/7*0./
DATA WEIGHT/11900./
DATA ITER,EPTRIM,FACTOR,PTRIM/10,.01,.5,0./
DATA DOF/4*0,4*1,26*0,4*1,23*0,6*1,13*0/
DATA CON,GUST/3*1,5*0,1,7*0,1,4*0,5*1,3*1/
DATA OPSYMM/0/
DATA DOFSYM,CONSYM,GUSSYM/61*0,-1,1,-1,1,-1,1,13*0,1,-1,1,5*0,-1,71*0,1,4*0,1,-1,1,-1,1,1,-1,1/
DATA OPFLOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREA1D,OPGRND,OPAERO,DEBUG,TONLY/-1,2*1,1,2*0,0,0,1,0,0,0,1,0,1/
DATA MPER,JCNST/24,12/
DATA NSYSAN,NSTEP,NFREQ,FREQ/0,0,0,100*0./
DATA LOCUS,PLOTMX,RESPON,NF,CONP,DOFP/-1,2.99,-1,100,29*0,80*0/
DATA K12H,K12F,K21H,K21F/4*0./
DATA KLTRIM/0./
DATA KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2/2*0.,1.8,2*0.,.7/
DATA RTRANS,IENG,KMAST1,KMAST2,KICS,KENG,ENGPOS/6.,6.,4*3.E4,1/
DATA KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV/6*0./
DATA THRTLC,NRPOWR,EDAMPC/1.,1275.,1./
DATA KPMC1,KPMS1,KPMC2,KPMS2/40*0./
DATA FSR1,WLR1,BLR1,FSR2,WLR2,BLR2/0.,8.,0.,33.,0.,0./
DATA FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT/4*0.,28.,0.,31.,0./
DATA IXX,IYY,IZZ,IXZ/5895.,27500.,23088.,0./
DATA KCCFE,KSCFE,PCCFE,PSCFE,K0CFE,KPCFE/2*1.,2*0.,2*1./
DATA KFCFE,KACFE,KECFE,KRCFE,KTCFE/5*0./
DATA CNTRLZ/11*0./
DATA QMASS,QFREQ,QDAMP/20*1.,10*.01/
DATA ZETAR1,GAMAR1,ZETAR2,GAMAR2/120*0./
DATA EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,1FXH,FTH,FYW,FFH,FYV,FFV,FPV/860*0./
DATA LFTAW,LFTDW,LFTFW,LFTEW/4*0./
DATA MOM0W,MOMAW,MOMDW,MOMFW/0.,570.,2*0./
DATA DRG0W,DRGDW,DRGFW,DRGEW,PIEL2/36.5,3*0.,140./
DATA AMAXW,DRGSW/20.,60./
DATA SIDEB,ROLLB,YAWB,KCLB/-170.,0.,-1700.,0./
DATA IWB,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW/6*0./
DATA LFLAPO,LFLAPI,LAILO,LAILI/4*0./
DATA LFTAH,LFTEH,AMAXH,IHT,SPANH/20.,0.,20.,5.,5.55/
DATA LFTAV,LFTEV,AMAXV,IVT,SPANV/56.,0.,20.,0.,6./
C
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA

TYPE1/4HMAIN/
OPAER1,OPMOD1,INCMP1/-1,1,1/
MRAER1,MRMOD1,NONRT1,NCOLB1,NCOLT1/9,40,0,2,2/
NTRM01,NTRM11,MPSI1,MREV1,EPMOD1/1,2,24,20,.001/
VTIPN1,ROTAT1/650.,1/
NBLAD1,GAMMA1,RADUS1,SIGMA1,IB1/4,9.94,28.,.0622,1145./
BBLD1,RROOT1/.97,.1/
OFFST1,ELAG1/2*.0357/
GGS1,GQS1,GPS1/0.,5*.01,3*.01/
WT01,WTCH1,WTSH1,WTRH1/5.8,5.3,5.3,6.5/
CONE1,DROOP1,SWEEP1,FDROP1,FSWEP1/5*0./
XFA1,ZFA1,RFA1/0.,0.,.08/
KPIN1,ATNKP1,DEL3G1/2,5*0.,0./
RPH1,XPH1,PHIPL1,PHIPH1/.05,.05,0.,0./

DATA HINGE1,KFLAP1,KLAG1,GIMBL1,HSPRC1,HSPRS1/0,0.,0.,0,0.,0./
DATA LDAMP1,EXDMP1/1700.,-.85/
DATA KHLMD1,KFLMD1/1.15,1.04/
DATA CHORD1,XA1,ITHET1,EIZZ1,EIXX1,GJ1,XI1,XC1/51*.0488,51*0.,51*.1016,51*102100.,51*1107000.,51*122200.,51*0.,51*0./
DATA MASS1/2*.1,3.36,2.05,1.77,.586,.172,.131,4*.138,36*.179,.317,1.347,.01/
DATA TWIST1/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
C
DATA TYPE2/4HTAIL/
DATA OPAER2,OPMOD2,INCMP2/-1,4,1/
DATA MRAER2,MRMOD2,NONRT2,NCOLB2,NCOLT2/9,40,0,2,2/
DATA NTRM02,NTRM12,MPSI2,MREV2,EPMOD2/2,2,24,20,.001/
DATA VTIPN2,ROTAT2/650.46,1/
DATA NBLAD2,GAMMA2,RADUS2,SIGMA2,IB2/4,4.54,4.67,.167,.87/
DATA BBLD2,RROOT2/.93,.1/
DATA OFFST2,ELAG2/0.,.1/
DATA GGS2,GQS2,GPS2/0.,5*.01,3*.01/
DATA WT02,WTCH2,WTSH2,WTRH2/5.8,5.3,5.3,6.5/
DATA CONE2,DROOP2,SWEEP2,FDROP2,FSWEP2/5*0./
DATA XFA2,ZFA2,RFA2/0.,0.,.15/
DATA KPIN2,ATNKP2,DEL3G2/2,0.,45.,3*0.,0./
DATA RPH2,XPH2,PHIPL2,PHIPH2/.05,.05,0.,0./
DATA HINGE2,KFLAP2,KLAG2,GIMBL2,HSPRC2,HSPRS2/0,0.,0.,0,0.,0./
DATA LDAMP2,EXDMP2/10.0,0./
DATA KHLMD2,KFLMD2/1.2,1.1/
DATA CHORD2,XA2,MASS2,ITHET2,EIZZ2,EIXX2,GJ2,XI2,XC2/51*.1312,51*01.,51*.0256,51*.00032,51*790.,51*23000.,51*2000.,51*0.,51*0./
DATA TWIST2/51*0./
END
BLOCK DATA
COMMON /R1DATA/TYPE1,NBLAD1,DEL3G1,NTRM01,NTRM11,MU1,LAMDA1,MTIP1,1NBM1,NTM1,GIMBL1,GAMMA1,OMEGA1,RADUS1,SIGMA1,OPAER1,MRAER1,BBLD1,R2ROOT1,INCMP1,T1C1,T1S1,T751,MPSI1,MREV1,OPMOD1,OFFST1,MRMOD1,NONRT31,NCOLB1,CONE1,IB1,DEL11,XFA1,KPG1,GGS1,GQS1(5),GPS1(3),NCOLT1,WT041,RFA1,HSPRC1,HSPRS1,KPB1(5),RPH1,XPH1,WTCH1,WTSH1,WTRH1,WTC1,WTS15,WTR1,NUGC1,NUGS1,KFLAP1,KLAG1,HINGE1,KPIN1,ATNKP1(5),DROOP1,SWEEP61,DEL21,DEL31,ZFA1,VTIPN1,DEBUG1,EPMOD1,GLAG1,LDAMP1,EXDMP1,PHIPH17,PHIPL1,LDMPC1,MUY1,MUZ1,CHORD1(51),XA1(51),MASS1(51),XI1(51),THET8A1(51),ITHET1(51),TWIST1(51),EIZZ1(51),EIXX1(51),GJ1(51),XC1(51),D9ENSE1,OPUNT1,ELAG1,VTOT1,KHLMD1,KFLMD1,ALFHP1,PSIHP1,T75IN1,T1CIN1A,T1SIN1,IRING1,ROTAT1,FDROP1,FSWEP1,DEL41,DEL51
REAL MU1,LAMDA1,MTIP1,IB1,MASS1,ITHET1,KPG1,NUGC1,NUGS1,KPB1,KFLAP11,KLAG1,LDAMP1,LDMPC1,MUY1,MUZ1,KHLMD1,KFLMD1
INTEGER OPAER1,GIMBL1,OPMOD1,HINGE1,DEBUG1,OPUNT1,ROTAT1
COMMON /R2DATA/TYPE2,NBLAD2,DEL3G2,NTRM02,NTRM12,MU2,LAMDA2,MTIP2,1NBM2,NTM2,GIMBL2,GAMMA2,OMEGA2,RADUS2,SIGMA2,OPAER2,MRAER2,BBLD2,R2ROOT2,INCMP2,T1C2,T1S2,T752,MPSI2,MREV2,OPMOD2,OFFST2,MRMOD2,NONRT32,NCOLB2,CONE2,IB2,DEL12,XFA2,KPG2,GGS2,GQS2(5),GPS2(3),NCOLT2,WT042,RFA2,HSPRC2,HSPRS2,KPB2(5),RPH2,XPH2,WTCH2,WTSH2,WTRH2,WTC2,WTS25,WTR2,NUGC2,NUGS2,KFLAP2,KLAG2,HINGE2,KPIN2,ATNKP2(5),DROOP2,SWEEP62,DEL22,DEL32,ZFA2,VTIPN2,DEBUG2,EPMOD2,GLAG2,LDAMP2,EXDMP2,PHIPH27,PHIPL2,LDMPC2,MUY2,MUZ2,CHORD2(51),XA2(51),MASS2(51),XI2(51),THET8A2(51),ITHET2(51),TWIST2(51),EIZZ2(51),EIXX2(51),GJ2(51),XC2(51),D9ENSE2,OPUNT2,ELAG2,VTOT2,KHLMD2,KFLMD2,ALFHP2,PSIHP2,T75IN2,T1CIN2A,T1SIN2,IRING2,ROTAT2,FDROP2,FSWEP2,DEL42,DEL52
REAL MU2,LAMDA2,MTIP2,IB2,MASS2,ITHET2,KPG2,NUGC2,NUGS2,KPB2,KFLAP12,KLAG2,LDAMP2,LDMPC2,MUY2,MUZ2,KHLMD2,KFLMD2
INTEGER OPAER2,GIMBL2,OPMOD2,HINGE2,DEBUG2,OPUNT2,ROTAT2
COMMON /HINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,OPUNIT,ALTMSL,TEMP,ALTA1GL,ACLIMB,AYAW,APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,WEIGHT,AFLAP,-

2ITER,EPTRIM,FACTOR,PTRIM,DOF(80),CON(26),GUST(3),OPSYMM,DOFSYM(80)3,CONSYM(26),GUSSYM(3),OPFLOW,OPTORS(2),OPMOTR,OPLMDA(2),OPINTR,OPQ4RTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,JCNST,NSYS5AN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(80),CONP(29),CNTRL0,CNT6RLC,CNTRLS,CNTRLP,THETFP,PSIFP,THETAF,PHIF,VELX,VELY,VELZ,GRAV,DEN7SE,CSOUND,FREQ(100)
COMMON /GINPUT/ORIENT,ASHAFT(2),ACANT(2),K12H,K12F,K12L,K21H,K21F,1K21L,KLTRIM,KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,CWRTR1,CWRTR22,CHRTR1,CHRTR2,CVRTR1,CVRTR2,TAUW,TAUH,TAUV,RTRANS,IENG,KMAST1,KM3AST2,KICS,KENG,ENGPOS,KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,TH4RTLC,NRPOWR,EDAMPC,KPMC1(10),KPMS1(10),KPMC2(10),KPMS2(10),FSR1,WL5R1,BLR1,FSR2,WLR2,BLR2,FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,IXX6,IYY,IZZ,IXZ,XRTR1,YRTR1,ZRTR1,XRTR2,YRTR2,ZRTR2,XWB,ZWB,XHT,ZHT,X7VT,ZVT,AMASS,AIXX,AIYY,AIZZ,AIXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCF8E,PRCCFE,PFPCFE,PRPCFE,KF0CFE,KR0CFE,KFSCFE,KRSCFE,KCCFE,KSCFE,PCC9FE,PSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,CNTRLZ(11),TCFEA(11,5),QMASS(10),QFREQ(10),QDAMP(10),ZETAR1(3,10),ZETAR2(3,10),GAMBAR1(3,10),GAMAR2(3,10),NEM,RSF1(3,3),RSF2(3,3),RFV(3,3),RSV1(3,3),CRSV2(3,3),ZAGL1,ZAGL2,CSEGE1,CSEGE2,LMDA01,LMDA02,EDAMP,QTHRTL
INTEGER OPUNIT,ORIENT,DOF,CON,GUST,OPSYMM,DOFSYM,CONSYM,GUSSYM,OPF1LOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND2,OPAERO,DEBUG,TONLY,RESPON,DOFP,CONP,ENGPOS
REAL LNGCYC,LATCYC,K12H,K12F,K12L,K21H,K21F,K21L,KLTRIM,KWRTR1,KWR1TR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,IENG,KMAST1,KMAST2,KICS,KENG,KPR1G2V,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,NRPOWR,KPMC1,KPMC2,KPMS1,KPMS32,IXX,IYY,IZZ,IXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,KF0CFE,KR0CFE,KFSCFE4,KRSCFE,KCCFE,KSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,LMDA501,LMDA02
COMMON /FINPUT/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ(10,10),EXX(140,10),EZX(10,10),EXZ(10,10),EZ(10),EX(10),ET(10),EZT(10,10),EXT(105,10),ETT(10,10),EZD(10),EXD(10),EZL(10),EXL(10),ETL(10),FZH(10),FX6H(10),FTH(10),FYW(10),FFH(10),FYV(10),FFV(10),FPV(10)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
C
DATA TITLE/72H
MEDIUM UTILITY HELICOPTER
1
/
DATA OPUNIT/1/
DATA ORIENT/1/
DATA ALTMSL,TEMP,ALTAGL/0.,1000.,0./
DATA ASHAFT,ACANT/5.,0.,2*0./
DATA ACLIMB,AYAW/0.,0./
DATA APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,AFLAP/7*0./
DATA WEIGHT/8000./
DATA ITER,EPTRIM,FACTOR,PTRIM/10,.01,.5,0./
DATA DOF/3*0,1,20*0,2*1,28*0,2*1,5*0,6*1,13*0/
DATA CON,GUST/3*1,5*0,1,7*0,1,4*0,5*1,3*1/
DATA OPSYMM/0/
DATA DOFSYM,CONSYM,GUSSYM/61*0,-1,1,-1,1,-1,1,13*0,1,-1,1,5*0,-1,71*0,1,4*0,1,-1,1,-1,1,1,-1,1/
DATA OPFLOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREA1D,OPGRND,OPAERO,DEBUG,TONLY/-1,2*1,1,2*0,0,4,1,0,0,0,1,0,1/
DATA MPER,JCNST/24,12/
DATA NSYSAN,NSTEP,NFREQ,FREQ/0,0,0,100*0./
DATA LOCUS,PLOTMX,RESPON,NF,CONP,DOFP/-1,2.99,-1,100,29*0,80*0/
DATA K12H,K12F,K21H,K21F/4*0./
DATA KLTRIM/0./
DATA KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2/2*0.,1.8,2*0.,.7/
DATA RTRANS,IENG,KMAST1,KMAST2,KICS,KENG,ENGPOS/5.122,10.,4*3.E4,11/
DATA KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV/6*0./
DATA THRTLC,NRPOWR,EDAMPC/1.,1100.,1./

DATA KPMC1,KPMS1,KPMC2,KPMS2/40*0./
DATA FSR1,WLR1,BLR1,FSR2,WLR2,BLR2/11.09,11.81,0.,33.94,11.46,1.261/
DATA FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT/11.42,4.79,11.42,4.791,30.25,5.4,37.53,8.73/
DATA IXX,IYY,IZZ,IXZ/2925.,10830.,9250.,1250./
DATA KCCFE,KSCFE,PCCFE,PSCFE,K0CFE,KPCFE/2*1.,2*0.,2*1./
DATA KFCFE,KACFE,KECFE,KRCFE,KTCFE/5*0./
DATA CNTRLZ/11*0./
DATA QMASS,QFREQ,QDAMP/20*1.,10*.01/
DATA ZETAR1,GAMAR1,ZETAR2,GAMAR2/120*0./
DATA EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,1FXH,FTH,FYW,FFH,FYV,FFV,FPV/860*0./
DATA LFTAW,LFTDW,LFTFW,LFTEW/40.1,3*0./
DATA MOM0W,MOMAW,MOMDW,MOMFW/120.,-1170.,2*0./
DATA DRG0W,DRGDW,DRGFW,DRGEW,PIEL2/11.65,3*0.,77./
DATA AMAXW,DRGSW/20.,45./
DATA SIDEB,ROLLB,YAWB,KCLB/-107.,0.,-599.,0./
DATA IWB,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW/-5.,5*0./
DATA LFLAPO,LFLAPI,LAILO,LAILI/4*0./
DATA LFTAH,LFTEH,AMAXH,IHT,SPANH/55.,0.,20.,-7.71,9.33/
DATA LFTAV,LFTEV,AMAXV,IVT,SPANV/20.,0.,20.,0.,4./
C
DATA TYPE1/4HMAIN/
DATA OPAER1,OPMOD1,INCMP1/-1,4,1/
DATA MRAER1,MRMOD1,NONRT1,NCOLB1,NCOLT1/9,40,0,2,2/
DATA NTRM01,NTRM11,MPSI1,MREV1,EPMOD1/1,2,24,20,.001/
DATA VTIPN1,ROTAT1/814.,1/
DATA NBLAD1,GAMMA1,RADUS1,SIGMA1,IB1/3,6.42,24.,.0464,817.65/
DATA BBLD1,RROOT1/.97,.1/
DATA OFFST1,ELAG1/2*.01/
DATA GGS1,GQS1,GPS1/0.,5*.01,3*.01/
DATA WT01,WTCH1,WTSH1,WTRH1/5.8,5.3,5.3,6.5/
DATA CONE1,DROOP1,SWEEP1,FDROP1,FSWEP1/2.75,4*0./
DATA XFA1,ZFA1,RFA1/0.,.0181,.06/
DATA KPIN1,ATNKP1,DEL3G1/2,5*0.,0./
DATA RPH1,XPH1,PHIPL1,PHIPH1/.05,.05,0.,0./
DATA HINGE1,KFLAP1,KLAG1,GIMBL1,HSPRC1,HSPRS1/1,0.,0.,1,0.,0./
DATA LDAMP1,EXDMP1/0.,0./
DATA KHLMD1,KFLMD1/1.15,1.04/
DATA CHORD1,XA1,MASS1,ITHET1,EIZZ1,EIXX1,GJ1,XI1,XC1/51*.0486,51*01.,51*.177,51*.0081,51*225000.,51*6.7E6,51*78000.,51*0.,51*0./
DATA TWIST1/8.18,7.96,7.74,7.52,7.3,7.09,6.87,6.65,6.43,6.21,6.,5.178,5.56,5.34,5.12,4.91,4.69,4.47,4.25,4.03,3.82,3.6,3.38,3.16,2.942,2.73,2.51,2.29,2.07,1.85,1.64,1.42,1.2,.98,.76,.55,.33,.11,-.11,-3.33,-.55,-.76,-.98,-1.2,-1.42,-1.64,-1.85,-2.07,-2.29,-2.51,-2.72/
C
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA

TYPE2/4HTAIL/
OPAER2,OPMOD2,INCMP2/-1,4,1/
MRAER2,MRMOD2,NONRT2,NCOLB2,NCOLT2/9,40,0,2,2/
NTRM02,NTRM12,MPSI2,MREV2,EPMOD2/2,2,24,20,.001/
VTIPN2,ROTAT2/738.3,-1/
NBLAD2,GAMMA2,RADUS2,SIGMA2,IB2/3,3.07,4.25,.1436,.92/
BBLD2,RROOT2/.93,.1/
OFFST2,ELAG2/2*.01/
GGS2,GQS2,GPS2/0.,5*.01,3*.01/
WT02,WTCH2,WTSH2,WTRH2/5.8,5.3,5.3,6.5/
CONE2,DROOP2,SWEEP2,FDROP2,FSWEP2/1.5,4*0./
XFA2,ZFA2,RFA2/0.,0.,.15/
KPIN2,ATNKP2,DEL3G2/2,5*0.,45./
RPH2,XPH2,PHIPL2,PHIPH2/.05,.05,0.,0./
HINGE2,KFLAP2,KLAG2,GIMBL2,HSPRC2,HSPRS2/1,0.,0.,1,0.,0./
LDAMP2,EXDMP2/0.,0./
KHLMD2,KFLMD2/1.2,1.1/
CHORD2,XA2,MASS2,ITHET2,EIZZ2,EIXX2,GJ2,XI2,XC2/51*.1503,51*0-

1.,51*.036,51*.0005,51*1200.,51*35000.,51*4000.,51*0.,51*0./
DATA TWIST2/51*0./
END
BLOCK DATA
COMMON /R1DATA/TYPE1,NBLAD1,DEL3G1,NTRM01,NTRM11,MU1,LAMDA1,MTIP1,1NBM1,NTM1,GIMBL1,GAMMA1,OMEGA1,RADUS1,SIGMA1,OPAER1,MRAER1,BBLD1,R2ROOT1,INCMP1,T1C1,T1S1,T751,MPSI1,MREV1,OPMOD1,OFFST1,MRMOD1,NONRT31,NCOLB1,CONE1,IB1,DEL11,XFA1,KPG1,GGS1,GQS1(5),GPS1(3),NCOLT1,WT041,RFA1,HSPRC1,HSPRS1,KPB1(5),RPH1,XPH1,WTCH1,WTSH1,WTRH1,WTC1,WTS15,WTR1,NUGC1,NUGS1,KFLAP1,KLAG1,HINGE1,KPIN1,ATNKP1(5),DROOP1,SWEEP61,DEL21,DEL31,ZFA1,VTIPN1,DEBUG1,EPMOD1,GLAG1,LDAMP1,EXDMP1,PHIPH17,PHIPL1,LDMPC1,MUY1,MUZ1,CHORD1(51),XA1(51),MASS1(51),XI1(51),THET8A1(51),ITHET1(51),TWIST1(51),EIZZ1(51),EIXX1(51),GJ1(51),XC1(51),D9ENSE1,OPUNT1,ELAG1,VTOT1,KHLMD1,KFLMD1,ALFHP1,PSIHP1,T75IN1,T1CIN1A,T1SIN1,IRING1,ROTAT1,FDROP1,FSWEP1,DEL41,DEL51
REAL MU1,LAMDA1,MTIP1,IB1,MASS1,ITHET1,KPG1,NUGC1,NUGS1,KPB1,KFLAP11,KLAG1,LDAMP1,LDMPC1,MUY1,MUZ1,KHLMD1,KFLMD1
INTEGER OPAER1,GIMBL1,OPMOD1,HINGE1,DEBUG1,OPUNT1,ROTAT1
COMMON /R2DATA/TYPE2,NBLAD2,DEL3G2,NTRM02,NTRM12,MU2,LAMDA2,MTIP2,1NBM2,NTM2,GIMBL2,GAMMA2,OMEGA2,RADUS2,SIGMA2,OPAER2,MRAER2,BBLD2,R2ROOT2,INCMP2,T1C2,T1S2,T752,MPSI2,MREV2,OPMOD2,OFFST2,MRMOD2,NONRT32,NCOLB2,CONE2,IB2,DEL12,XFA2,KPG2,GGS2,GQS2(5),GPS2(3),NCOLT2,WT042,RFA2,HSPRC2,HSPRS2,KPB2(5),RPH2,XPH2,WTCH2,WTSH2,WTRH2,WTC2,WTS25,WTR2,NUGC2,NUGS2,KFLAP2,KLAG2,HINGE2,KPIN2,ATNKP2(5),DROOP2,SWEEP62,DEL22,DEL32,ZFA2,VTIPN2,DEBUG2,EPMOD2,GLAG2,LDAMP2,EXDMP2,PHIPH27,PHIPL2,LDMPC2,MUY2,MUZ2,CHORD2(51),XA2(51),MASS2(51),XI2(51),THET8A2(51),ITHET2(51),TWIST2(51),EIZZ2(51),EIXX2(51),GJ2(51),XC2(51),D9ENSE2,OPUNT2,ELAG2,VTOT2,KHLMD2,KFLMD2,ALFHP2,PSIHP2,T75IN2,T1CIN2A,T1SIN2,IRING2,ROTAT2,FDROP2,FSWEP2,DEL42,DEL52
REAL MU2,LAMDA2,MTIP2,IB2,MASS2,ITHET2,KPG2,NUGC2,NUGS2,KPB2,KFLAP12,KLAG2,LDAMP2,LDMPC2,MUY2,MUZ2,KHLMD2,KFLMD2
INTEGER OPAER2,GIMBL2,OPMOD2,HINGE2,DEBUG2,OPUNT2,ROTAT2
COMMON /HINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,OPUNIT,ALTMSL,TEMP,ALTA1GL,ACLIMB,AYAW,APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,WEIGHT,AFLAP,2ITER,EPTRIM,FACTOR,PTRIM,DOF(80),CON(26),GUST(3),OPSYMM,DOFSYM(80)3,CONSYM(26),GUSSYM(3),OPFLOW,OPTORS(2),OPMOTR,OPLMDA(2),OPINTR,OPQ4RTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,JCNST,NSYS5AN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(80),CONP(29),CNTRL0,CNT6RLC,CNTRLS,CNTRLP,THETFP,PSIFP,THETAF,PHIF,VELX,VELY,VELZ,GRAV,DEN7SE,CSOUND,FREQ(100)
COMMON /GINPUT/ORIENT,ASHAFT(2),ACANT(2),K12H,K12F,K12L,K21H,K21F,1K21L,KLTRIM,KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,CWRTR1,CWRTR22,CHRTR1,CHRTR2,CVRTR1,CVRTR2,TAUW,TAUH,TAUV,RTRANS,IENG,KMAST1,KM3AST2,KICS,KENG,ENGPOS,KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,TH4RTLC,NRPOWR,EDAMPC,KPMC1(10),KPMS1(10),KPMC2(10),KPMS2(10),FSR1,WL5R1,BLR1,FSR2,WLR2,BLR2,FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,IXX6,IYY,IZZ,IXZ,XRTR1,YRTR1,ZRTR1,XRTR2,YRTR2,ZRTR2,XWB,ZWB,XHT,ZHT,X7VT,ZVT,AMASS,AIXX,AIYY,AIZZ,AIXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCF8E,PRCCFE,PFPCFE,PRPCFE,KF0CFE,KR0CFE,KFSCFE,KRSCFE,KCCFE,KSCFE,PCC9FE,PSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,CNTRLZ(11),TCFEA(11,5),QMASS(10),QFREQ(10),QDAMP(10),ZETAR1(3,10),ZETAR2(3,10),GAMBAR1(3,10),GAMAR2(3,10),NEM,RSF1(3,3),RSF2(3,3),RFV(3,3),RSV1(3,3),CRSV2(3,3),ZAGL1,ZAGL2,CSEGE1,CSEGE2,LMDA01,LMDA02,EDAMP,QTHRTL
INTEGER OPUNIT,ORIENT,DOF,CON,GUST,OPSYMM,DOFSYM,CONSYM,GUSSYM,OPF1LOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND2,OPAERO,DEBUG,TONLY,RESPON,DOFP,CONP,ENGPOS
REAL LNGCYC,LATCYC,K12H,K12F,K12L,K21H,K21F,K21L,KLTRIM,KWRTR1,KWR1TR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,IENG,KMAST1,KMAST2,KICS,KENG,KPR1G2V,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,NRPOWR,KPMC1,KPMC2,KPMS1,KPMS32,IXX,IYY,IZZ,IXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,KF0CFE,KR0CFE,KFSCFE4,KRSCFE,KCCFE,KSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,LMDA501,LMDA02
COMMON /FINPUT/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,-

3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ(10,10),EXX(140,10),EZX(10,10),EXZ(10,10),EZ(10),EX(10),ET(10),EZT(10,10),EXT(105,10),ETT(10,10),EZD(10),EXD(10),EZL(10),EXL(10),ETL(10),FZH(10),FX6H(10),FTH(10),FYW(10),FFH(10),FYV(10),FFV(10),FPV(10)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
C
DATA TITLE/72H
LIGHT HELICOPTER
1
/
DATA OPUNIT/1/
DATA ORIENT/1/
DATA ALTMSL,TEMP,ALTAGL/0.,1000.,0./
DATA ASHAFT,ACANT/3.,0.,2*0./
DATA ACLIMB,AYAW/0.,0./
DATA APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,AFLAP/7*0./
DATA WEIGHT/2550./
DATA ITER,EPTRIM,FACTOR,PTRIM/10,.01,.5,0./
DATA DOF/4*0,4*1,46*0,2*1,5*0,6*1,13*0/
DATA CON,GUST/3*1,5*0,1,7*0,1,4*0,5*1,3*1/
DATA OPSYMM/0/
DATA DOFSYM,CONSYM,GUSSYM/61*0,-1,1,-1,1,-1,1,13*0,1,-1,1,5*0,-1,71*0,1,4*0,1,-1,1,-1,1,1,-1,1/
DATA OPFLOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREA1D,OPGRND,OPAERO,DEBUG,TONLY/-1,2*1,1,2*0,0,0,1,0,0,0,1,0,1/
DATA MPER,JCNST/24,12/
DATA NSYSAN,NSTEP,NFREQ,FREQ/0,0,0,100*0./
DATA LOCUS,PLOTMX,RESPON,NF,CONP,DOFP/-1,2.99,-1,100,29*0,80*0/
DATA K12H,K12F,K21H,K21F/4*0./
DATA KLTRIM/0./
DATA KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2/2*0.,1.8,2*0.,.7/
DATA RTRANS,IENG,KMAST1,KMAST2,KICS,KENG,ENGPOS/6.445,2.46,28000.,1100000.,584.,420000.,1/
DATA KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV/6*0./
DATA THRTLC,NRPOWR,EDAMPC/1.,320.,1./
DATA KPMC1,KPMS1,KPMC2,KPMS2/40*0./
DATA FSR1,WLR1,BLR1,FSR2,WLR2,BLR2/8.333,8.333,0.,23.5,5.942,-.97/
DATA FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT/8.38,4.133,8.38,4.1331,23.358,6.833,23.88,6.847/
DATA IXX,IYY,IZZ,IXZ/329.,894.,722.,94.5/
DATA KCCFE,KSCFE,PCCFE,PSCFE,K0CFE,KPCFE/2*1.,2*0.,2*1./
DATA KFCFE,KACFE,KECFE,KRCFE,KTCFE/5*0./
DATA CNTRLZ/11*0./
DATA QMASS,QFREQ,QDAMP/20*1.,10*.01/
DATA ZETAR1,GAMAR1,ZETAR2,GAMAR2/120*0./
DATA EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,1FXH,FTH,FYW,FFH,FYV,FFV,FPV/860*0./
DATA LFTAW,LFTDW,LFTFW,LFTEW/6.3,0.,0.,0./
DATA MOM0W,MOMAW,MOMDW,MOMFW/-3.,145.4,0.,0./
DATA DRG0W,DRGDW,DRGFW,DRGEW,PIEL2/5.4,3*0.,50./
DATA AMAXW,DRGSW/20.,45./
DATA SIDEB,ROLLB,YAWB,KCLB/-10.6,0.,-149.,0./
DATA IWB,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW/6*0./
DATA LFLAPO,LFLAPI,LAILO,LAILI/4*0./
DATA LFTAH,LFTEH,AMAXH,IHT,SPANH/19.0,0.,14.,.8,5.3/
DATA LFTAV,LFTEV,AMAXV,IVT,SPANV/22.0,0.,12.,0.,6.4/
C
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA
DATA

TYPE1/4HMAIN/
OPAER1,OPMOD1,INCMP1/-1,4,1/
MRAER1,MRMOD1,NONRT1,NCOLB1,NCOLT1/9,40,0,2,2/
NTRM01,NTRM11,MPSI1,MREV1,EPMOD1/1,2,24,20,.001/
VTIPN1,ROTAT1/646.,1/
NBLAD1,GAMMA1,RADUS1,SIGMA1,IB1/4,4.88,13.167,.0544,47./
BBLD1,RROOT1/.97,.1/
OFFST1,ELAG1/.035,.09/
GGS1,GQS1,GPS1/0.,5*.01,3*.01/

DATA WT01,WTCH1,WTSH1,WTRH1/5.8,5.3,5.3,6.5/
DATA CONE1,DROOP1,SWEEP1,FDROP1,FSWEP1/5*0./
DATA XFA1,ZFA1,RFA1/0.,0.,.11/
DATA KPIN1,ATNKP1,DEL3G1/2,5*0.,0./
DATA RPH1,XPH1,PHIPL1,PHIPH1/.05,.05,0.,0./
DATA HINGE1,KFLAP1,KLAG1,GIMBL1,HSPRC1,HSPRS1/0,0.,0.,0,0.,0./
DATA LDAMP1,EXDMP1/115.,0./
DATA KHLMD1,KFLMD1/1.15,1.04/
DATA CHORD1,XA1,MASS1,ITHET1,EIZZ1,EIXX1,GJ1,XI1,XC1/51*.0427,51*01.,51*.062,51*.0007,51*15000.,51*450000.,51*4300.,51*0.,51*0./
DATA TWIST1/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
C
DATA TYPE2/4HTAIL/
DATA OPAER2,OPMOD2,INCMP2/-1,4,1/
DATA MRAER2,MRMOD2,NONRT2,NCOLB2,NCOLT2/9,40,0,2,2/
DATA NTRM02,NTRM12,MPSI2,MREV2,EPMOD2/2,2,24,20,.001/
DATA VTIPN2,ROTAT2/671.94,1/
DATA NBLAD2,GAMMA2,RADUS2,SIGMA2,IB2/3,1.11,2.125,.12,.0667/
DATA BBLD2,RROOT2/.93,.1/
DATA OFFST2,ELAG2/0.,.1/
DATA GGS2,GQS2,GPS2/0.,5*.01,3*.01/
DATA WT02,WTCH2,WTSH2,WTRH2/5.8,5.3,5.3,6.5/
DATA CONE2,DROOP2,SWEEP2,FDROP2,FSWEP2/5*0./
DATA XFA2,ZFA2,RFA2/0.,0.,.15/
DATA KPIN2,ATNKP2,DEL3G2/2,5*0.,30./
DATA RPH2,XPH2,PHIPL2,PHIPH2/.05,.05,0.,0./
DATA HINGE2,KFLAP2,KLAG2,GIMBL2,HSPRC2,HSPRS2/1,0.,0.,1,0.,0./
DATA LDAMP2,EXDMP2/0.,0./
DATA KHLMD2,KFLMD2/1.2,1.1/
DATA CHORD2,XA2,MASS2,ITHET2,EIZZ2,EIXX2,GJ2,XI2,XC2/51*.1255,51*01.,51*.0208,51*.00005,51*140.,51*4200.,51*330.,51*0.,51*0./
DATA TWIST2/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
END
BLOCK DATA
COMMON /R1DATA/TYPE1,NBLAD1,DEL3G1,NTRM01,NTRM11,MU1,LAMDA1,MTIP1,1NBM1,NTM1,GIMBL1,GAMMA1,OMEGA1,RADUS1,SIGMA1,OPAER1,MRAER1,BBLD1,R2ROOT1,INCMP1,T1C1,T1S1,T751,MPSI1,MREV1,OPMOD1,OFFST1,MRMOD1,NONRT31,NCOLB1,CONE1,IB1,DEL11,XFA1,KPG1,GGS1,GQS1(5),GPS1(3),NCOLT1,WT041,RFA1,HSPRC1,HSPRS1,KPB1(5),RPH1,XPH1,WTCH1,WTSH1,WTRH1,WTC1,WTS15,WTR1,NUGC1,NUGS1,KFLAP1,KLAG1,HINGE1,KPIN1,ATNKP1(5),DROOP1,SWEEP61,DEL21,DEL31,ZFA1,VTIPN1,DEBUG1,EPMOD1,GLAG1,LDAMP1,EXDMP1,PHIPH17,PHIPL1,LDMPC1,MUY1,MUZ1,CHORD1(51),XA1(51),MASS1(51),XI1(51),THET8A1(51),ITHET1(51),TWIST1(51),EIZZ1(51),EIXX1(51),GJ1(51),XC1(51),D9ENSE1,OPUNT1,ELAG1,VTOT1,KHLMD1,KFLMD1,ALFHP1,PSIHP1,T75IN1,T1CIN1A,T1SIN1,IRING1,ROTAT1,FDROP1,FSWEP1,DEL41,DEL51
REAL MU1,LAMDA1,MTIP1,IB1,MASS1,ITHET1,KPG1,NUGC1,NUGS1,KPB1,KFLAP11,KLAG1,LDAMP1,LDMPC1,MUY1,MUZ1,KHLMD1,KFLMD1
INTEGER OPAER1,GIMBL1,OPMOD1,HINGE1,DEBUG1,OPUNT1,ROTAT1
COMMON /R2DATA/TYPE2,NBLAD2,DEL3G2,NTRM02,NTRM12,MU2,LAMDA2,MTIP2,1NBM2,NTM2,GIMBL2,GAMMA2,OMEGA2,RADUS2,SIGMA2,OPAER2,MRAER2,BBLD2,R2ROOT2,INCMP2,T1C2,T1S2,T752,MPSI2,MREV2,OPMOD2,OFFST2,MRMOD2,NONRT32,NCOLB2,CONE2,IB2,DEL12,XFA2,KPG2,GGS2,GQS2(5),GPS2(3),NCOLT2,WT042,RFA2,HSPRC2,HSPRS2,KPB2(5),RPH2,XPH2,WTCH2,WTSH2,WTRH2,WTC2,WTS25,WTR2,NUGC2,NUGS2,KFLAP2,KLAG2,HINGE2,KPIN2,ATNKP2(5),DROOP2,SWEEP62,DEL22,DEL32,ZFA2,VTIPN2,DEBUG2,EPMOD2,GLAG2,LDAMP2,EXDMP2,PHIPH27,PHIPL2,LDMPC2,MUY2,MUZ2,CHORD2(51),XA2(51),MASS2(51),XI2(51),THET8A2(51),ITHET2(51),TWIST2(51),EIZZ2(51),EIXX2(51),GJ2(51),XC2(51),D-

9ENSE2,OPUNT2,ELAG2,VTOT2,KHLMD2,KFLMD2,ALFHP2,PSIHP2,T75IN2,T1CIN2A,T1SIN2,IRING2,ROTAT2,FDROP2,FSWEP2,DEL42,DEL52
REAL MU2,LAMDA2,MTIP2,IB2,MASS2,ITHET2,KPG2,NUGC2,NUGS2,KPB2,KFLAP12,KLAG2,LDAMP2,LDMPC2,MUY2,MUZ2,KHLMD2,KFLMD2
INTEGER OPAER2,GIMBL2,OPMOD2,HINGE2,DEBUG2,OPUNT2,ROTAT2
COMMON /HINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,OPUNIT,ALTMSL,TEMP,ALTA1GL,ACLIMB,AYAW,APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,WEIGHT,AFLAP,2ITER,EPTRIM,FACTOR,PTRIM,DOF(80),CON(26),GUST(3),OPSYMM,DOFSYM(80)3,CONSYM(26),GUSSYM(3),OPFLOW,OPTORS(2),OPMOTR,OPLMDA(2),OPINTR,OPQ4RTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,JCNST,NSYS5AN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(80),CONP(29),CNTRL0,CNT6RLC,CNTRLS,CNTRLP,THETFP,PSIFP,THETAF,PHIF,VELX,VELY,VELZ,GRAV,DEN7SE,CSOUND,FREQ(100)
COMMON /GINPUT/ORIENT,ASHAFT(2),ACANT(2),K12H,K12F,K12L,K21H,K21F,1K21L,KLTRIM,KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,CWRTR1,CWRTR22,CHRTR1,CHRTR2,CVRTR1,CVRTR2,TAUW,TAUH,TAUV,RTRANS,IENG,KMAST1,KM3AST2,KICS,KENG,ENGPOS,KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,TH4RTLC,NRPOWR,EDAMPC,KPMC1(10),KPMS1(10),KPMC2(10),KPMS2(10),FSR1,WL5R1,BLR1,FSR2,WLR2,BLR2,FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,IXX6,IYY,IZZ,IXZ,XRTR1,YRTR1,ZRTR1,XRTR2,YRTR2,ZRTR2,XWB,ZWB,XHT,ZHT,X7VT,ZVT,AMASS,AIXX,AIYY,AIZZ,AIXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCF8E,PRCCFE,PFPCFE,PRPCFE,KF0CFE,KR0CFE,KFSCFE,KRSCFE,KCCFE,KSCFE,PCC9FE,PSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,CNTRLZ(11),TCFEA(11,5),QMASS(10),QFREQ(10),QDAMP(10),ZETAR1(3,10),ZETAR2(3,10),GAMBAR1(3,10),GAMAR2(3,10),NEM,RSF1(3,3),RSF2(3,3),RFV(3,3),RSV1(3,3),CRSV2(3,3),ZAGL1,ZAGL2,CSEGE1,CSEGE2,LMDA01,LMDA02,EDAMP,QTHRTL
INTEGER OPUNIT,ORIENT,DOF,CON,GUST,OPSYMM,DOFSYM,CONSYM,GUSSYM,OPF1LOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND2,OPAERO,DEBUG,TONLY,RESPON,DOFP,CONP,ENGPOS
REAL LNGCYC,LATCYC,K12H,K12F,K12L,K21H,K21F,K21L,KLTRIM,KWRTR1,KWR1TR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,IENG,KMAST1,KMAST2,KICS,KENG,KPR1G2V,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,NRPOWR,KPMC1,KPMC2,KPMS1,KPMS32,IXX,IYY,IZZ,IXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,KF0CFE,KR0CFE,KFSCFE4,KRSCFE,KCCFE,KSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,LMDA501,LMDA02
COMMON /FINPUT/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ(10,10),EXX(140,10),EZX(10,10),EXZ(10,10),EZ(10),EX(10),ET(10),EZT(10,10),EXT(105,10),ETT(10,10),EZD(10),EXD(10),EZL(10),EXL(10),ETL(10),FZH(10),FX6H(10),FTH(10),FYW(10),FFH(10),FYV(10),FFV(10),FPV(10)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
C
DATA TITLE/72H
SOFT-INPLANE HINGELESS ROTOR HELICOPTER (BO-1105)
/
DATA OPUNIT/1/
DATA ORIENT/1/
DATA ALTMSL,TEMP,ALTAGL/0.,1000.,0./
DATA ASHAFT,ACANT/3.,-4.,2*0./
DATA ACLIMB,AYAW/0.,0./
DATA APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,AFLAP/7*0./
DATA WEIGHT/4640./
DATA ITER,EPTRIM,FACTOR,PTRIM/10,.01,.5,0./
DATA DOF/8*1,7*0,8*1,4*0,1,2*0,24*0,2*1,0,1,3*0,6*1,13*0/
DATA CON,GUST/3*1,5*0,1,7*0,1,4*0,5*1,3*1/
DATA OPSYMM/0/
DATA DOFSYM,CONSYM,GUSSYM/61*0,-1,1,-1,1,-1,1,13*0,1,-1,1,5*0,-1,71*0,1,4*0,1,-1,1,-1,1,1,-1,1/
DATA OPFLOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREA1D,OPGRND,OPAERO,DEBUG,TONLY/-1,0,1,1,2*0,0,2,1,0,0,0,1,0,1/
DATA MPER,JCNST/24,12/
DATA NSYSAN,NSTEP,NFREQ,FREQ/0,0,0,100*0./
DATA LOCUS,PLOTMX,RESPON,NF,CONP,DOFP/-1,2.99,-1,100,29*0,80*0/

DATA K12H,K12F,K21H,K21F/4*0./
DATA KLTRIM/0./
DATA KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2/2*0.,1.8,2*0.,.7/
DATA RTRANS,IENG,KMAST1,KMAST2,KICS,KENG,ENGPOS/5.527,4.,4*3.E4,1/
DATA KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV/6*0./
DATA THRTLC,NRPOWR,EDAMPC/1.,540.,1./
DATA KPMC1,KPMS1,KPMC2,KPMS2/40*0./
DATA FSR1,WLR1,BLR1,FSR2,WLR2,BLR2/8.2,5.1,0.,27.92,5.73,-1.04/
DATA FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT/8.37,.58,8.37,-.16,231.12,2.15,25.52,3.94/
DATA IXX,IYY,IZZ,IXZ/1330.,3600.,3260.,250./
DATA KCCFE,KSCFE,PCCFE,PSCFE,K0CFE,KPCFE/2*1.,2*10.,2*1./
DATA KFCFE,KACFE,KECFE,KRCFE,KTCFE/5*0./
DATA CNTRLZ/11*0./
DATA QMASS,QFREQ,QDAMP/20*1.,10*.01/
DATA ZETAR1,GAMAR1,ZETAR2,GAMAR2/120*0./
DATA EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,1FXH,FTH,FYW,FFH,FYV,FFV,FPV/860*0./
DATA LFTAW,LFTDW,LFTFW,LFTEW/32.6,3*0./
DATA MOM0W,MOMAW,MOMDW,MOMFW/-29.3,90.28,2*0./
DATA DRG0W,DRGDW,DRGFW,DRGEW,PIEL2/9.1,3*0.,75./
DATA AMAXW,DRGSW/20.,45./
DATA SIDEB,ROLLB,YAWB,KCLB/-47.,97.,-57.,0./
DATA IWB,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW/-4.94,5*0./
DATA LFLAPO,LFLAPI,LAILO,LAILI/4*0./
DATA LFTAH,LFTEH,AMAXH,IHT,SPANH/32.,0.,15.,0.,8.4/
DATA LFTAV,LFTEV,AMAXV,IVT,SPANV/30.,0.,20.,0.,6./
C
DATA TYPE1/4HMAIN/
DATA OPAER1,OPMOD1,INCMP1/1,1,1/
DATA MRAER1,MRMOD1,NONRT1,NCOLB1,NCOLT1/9,40,0,2,2/
DATA NTRM01,NTRM11,MPSI1,MREV1,EPMOD1/1,2,24,20,.001/
DATA VTIPN1,ROTAT1/715.,1/
DATA NBLAD1,GAMMA1,RADUS1,SIGMA1,IB1/4,5.,16.11,.0701,161.9/
DATA BBLD1,RROOT1/.97,.1/
DATA OFFST1,ELAG1/2*.04/
DATA GGS1,GQS1,GPS1/0.,5*.06,3*.06/
DATA WT01,WTCH1,WTSH1,WTRH1/5.8,5.3,5.3,6.5/
DATA CONE1,DROOP1,SWEEP1,FDROP1,FSWEP1/2.5,4*0./
DATA XFA1,ZFA1,RFA1/0.,0.,.035/
DATA KPIN1,ATNKP1,DEL3G1/1,5*0.,0./
DATA RPH1,XPH1,PHIPL1,PHIPH1/.0345,.0335,0.,0./
DATA HINGE1,KFLAP1,KLAG1,GIMBL1,HSPRC1,HSPRS1/1,0.,0.,0,0.,0./
DATA LDAMP1,EXDMP1/0.,0./
DATA KHLMD1,KFLMD1/1.15,1.04/
DATA CHORD1,XA1/51*.0550,51*0./
DATA XI1/0.,.00347,2*.00185,3*.00044,2*-.00307,3*.00307,39*.00006/
DATA MASS1/3*1.25,3*.778,2*.361,5*.1084,38*.1153/
DATA EIZZ1/5.9E6,1.396E6,2*9.44E4,3*2.54E4,2*2.24E4,3*1.82E4,39*1.165E4/
DATA EIXX1/6.14E6,1.417E6,2*4.76E5,3*2.65E5,2*2.64E5,3*4.02E5,39*41.125E5/
DATA GJ1/2.85E4,2.85E4,2*2.84E4,3*2.23E4,2*1.49E4,3*1.13E4,39*1.061E4/
DATA XC1/4*0.,3*-.00238,2*-.00057,3*-.0038,39*-.00385/
DATA ITHET1/51*.00501/
DATA TWIST1/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
C
DATA TYPE2/4HTAIL/
DATA OPAER2,OPMOD2,INCMP2/-1,4,1/
DATA MRAER2,MRMOD2,NONRT2,NCOLB2,NCOLT2/9,40,0,2,2/

DATA NTRM02,NTRM12,MPSI2,MREV2,EPMOD2/2,2,24,20,.001/
DATA VTIPN2,ROTAT2/764.6,1/
DATA NBLAD2,GAMMA2,RADUS2,SIGMA2,IB2/3,1.544,3.117,.12,.325/
DATA BBLD2,RROOT2/.93,.1/
DATA OFFST2,ELAG2/0.,.1/
DATA GGS2,GQS2,GPS2/0.,5*.01,3*.01/
DATA WT02,WTCH2,WTSH2,WTRH2/5.8,5.3,5.3,6.5/
DATA CONE2,DROOP2,SWEEP2,FDROP2,FSWEP2/5*0./
DATA XFA2,ZFA2,RFA2/0.,0.,.15/
DATA KPIN2,ATNKP2,DEL3G2/2,5*0.,45./
DATA RPH2,XPH2,PHIPL2,PHIPH2/.05,.05,0.,0./
DATA HINGE2,KFLAP2,KLAG2,GIMBL2,HSPRC2,HSPRS2/1,0.,0.,1,0.,0./
DATA LDAMP2,EXDMP2/0.,0./
DATA KHLMD2,KFLMD2/1.2,1.1/
DATA CHORD2,XA2,MASS2,ITHET2,EIZZ2,EIXX2,GJ2,XI2,XC2/51*.1257,51*01.,51*.0322,51*.00017,51*610.,51*18000.,51*1400.,51*0.,51*0./
DATA TWIST2/51*0./
END
!!! ERASE SOURCE.BO105
BLOCK DATA
COMMON /R1DATA/TYPE1,NBLAD1,DEL3G1,NTRM01,NTRM11,MU1,LAMDA1,MTIP1,1NBM1,NTM1,GIMBL1,GAMMA1,OMEGA1,RADUS1,SIGMA1,OPAER1,MRAER1,BBLD1,R2ROOT1,INCMP1,T1C1,T1S1,T751,MPSI1,MREV1,OPMOD1,OFFST1,MRMOD1,NONRT31,NCOLB1,CONE1,IB1,DEL11,XFA1,KPG1,GGS1,GQS1(5),GPS1(3),NCOLT1,WT041,RFA1,HSPRC1,HSPRS1,KPB1(5),RPH1,XPH1,WTCH1,WTSH1,WTRH1,WTC1,WTS15,WTR1,NUGC1,NUGS1,KFLAP1,KLAG1,HINGE1,KPIN1,ATNKP1(5),DROOP1,SWEEP61,DEL21,DEL31,ZFA1,VTIPN1,DEBUG1,EPMOD1,GLAG1,LDAMP1,EXDMP1,PHIPH17,PHIPL1,LDMPC1,MUY1,MUZ1,CHORD1(51),XA1(51),MASS1(51),XI1(51),THET8A1(51),ITHET1(51),TWIST1(51),EIZZ1(51),EIXX1(51),GJ1(51),XC1(51),D9ENSE1,OPUNT1,ELAG1,VTOT1,KHLMD1,KFLMD1,ALFHP1,PSIHP1,T75IN1,T1CIN1A,T1SIN1,IRING1,ROTAT1,FDROP1,FSWEP1,DEL41,DEL51
REAL MU1,LAMDA1,MTIP1,IB1,MASS1,ITHET1,KPG1,NUGC1,NUGS1,KPB1,KFLAP11,KLAG1,LDAMP1,LDMPC1,MUY1,MUZ1,KHLMD1,KFLMD1
INTEGER OPAER1,GIMBL1,OPMOD1,HINGE1,DEBUG1,OPUNT1,ROTAT1
COMMON /R2DATA/TYPE2,NBLAD2,DEL3G2,NTRM02,NTRM12,MU2,LAMDA2,MTIP2,1NBM2,NTM2,GIMBL2,GAMMA2,OMEGA2,RADUS2,SIGMA2,OPAER2,MRAER2,BBLD2,R2ROOT2,INCMP2,T1C2,T1S2,T752,MPSI2,MREV2,OPMOD2,OFFST2,MRMOD2,NONRT32,NCOLB2,CONE2,IB2,DEL12,XFA2,KPG2,GGS2,GQS2(5),GPS2(3),NCOLT2,WT042,RFA2,HSPRC2,HSPRS2,KPB2(5),RPH2,XPH2,WTCH2,WTSH2,WTRH2,WTC2,WTS25,WTR2,NUGC2,NUGS2,KFLAP2,KLAG2,HINGE2,KPIN2,ATNKP2(5),DROOP2,SWEEP62,DEL22,DEL32,ZFA2,VTIPN2,DEBUG2,EPMOD2,GLAG2,LDAMP2,EXDMP2,PHIPH27,PHIPL2,LDMPC2,MUY2,MUZ2,CHORD2(51),XA2(51),MASS2(51),XI2(51),THET8A2(51),ITHET2(51),TWIST2(51),EIZZ2(51),EIXX2(51),GJ2(51),XC2(51),D9ENSE2,OPUNT2,ELAG2,VTOT2,KHLMD2,KFLMD2,ALFHP2,PSIHP2,T75IN2,T1CIN2A,T1SIN2,IRING2,ROTAT2,FDROP2,FSWEP2,DEL42,DEL52
REAL MU2,LAMDA2,MTIP2,IB2,MASS2,ITHET2,KPG2,NUGC2,NUGS2,KPB2,KFLAP12,KLAG2,LDAMP2,LDMPC2,MUY2,MUZ2,KHLMD2,KFLMD2
INTEGER OPAER2,GIMBL2,OPMOD2,HINGE2,DEBUG2,OPUNT2,ROTAT2
COMMON /HINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,OPUNIT,ALTMSL,TEMP,ALTA1GL,ACLIMB,AYAW,APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,WEIGHT,AFLAP,2ITER,EPTRIM,FACTOR,PTRIM,DOF(80),CON(26),GUST(3),OPSYMM,DOFSYM(80)3,CONSYM(26),GUSSYM(3),OPFLOW,OPTORS(2),OPMOTR,OPLMDA(2),OPINTR,OPQ4RTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,JCNST,NSYS5AN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(80),CONP(29),CNTRL0,CNT6RLC,CNTRLS,CNTRLP,THETFP,PSIFP,THETAF,PHIF,VELX,VELY,VELZ,GRAV,DEN7SE,CSOUND,FREQ(100)
COMMON /GINPUT/ORIENT,ASHAFT(2),ACANT(2),K12H,K12F,K12L,K21H,K21F,1K21L,KLTRIM,KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,CWRTR1,CWRTR22,CHRTR1,CHRTR2,CVRTR1,CVRTR2,TAUW,TAUH,TAUV,RTRANS,IENG,KMAST1,KM3AST2,KICS,KENG,ENGPOS,KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,TH4RTLC,NRPOWR,EDAMPC,KPMC1(10),KPMS1(10),KPMC2(10),KPMS2(10),FSR1,WL5R1,BLR1,FSR2,WLR2,BLR2,FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,IXX6,IYY,IZZ,IXZ,XRTR1,YRTR1,ZRTR1,XRTR2,YRTR2,ZRTR2,XWB,ZWB,XHT,ZHT,X7VT,ZVT,AMASS,AIXX,AIYY,AIZZ,AIXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCF8E,PRCCFE,PFPCFE,PRPCFE,KF0CFE,KR0CFE,KFSCFE,KRSCFE,KCCFE,KSCFE,PCC-

9FE,PSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,CNTRLZ(11),TCFEA(11,5),QMASS(10),QFREQ(10),QDAMP(10),ZETAR1(3,10),ZETAR2(3,10),GAMBAR1(3,10),GAMAR2(3,10),NEM,RSF1(3,3),RSF2(3,3),RFV(3,3),RSV1(3,3),CRSV2(3,3),ZAGL1,ZAGL2,CSEGE1,CSEGE2,LMDA01,LMDA02,EDAMP,QTHRTL
INTEGER OPUNIT,ORIENT,DOF,CON,GUST,OPSYMM,DOFSYM,CONSYM,GUSSYM,OPF1LOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND2,OPAERO,DEBUG,TONLY,RESPON,DOFP,CONP,ENGPOS
REAL LNGCYC,LATCYC,K12H,K12F,K12L,K21H,K21F,K21L,KLTRIM,KWRTR1,KWR1TR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,IENG,KMAST1,KMAST2,KICS,KENG,KPR1G2V,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,NRPOWR,KPMC1,KPMC2,KPMS1,KPMS32,IXX,IYY,IZZ,IXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,KF0CFE,KR0CFE,KFSCFE4,KRSCFE,KCCFE,KSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,LMDA501,LMDA02
COMMON /FINPUT/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ(10,10),EXX(140,10),EZX(10,10),EXZ(10,10),EZ(10),EX(10),ET(10),EZT(10,10),EXT(105,10),ETT(10,10),EZD(10),EXD(10),EZL(10),EXL(10),ETL(10),FZH(10),FX6H(10),FTH(10),FYW(10),FFH(10),FYV(10),FFV(10),FPV(10)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
C
DATA TITLE/72H
TANDEM-ROTOR HELICOPTER
1
/
DATA OPUNIT/1/
DATA ORIENT/2/
DATA ALTMSL,TEMP,ALTAGL/0.,1000.,0./
DATA ASHAFT,ACANT/9.5,7.,2*0./
DATA ACLIMB,AYAW/0.,0./
DATA APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,AFLAP/7*0./
DATA WEIGHT/19000./
DATA ITER,EPTRIM,FACTOR,PTRIM/10,.01,.5,0./
DATA DOF/3*0,3*1,21*0,1,5*0,3*1,21*0,1,3*0,6*1,13*0/
DATA CON,GUST/3*1,5*0,3*1,5*0,1,4*0,5*1,3*1/
DATA OPSYMM/0/
DATA DOFSYM,CONSYM,GUSSYM/61*0,-1,1,-1,1,-1,1,13*0,1,-1,1,5*0,1,-11,1,5*0,1,4*0,1,-1,1,-1,1,1,-1,1/
DATA OPFLOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREA1D,OPGRND,OPAERO,DEBUG,TONLY/-1,2*1,1,2*0,0,0,1,0,0,0,1,0,1/
DATA MPER,JCNST/24,12/
DATA NSYSAN,NSTEP,NFREQ,FREQ/0,0,0,100*0./
DATA LOCUS,PLOTMX,RESPON,NF,CONP,DOFP/-1,2.99,-1,100,29*0,80*0/
DATA K12H,K12F,K21H,K21F/.062,0.,.062,1.76/
DATA KLTRIM/0./
DATA KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2/6*0./
DATA RTRANS,IENG,KMAST1,KMAST2,KICS,KENG,ENGPOS/1.,10.,4*3.E4,2/
DATA KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV/6*0./
DATA THRTLC,NRPOWR,EDAMPC/1.,2500.,1./
DATA KPMC1,KPMS1,KPMC2,KPMS2/40*0./
DATA FSR1,WLR1,BLR1,FSR2,WLR2,BLR2/-17.09,6.24,0.,15.93,10.93,0./
DATA FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT/8*0./
DATA IXX,IYY,IZZ,IXZ/14000.,120000.,110000.,0./
DATA KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCFE,PRCCFE,PFPCFE,PRPCFE,KF0CF1E,KR0CFE,KFSCFE,KRSCFE/1.09,.91,2*1.,4*0.,4*1./
DATA KFCFE,KACFE,KECFE,KRCFE,KTCFE/5*0./
DATA CNTRLZ/11*0./
DATA QMASS,QFREQ,QDAMP/20*1.,10*.01/
DATA ZETAR1,GAMAR1,ZETAR2,GAMAR2/120*0./
DATA EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,1FXH,FTH,FYW,FFH,FYV,FFV,FPV/860*0./
DATA LFTAW,LFTDW,LFTFW,LFTEW/253.,3*0./
DATA MOM0W,MOMAW,MOMDW,MOMFW/-100.,1150.,2*0./
DATA DRG0W,DRGDW,DRGFW,DRGEW,PIEL2/37.,3*0.,190./
DATA AMAXW,DRGSW/20.,90./

DATA
DATA
DATA
DATA
DATA

SIDEB,ROLLB,YAWB,KCLB/-382.,-1063.,-573.,0./
IWB,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW/.9,5*0./
LFLAPO,LFLAPI,LAILO,LAILI/4*0./
LFTAH,LFTEH,AMAXH,IHT,SPANH/5*0./
LFTAV,LFTEV,AMAXV,IVT,SPANV/5*0./

C
DATA TYPE1/4HFRNT/
DATA OPAER1,OPMOD1,INCMP1/-1,4,1/
DATA MRAER1,MRMOD1,NONRT1,NCOLB1,NCOLT1/8,40,0,2,2/
DATA NTRM01,NTRM11,MPSI1,MREV1,EPMOD1/1,2,24,20,.001/
DATA VTIPN1,ROTAT1/690.,-1/
DATA NBLAD1,GAMMA1,RADUS1,SIGMA1,IB1/3,8.66,25.,.0573,895.9/
DATA BBLD1,RROOT1/.97,.2/
DATA OFFST1,ELAG1/2*.017/
DATA GGS1,GQS1,GPS1/0.,5*.01,3*.01/
DATA WT01,WTCH1,WTSH1,WTRH1/5.8,5.3,5.3,6.5/
DATA CONE1,DROOP1,SWEEP1,FDROP1,FSWEP1/5*0./
DATA XFA1,ZFA1,RFA1/0.,0.,.05/
DATA KPIN1,ATNKP1,DEL3G1/2,5*0.,0./
DATA RPH1,XPH1,PHIPL1,PHIPH1/.05,.05,0.,0./
DATA HINGE1,KFLAP1,KLAG1,GIMBL1,HSPRC1,HSPRS1/0,0.,0.,0,0.,0./
DATA LDAMP1,EXDMP1/1200.,0./
DATA KHLMD1,KFLMD1/1.15,1.04/
DATA CHORD1,XA1,MASS1,ITHET1,EIZZ1,EIXX1,GJ1,XI1,XC1/51*.0600,51*01.,51*.172,51*.0131,51*170000.,51*5.E6,51*91000.,51*0.,51*0./
DATA TWIST1/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
C
DATA TYPE2/4HREAR/
DATA OPAER2,OPMOD2,INCMP2/-1,4,1/
DATA MRAER2,MRMOD2,NONRT2,NCOLB2,NCOLT2/8,40,0,2,2/
DATA NTRM02,NTRM12,MPSI2,MREV2,EPMOD2/1,2,24,20,.001/
DATA VTIPN2,ROTAT2/690.,1/
DATA NBLAD2,GAMMA2,RADUS2,SIGMA2,IB2/3,8.66,25.,.0573,895.9/
DATA BBLD2,RROOT2/.97,.2/
DATA OFFST2,ELAG2/2*.017/
DATA GGS2,GQS2,GPS2/0.,5*.01,3*.01/
DATA WT02,WTCH2,WTSH2,WTRH2/5.8,5.3,5.3,6.5/
DATA CONE2,DROOP2,SWEEP2,FDROP2,FSWEP2/5*0./
DATA XFA2,ZFA2,RFA2/0.,0.,.05/
DATA KPIN2,ATNKP2,DEL3G2/2,5*0.,0./
DATA RPH2,XPH2,PHIPL2,PHIPH2/.05,.05,0.,0./
DATA HINGE2,KFLAP2,KLAG2,GIMBL2,HSPRC2,HSPRS2/0,0.,0.,0,0.,0./
DATA LDAMP2,EXDMP2/1200.,0./
DATA KHLMD2,KFLMD2/1.15,1.04/
DATA CHORD2,XA2,MASS2,ITHET2,EIZZ2,EIXX2,GJ2,XI2,XC2/51*.0600,51*01.,51*.172,51*.0131,51*170000.,51*5.E6,51*91000.,51*0.,51*0./
DATA TWIST2/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
END
BLOCK DATA
COMMON /R1DATA/TYPE1,NBLAD1,DEL3G1,NTRM01,NTRM11,MU1,LAMDA1,MTIP1,1NBM1,NTM1,GIMBL1,GAMMA1,OMEGA1,RADUS1,SIGMA1,OPAER1,MRAER1,BBLD1,R2ROOT1,INCMP1,T1C1,T1S1,T751,MPSI1,MREV1,OPMOD1,OFFST1,MRMOD1,NONRT31,NCOLB1,CONE1,IB1,DEL11,XFA1,KPG1,GGS1,GQS1(5),GPS1(3),NCOLT1,WT041,RFA1,HSPRC1,HSPRS1,KPB1(5),RPH1,XPH1,WTCH1,WTSH1,WTRH1,WTC1,WTS15,WTR1,NUGC1,NUGS1,KFLAP1,KLAG1,HINGE1,KPIN1,ATNKP1(5),DROOP1,SWEEP61,DEL21,DEL31,ZFA1,VTIPN1,DEBUG1,EPMOD1,GLAG1,LDAMP1,EXDMP1,PHIPH17,PHIPL1,LDMPC1,MUY1,MUZ1,CHORD1(51),XA1(51),MASS1(51),XI1(51),THET-

8A1(51),ITHET1(51),TWIST1(51),EIZZ1(51),EIXX1(51),GJ1(51),XC1(51),D9ENSE1,OPUNT1,ELAG1,VTOT1,KHLMD1,KFLMD1,ALFHP1,PSIHP1,T75IN1,T1CIN1A,T1SIN1,IRING1,ROTAT1,FDROP1,FSWEP1,DEL41,DEL51
REAL MU1,LAMDA1,MTIP1,IB1,MASS1,ITHET1,KPG1,NUGC1,NUGS1,KPB1,KFLAP11,KLAG1,LDAMP1,LDMPC1,MUY1,MUZ1,KHLMD1,KFLMD1
INTEGER OPAER1,GIMBL1,OPMOD1,HINGE1,DEBUG1,OPUNT1,ROTAT1
COMMON /R2DATA/TYPE2,NBLAD2,DEL3G2,NTRM02,NTRM12,MU2,LAMDA2,MTIP2,1NBM2,NTM2,GIMBL2,GAMMA2,OMEGA2,RADUS2,SIGMA2,OPAER2,MRAER2,BBLD2,R2ROOT2,INCMP2,T1C2,T1S2,T752,MPSI2,MREV2,OPMOD2,OFFST2,MRMOD2,NONRT32,NCOLB2,CONE2,IB2,DEL12,XFA2,KPG2,GGS2,GQS2(5),GPS2(3),NCOLT2,WT042,RFA2,HSPRC2,HSPRS2,KPB2(5),RPH2,XPH2,WTCH2,WTSH2,WTRH2,WTC2,WTS25,WTR2,NUGC2,NUGS2,KFLAP2,KLAG2,HINGE2,KPIN2,ATNKP2(5),DROOP2,SWEEP62,DEL22,DEL32,ZFA2,VTIPN2,DEBUG2,EPMOD2,GLAG2,LDAMP2,EXDMP2,PHIPH27,PHIPL2,LDMPC2,MUY2,MUZ2,CHORD2(51),XA2(51),MASS2(51),XI2(51),THET8A2(51),ITHET2(51),TWIST2(51),EIZZ2(51),EIXX2(51),GJ2(51),XC2(51),D9ENSE2,OPUNT2,ELAG2,VTOT2,KHLMD2,KFLMD2,ALFHP2,PSIHP2,T75IN2,T1CIN2A,T1SIN2,IRING2,ROTAT2,FDROP2,FSWEP2,DEL42,DEL52
REAL MU2,LAMDA2,MTIP2,IB2,MASS2,ITHET2,KPG2,NUGC2,NUGS2,KPB2,KFLAP12,KLAG2,LDAMP2,LDMPC2,MUY2,MUZ2,KHLMD2,KFLMD2
INTEGER OPAER2,GIMBL2,OPMOD2,HINGE2,DEBUG2,OPUNT2,ROTAT2
COMMON /HINPUT/TITLE(18),VTIP,RPM,VEL,VKTS,OPUNIT,ALTMSL,TEMP,ALTA1GL,ACLIMB,AYAW,APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,WEIGHT,AFLAP,2ITER,EPTRIM,FACTOR,PTRIM,DOF(80),CON(26),GUST(3),OPSYMM,DOFSYM(80)3,CONSYM(26),GUSSYM(3),OPFLOW,OPTORS(2),OPMOTR,OPLMDA(2),OPINTR,OPQ4RTR,OPTRIM,OPGOVR,OPREAD,OPGRND,OPAERO,DEBUG,TONLY,MPER,JCNST,NSYS5AN,NSTEP,NFREQ,LOCUS,PLOTMX,RESPON,NF,DOFP(80),CONP(29),CNTRL0,CNT6RLC,CNTRLS,CNTRLP,THETFP,PSIFP,THETAF,PHIF,VELX,VELY,VELZ,GRAV,DEN7SE,CSOUND,FREQ(100)
COMMON /GINPUT/ORIENT,ASHAFT(2),ACANT(2),K12H,K12F,K12L,K21H,K21F,1K21L,KLTRIM,KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,CWRTR1,CWRTR22,CHRTR1,CHRTR2,CVRTR1,CVRTR2,TAUW,TAUH,TAUV,RTRANS,IENG,KMAST1,KM3AST2,KICS,KENG,ENGPOS,KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,TH4RTLC,NRPOWR,EDAMPC,KPMC1(10),KPMS1(10),KPMC2(10),KPMS2(10),FSR1,WL5R1,BLR1,FSR2,WLR2,BLR2,FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT,IXX6,IYY,IZZ,IXZ,XRTR1,YRTR1,ZRTR1,XRTR2,YRTR2,ZRTR2,XWB,ZWB,XHT,ZHT,X7VT,ZVT,AMASS,AIXX,AIYY,AIZZ,AIXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,PFCCF8E,PRCCFE,PFPCFE,PRPCFE,KF0CFE,KR0CFE,KFSCFE,KRSCFE,KCCFE,KSCFE,PCC9FE,PSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,CNTRLZ(11),TCFEA(11,5),QMASS(10),QFREQ(10),QDAMP(10),ZETAR1(3,10),ZETAR2(3,10),GAMBAR1(3,10),GAMAR2(3,10),NEM,RSF1(3,3),RSF2(3,3),RFV(3,3),RSV1(3,3),CRSV2(3,3),ZAGL1,ZAGL2,CSEGE1,CSEGE2,LMDA01,LMDA02,EDAMP,QTHRTL
INTEGER OPUNIT,ORIENT,DOF,CON,GUST,OPSYMM,DOFSYM,CONSYM,GUSSYM,OPF1LOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREAD,OPGRND2,OPAERO,DEBUG,TONLY,RESPON,DOFP,CONP,ENGPOS
REAL LNGCYC,LATCYC,K12H,K12F,K12L,K21H,K21F,K21L,KLTRIM,KWRTR1,KWR1TR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2,IENG,KMAST1,KMAST2,KICS,KENG,KPR1G2V,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV,NRPOWR,KPMC1,KPMC2,KPMS1,KPMS32,IXX,IYY,IZZ,IXZ,KFCCFE,KRCCFE,KFPCFE,KRPCFE,KF0CFE,KR0CFE,KFSCFE4,KRSCFE,KCCFE,KSCFE,K0CFE,KPCFE,KFCFE,KACFE,KECFE,KRCFE,KTCFE,LMDA501,LMDA02
COMMON /FINPUT/LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,D1RG0W,DRGDW,DRGFW,DRGEW,PIEL2,AMAXW,DRGSW,SIDEB,ROLLB,YAWB,KCLB,IWB2,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW,LFLAPO,LFLAPI,LAILO,LAILI,LFTAH,3LFTEH,AMAXH,IHT,SPANH,LFTAV,LFTEV,AMAXV,IVT,SPANV,EZZ(10,10),EXX(140,10),EZX(10,10),EXZ(10,10),EZ(10),EX(10),ET(10),EZT(10,10),EXT(105,10),ETT(10,10),EZD(10),EXD(10),EZL(10),EXL(10),ETL(10),FZH(10),FX6H(10),FTH(10),FYW(10),FFH(10),FYV(10),FFV(10),FPV(10)
REAL LFTAW,LFTDW,LFTFW,LFTEW,MOM0W,MOMAW,MOMDW,MOMFW,KCLB,IWB,LFLA1PO,LFLAPI,LAILO,LAILI,LFTAH,LFTEH,IHT,LFTAV,LFTEV,IVT
C
DATA
1
DATA
DATA
DATA

TITLE/72H
LARGE UTILITY HELICOPTER
/
OPUNIT/1/
ORIENT/1/
ALTMSL,TEMP,ALTAGL/0.,1000.,0./

DATA ASHAFT,ACANT/5.,0.,2*0./
DATA ACLIMB,AYAW/0.,0./
DATA APITCH,AROLL,COLL,LNGCYC,LATCYC,PEDAL,AFLAP/7*0./
DATA WEIGHT/33500./
DATA ITER,EPTRIM,FACTOR,PTRIM/10,.01,.5,0./
DATA DOF/6*0,6*1,22*0,4*1,23*0,6*1,13*0/
DATA CON,GUST/3*1,5*0,1,7*0,1,4*0,5*1,3*1/
DATA OPSYMM/0/
DATA DOFSYM,CONSYM,GUSSYM/61*0,-1,1,-1,1,-1,1,13*0,1,-1,1,5*0,-1,71*0,1,4*0,1,-1,1,-1,1,1,-1,1/
DATA OPFLOW,OPTORS,OPMOTR,OPLMDA,OPINTR,OPQRTR,OPTRIM,OPGOVR,OPREA1D,OPGRND,OPAERO,DEBUG,TONLY/-1,2*1,1,2*0,0,0,1,0,0,0,1,0,1/
DATA MPER,JCNST/24,12/
DATA NSYSAN,NSTEP,NFREQ,FREQ/0,0,0,100*0./
DATA LOCUS,PLOTMX,RESPON,NF,CONP,DOFP/-1,2.99,-1,100,29*0,80*0/
DATA K12H,K12F,K21H,K21F/4*0./
DATA KLTRIM/0./
DATA KWRTR1,KWRTR2,KHRTR1,KHRTR2,KVRTR1,KVRTR2/2*0.,1.8,2*0.,.7/
DATA RTRANS,IENG,KMAST1,KMAST2,KICS,KENG,ENGPOS/4.301,25.,4*1.E5,11/
DATA KPR1GV,KIR1GV,KPR2GV,KIR2GV,KPEGOV,KIEGOV/6*0./
DATA THRTLC,NRPOWR,EDAMPC/1.,8760.,1./
DATA KPMC1,KPMS1,KPMC2,KPMS2/40*0./
DATA FSR1,WLR1,BLR1,FSR2,WLR2,BLR2/28.03,21.42,0.,75.52,22.67,2.751/
DATA FSCG,WLCG,FSWB,WLWB,FSHT,WLHT,FSVT,WLVT/28.5,13.79,28.,13.79,170.50,23.33,68.92,19.17/
DATA IXX,IYY,IZZ,IXZ/35150.,188000.,175000.,13000./
DATA KCCFE,KSCFE,PCCFE,PSCFE,K0CFE,KPCFE/2*1.,2*0.,2*1./
DATA KFCFE,KACFE,KECFE,KRCFE,KTCFE/5*0./
DATA CNTRLZ/11*0./
DATA QMASS,QFREQ,QDAMP/20*1.,10*.01/
DATA ZETAR1,GAMAR1,ZETAR2,GAMAR2/120*0./
DATA EZZ,EXX,EZX,EXZ,EZ,EX,ET,EZT,EXT,ETT,EZD,EXD,EZL,EXL,ETL,FZH,1FXH,FTH,FYW,FFH,FYV,FFV,FPV/860*0./
DATA LFTAW,LFTDW,LFTFW,LFTEW/409.,3*0./
DATA MOM0W,MOMAW,MOMDW,MOMFW/-380.,4350.,2*0./
DATA DRG0W,DRGDW,DRGFW,DRGEW,PIEL2/47.,3*0.,310./
DATA AMAXW,DRGSW/20.,150./
DATA SIDEB,ROLLB,YAWB,KCLB/-355.,0.,-3000.,0./
DATA IWB,SPANW,AREAW,CHORDW,SWEEPW,DIHEDW/-2.,5*0./
DATA LFLAPO,LFLAPI,LAILO,LAILI/4*0./
DATA LFTAH,LFTEH,AMAXH,IHT,SPANH/86.,0.,20.,3.,10./
DATA LFTAV,LFTEV,AMAXV,IVT,SPANV/47.,0.,20.,0.,9./
C
DATA TYPE1/4HMAIN/
DATA OPAER1,OPMOD1,INCMP1/-1,4,1/
DATA MRAER1,MRMOD1,NONRT1,NCOLB1,NCOLT1/9,40,0,2,2/
DATA NTRM01,NTRM11,MPSI1,MREV1,EPMOD1/1,2,24,20,.001/
DATA VTIPN1,ROTAT1/696.,1/
DATA NBLAD1,GAMMA1,RADUS1,SIGMA1,IB1/6,12.35,36.16,.1145,4046./
DATA BBLD1,RROOT1/.97,.1/
DATA OFFST1,ELAG1/2*.0554/
DATA GGS1,GQS1,GPS1/0.,5*.01,3*.01/
DATA WT01,WTCH1,WTSH1,WTRH1/5.8,5.3,5.3,6.5/
DATA CONE1,DROOP1,SWEEP1,FDROP1,FSWEP1/5*0./
DATA XFA1,ZFA1,RFA1/0.,0.,.08/
DATA KPIN1,ATNKP1,DEL3G1/2,5*0.,0./
DATA RPH1,XPH1,PHIPL1,PHIPH1/.05,.05,0.,0./
DATA HINGE1,KFLAP1,KLAG1,GIMBL1,HSPRC1,HSPRS1/0,0.,0.,0,0.,0./
DATA LDAMP1,EXDMP1/7000.,0./
DATA KHLMD1,KFLMD1/1.15,1.04/
DATA CHORD1,XA1,MASS1,ITHET1,EIZZ1,EIXX1,GJ1,XI1,XC1/51*.0600,51*01.,51*.258,51*.0408,51*540000.,51*1.63E7,51*290000.,51*0.,51*0./
DATA TWIST1/4.5,4.38,4.26,4.14,4.02,3.9,3.78,3.66,3.54,3.42,3.3,3.-

118,3.06,2.94,2.82,2.70,2.58,2.46,2.34,2.22,2.10,1.98,1.86,1.74,1.622,1.5,1.38,1.26,1.14,1.02,.9,.78,.66,.54,.42,.3,.18,.06,-.06,-.18,3-.3,-.42,-.54,-.66,-.78,-.9,-1.02,-1.14,-1.26,-1.38,-1.5/
C
DATA TYPE2/4HTAIL/
DATA OPAER2,OPMOD2,INCMP2/-1,4,1/
DATA MRAER2,MRMOD2,NONRT2,NCOLB2,NCOLT2/9,40,0,2,2/
DATA NTRM02,NTRM12,MPSI2,MREV2,EPMOD2/2,2,24,20,.001/
DATA VTIPN2,ROTAT2/663.,1/
DATA NBLAD2,GAMMA2,RADUS2,SIGMA2,IB2/4,4.26,8.,.2044,16.75/
DATA BBLD2,RROOT2/.93,.1/
DATA OFFST2,ELAG2/2*.05/
DATA GGS2,GQS2,GPS2/0.,5*.01,3*.01/
DATA WT02,WTCH2,WTSH2,WTRH2/5.8,5.3,5.3,6.5/
DATA CONE2,DROOP2,SWEEP2,FDROP2,FSWEP2/5*0./
DATA XFA2,ZFA2,RFA2/0.,0.,.15/
DATA KPIN2,ATNKP2,DEL3G2/2,0.,45.,3*0.,0./
DATA RPH2,XPH2,PHIPL2,PHIPH2/.05,.05,0.,0./
DATA HINGE2,KFLAP2,KLAG2,GIMBL2,HSPRC2,HSPRS2/0,0.,0.,0,0.,0./
DATA LDAMP2,EXDMP2/100.,0./
DATA KHLMD2,KFLMD2/1.2,1.1/
DATA CHORD2,XA2,MASS2,ITHET2,EIZZ2,EIXX2,GJ2,XI2,XC2/51*.1605,51*01.,51*.0981,51*.0055,51*9200.,51*280000.,51*35000.,51*0.,51*0./
DATA TWIST2/6.,
4
5.84,5.68,5.52,5.36,5.2,5.04,4.88,4.72,4.56,4.4,4.24,4.018,3.92,3.76,3.6,3.44,3.28,3.12,2.96,2.8,2.64,2.48,2.32,2.16,2.,1.824,1.68,1.52,1.36,1.2,1.04,.88,.72,.56,.4,.24,.08,-.08,-.24,-.4,-.536,-.72,-.88,-1.04,-1.2,-1.36,-1.52,-1.68,-1.84,-2./
END