1      SUBROUTINE M31O37
2C
3C***  EXEC FOR OVERLAY 31, SUBSONIC WING CM, BODY CA,CN
4C
5      COMMON /OVERLY/ NLOG,NMACH,I,NALPHA,IG,IJKDUM(3),NOVLY
6      COMMON /CONSNT/ PI,DEG,UNUSED,RAD
7      COMMON /IWING/  PWING, WING(400)
8      COMMON /WINGI/  WINGIN(100)
9      COMMON /WINGD/  A(195),B(49)
10      COMMON /WHAERO/ C(51)
11      COMMON /FLGTCD/ FLC(95)
12      COMMON /EXPER/  KK(105),KWING
13      DIMENSION C1(6),C3(6),ROUTID(2)
14      LOGICAL FLAG,KWING
15      DATA ROUTID/ 4HM31O, 4H37  /
16      DATA STRA/ 4HSTRA /
17      NOVLY=31
18      CALL CMALPH(A,B,C,WINGIN,WING)
19      CALL CACALC(B,WING)
20C
21C     SET IOM CN, CA, CLA, AND CMA
22C
23      CALL EXSUBT
24      IN = 0
25      IM = 0
26      CLA=WING(101)
27      CMA=WING(121)
28      DO 1000 J=1,NALPHA
29         CA = COS(FLC(J+22)/RAD)
30         SA = SIN(FLC(J+22)/RAD)
31         WING(J+60) = WING(J+20)*CA + WING(J)*SA
32         WING(J+80) = WING(J)*CA - WING(J+20)*SA
33         CALL TBFUNX(FLC(J+22),X,WING(J+100),NALPHA,FLC(23),WING(21),
34     1               C1,IN,MI,NG,0,0,4HCLAW,1,ROUTID)
35         CALL TBFUNX(FLC(J+22),X,WING(J+120),NALPHA,FLC(23),WING(41),
36     1               C3,IM,MI,NG,0,0,4HCMAW,1,ROUTID)
37 1000 CONTINUE
38      IN=0
39      IM=0
40      IF(KWING)CALL TBFUNX(0.,X,CLA,NALPHA,FLC(23),
41     1              WING(21),C1,IN,MI,NG,0,0,4HCLAW,1,ROUTID)
42      IF(KWING)CALL TBFUNX(0.,X,CMA,NALPHA,FLC(23),
43     1              WING(41),C3,IM,MI,NG,0,0,4HCMAW,1,ROUTID)
44C
45C***  IF THE WING LIFT DEVIATES FROM THE LINEAR VALUE BY
46C***  15 PERCENT OR MORE SET CM AND CMA TO NA (2*UNUSED)
47C
48      IF(A(7) .LE. (6./A(124)) .AND. WINGIN(15) .EQ. STRA) GO TO 1020
49      IF(KWING) GO TO 1020
50      FLAG=.FALSE.
51      DO 1010 J=2,NALPHA
52         DEL = 100.*ABS((WING(J+100)/CLA-1.0))
53      IF(DEL.GT.90.0) FLAG=.TRUE.
54         IF(FLAG) WING(J+40)  = 2.0*UNUSED
55         IF(FLAG) WING(J+120) = 2.0*UNUSED
56 1010 CONTINUE
57 1020 CONTINUE
58      WING(101)=CLA
59      WING(121)=CMA
60      RETURN
61      END
62