100 SUBROUTINE CORRI ( YSBEG, YSEND, Q, B, YCOLD, E, F, P, R, Y, YS, N) 



110 C* CORRI recalculates transport rates <Q) due to limited sand 



120 C* volume in front of a seawall and adjusts the shoreline 



130 C* position as necessary. Implicit calculation scheme. 



140 C* 



150 INTEGER YSBEG, YSEND 



160 REAL Q(41) ,Y(40) , YS(40) ,YC0LD(40) 



170 REAL E(40) ,F(40) ,P(41) ,R(41) 



180 C* 



190 I=YSBEG 



200 IF(Q(I) .GT.0) THEN 



210 C* 



220 C* Q positive: Calc of shoreline with correction of Q and Y 



230 C# as necessary. 



240 



C* 





250 



10 



IF(Q(I+1) .GE.0) THEN 



260 





Qd + l)=Pd + l)*9d)+Rd + l) 



270 





YC=2*B#(Qd> -9(1+1) )+YC0LD(I) 



280 





IF(YC.LT.YSd) ) THEN 



290 





DQ=(YS(I) -YCOLD (I) ) / (2*B) 



300 





Q(I+1)=Q(I) -DQ 



310 





ENDIF 



320 





Yd)=B*(Qd)-Qd + l> >+YC0LDd> 



330 





1 = 1 + 1 



340 





IFd.EQ. YSEND + 1) GOTO 100 



350 





GOTO 10 



360 





ENDIF 



370 





K=I 



380 





1 = 1 + 1 



390 





IFII.EQ. YSEND+1) THEN 



400 





Yd-l)=B*(Qd-l)-Qd) ) +YC0LDII-1) 



410 





GOTO 100 



420 





ENDIF 



430 





IFd.EQ. YSEND) THEN 



440 





1 = 1 + 1 



450 





GOTO 30 



460 





ENDIF 



470 





ELSE 



480 





K=YSBEG 



490 



C* 





500 



C* 



Q negative: Search for minus piont. 1+ absent, 



510 



C# 



for right end element. Correct Q as necessary. 



520 



C* 





530 





ENDIF 



540 



20 



IF(Q(I+1) .LT.O) THEN 



550 





1 = 1 + 1 



560 





IFd.EQ. YSEND) THEN 



570 





IF(Q(I+1) .LE.O) THEN 



580 





YC=2*B*(Qd) -Qd + 1) ) +YC0LD(I) 



590 





IFIYC.LT. YS(I) ) THEN 



600 





DQ= ( YS ( I ) -YCOLD (I ) ) / <2*B) 



610 





Q(I)=Q(I+1) +DQ 



620 





ENDIF 



630 





Yd)=B*(Q(I) -Qd + 1) > +YC0LD(I) 



640 





GOTO 30 



650 





ENDIF 



660 





ENDIF 



670 





GOTO 20 



680 





ENDIF 



690 



C* 





700 



c* 



Minus point: Corr of Q out o-f the element if s 



710 



c* 



behi nd seawal 1 . 



720 



c* 





730 





YC=2*B*(Q(I)-Q(I+1) ) +YC0LDII) 



740 





IF<YC.LT. YSd) ) THEN 



shoreline moves 



A9 



