APPENDIX E 

 USER DOCUMENTATION AND INPUT AND OUTPUT FOR PROGRAM VERIFICATION 



The computer program presented in Appendix B was run on a Burroughs 

 B-7700 computer. The B7000/B6000 series FORTRAN language was designed so 

 several existing programs written in FORTRAN would be compatible with minimal 

 changes. It was designed to be compatible with Fortram IV, H level and to 

 contain ANSI X3. 9-1966 Standard FORTRAN as a subset. 



Line 37,200 of the coding (see App. B) requires a subroutine from the 

 IMSL subroutine package, LEQTIB and its associated subroutines. If the 

 user's computing center has access to this package of subroutine programs 

 they need only bind them to the program (note: copyright laws prohibited the 

 inclusion of the IMSL coding). If not, a substitute subroutine must be user 

 supplied. It must facilitate the solution of a banded storage mode matrix. 



The program input will be described here using a card deck set-up, 

 however, the use of diskpack or magnetic tape input follows directly. Lines 

 3100, 4100, 5500, 5900, 6800, 7500, and 12,900 are read statements. The 

 cards used for the simulation presented in this appendix are shown in Figure 

 E-1. The first card contains the value of WDEPTH, the depth of water (in 

 meters) to which the input wave conditions are to be transformed (a partial 

 list of variables used in the program is presented beginning on page A-8 of 

 Appendix A). The format statements are obviously in the program coding. 



The second data input card is read by line 4100 where the variables 

 SJETTY, BERM, SFACE, and DIAM are required (length of the structure, berm 

 height, shore face slope, and sediment diameter, respectively). 



Lines 5500 reads MMAX, the number of structures to be simulated (as 

 set-up here, a maximum of 10 structures can be modeled, however, appropriate 

 changes in array dimensions would allow additions (structures). Line 5900, 

 which is in a "DO" loop reads the lesser I grid value adjacent to where the 

 structure is desired. The number of structures, MMAX, determines the number 

 of data cards required here; 3 structures require 3 cards with the 3 I grid 

 locations (note, the present configuration of the refraction and diffraction 

 subroutines requires evenly spaced structures, however this can be altered if 

 necessary). 



The parameter ADEAN, which represents the value of A in the equilibrium 

 profile used is the next value input (line 6800). As mentioned previously, 

 whenever possible a site-specific value should be used. The space-step and 

 time-step (DX and DELT in the coding) are input next (line 7500). 



The last input values are the wave data, HS, T, ALPWIS read by line 

 12,900. This statement is in a loop made by the unconditional GO TO statement 

 (line 16,400) and the read statement. There is an action specifier included 

 in the read statement to transfer the program to statement 1000, thereby stop- 

 ping execution of the program once all the wave climate data have been used. 

 The number of data cards required for this read statement is dictated by the 

 length of the simulation and the time-step used. 



The input file and output for program verification follow. 



Ill 



