|
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />f ILF.
<br />fILL
<br />FILE
<br />
<br />C
<br />C
<br />
<br />C
<br />C
<br />
<br />101
<br />
<br />200
<br />C
<br />
<br />23
<br />C
<br />
<br />Appendix C
<br />Program FLOOD
<br />
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />C
<br />
<br />THIS PROG~A~ FI~08 ESTJ~AT~S FUP THt PAPAMfTEPS ALPH4Cl),
<br />ALP~A(2), THETACI), AND THLTA(2)
<br />qfOUIREO !NPlJT IpJCLUDlS ", TH~ ~UM8ER of
<br />YEARS UF THE OA'_ RFCO~~, ,B8, THE MAXIHllM POSSIBLE FLOOD
<br />HEIGHT, AND ee, THE SCAlt FAtTOR. THE 0RDfPFD FLOOD D414,
<br />THE tL(l), AND THf l(l) ^~t READ INTO AR~AYS X, EL, AND ~
<br />RESPECTIV~LY (T~F WCIl .kE C~MPUTED AND 51DOEO I, A~R'Y .
<br />DYPI'G E~lCUTION). THE HI'I"IZATION POOCfSS IS 'CHIEVLD WITH
<br />A SUBROUTINE FROM THE I~SL CIQ77) LIBR,Ry C'LLED ZSYST"
<br />
<br />THIS SUB~uUrIN~ RfGuJPES AN fXTfRN~L FUNCTION (F),
<br />Two CONVERGENCE CRITERI' (EPS AND 'SIG), THE NU"bER ~F
<br />UNKNOWNS (N), IHE HIXIHUh NUMBER OF IlfAATIO'S 0' THE
<br />EXTfRNAL FUNCTION' (IlM'n, ^ .ORK AAEI Of' CDMPlln"
<br />SlOR'GE (W~), AN ARRAV Fl'~ PASSING PAHAHETFkS (PAP, ~HICH
<br />IS NUT USED IN THIS ~TUOY), AN EHPO~ MESS.GE V'.I'~LE (IER),
<br />AND SHRTI"G VALUES ,OR Tllf ALPHAS. THE STAHTI'G VALL'ES
<br />FOR AlPHA(1) ANO ALPH~(~) AR~ C0MPUTED FPO~ THE O~D[QfD
<br />(JATA. OUTPUT CONSrSTS OF ALMPA(l), ALPHA(C'), THF1A.(1),
<br />THETA(2), ITMAX, ANn IfR, TH~ ERROR MESS'GE, IE..O MEANS
<br />THERE ARf NO EPRORS AN" "I"IMIZATION .AS L~HPL~T~O TO T"~
<br />ACCUkACY SPECIFlfD 6Y IHt eONV'RGfNCE CRllfPI',
<br />FDR MORE DETAILED INFO...TION ON THE SUBROUTIN~ ZSYSTEM,
<br />SEE THE IHSL (1077) LI..AR.,
<br />I CKIND'0ISK,TITLE'"(87~073)SAGUENAY/OATA")
<br />2CKIND.DISK,TITLE."(878073)SAGUENAY/EL")
<br />3(KINO'OISK,TITLE."(R70"73)SAGDENAY/.")
<br />
<br />ExTfRNAL ,
<br />DIMENSION ALPHAC2J,.A(Z~),rARI2),XR~G(2Q"),YkEG(200)
<br />CU~MON M,B8,CC,X(C'OO),~(200),EL(20u),ThET~(2),Y(2PO)
<br />tP5=t .O(-q
<br />NSIG'S
<br />N'2
<br />ITMIX'I~O
<br />IEH'O
<br />H--TH~ "UMbER OF YEAR~ OF OAT', Be--THE MAXIMUM POSSIBL~ 'LDun
<br />HEIGHT, AND ce--Thf SCALE 'ACTUH ARE READ.
<br />~EAnIS,Il"
<br />
<br />RfAOIS,/)BB,ee
<br />THE ORDERED OAIA, THE fLII), AND THE Zll) AR~ RfAO INTO ARR..S
<br />X, EL, AND W ~EsPEeTIVfLY,
<br />kEAOll,IOI)(X(I),lol,M)
<br />fOHMATlIX,f12.2)
<br />kFAD(Z,200)(EL(I),I=I,M)
<br />PEAOt3,?O(f) 010 ),l=l,M)
<br />'OA~AT(IX,'40,IS)
<br />THE .(1) IRE CALCULATED.
<br />DO 23 I=I,M
<br />.(I)'I<CIl/.(H)
<br />CON 11 NUE
<br />STARTING VALUES ARE DETfkMI"EO fOH ALPH.C!) ANO ALPHAI2J.
<br />SUMXl=€,.
<br />SUMV1=O.
<br />SUMXY1=O.
<br />SUMXXlcQ.
<br />DO IS 1:1,4
<br />XRfGII)'ALOG(SP_X(t))
<br />YREGII).ALOG(_'LDG(I/(M+I,)))
<br />
<br />31
<br />
|