DIMENSION DY(7),PHID(7),GAM(7),XX(2),CC(2) C------DATA VALUES PI=3.14159 DTR=PI/180. RTD=1./DTR AL=3.289 HEL=AL*.5 WD=.014 PSI=PI/2. GS=.5*(SQRT(5.0)-1.0) DO 45 I=1,7 45 DY(I)=.045-(I-1)*.015 GAM(1)=31.5781 GAM(2)=330.9301 GAM(3)=63.4400 GAM(4)=0.0 GAM(5)=299.1013 GAM(6)=91.0496 GAM(7)=270. PHID(1)=14.2601 PHID(2)=PHID(1) PHID(3)=24.2763 PHID(4)=0.0 PHID(5)=PHID(3) PHID(6)=PHID(1) PHID(7)=PHID(1) C------FIND THE X WHERE THE SPACING IS A MINIMUM DO 16 II=1,6 JJ=II+1 KNT=1 XS=0.0 XF=1.5 F1=PHID(II)*DTR*.5 F2=PHID(JJ)*DTR*.5 SN1=SIN(GAM(II)*DTR) SN2=SIN(GAM(JJ)*DTR) CS1=COS(GAM(II)*DTR) CS2=COS(GAM(JJ)*DTR) C------GOLDEN SECTION SEARCH JN=3 11 DX=GS*(XF-XS) XX(1)=XF-DX XX(2)=XS+DX DO 7 J=1,2 IF(J.EQ.JN) GO TO 7 X=XX(J) CALL BN2(X,Y1,A1,EL,F1,HEL,PSI) CALL BN2(X,Y2,A2,EL,F2,HEL,PSI) DYY=Y1*CS1+DY(II)-Y2*CS2-DY(JJ) DZZ=Y1*SN1-Y2*SN2 CC(J)=SQRT(DYY**2+DZZ**2) 7 CONTINUE C-----INCREMENT COUNT AND QUIT IF OVER 20 KNT=KNT+1 IF(KNT.GT.20) GO TO 1 IF(CC(1)-CC(2)) 8,9,10 C------C(1)C(2) : RIGHT SIDE OF INTERVAL 10 CC(1)=CC(2) JN=1 XS=XX(1) GO TO 11 1 C=.5*(CC(1)+CC(2))-.007*(1.0/COS(A1)+1.0/COS(A2)) 16 WRITE(5,4) II,JJ,X,Y1,Y2,C 4 FORMAT(//,' WIRES ',2I1,' X : ',F4.2,' Y1 : ',F8.6, 2' Y2 : ',F8.6,' C : ',F8.6) END