XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX SOURCE PROGRAM XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



• 88 



C 



c 

 c 



10 

 40 



20 

 100 



30 



1 



1000 



2 



7 

 2000 



200 



5000 

 3000 



400 



1620. FORTRAN II. DISTRIBUTION OF WAVE VEL0CIT1ES(GRID UNITS/ DISTV 01 



SECOND) OVER A GRID OF DEPTHS(FEET OR FATHOMS) AS A FUNCTION DISTV 02 



OF WAVE PERIOD(SEC.) .PROGRAMED BY W.S.WILSON. MAY 19.1964. DISTV 03 



DIMENSION CXX{1025) .CMAT(200) DISTV 04 



READ 40, MM.NN.FMOP.GRID.TT DISTV 05 



FORMAT ( 14.14. F2.0.F7.0.F5. 1) DISTV 06 



XL = 0.5»5.118»(TT**2.0) DISTV OT 



^ = XL DISTV 08 



READ 20. (CXX(K) .<=1.L) DISTV 09 



FORMAT (5<5X.F7.2.3X)) DISTV 10 



PUNCH 100, TT. GRID DISTV 11 



FORMAT (8HPERI0D =.F5.1.17H SEC. .GRID SIZE ■♦F7.0.6H FEET./) DISTV 12 



DO 3000 J=1.NN DISTV 13 



READ 30. (CMAT( I) .I=1.MM> DISTV 14 



FORMAT (10F4.1) DISTV 15 



IF (FMOP) 2.2.1 DISTV 16 



DO 1000 I«1.MM DISTV 17 



CMAT(I) = 6.0*CMAT<I) DISTV 18 



DO 2000 I-l.MM DISTV 19 



NVALUE «= 1 DISTV 20 



IF (CMAT(I)) 3.2000,4 DISTV 21 



CMAT(I) = -(CMATd)) DISTV 22 



NVALUE = 2 DISTV 23 



IF (CMAT(I)-XL) 6.6.5 DISTV 24 



CMATd) = XL DISTV 25 



XK = CMATd) DISTV 26 



< = X< DISTV 27 



CMATd) « CXX(K)/GRID DISTV 28 



CALL ROUND <CMATd ) .l.ES ) DISTV 29 



GO TO (2000. 7), NVALUE DISTV 30 



CMATd) « -(CMATd)) DISTV 31 



CONTINUE DISTV 32 



MM5 = MM/5 DISTV 33 



II » 1 DISTV 34 

 DO 5000 JMal.MMS DISTV 35 

 JX = Il-l DISTV 36 

 JY " J-1 DISTV 37 



III « II+4 DISTV 38 

 PUNCH 200. (CMATd). I«II.Iin»TT.JX»JY DISTV 39 

 FORMAT (5(F10.8.3X) »2X.F5.1.I4,I4) DISTV 40 

 II = III+l DISTV 41 

 CONTINUE DISTV 42 

 CONTINUE _^ DISTV 43 

 TYPE 400. TT DISTV 44 

 FORMAT (F5.1t20H SEC. GRID COMPLETED.) DISTV 45 

 PAUSE DISTV 46 

 PUNCH 300 DISTV 47 



36 



