21700 1000 CONTINUE 



21800 RH0= 1 . 99 



2 1900 RH0S=5 . 14 



22000 POROS=0.40 



22100 C0NST=O.77 



22200 CAPPA=0.78 



22300 TAU=0.25 



22400 TKSI = (CONST*RHO*SORT(G))/((RHOS-RHO)*( 1 . 0-POROS ) * 16 . 0*SORT ( CAPPA ) 



22500 C* OXd.J) IS THE TRANSPORT BETWEEN THE (I,J+1) AND (I, J) CONTOURS. 



22600 C*THE 'DO 1 LOOP' SIMULATES TIME T IME = DELT*NT IMES . 



22700 COFF =0.00001 



22800 GAMMA=RHO*G 



22900 DO 1 NTIME=1 ,NTIMES 



23000 NUNIV=NUNIV+1 



23100 C*THE MATRICES ABAND AND BMATRX MUST BE "ZEROED OUT" 



23200 K = 



23300 DO 26 I=2.IMAX-1 



23400 DO 26 d= 1 , JMAX 



23500 K=K+1 



23600 BMATRX(K)=0.0 



23700 DO 26 L= 1 , JMAX+ 1+JMAX 



23800 26 ABAND(K,L)=0.0 



23900 XNTIME=1 .0*(NTIME) 



24000 CALL PREDIF 



24100 C*SM00THING OF THE WAVE ANGLE. THETA, IS RE'D TO ACCT FOR DIFF EFFECTS. 



24200 CALL SMOOTH( THETA , IMAX , JMAX , I JET . SJETTY . MMAX , Y ) 



24300 CALL OTRAN 



24400 C*FIRST THE LONGSHORE SEDIMENT TRANSPORT WILL BE DISTRIBUTED 



24500 C****ACR05S THE SURF ZONE.... 



24600 CC=1.25 



24700 C***OX(I.J) WILL BE DETERMINED BY SUBTRACTING FROM THE INTEGRAL 



24800 C**0F OX FROM DEEP(I.J-I) TO INFINITY, THE INTEGRAL OF OX FROM DEEP(I,J) 



24900 C***TO INFINITY. IN THIS WAY THE SEDIMENT TRANS FROM JMAX OUT GETS 



25000 C***INCLUDED IN OX(I,JMAX). TO INCLUDE THE SWASH TRANS, WHEN J=1 



25100 C*WE WILL SUBTRACT FROM 2 TO INFINITY FROM 1.0 



25200 C*LOOP FOR VALUES WHICH ARE HELD CONST AND STORED. 



25300 THETAB( 1 , 1 ) =0 . 5* ( THETA( 1 , 1 )+0.0) 



25400 R(1,1)=0.5/(DX*(DEEP(1,1 )+BERM/2. ) ) 



25500 DO 290 1=2, IMAX 



25600 R(I , 1 )=0.5/(DX*(DEEP(I , 1 )+B.ERM/2. ) ) 



25700 C* THETAB( I , 1 )=0.25*(THETA( I , 1 )+THETA( 1-1 , 1 )+0.+0. ) 



25800 THETABd , 1 ) =0 . 5* ( THETA (I , 1 ) + THETA(I-1 , 1 )) 



25900 C*NO NEED TO COMPUTE PROP ANGLE AT STRUCTS BECAUSE OX =0.0 AT IJET(M)+1 



26000 ANGLOCd , 1 )=ATAN( (Y(I,1)-Y(I-1,1 ))/DX) 



26100 C*HBO( IJET(M)+1 ) IS PROPERLY SET IN THE SUBROUTINE OTRAN. 



26200 DISTRd , 1 ) = 1 . 0-EXP( -( (DEEP( I , 1 ) ** 1 .5+HB0(I )*ADEAN**1 .5)/ 



26300 * (CC*DEEPB(I )**1 .5))**3) 



26400 DISTR( I , 1 )=DISTRd , 1 ) *TKSI *HB0( I ) **2 . 5 



26500 DO 290 J=2,JMAX 



26600 Rd,J)=0.5/(DX'(DEEP(I,J)-DEEPd,J-1))) 



26700 THETABd , J ) =0 . 5* ( THETA ( I ,d)+THETA(I-1 , J)) 



26800 ANGLOC( I,J)=ATAN((Y(I,J)-Yd-1.J) )/DX) 



26900 DISTRd , J) = EXP(-((DEEP(I ,J-1)**1 . 5+HB0( I ) *ADEAN* * 1 . 5 )/(CC*DEEPB( I ) 



27000 * **1 .5))**3)-EXP(-((DEEP(I , J)**1 . 5+HB0( I ) *ADEAN** 1 .5)/(CC* 



27100 * DEEPBd )**1 .5) )**3) 



27200 DISTRd , J)=DISTR(I , J ) *TKSI *HBO(I ) **2 . 5 



27300 290 CONTINUE 



27400 DO 301 J=1,JMAX 



27500 DO 301 1=2, IMAX 



27600 AWARE (I , J ) =DE LT*R (I , J ) * ( OX ( I . J ) -QX ( 1+ 1 , J)+OYd , J)-OY(I , J+1))+Y(I , J 



27700 * ) 



27800 S 1 =2. *SIN( THETABd , J ) ) *COS ( THETAB ( I ,J))*(-1.+2. *(COS( 



27900 • ANGLOCd , J) ) )**2) 



28000 S2=C0S(2 . *THETAB( I , J) ) *COS( ANGL0C( I , J ) )/( SORT (DX**2+ 



28100 * (Yd , J)-Y(I-1 , J))**2) ) 



28200 S3(I , J)=S2*DISTR(I , J) 



28300 IF(SJETTY . EO.0.0) GO TO 302 



28400 DO 325 M=1,MMAX 



28500 IF(I NE . IJET(M)+1) GO TO 325 



28600 IF(THETAO(M) GE .0.0) I SIDE = I JET (M) 



28700 IF(THETAO(M) .LT.0.0) I SIDE=I JET(M)+ 1 



28800 YSEA=0.5*(Y(ISIDE,J)+Y(ISIDE,J+1)) 



28900 YSH0RE=0.5*(Y( ISIDE,J)+Y(ISIDE,J-1)) 



69 



