|
<br />15
<br />
<br />SUMXl=SUMXI+XREGII)
<br />SUMYlcSUMYl+YRf.G(I)
<br />SU"XYI=SUMX,,+XRlG(I)"'Er.lI)
<br />SUMX~1=SUMtX'+1PfG(I)**?
<br />CONTINUE
<br />X16AR=SUMXI/4.
<br />YIBAR=SlJMYl/..
<br />ALPH.(1)=(SUMXYI-~.*~1~~R*YIBAR)/lSUMXXI_U..~IB~R.*2)
<br />SUMx<=n.
<br />SUMY2=O.
<br />SUM:lCY2=O.
<br />SUM02=0.
<br />DU Ih J=~-~,M
<br />X"FGIJ) =A'LuG lRb.. lJ))
<br />Yk!GlJ)=ALuG(.ALOGIJ/("+l.)))
<br />SUMX2=SUMX2+XkFG(J)
<br />SUM'2=SUMY2+YReGlJ)
<br />SUMXY2=SUMXY2+XRFGlJ)',R,GIJ)
<br />SU~XX2:SUMxX2+XReG(J)**2
<br />CONTINUE
<br />X28AR=SUMX21..
<br />Yi'8A.=SUMY21..
<br />ALPHA(2)=(SUMXY2-4.*~2bAK*Y2BAR)/(SUM~A2_~..X2~AR**2)
<br />wRITt(~,50)M,B~,CC
<br />fVRMAT(IX,'THE NUhBfR OF YEARS Of THE DATA R~CORD=',115,11,
<br />.IX,'THE MAXIMU~ POSSl8Lf FLOOD HfIGHT=',I15,II,t_,
<br />* 'THE SCALe FACTor=',I15,11)
<br />.RIIEIb,bO)ALPHAr!),.LPMAI2)
<br />fOR~AT(lX,ITHE STAPTl~G VALUES Akfl,l,lX,'ALPHAC1J:t,F15.S,
<br />.. sx, 'ALPHA(2)=',FI5.5,11111)
<br />ZSYS!fl IS CALLED TO MINIMIZE FQUATIQN (20) AND OuTpUT THE
<br />FSTIMATED PARA"'ET!RS
<br />CALL ZSYSTM(F,EP5,N~iG,N,ALPHA,IT~AX/~A,PAR,IE~)
<br />,RITElb,70)ITM'X,I!R
<br />FOR~AT(lX,fNUM6EP or JTFRATIu~s OF [XTERt~.L fU~CTIO~=',J5,
<br />. I/,IX,IERPOR H[,SSAGE=',IS,I/II/)
<br />.R!TUb,Bo)
<br />FOR"'ATI1X,'PARAMEH'R !STI"ATES ARE',11l
<br />.RITE(b,QO)ALPHA(ll,ALFHA(2),THETA(I),THETAI?)
<br />FORH~T(lx,'ALPHA(1)=I,F20.IO,5x,'~LP~A(2)=r,F2Q.I~,II,Ix,
<br />* tTHtTA(I)~I,F20.IO,5X,ITHF.TA(2)=I,F20.JO)
<br />STOP
<br />END
<br />
<br />1 ~
<br />
<br />50
<br />
<br />bO
<br />
<br />C
<br />C
<br />
<br />70
<br />
<br />BO
<br />
<br />qO
<br />
<br />IIHgHOIOHIOOO#HHOHOOHOOHO#H##OO###OHO######O#U#O##########I##HOO#####O######
<br />
<br />fUNCTION FlALPHA,K',P.k)
<br />C THIS fUNCTION EVALUATES l"'E T"O [QUATIONS IN T"O U~KNO.NS.
<br />DIMENSION .EL(200),.L~12no),~ELXI200)"LPHA(2),PAR(2)
<br />tO~MON M,~O,CC'X(200),~(20P),EL(20~),THET'(2),Y(2QO)
<br />00 10 1=1,'"
<br />.ELll)=.II)'ELII)
<br />YlI)=I~B-~II))/CC
<br />.LXII)="II)*ALD~IY(I))
<br />10 .ELXlI)="LXCI)*ElIJI
<br />A3=2.*ALPMA(I)IA.=?.ALP.AI2)/'5=ALPHAII)t'LPHAlZ)
<br />Zl=O.IZ2=O.;Vl=O.'Y2=o.;Y3=".,y~c~,'YS=o.;Y~=o.;Y7=O.'Y~=o.
<br />Bl=0./BZ=0./B3=0./B"=O./P5=0.
<br />00 20 I=t,M
<br />YAI=Ylll"ALPHA(I)
<br />YAZ=YII)..ALFHA(ZI
<br />
<br />32
<br />
|