<br />1/100.178.,S8.147.,34.1Z4./
<br />DATA ~(IZ,l),"(IZ,Z),~(IZ,3),F(IZ,4),F(lZ,S),F(lZ,6)
<br />1/100.,91., 7S. ,6Z. ,48 '134.
<br />DATA F(13,1),F(13,Z),"{13,3),F(13,4),F(13,S),F(13,6)
<br />1/100.194.,B3.174.,60.144.
<br />DATA ~(14,1),"(14,Z),~(14,3),F(14,4),F(14,S),F(14,6)
<br />1/100.,9S.,B7.L81.,69.,S3.
<br />DATA F(IS,I),~(IS,Z),F(IS,3),F(15,4),F(IS,S),F(IS,6)
<br />1/100.,97.,91.,86.,76.,60.
<br />END
<br />c*********************************************************************
<br />C **************************** SUBROUTINES ***************************
<br />c*********************************************************************
<br />C
<br />C
<br />C
<br />C
<br />C********************************************** DEPTH CALCULATION ROUTINE
<br />SUBROUTINE STAGE(K,M,N,S2,Q,H3,Q3,Hl)
<br />C INPUT PARMS ARE:
<br />C K,M =CHNL. FIT PARMS, Q3=QV, Q-FLOW, SZ=SLOPE, N=MANNING N
<br />C
<br />REAL K KO K1,M MO,Ml,N
<br />COMMON'ARAY/ F(lS,6),X(6),D(8),E(B,8),T(8,B),C(3,B),H(Z,B),B(4,8),
<br />IP(8,S)
<br />CO~~ON/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,Nl,SO
<br />Al = (1.49/N)*(SZ**0.S)*K/\(M+l.J**1.6b7)
<br />Bl=1.667+M
<br />IF(Q .GT. Q3) GO TO SO
<br />Hl=(Q/A1)**(1./Bl)
<br />GO TO 100
<br />R1=Al*((M+1.)**1.667)*(H3**M)
<br />Rl=(1/Rl)**.6
<br />Hl=(Q**.6)*Rl+H3*(M/(M+l.))
<br />100 RETURN
<br />END
<br />
<br />C
<br />
<br />CASE 1 STG < HV
<br />
<br />C
<br />SO
<br />
<br />CA1SE Z STG > HV
<br />
<br />C
<br />C********************************************** HV CALCULATION ROUTINE
<br />SUBROUTINE COMPHV(I,H3)
<br />C THIS ROUTINE COMPUTES THE WEIGHTED AVERAGE HV ELEVATION OF REACH
<br />REAL K KO Kl,M MO,Ml,N
<br />CO~~ION'AiAy/ F(IS,6),X(6),D(B),E(B,8),T(B,8),C(3,8),H(Z,B),B(4,8),
<br />IP(8,S)
<br />CO~~ON/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,N1,SO
<br />H3=0
<br />DO 100 J1 = 2,1 ( ,
<br />H3 = H3+((C(Z,Jl-1)+C(Z,J1))/Z.)*(D(J1)-D(Jl-1))
<br />100 CONTINUE
<br />H3=H3/(DlI)-D(l))
<br />RETURN
<br />END
<br />
<br />C
<br />C********************************************** QV CALCULATION ROUTINE
<br />SUBROUTINE COMPQV(K M N SZ H3 Q3)
<br />C THIS ROUTINE COMPUtES tHt FLOW QV ASSOC. WITH HV
<br />REAL K KO Kl,M MO,Ml,N
<br />CO~~ON'ARAY/ F(lS,6),X(6),D(8),E(8,8),T(8,8),C(3,B),H(Z,8),B(4,B),
<br />IP(B,S)
<br />COMMON/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,N1,SO
<br />C INPUT PARMS ARE:
<br />C K,M =CHNL. FIT PARMS, SZ=SLOPE, N=MANNING N, H3=IlV
<br />Al=(1.49/N)*(SZ**.S)*K/((M+1.)**1.667)
<br />B1=1.667+M '
<br />Q3=A1 *(H3**Bl)
<br />RETURN
<br />END
<br />
<br />C
<br />C****************************************************BREACH ROUTINE
<br />SUBROUTINE DAMBRK(K,M,N,SZ,HZ,H3,H4,Q3,Q,Hl)
<br />C THIS ROUTINE DET. THE MAX OUTFLOW AT DAM, ACCTNG FOR SUBMERGENCE
<br />REAL K KO Kl,M MO,Ml,N
<br />CO~~ON'ARAY/ F(lS,6),X(6),D(B),E(B,B),T(8,8),C(3,8),H(Z,8),B(4,81,
<br />II-I0
<br />
|