ifwindo. .windowing flag: [0]=no windowing of data segments, 



[1]=segments windowed (Kaiser-Bessel window) 

 ifdtrnd..det rending flag: [0]=no detrending, [1] = linear trend 

 removed from data segments 

 nfft..# of points in each data ensemble 

 nensb..# of half- lapped segments of cross-spectral computations 

 nband..# of raw frequency bands averaged in frequency smoothing 

 idgfr. .degrees of freedom in cross-spectral computations 

 (based on contiguous segments only) 



nfrq.. number of output frequency bands, equals range of index 

 nf 

 detfs..[Hz] output frequency bandwidth 

 nang. .number of output angle bins, equals range of index na 

 delang. . [deg] output angle bin width 



dmin..[m] minimum ensemble segment water depth at reference 



gauge 'rname' during collection 

 dbar..[m] mean water depth at gauge 'rname' during collection 

 dmax..[m] maximum ensemble segment water depth at reference 



gauge 'rname' during collection 

 rname. . [character*5] reference gauge id for depth computations 



character*5 rname 

 character*6 gpat(13) 

 character*10 indattim, datetime 

 character*80 infile 



dimension f(13), sf(13), iter(13) 

 dimension angle(181), ddf(13,181), fds(13,181) 

 c 



c get file-naming date/time group from user 

 c 



writer.*, ' (2x, ' 'Enter date/time group (yymmddhhmm). ..'')' ) 

 read(*, '(a10)') indattim 'date/time string 



c 



c define input data file 

 c 



infile='HP'//indattim(1:10)//'.ASC 

 c 



c open, read, and close data file 

 c 



open(10,f i le=inf i le,status='old' ,form=' formatted' ) 

 c 



readCIO, 

 & '( a10, f10.2, f10.5, f10.1, i10, i10, 

 & f10.5, ilO,/, i10, 1 10. i"10, i 10, 

 & i10, i 10, f 10.5, i 10,/, f 10. 1 , f10.2 r 

 & f 10.2, f 10.2, 5x,a5)') 



& datetime, Hmo, fp, thp, ifimle, istot, 

 & sfrq, ifwindo, ifdtrnd, nfft, nensb, nband, 

 & idgfr, nfrq, delfs, nang, delang, dmin, 

 & dbar, dmax, rname 

 c 



read(10,'(10f8.1)') (angle(na),na=1,nang) 

 c 



do 10 nf=1,nfrq 

 readOO, 

 & '( i 10, f10.5, f10.6, 4x,a6, i10)') 



& if, f(nf), sf(nf), gpat(nf), iter(nf) 



read(10,'(8f10.7)') (ddf (nf ,na),na=1,nang) 

 1 cont i nue 

 c 



close(10) 



Figure C1. (Sheet 2 of 3) 



C2 



Appendix C Listing of FORTRAN Computer Program 



