DIMENSION X(2,325),Y(2,325) PI=3.1415926 PHID=14.2601 DIA=.032 PSI=.5*PI ELS=1.15334 HEL=.5*3.289 SCLF=10. NPT=50 FI=PI*PHID/360. DO 1 IL=1,NPT ELD=ELS+(IL-1)*1.625/(NPT-1) CALL BN3(XCL,YCL,ANG,ELD,FI,HEL,PSI) DX=.5*DIA*SIN(ANG) DY=.5*DIA*COS(ANG) ANG=ANG*180./PI WRITE(3,48) IL,ELD,ANG,XCL,YCL,DX,DY 48 FORMAT(1X,I3,6(4X,F10.5)) X(1,IL)=XCL+DX X(2,IL)=XCL-DX Y(1,IL)=YCL-DY 1 Y(2,IL)=YCL+DY C WRITE(3,10) 10 FORMAT(//,' CONCAVE SIDE',//) DO 2 I=2,NPT DX=SCLF*(X(1,I)-X(1,I-1)) DY=SCLF*(Y(1,IL)-Y(1,I-1)) C2 WRITE(3,3) I,DY,DX 2 CONTINUE 3 FORMAT(2X,I5,2(3X,F10.5)) C WRITE(3,11) 11 FORMAT(//,' CONVEX SIDE',//) DO 4 II=2,NPT I=NPT+2-II DX=SCLF*(X(2,I)-X(2,I-1)) DY=SCLF*(Y(2,I)-Y(2,I-1)) C4 WRITE(3,3) I,DX,DY 4 CONTINUE END