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