RUNT VERilON FEB 74 9 17S12 04/23/76 



6 

 6 

 6 



7 

 23 

 42 

 47 

 53 

 55 

 56 

 57 

 60 

 61 

 63 



66 



70 



104 



123* 



126 

 136 

 145 



146 

 166 

 200 

 205 

 207 

 211 

 213 



214 

 215 

 216 

 224 

 255 

 306 

 311 

 326 



344 

 361 



SUBROUTINF PHYSCL(K) 

 COMMON Rn2(25.25)»RK56(25.25)» PO 

 1FI(25,6), RI(25.25)f RJ(25»25), R 

 2RMU(3»3»10). RLAM(3,3»lu)» F6(3»1C 

 3 0ELW(25,6»101 . xaL(25.10) 



COMMON/ONE/X(25)fY(25»,X8(25).Y8(2 

 1»FE1N(25). FirN(25), RNQR^"( 25, 3 ) » 

 COMMOM /TWO/ N,NNW» NWAVEL, ISYM, 

 COMMON/THREE/ W»VEL(10), WN(IO). B 

 COMMON /SEVEN/ ARtAf B, D, ROE» G 

 REAL K 



DO 3 I = l.N 

 C***«*MODE 6 = INCIDENT ♦ 

 FEtI»6) « Fe(l,4) ♦ 

 3 FI<I,6) » FI(I>4) ♦ 

 FACM • (B**2)/*REA 

 FACM+SORT(K) 

 FACM * K 



T(»5f?5), H0W(25,25)»FE(25,f ), 



<(?•!, 4), PL(25,4). 



). DEL'^'^O.IO). HWB(25,6»10). 



5)»ANr,(25).DEL(25).VV(25) 



JC«5) 



ISKIB, NC. PIE.GA^I-^A.M.TK.TP 



01(10), TL 



EE. BTTTL£(3), TITLF(e) 



OIFRACTED POTENTIAL "^ 



FE(I.5) 



FI(I,5) 



FACL 

 FACF 



DO 1 L ' 1,3 

 DDl ri ' 1,6 

 RA • 0.0 

 RM = 0.0 



IF(M1 .EO, 4) 60 TO 1 

 IF(M1 .EO. 5» GO TO 1 

 C+*#**INTIGRAT£ PRESSURE COMPONENTS OVER 

 DO 5 I = 1,N 



RM « RM + FF(I,M1I ♦ RNORMd.L) 

 5 RA - RA + F1(I,M1» ♦ Pf40«K{I,L) 

 1F(M1 .GT. 4) GO TO 8 

 C*****ADDED MASS AND DAMPING IN DIRECTION L ODE TO tlOTIDN Ml AT 

 C*****WAVELENGTH IL. 



RMU(L,M1, ID = R1»FACM 

 RLAM(L,M1, IL) » RA*FACL 

 GO TO 1 

 C*****HAVE FORCE AMPLITUDE AND PHASE IN 

 C*****1NCI0ENT WAVE AT WAVELENGTH IL 

 8 FD(L, IL) = S0RT(RW**2 -f RA**2) 



DELFB(L,IL) • ATAN2(-PA,RM) 

 1 CONTINUE 



IW - N + 1 

 IMAX « N + NNW 

 DO 30 I • IW,IMAX 

 DO 6 L ■ 1,4 

 C+**#*C0MPUTE POTENTIAL AT FREE SURFACE POINTS USING GREENS TMEORUH 

 DO 4 J ' 1, N 

 ■ 1,5 



JCdO) GO TO 4 

 FE(I,L) 

 FI(I,L) 



DFL(T) 



nPL(i) 



ntRF^TION Ml DUE TO 

 * P4CF 



10 



DO 10 IC 



IF(J .EO 



FEd.L) 



FI(I,L) 



4 CONTINUE 

 F£(I,L) 



6 FI(I,L) 

 C*****MO0E 6 ' 

 FE(I,6) 

 FI{I,6) 



FE( J,L)'«>RI(I,J 

 FE{J,L)*RJ(I,J 



) -'=T(J,L)*RJ(I,J) 

 ) ♦FT(J,L)*RI(I,J) 



(FE(I,L) - RK(I,L))/TP 

 (FI(I,L) - RLd.L) )/TP 

 INCIDtNT + OIFRACTED POTPNTIAL*: 

 FECI, 4) + FE(I,5) 

 FI(I.4) ♦ Fid, 5) 



Table D-2. Continued 



125 



