100 C* Program YSEXP calculates shoreline change according to one line 



110 C* theory, taking into account the e-f-fects o-f a seawall. 



120 INTEGER YSBEG, YSEND 



130 REAL K1.KAP1 



140 DIMENSION Y < 40) , YS ( 40) , Q (41 > , Z (40) , H (40) , Y0 < 40) 



150 DATA YSBEG/26/, YSEND/40/ , DX/50 . /,DT/6. / 



160 DATA DENOM/2. 362 / , NTIMES/44/ , N/40/ , IT1 / 15/ , IT2/31 / 



170 DATA Kl/0. 12/, T/8. 0/ , G/9. 806/ , GAMMA/0. 78/, RADIUS/12000. / 



180 WRITE!*,*) '*******EXPLICIT CALCULATION*******' 



190 WRITE!*,*) 'YSBEG=* ,YSBEG, ' YSEND= ' , YSEND 



200 C* 



210 C* 



220 C* Initialize arrays 



230 C* Straight shoreline 



240 DO 100 1=1, N 



250 Q(I)=0. 



260 Y(I)=0. 



270 100 CONTINUE 



280 DO 105 1=1, N 



290 YS(I)=-7. 



300 105 CONTINUE 



310 Q(N+1)=0. 



320 DCL0S=0. 



330 GOTO 120 



340 C* Curved shoreline 



350 DO 110 1=1, N 



360 BET=ASIN(FL0AT(21-I)*DX/RADIUS) 



370 Y(I)=RADIUS*(1. -COS(BET) ) 



380 Y0(I)=Y(I) 



390 110 CONTINUE 



400 DO 115 I=YSBEG, YSEND 



410 YS(I)=Y(I) -4. 



420 115 CONTINUE 



430 120 CONTINUE 



440 C* 



450 WRITE(*,10) (YS(I) , 1=1, N) 



460 KAP1=K1/ (16.*DEN0M) 



470 C* QL=longshore transport rate over open boundary 



480 C* QL=0 



490 DO 200 IT=1,NTIMES+1 



500 IFdT.EQ.l.OR. IT.EQ.IT1.0R. IT.EQ. IT2) IC = 1 



510 C* Subroutine INDATA computes relevant input wave data 



520 C* at any desired time step. 



530 IF(IC.EQ.l) CALL INDATA ( IT, IT1 , IT2, H, Z, N, DT) 



540 IF(IT.EQ.NTIMES+1) THEN 



550 IC=2 



560 IH0URS=(IT-1)*INT(DT) 



570 WRITE!*, 40) 



580 WRITE(*,*) 'FINAL CONDITIONS (after ' , IHOURS, ' hours)' 



590 ENDIF 



600 IF(IC.GE.l) THEN 



610 WRITE!*, 40) 



620 WRITE!*, 30) (Yd), 1 = 1, N) 



630 WRITE!*, 40) 



640 WRITE!*, 20) (Q(I),I=1,N) 



650 ENDIF 



660 IF(IC.E9.2) GOTO 999 



670 IC=0 



680 DCL0S=2.28*H(1) -68.5*(H(1) /T)**2/G 



690 B=DT*3600. / (2 . *DCL0S*DX) 



700 B2=2.*B 



710 C* 



720 DO 300 1=2, N 



730 ZBS=Z(I) -ATAN! (Yd) -Y(I-l) ) /DX) 



740 Q d ) =H ( I ) **2*SQRT ( G/GAMMA*H d ) ) *KAP1*SIN ( 2*ZBS ) 



A2 



