RUNT VERSION FES T^ 8 17«12 04/23/76 



320 IF(RKMQR( I.J) .60. -0.0) RKMOR(I»J) ' 3.0 



331 RKTB(I.J)- (RKHYD(I.J) + RKMOR(I.J))/0 



3<t6 130 IFd.EO. 3 .OR. J .EO. 3)RKTB(I.J) -RKTB ( I, J) /B 



372 PKTB(3»3) • RKTa(3»3)/9 



^02 PM4SSB • RM4SS/(AREA*B0F) 



<i06 RINEPB » RINERT/( A3EA*R0E*B*e) 



'»12 XG8 = XG/a 



413 YGB • rG/B 



C*«***START WAVELENGTH LOOP 



415 00 150 IL • 1.10 



415 IF (IL .GT. NWAVEL) GO TO 160 



C*#***SET VALUES IN NQNDIMENTION ALI2ED ALGE^RATC EQUATIONS OF MOTTTN 



^21 Atl,l) • RKTB(i»i) - WN(IL) ♦ (RMA'n'^B + RMU(1,1,IL)) 



i,<,5 A(2»2)» RKT8(2.2) - WN(IL) * (RHAS'^R ♦ ' ^'U( 2. ?» IL ) ) 



^71 A(3,3)- RKTB(3»3) -YGB + RMASSB -WN(TL) ♦ (i?TMERB + RMU(3.3.TL) + 



1 (X68**2 •fyGB**2) ♦ Rf^ASSB) 



531 A(l*2) ' RKTe(l,2) - WN(IL) ♦ RMU(1.?.TL) 



552 A(l,3) » RKTB(1»3) - WM(IL) ♦ (RMU(1.3,TL) -YG8*RM*SSB) 



576 A(2»3) = RKTB(2.3) - WN(IL) ♦ (RHU(2.3.IL) +XGB*RmaSSB) 



623 A(2>1) • A(1.2) 



633 A(3.1) • A(l,3) 



643 A(3.2) = A(2.3) 



652 DO 20 I = 1.3 



654 DC 10 J = 1»3 



655 A(I*3# J*3) • A(I,J) 



666 A(I,J+3) «RLAMti, J.IL)*SORT(WN(IL) ) 



C ADO CORRECTION FOR VISCOUS DAMPING 



705 IF (I .£0. J) A(I.J + 3) » (1.0 ♦ DA*IP( I) ) *A ( I. J + 3 ) 



723 10 A(I+3#J) • - A(I,J+3) 



736 C(I) - FBd.ID* SIN(DELFB(I,IL)) 



753 20 C(I+3) « FBd.ID ♦ COS ( DELF3 ( I . IL ) ) 



771 SCALE » 1. 



773 NN « 6 



C*«**«S0LVE ALGEBRAIC EQS OF MOTION. 8(1). B(?). B(3) « AMOLITUOFS 

 C««*4i*0F COS(WT). 6(4). 3(5). 9(6) » AMPLITUDES OP STN(WT) FOR SWAY, HEA 

 C»****AND ROLL AT CENTER FO USERS COORDINATE SYST?'!. 



774 LL • LNEQF(6.NN,1, A. C. SCALE. ERASE) 



1005 DO 30 I » 1.3 

 C*****AMPLITU0E AND PHASE OF RESPONSE 



1006 PARd.IL) = S0RT(C(I)**2 ♦ C(I*3)**2) 

 1032 30 DELRd.IL) • ATAN2( C ( I ) .C d + 3 ) ) 



1047 DC 41^ I « 1. NNW 



1050 AW =• 0. 



1051 BW • 0. 



1052 90 DO 50 J ■ 1.3 



C****«PESULTANT WAVE AMPLITUDE AND PHASF FHR SWAY, HEAVE ROLL. 

 1054 HwRd.J.IL) = HWBd.J.lL) * RAR(J.TL) 



1072 OELWRd.J.IL) ' OELWd,J,IL) ♦ DELR(J.tl) 



1111 AW • AW ♦ HWR(I.J» IL) * SIN(DELWR(T. J.TD) 



1132 50 BW » 8W + HWRd.J.IL) * COS ( DE LWR d . J . IL ) ) 

 1155 IF(XOLd.IL).LT. 0.) GO TO 70 



C****+T0TaL REFLcCTEO wave (VECTOR ADOITTT^) 

 1163 AW • AW + rlWB(I.6.IL) ♦ S IN ( DELW( T . «•. TL ) ) 



1204 BW « BW + HWBd,6,IL) * COS ( DELW( I . ft. ^L ) ) 



1225 GO TO 45 



C«****TOTAL TRANSMITTED WAVE(VECTOR ADOTTTTN) 



Table D-2. Continued 



117 



