|
<br />IF(FB .EQ. 3) DO=DO/3.
<br />
<br />C
<br />C CALL ROUTINE FOR DOWNSTREAM PEAK
<br />FO=F8
<br />VS=V8
<br />CALL PREDIC(X2,FO,VS,Q,QB)
<br />n:~~
<br />CALL PREDIC(X2,FO,VS,Q,QB)
<br />QZ=QB
<br />C INTERPOL FOR SPEC FC AND LOWER V*
<br />Q4=Q1+lQ2-Ql)*DO
<br />FO=F8
<br />VS =V9
<br />CALL PREDIC(XZ,FO,VS,Q,QB)
<br />Q1=Q8
<br />FO=F9
<br />CALL PREDIC(X2,FO,VS,Q,QB)
<br />C INTE~~-~gR SPEC FC AND UPPER V*
<br />QS=Q1+(QZ-Ql)*DO
<br />C INTERP FOR SPECIFIC V*
<br />Q7=Q4+(QS-Q4)*(V7-FLOAT(VS))
<br />IF(V7 .LT. 1.0) Q7=Q4
<br />
<br />FLOWS
<br />
<br />C
<br />
<br />RETURN
<br />END
<br />
<br />C
<br />C****************************************SUBROUTINE PREDIC
<br />C
<br />SUBROUTINE PREDIC(X2,FO,VS,0,QB)
<br />C THIS ROUTINE INTERPOLATES WITH!N THE ROUTING CURVES TO
<br />C FIND THE RATIO OF QP TO QBMAX FOR A SPECIFIC (INTEGER)
<br />C VALUE OF V* AND AN FC OF O.ZS, O.SO, OR 0.7S.
<br />C
<br />C INPUT VARIABLES ARE:
<br />C FO = Fe POINTER, VS = V* POINTER, X2 = X*, Q a PEAK FLOW AT DAM
<br />C RETURN VALUE IS Q8
<br />C
<br />INTEGER VS FO
<br />CO~illON/ARAY/ F(lS,6),X(6),D(8),E(8,B),T(8,B),C(3,B),H(Z,S),B(4,8),
<br />IP(B,S)
<br />CO~~IDN/SCALAR/ CO,BO,VO,AO,WO,TO,QO,NO,N1,SO
<br />DO SO 11 = 1,6
<br />IZ=11
<br />X3=X(11)/10.
<br />'IF((XZ-X3) .LE. 0.00) GO TO 100
<br />SO CONTINUE
<br />C BRACKETING X'S FOUND
<br />100 13=IZ-1
<br />JO=O
<br />IF(FO .GT. 1) JO=JO+S
<br />IF(FO .GT. Z) JO=JO+S
<br />J3=JO+VS
<br />RL1=ALOG!F~J3'I3))
<br />RL2=ALOG F ,J3lIZ))
<br />RL3=ALOG X 13 )
<br />RL4=ALOG X IZ )
<br />RLS=ALOG X2*10.)
<br />D3=(RLZ-RL1)/(RL4-RL3)
<br />RL6=(RLS-RL3)*D3 + RL1
<br />IFl(XZ*10.) .GT. X(6)) RL6 = (RLS-RL4)*D3 + RLZ
<br />IF(RL6 .LT. -3.0) RL6 = -3.0
<br />Q8 = (EXP(RL6)/100.) * Q
<br />RETURN
<br />END
<br />
<br />C
<br />
<br />1I-13
<br />
|