7300 



7400 



7500 



7600 



7700 



7800 



7900 



8000 



8100 



8200 



8300 



8400 



8500 



8600 



8700 



8800 



8900 



9000 



9100 



9200 



9300 



9400 



9500 



9600 



9700 



9800 



9900 



10000 



10100 



10200 



10300 



10400 



10500 



10600 



10700 



10800 



10900 



11000 



11100 



1 1200 



1 1300 



1 1400 



1 1500 



1 1600 



1 1700 



1 1800 



1 1900 



12000 



12100 



12200 



12300 



12400 



12500 



12600 



12700 



12800 



12900 



13000 



13100 



13200 



13300 



13400 



13500 



13600 



13700 



13800 



13900 



14000 



14 100 



14200 



14 300 



14400 



WRITE(6,772) 

 772 F0RMAT(2X, "READ IN THE SPACE STEP , TIMESTEP " , / ) 



READ(5,775) DX.DELT 

 775 F0RMAT(2(F10.3) ) 



WRITE(6,737) DX 

 737 F0RMAT(2X, 'THE VALUE OF THE LONGSHORE SPACE-STEP. DX= 



WRITE(6,736) DELT 

 736 F0RMAT(2X, 'THE TIME-STEP IN SECONDS, DELT= ',F10.3) 



,F10.3) 



1 1 . . 14. , 17 . .25. .32.808. 10*0.0/ 



0.35. .25. 

 6 . .8. , 10. 



21 , 



.20. . 19. 

 ,14./ 



19/ 



'THE DEPTHS BETWEEN CONTOURS ARE AS FOLLOWS'./) 

 (DEEP( 1 . J). J=1 .dMAX+2) 



DATA CHANGE/1 . .2. . 3. .5. ,7. , 



DO 220 d-1 . JMAX+2 



DO 220 1 = 1 . IMAX 

 220 DEEP( I , J)=CHANGE( J) 



DATA(HC(I ).I=1.8)/1.87.0.5 



DATA(TC(I ) . 1 = 1 .8)/2. .3. .4. . 



DO 200 J=1 .dMAX+2 



DO 200 1=1 . IMAX 

 200 Y(I,J+1) = (0.5*(DEEP(I.J+1)+DEEP(I . d ) )/ADEAN ) ** 1 .5+Y(I . 1) 



WRITE(6.732) 



C*WE WILL ALWAYS REQUIRE Y(I. dMAX+2) TO COMPUTE DY AND YEAR. 



C*WE WILL ALWAYS REQUIRE DEEP( I . dMAX + 2 ) TO COMP SEDIMENT TRANSPORT. 



WRITE(6,734) 



734 F0RMAT(2X. 'THE BOUNDARY Y-VALUES. 1=1. IMAX ARE AS FOLLOWS'./) 

 WRITE(6,801) (Y(1 .d),d=1 , dMAX + 2) 

 WRITE(6.801 ) ( Y( IMAX.d) . d= 1 , dMAX+2 ) 

 WRITE(6,732) 

 WRITE(6,735) 



735 F0RMAT(/,2X. 

 WRITE(6,801 ) 

 WRITE(6. 732) 



801 F0RMAT(2X, 10(F8.2) ) 

 DO 2 1 = 1 , IMAX 

 2 YZEROd )=Y(I , 1 )-(BERM/SFACE) 

 C*WILL COMPUTE THE EQUIL WIDTH BETWEEN CONTOURS. HERE. 

 DO 1 1=1. IMAX 

 WEQd . 1 ) = Y(I . 1 )-YZERO(I) 

 DO 1 d=1,dMAX 

 IF(d NE. 1 ) QO TO 32 

 YTEMP1=0.0 

 GO TO 33 



32 YTEMP1=((0.5*(DEEP(I . d- 1 )+DEEP( I .d) ) )/ADEAN)**1 .5 



33 YTEMP2=( (0.5*(DEEP(I,d)+DEEP(I.d+1)) )/ADEAN)**1 .5 

 WEQd .d+1 )=YTEMP2-YTEMP1 



1 CONTINUE 

 C*LET'S STORE THE ORIG VALUES TO COMPUTE VOL CHANGES OVER CONTOURS . LATER 

 DO 79G 1=1 , IMAX+1 

 YZEROOd )=YZERO(I ) 

 DO 796 d=1, dMAX+2 



796 YORIGd .d) = Y(I ,d) 



Q*********************************************************************** 



C*READ THE DISK FILE AND TRANSFORM PARAMETERS INTO MY UNITS. 

 C*t************«****i I I j ] I 1 I I I I ! I I I I I I I I I ! I I I I ! ! ] ! I 1**** + * + ************* 



C*ALL ADdUSTMENTS TO WAVE ANGLE . HE I GHT . CELERI TY . GROUP VEL. WILL BE MADE 

 C**HERE. AND THRUOUT THE REST OF THE PROG, THEY WILL BE AS IF OCCURRED 

 C***AT WDEPTH! 



798 READ(5.799.END=1000) HS.T.ALPWIS 



799 F0RMAT( 10X.3F6. 1 ) 

 NTIMES=1 



NCHECK=NUNIV+NTIMES 

 HGEN=0.707107*HS 

 SIGMA=TW0PI/T 



G = 32. 17 



C0=G*T/TW0PI 



ELO=CO*T 



IF(T.LE .2 .0) GO TO 797 



HCC=0.23 



DO 444 1=2.7 



T2 = TC(I ) 



IF(T .GT.T2) GO TO 444 



T1=TC( 1-1 ) 



DELTAT=T2-T1 



67 



