100 C* Program YSIMP is an implicit version of program YSEXP and 



110 C* calculates shoreline change according to one line theory, 



120 C* taking into account the e-f-fects o-f a seawall. 



130 INTEGER YSBEG, YSEND 



140 REAL K1.KAP1 



150 DIMENSION Z (40) , Y (40) , YS (40) , Y0 (40) , Q (41 ) , YC0LD<40> ,E(40) ,F(40) 



160 DIMENSION EP (40) , FP (40) , BP (40) , P (41 ) ,R (41 ) , QQ (41 ) , H (40) 



170 DATA YSBEG/l/,YSEND/40/,DX/50. /,DT/6. / 



180 DATA DENOM/2.362 / , NTIMES/56/ , N/40/ , IT1 /22/ , IT2/45/ 



190 DATA Kl/0. 12/, T/8.0/.G/9. 806/, GAMMA/0. 78/, RADIUS/ 12000. / 



200 WRITE!*,*) '******* IMPLICIT CALCULATION*******' 



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



220 C* 



230 C* Initialize arrays 



240 C* Straight shoreline 



250 DO 100 1=1, N 



260 Q(I)=0. 



270 Y(I)=0. 



280 100 CONTINUE 



290 DO 105 I=YSBEG, YSEND 



300 YS(I)=-7. 



310 105 CONTINUE 



320 Q(N+1)=0. 



330 D0LD=0. 



340 DCL0S=0. 



350 C* GOTO 120 



360 C* Curved shoreline 



370 DO 110 1=1, N 



380 BET=ASIN (FLOAT (21-1 )*DX/RADIUS) 



390 Y(I)=RADIUS*(1.-C0S(BET) ) 



400 Y0(I)=Y(I) 



410 110 CONTINUE 



420 DO 115 I=YSBEG, YSEND 



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



440 115 CONTINUE 



450 120 CONTINUE 



460 C* 



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



480 KAP1=K1/ (16.*DEN0M) 



490 C* QL=longshore transport rate over open boundary 



500 C* QL=0. 



510 C* C=correction term in continuity calculation 



520 C* C=1.0 



530 DO 200 IT=1,NTIMES+1 



540 IFdT.EQ. 1. OR. IT.EQ. IT1. OR. IT.EQ. IT2) IC = 1 



550 C* Subroutine INDATA computes relevant input wave data 



560 C* at any desired time step. 



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



580 IFdT.EQ. NTIMES+1) THEN 



590 IC=2 



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



610 WRITE(*,40> 



620 WRITEf*,*) 'FINAL CONDITIONS (a-fter ', I HOURS,' hours)' 



630 ENDIF 



640 IF(IC.GE.l) THEN 



650 WRITE(*,40) 



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



670 WRITE!*, 40) 



680 WRITEf*, 20) (Q(I),I=1,N) 



690 ENDIF 



700 IF(IC.EQ.2) GOTO 999 



710 IC=0 



720 D0LD=DCL0S 



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



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



A6 



