1570 B(l) = (SX)(tSLLY-S)(LLY*S)l)/(NtSl(X-Sl(«»2) 



1580 B(2) = (SLl(l(tSLLQY-T00BlG»SL){)/(NtSLn-SL)(*»2) 



1590 B(3) = (SL)ll(tSLLY-SLl(LLYtSLK)/(N*SLn-SL)(M2) 



16B0C CftLCULATE PARAMETERS OF EACH DISTRIBUTION FROM SLOPE AND INTERCEPT DATA 



1610 PHI=1.0/A(1) 



1620 EPSI=-B(1)/A(1) 



1630 C=A(2) 



1640 SI6«A=E1(P<-B(2)/A(2)) 



1650 U=A(3) 



1660 SI6HA2=E)!P(-B(3)/A(3)) 



1670 



16B0C ASSIGN ARRAYS ALPHA AND BETA THE PARAMETERS OF EACH DISTRIBUTION 



1690C FOR EASY PRINTOUT OF DATA 



1700 ALPHA! 1)=EPSI 



1710 BETA(1)=PHI 



1720 ALPHA (2) =C 



1730 BETA(2)=SI6HA 



1740 ALPHA! 3) =U 



1750 BETA(3)=SIG«A2 



1760C CALCULATE APPROXIMATE PROBABILITY AS ESTIMATED BY DISTRIBUTION 



1770 DO 100 J=1,N 



1780 YEST(J,1)=F1(HS(J)) 



1790 YEST(J,2)=F2(HS(J)) 



1800 YEST(J,3)=F3(HS(J)) 



1810 100 CONTINUE 



1820 



1830C CALCULATE AVERAGE PROBABILITY AND CORRELATION COEFFICIENTS 



1840 DO 110 K=l,3 



1850 YAVG(K)=SY/FLOAT(N) 



1860 ST(K)=0 



1870 S6(K)=0 



1880 110 CONTINUE 



1890 



1900 DO 120 K=l,3 



1910 DO 130 1=1, N 



1920 ST(K)=ST(K)+(YACT(I,K)-YEST(I,K))»«2 



1930 SB(K)=SB(K) + (YACT(I,K)-YAV6(K))H2 



1940 130 CONTINUE 



1958 CORR(K)=1.0-ST(K)/SB(K) 



1960 120 CONTINUE 



1970 



1980C CALCULATE DATA FOR RETURN PERIOD TABLES 



1990 DO 57 J=l,5 



2000 PROB=1.0-1.0/(LAMBDA*RET(J)) 



2010 IF(PROB .LE. 0) PROB=. 0000001 



2020 CHS(J,1)=-AL0G(-AL06(PR0B))»PHI+EPSI 



2030 CHS(J,2)=(-ALOG(1.0-PROB) )♦»(!. 0/C)tSIGMA 



2040 CHS(J,3)=SI6MA2/((-ALDS(PROB))t»(1.0/U)) 



2050 57 CONTINUE 



2060 



2070C CALCULATE MEAN AND VARIANCE FOR EACH DISTRIBUTION 



2080 MEAN(l)=EPSI+EULERtPHI 



2090 VAR(1)=1.6449341«PHI*»2 



C9 



