750 300 CONTINUE 

 760 C* Boundary conditions: 

 770 C# Pinned beach 

 780 Q41)=Q(2) 



790 C* Groin(s) 



800 



C* 



QU)=0. 



810 





Q(N+1)=0. 



820 



C* 





830 





IF (YSBEG.GE.3) THEN 



840 





DO 400 I=l,YSBEG-2 



850 





Y(I)=Y(I) -B2*(Q(I+1) 



860 



400 



CONTINUE 



870 





ENDIF 



880 C* Correction o-f shoreline in -front o-f seawall i -f necessary 



890 CALL C0RRE(YSBEG,YSEND,Q,B2,Y,YS) 



900 IF(YSEND.NE.N) THEN 



910 DO 500 I=YSEND+1,N 



920 Y(I)=Y(I> -B2*(Q(I+1) -Q(I> > 



930 500 CONTINUE 



940 ENDIF 



950 C* 



960 C* 



970 C* Error calculation (DIFF: closed boundaries, AR0UT: open boundary) 



980 C* QL=QL+Q(1) 



990 200 CONTINUE 



1000 999 CONTINUE 



1010 DIFF=0. 



1020 AAREA=0. 



1030 DO 600 1 = 1, N 



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



1050 AAREA=AAREA+ABS(YO(I)-Y(I) ) 



1060 600 CONTINUE 



1070 ERR0R=DIFF/AAREA 



1080 C* AR0UT=QL*DT*3600. /DCL0S-DIFF*DX 



1090 C* ERR0R=AR0UT/AAREA 



1100 C* 



1110 C*** Output *** 



1120 WRITE(*,*) 



1130 WRITE!*,*) 'LOST SAND V0LUME=' , ERR0R*100, ' %' 



1140 C* WRITE!*,*) ' (QL*DT/D-AREA) / ABSAREA*100=' ,ERR0R*100, '%' 



1150 10 FORMAT (IX, 'SEAWALL POSITION '/( IX, 10F8 . 2) ) 



1160 20 FORMAT (IX, 'LONGSHORE TRANSPORT */( IX , 10F8. 4) ) 



1170 30 FORMAT (IX, 'SHORELINE POSITION' /( IX, 10F8. 2) ) 



1180 40 F0RMAT(//> 



1190 STOP 



1200 END 



A3 



