1      SUBROUTINE XYCORD(IVAL,INT)
2C
3C  CALCULATE MEAN LINE AND THICKNESS DISTRIBUTION FROM X,Y COORDINATES
4C
5      COMMON /IBODY/ PB,NACA(80)
6      COMMON /IWING/ PW, X(60)
7      COMMON / IHT / PHT, XU(60),XL(60),YU(60),YL(60)
8      COMMON / IBW / PBW, N
9      COMMON / IVT / PVT, YUN(60),YLN(60)
10      COMMON / IBH / PBH, THN(60),CAM(60)
11      COMMON /IBWHV/ PBWHV, RHO,TMAX
12      DIMENSION ROUT(2),C(6)
13      REAL NACA
14      DATA ROUT /4HXYCO, 4HRD  /
15      DATA A4W,A4H,A4V,A4F/4HW   ,4HH   ,4HV   ,4HF   /
16      IF(IVAL.NE.0) GO TO 1010
17      DO 1000 I=1,N
18         THN(I)=0.5*(YU(I)-YL(I))
19         CAM(I)=0.5*(YU(I)+YL(I))
20 1000 CONTINUE
21 1010 THN(1)=0.0
22      THN(N)=0.0
23      CAM(1)=0.0
24      CAM(N)=0.0
25      IN=0
26      DO 1020 I=1,N
27         CALL TBFUNX(X(I),Y,DYDX,N,X,CAM,C,IN,MI,NG,0,0,4HDCAM,1,ROUT)
28         THETA=ATAN(DYDX)
29         SA=SIN(THETA)
30         CA=COS(THETA)
31         XU(I)=X(I)-THN(I)*SA
32         XL(I)=X(I)+THN(I)*SA
33         YUN(I)=CAM(I)+THN(I)*CA
34         YLN(I)=CAM(I)-THN(I)*CA
35 1020 CONTINUE
36      XU(1)=0.0
37      XL(1)=0.0
38      XU(N)=1.0
39      XL(N)=1.0
40      YUN(1)=0.0
41      YUN(N)=0.0
42      YLN(1)=0.0
43      YLN(N)=0.0
44      IF(INT.EQ.0)RETURN
45C
46C  DOUBLE THICKNESS TO REFLECT USER INPUTS
47C
48      DO 1030 I=1,N
49 1030 THN(I)=2.*THN(I)
50      WRITE(6,1050)
51      IF(INT.EQ.1)WRITE(6,1060)
52      IF(INT.EQ.2)WRITE(6,1061)
53      IF(INT.EQ.3)WRITE(6,1062)
54      IF(INT.EQ.4)WRITE(6,1063)
55      IF(INT.EQ.1)NACA(6)=A4W
56      IF(INT.EQ.2)NACA(6)=A4H
57      IF(INT.EQ.3)NACA(6)=A4V
58      IF(INT.EQ.4)NACA(6)=A4F
59      WRITE(6,1070)
60      WRITE(6,1080)(XU(M),YUN(M),XL(M),YLN(M),X(M),CAM(M),THN(M),
61     1              M=1,N)
62C
63C  NOW HALVE THICKNESS FOR PROGRAM COMPUTATIONS
64C
65      DO 1040 I=1,N
66 1040 THN(I)=.5*THN(I)
67 1050 FORMAT(1H1,29X,75HAUTOMATED STABILITY AND CONTROL METHODS PER APRI
68     1L 1976 VERSION OF DATCOM     )
69 1060 FORMAT(1H ,55X,25HUSER DEFINED WING SECTION)
70 1061 FORMAT(1H ,50X,36HUSER DEFINED HORIZONTAL TAIL SECTION)
71 1062 FORMAT(1H ,50X,35HUSER DEFINED VERTICAL TAIL SECTION )
72 1063 FORMAT(1H ,51X,34HUSER DEFINED VENTRAL FIN SECTION  )
73 1070 FORMAT(1H ,3X,15H UPPER ABSCISSA,4X,15H UPPER ORDINATE,4X,
74     1              15H LOWER ABSCISSA,4X,15H LOWER ORDINATE,6X,
75     2              16HX-FRACTION CHORD,4X,9HMEAN LINE,4X,9HTHICKNESS)
76 1080 FORMAT(8X,F8.5,11X,F8.5,10X,F8.5,10X,F8.5,13X,F8.5,10X,
77     1       F8.5,6X,F8.5)
78      RETURN
79      END
80