750 B0LD=DT#3600. / <2.*DOLD*DX) 



760 YC0LD(1)=Y(1) +B0LD*(Q<1) -Q(2) ) 



770 C# Boundary conditions: 1 



780 C* Groin causing 8(1 )=0. 



790 E(1)=0. 



800 F(1)=0. 



810 C* Pinned beach as S(1)=Q<2) 



820 C* E(l)=l. 



830 C* F(1)=0. 



840 C* 



850 DO 300 1=2, N 



860 YC0LD(I)=Y(I>+B0LD*(Q<I)-Q(I+1) ) 



870 ZS=ATAN( (Y(I)-Y(I-l) WDX) 



880 Z2=2.*Z(I) 



890 PWR=H(I)**2*SQRT(G/GAMMA*H(I) ) 



900 EP(I)=PWR*KAP1*2*C0S(Z2)*<C0S<ZS) )**2/DX 



910 FP(I)=PWR#KAP1*SIN<Z2>*<2*(C0S<ZS) )**2-l. ) 



920 BP<I)=B*EP(I> 



930 DEN=1.+BP(I)*<2. -E(I-l) ) 



940 E(I)=BP(I) /DEN 



950 F<I)=<FP<I)+EP<I)*(YC0LD(I-1>-YC0LD<I) ) +BP ( I ) *F ( I - 1 ) ) /DEN 



960 300 CONTINUE 



970 C* Boundary condition 2: groin 



980 Q(N+1)=0. 



990 C* 



1000 DO 400 I=N,1,-1 



1010 Q(I)=E(I)*Q(I+1) +F(I) 



1020 IF(IT.EQ.l) THEN 



1030 YC0LD(I)=Y(I) +B* (9 ( I ) -Q ( I +1 ) ) 



1040 ENDIF 



1050 400 CONTINUE 



1060 C*** Reversed double sweep *** 



1070 C* Boundary conditions 3: groin 



1080 P(N+1)=0. 



1090 R(N+1)=0. 



1100 C* 



1110 DO 500 I=N,2,-1 



1120 P(I)=BP(I) / <1.+BP(I>*<2. -PII+1 > ) > 



1130 R(I)=(FP(I) +EP<I)*(YC0LD(I-1> -YCOLD(I) ) +BP ( I ) *R ( I +1 ) ) / 



1140 & (l.+BP(I)*(2.-P(I+l) ) ) 



1150 500 CONTINUE 



1160 C* Boundary condition 4 (alt 1: closed boundary, alt 2: open) 



1170 QQ(1)=0. 



1180 C* QQ(1)=R(2) / (1. -P(2) ) 



1190 C* 



1200 DO 550 1=2, N+l 



1210 9Q(I)=P(I)*QQ(I-1) +R(I) 



1220 CHECK=ABS(QQ(I) -Q ( I ) ) 



1230 C# IFfCHECK.GT. 0.0005) WRITE(*,*) 'TRANSPORT CALC. DIFFER' 



1240 550 CONTINUE 



1250 C* 



1260 C*Correction of shoreline in front of seawall if necessary 



1270 CALL C0RRI (YSBEG, YSEND,Q,B, YCOLD, E , F, P, R, Y, YS,N) 



1280 C* 



1290 C* 



1300 C* Error calculation (DIFF: closed boundaries, AROUT: open boundary) 



1310 C* IF(IT.EQ.NTIMES) C=0.5 



1320 C* QL=QL+C*Q(1> 



1330 200 CONTINUE 



1340 



999 



CONTINUE 



1350 





DIFF=0. 



1360 





AAREA=0. 



1370 





DO 600 1 = 1, N 



1380 





DIFF=DIFF+YO(I) -Y(I) 



1390 





AAREA=AAREA+ ABS ( YO ( I ) 



1400 



600 



CONTINUE 



A7 



