Program Listing 



7. Below is a program listing for computer program BWL0SS1 (FORTRAN 



version) . 



laC PROGRAM "BWLQSSl". 11/B5 VERSION 



20C DESIGN BRANCH-COASTAL ENGINEERING RESEARCH CENTER 



30C U.S. ARMY ENGINEERS WATERWAY EXPERIMENT STATION 



40C P. 0. BOX 631 



SBC VICKSBUR6, MS 39180-0631 



60C FOR FURTHER INFORMATION CONCERNING THE APPLICATION 



70C OF "BWLOSSl", CALL.. 



80C ORSON P. SMITH (601 ) -634-2013 FTS:542-2013 OR 



90C ROBERT B. LUND (601 ) -634-2068 FTS:542-2068 OR 



100C DOYLE L. JONES (601 ) -634-2069 FTS:542-2069 



110C 



120C FORTRAN 4 HONEYWELL DPS-8 



130C REF: "COMPUTER PROGRAM WAVDIST" CETN-I- 



140C REF: "PROBABILITY AND STATISTICS" BY MORRIS DE6R00T 



150C REF: "COST EFFECTIVE OPTIMIZATION OF RUBBLE-MOUND BREAKWATER 



160C CROSS-SECTIONS" BY ORSON P. SMITH 



170C REF: "EXTREMAL STATISTICS IN WAVE CLIMATOLOGY" BY BORGMAN AND RESIO 



180C 



190C N = THE NUMBER OF DATA POINTS 



200C X = THE ARRAY OF SIGNIFICANT WAVE HEIGHTS 



210C YH = THE ARRAY OF LOSSES CORRESPONDING TO EACH SIGNIFICANT WAVE HEIGHT 



220C Y » THE TRANSFORMED Y ARRAY USED IN THE METHOD OF LEAST SQUARES 



23BC VI = Hlo, THE MAXIMUM WAVE HEIGHT FOR WHICH LOSSES ARE NEGLI6BLE 



240c W = $Lmax, THE MAXIMUM CONCEIVABLE LOSS IN MILLONS OF DOLLARS 



250C V2 = A, THE REGRESSION COEFFICIENT A<0 



260C CORR = THE NON-LINEAR CORRELATION OF THE LOSS FUNCTION 



270C ST = THE SUM OF THE SQUARE RESIDUALS 



280C Dl = THE ARRAY THAT CONTAINS THE RESIDUAL FOR EACH DATA POINT 



290C Zl = THE LOSSES AS ESTIMATED BY THE LOSS CURVE 



3000 PDF(X) = THE PROBABILITY DENSITY FUNCTION OF EXTREMAL WAVES 



310C CDF(X) = THE CUMULATIVE DISTRIBUTION FUNCTION OF EXTREMAL WAVES 



320C G(X) = THE LOSS FUNCTION 



330C C = SIMPSON'S RULE COEFFICIENTS 



340 



350C INITIALIZE VARIABLES , STRI NGS , AND FUNCTIONS 



360 DIMENSION X ( 101 ) , Y (10 1 ) , YH ( 101 ) 



370 COMMON X,Y,YH 



380 CHARACTER»1 L 



390 CHARACTER»60 ST(20) 



400 G(X)=W«(1-EXP(V2*(X-V1) ) ) 



4J0 ST(1)='«*»**«***»»****»****»*»******«**** *♦♦*»♦♦♦** *♦****♦*«♦»♦ 



420 ST(2)='* "BWLOSSl" IS A PROGRAM WHICH FITS AN EXPONENTIAL ♦ 



430 ST(3)='* CURVE TO HISTORICAL INFORMATION ON ECONOMIC LOS- ♦ 



440 ST(4)='» SES CAUSED BY WAVE ATTACK, EACH STORM CAUSING ♦ 



450 ST(5)='* LOSSES IS ASSUMED TO BE CHARACTERIZED BY A SIN- ♦ 



460 ST(6)='* GLE SIGNIFICANT WAVE HEIGHT. THE PROGRAM RE- » 



470 ST(7)='* QUIRES ESTIMATES OF THE MAXIMUM LOSS SUSTAINABLE ♦ 



480 ST(8)='* FROM WAVE ATTACK, THE MAXIMUM SIGNIFICANT WAVE ♦ 



490 ST(9)='» HEIGHT FOR WHICH LOSSES CAN BE NEGLECTED AND AT » 



500 ST(10)='» LEAST ONE HISTORICAL LOSS WITH ASSOCIATED WAVE ♦ 



510 ST(11)='* HEIGHT. THE EXPONENTIAL CURVE IS COMPUTED BY ♦ 



A5 



