SUBROUTINE WEDGE (F , FABS , FPHA , XRL , XTH , WEDGEA , WAVEA , IDX ) WEGDE1 

 C ********************************** *NOTICEl 



C * THIS COMPUTER PROGRAM WAS WRITTEN BY H.S. CHEN OF CERC IN #N0TICE2 



C * 1986 UNDER CORP' CIVIL WORK R&D PROGRAM. NEITHER ANY OF *N0TICE3 



C * AGENCIES NOR ANY INDIVIDUAL ASSUMES ANY LEGAL LIABILITY OR *N0TICE4 



C * RESPONSIBILITY FOR THE ACCURACY OF THE PROGRAM. *N0TICE5 

 C ***#♦*******♦****#***#*###****#*#* *N0TICE6 



c ________ __ _ -WEDGE2 



C WAVE REFLECTION AND DIFFRACTION BY A VERTICAL WEDGE. WEDGE3 



C WEDGE4 



C INPUT : XRL = R/L, (RADIUS VECTOR) / (WAVE LENGTH) WEDGES 



C XTH = VECTORIAL ANGLE IN DEGREE WEDGE6 



C WEDGEA = WEDGE ANGLE IN DEGREE WEDGE7 



C WAVEA = INCIDENT WAVE ANGLE IN DEGREE WEDGEB 



C OUTPUT : F = PHI, WAVE RESPONSE NORMALIZED BY INCIDENT WEDGE9 



C WAVE AMPLITUDE WEDGEIO 



C FABS = ABSOLUTE VALUE OF PHI, THE AMPLIFICATION WEDGE11 



C FACTOR WEDGE12 



C FPHA = PHASE OF PHI IN RADIAN. WEDGE13 



c _____________ ____ _______ -WEDGE14 



PARAMETER(NN=200) WEDGE15 



DIMENSION BJ(NN) ,W(NN) ,XM(NN) WEDGE16 



COMPLEX F,TM WEDGE17 



DATA T0LR/1.E-8/, ITER/8/ WEDGE18 



PI=3, 141592654 WEDGE19 



C P I = P I / 1 8 . WEDGE20 



XKR=XRL*2.0*PI WEDGE21 



TH=XTH*CPI WEDGE22 



WA=WAVEA#CPI WEDGE23 



XNU= (360. -WEDGEA)/ 180. WEDGE24 



X M ( 1 ) = 1 . WEDGE25 



IF ( IDX. NE. 0) GOTO 14 WEDGE26 



CALL BESJ(XKR,0.0,1,W,NZ) WEDGE27 



BJ(1)=W(1) WEDGE28 



ICOUNT=0 WEDGE29 



DO 10 N=l ,NN WEDGE30 



N1=N+1 WEDGE31 



ICQUNT=IC0UNT+1 WEDGE32 



IFdCOUNT.LE. ITER) GOTO 8 WEDGE33 



NNN=N WEDGE34 



GOTO 14 WEDGE35 



8 XM(N1)=FL0AT(N)/XNU WEDGE36 

 M=INT(XM(N1) ) WEDGE37 

 ALPHA=XM(N1)-M WEDGE38 

 M1=M+1 WEDGE39 

 CALL BESJUKR, ALPHA, Ml, W,NZ) WEDGE40 

 IF(Nl.EQ.NN) WRITE(6,9) XKR , ALPHA ,M1 , W (Ml ) , NZ WEDGE41 



9 FORMAT(/' **** NO. OF TERMS FOR SUMMATION IS INSUFFICIENT ****', WEDGE42 

 1 /' XKR, ALPHA, Ml, W(M1) ,NZ = ' , 2F10. 4 , 15 ,E15. 6 , 15/ ) WEDGE43 



BJ(N1)=W(M1) WEDGE44 



Figure A3. Computer program list 2 (Sheet 1 of 25) 

 A5 



