131 

 132 

 133 

 134 

 135 

 136 

 137 

 138 

 139 

 140 

 141 

 142 

 143 

 144 

 145 

 146 

 147 

 148 

 149 

 150 

 151 

 152 

 153 

 154 

 155 

 156 

 157 

 158 

 159 

 160 

 161 

 162 

 163 

 164 

 165 



1 

 i 



3 

 4 

 5 

 6 



9 



10 

 11 



12 

 13 

 14 

 15 

 16 

 17 

 18 

 19 

 20 

 21 

 22 

 23 

 24 



DO 219 J=1,N 

 DUM1(1,J)=DUH1(2,J) 



DUHHM,J>=0.0 



DUMHM1,J)=0.0 



DUM1(N2,J)=0.0 



DO 219 1=1, M 



IF(DUhl(I,J).EQ.1.0)DUMHI,J)=11.0 

 219 CONTINUE 



CALL P0UT<IP1,IP2,INC,JP1.JP2.' BREAKER INDEX '.DUH1.1.0) 

 CAAAAAAAA ANGLES, HEIGHTS, AND HAVE NUMBERS ARE MULTIPLIED BY ZMULT , 

 CAAAAAAAA HMULT. AND RKHULT RESPECTIVELY IN IHE PRINT OUT. IHESE MUSI 

 CAAAAAAAA BE CHANGED INTERNALLY WITHIN IHE MAIN PROGRAM IE OTHER VALUES 

 CAAAAAAAA ARE DESIRED 

 C 



DO 23 1=1, H 



DO 23 J-l.N 



DUM3=Z(I,J)-PI 



DUM1(I.J)=DUM3ARAD 

 23 CONTINUE 



CALL P0UKIP1.IP2, INC, JP1,JP2,' WAVE ANGLES (DEG)',DUM1, 

 AZMULI) 

 C 



CAAAAAAAA CONVERI EROM AMPLITUDE EUNCTION TO WAVE HEIGHT 

 C 



DO 81 1=1. M 



DO 81 J=1,N 



DUM1(I,J)=H(I,J)AHEACI 

 81 CONTINUE 



CALL P0UI(IP1,IP2,INC.JP1.JP2,' WAVE HEIGHTS ,DUM1, 



lHHULI) 



CALL POUT ( IP1 , IP2 , INC , JP1 . JP2 , ' WAVE NUMBER ' . GRDK , 



ARKMULI) 

 99 REIURN 



END 

 CAAAAAAAAAAAAAAAAAAA*AAAAAAAA*AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

 C 



SUBROUT INE GRADK ( ISIARI , IEND ) 

 C 



CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

 C 



C THIS SUBROUTINE COMPUTES THE UPDAIED VALUES OE IHE GRADIENT 

 C OE THE WAVE PHASE EUNCTION ALONG A GIVEN ROW 

 C 



CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

 C 



PARAMETER (IQ=95,JQ=95) 

 C 



COMHON/ANGLES/Z ( IQ , JQ ) . S I ( IQ . JQ ) , CO ( IQ , JQ > 



COMMON/DEPIHS/D(IQ,JQ> 



COMMON/WVH/H( IQ, JQ) ,CCG< IQ, JQ) 



COMM0N/WAVNUM/RKA( IQ, JQ) ,GRDK( IQ, JQ) 



COMHON/CONSI/ G , P I , P 12 , RAD , HCONVR , SCONVR , DX , DY , DX2 , DY2, I , OMEG , 



AM,N.M1,M2,NM1.NM2,IWEI.IDRY.IWEIP1,CNTRANG,HFACI,DLEVEL 



C0MM0N/C0NS2/ IGR ID , IREE , IIAMX , IIHMX, ID IEE 



COMHON/MUSC/XMUC ( IQ ) , XMUS ( IQ ) , YMUC ( JQ ) . YMUS ( JQ ) 



COMMON /TRANE/ DECAY,STABL.IBRK(JQ),IBRKM(JQ) 

 DIMENSION DUMl ( IQ.JQ), X(JQ) 



REED IE 



REED IE 



REED IE 



REFDIE 



REEDIF 



REED IF 



REEDIE 



REFDIE 



REEDIF 



REFDIF 



REEDIF 



REFDIF 



REEDIF 



REFDIE 



REFDIE 



REFDIF 



REEDIF 



REFDIE 



REFDIF 



REFDIF 



REEDIF 



REFDIE 



REFDIF 



REEDIF 



REEDIF 



REFDIF 



REEDIE 



REFDIE 



REFDIF 



REEDIF 



REFDIF 



REFDIE 



REEDIE 



REFDIF 



REFDIF 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



PARAM 



PARAM 



PARAM 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



GRADK 



130 

 131 

 132 

 133 

 134 

 135 

 136 

 137 

 138 

 139 

 140 

 141 

 142 

 143 

 144 

 145 

 146 

 147 

 148 

 149 

 150 

 151 

 152 

 153 

 154 

 155 

 156 

 157 

 158 

 159 

 160 

 161 

 162 

 163 

 164 



3 

 4 



F9 



