14500 

 14G00 

 14700 

 14800 

 14900 

 15000 

 15100 

 15200 

 15300 

 15400 

 15500 

 15600 

 15700 

 15800 

 15900 

 16000 

 16100 

 16200 

 16300 

 16400 

 16500 

 16600 

 16700 

 16800 

 16900 

 17000 

 17100 

 17200 

 17300 

 17400 

 17500 

 17600 

 17700 

 17800 

 17900 

 18000 

 18100 

 18200 

 18300 

 18400 

 18500 

 18600 

 18700 

 18800 

 18900 

 19000 

 19100 

 19200 

 19300 

 19400 

 19500 

 19600 

 19700 

 19800 

 19900 

 20000 

 20100 

 20200 

 20300 

 20400 

 20500 

 20600 

 20700 

 20800 

 20900 

 21000 

 21 100 

 21200 

 21300 

 21400 

 21500 

 21600 



DT = (T-T1 )/DELTAT 



DTT=(T2-T)/DELT 



HCC=HC(I )*DT+HC(I-1 )*DTT 



GO TO 446 

 444 CONTINUE 

 446 CONTINUE 



IF(HGEN.LT.HCC) GO TO 797 



ANGGEN=ALPWIS*TWOP 1/360. 



CALL WVNUM(WDEPTH,T,DUMKK) 

 "ANGGEN.HGEN.CGEN.CGGEN REPRESENT THE WAVE ANGLE , HEIGHT , CELER 

 ►*GROUP VEL(RESPECT. ) OF THE SPECIFIED WAVE INPUT AT A DEPTH, 



CALL WVNUM( 1 1 .O.T.DUMKKK) 



CI 1=TW0PI/(T*DUMKKK) 



CGI 1=0.5*C11*( 1 . + (2. *DUMKKK*1 1 .0/SINH(2. *DUMKKK* 1 1.0))) 



CGEN=TWOPI/(T*DUMKK) 



CGGEN=0.5*CGEN*( 1 .+(2. *DUMKK*WDEPTH/SINH( 2 . *DUMKK*WDEPTH 



CALL TRANS 

 797 IF(NCHECK.NE.NUNIV) 

 709 GO TO 798 

 000 CONTINUE 



STOP 



END 



ITY AND 

 WDEPTH 



))) 



NUNIV=NCHECK 



k ****** 4 



SUBROUTINE TRANS 

 THIS SUBROUTINE WILL COMPUTE SED-IMENT TRANSPORT 



COMMON/ A/ C(G0,2C) ,RK(60,20),Y(60.20) . DEEP ( 60, 20) , ALPHAS ( 60 , 20) 



COMMON/ A A/ Y ZERO (GO) 



COMMON/BB/WE0(6O.2O) 



COMMON/B/ THETA(60,20) ,0XT0T(60) , OLDANG( 60, 20) , DY(60,20) 



COMMON/C/ H(60,20),CG(60.20) , HOLD (60,20 ) ,HB( 60,20 ) ,YB( 60) 



COMMON/N USED/dUSE,T,CO,CGEN,CGGEN, ANGGEN,DX,BERM.THETAO( 10) ,MMAX 



COMMON/D/SIGMA,G,ELO,dMAX, I MAX, PI ,TWOPI .PI02.HGEN. I JET ( 10) , S JETTY 



COMMON/ E/RHO , RHOS , POROS , CONST , TKS I 



COMMON/F/ADEAN,REPOSE,DIAM 



C0MM0N/G/IBREAK(60) ,HN0NBR(2O) 



COMMON/P/HB0(6O) ,DEEPB(60) 



COMMON/ZZZ/NTIME 



COMMON/AAA/DELT,NTIMES 



COMMON/COUNT/NUN IV 



DIMENSION Y0LD(60,20) , R(60, 20) . S{60, 20) , HC( 60 , 20) , 0Y(6O, 20) , YDI SS( 



60,20) 

 DIMENSION RHS1(G0,20) , S3 (60, 20) , THETAB( 60 , 20) , ANGL0C( 60 , 20) 

 DIMENSION DISTR(60,20) , AWARE (60, 20 



K*******+^ 



It ****** :( 



********** NOTE : SIZE OF ABAND AND XL HAVE TO BE CHANGED 

 ********** ACCORDING TO JMAX+1 + dMAX AND dMAX+ 1 , RESPECT . 



********** CHANGE REO'D AT 7040 AND 18650 



********************+****+*****************************+**********i 



* ) ,BMATRX(4 32) , ABAND (432, 1 9 ) , OX ( GO , 20) , XL ( 432 , 10) , CONSTG ( GO , 20 ) 



COMMON/MP/ RKB(60),HBI(60) ,DEEPBI(60) 



C0MM0N/EXPL/0YEXP(6O,2O), YIMP(60,20) 



DIMENSION SANGLE(20) 

 'S ZERO-OUT ALL OF THE DIMENSIONED MATRICES. 



DO 1000 U=1 ,0MAX+2 



SANGLE(d)=0.0 



DO 1000 1 = 1 , IMAX + 2 



YOLD( I ,d)=0.0 



R(I, J)=0.0 



S(I , J)=0.0 



HC(I,d)=0.0 



0Y( I , J)=0.0 



YDISSd, J)=0.0 



RHSKI ,d)=0.0 



S3(I, J)=0.0 



THETAB( I ,d)=0.0 



ANGLOCd ,d)=0.0 



DISTR( I ,d)=0.0 



AWARE( I ,d)=0.0 



OX(I ,d)=0.0 



CONSTGd , J)=0.0 



68 



