|
<br />.
<br />
<br />O(liF;'~5
<br />i.J ,~, '..1 V
<br />
<br />$STORAGE:2
<br />$URGE
<br />$DEBUG
<br />SUBROUTINE NATFLOW(PATHTR)
<br />CHARACTER*lO TlTLE,ATHEAD,AGHEAD,CNAME,ANAME,RDATE
<br />CHARACTER*8 ALHEAD
<br />CHARACTER*s ISTUDY
<br />REAL MM(12)
<br />COMMON CDATA(8s) ,TABLE(10,42,13) ,FLOW(224,13) ,HDATA(42 ,12) ,
<br />+IFRST,IGROUP(40) ,HOLD(224,13) ,ADATA(3s,12) ,ISTART,IEND,J, I,
<br />+NT,NCL,IE,NG,NL,NC,NA,NH,NYI,NRES,TITLE(2,8),ATHEAD(10,8),
<br />+ALHEAD(224 ,4) ,AGHEAD(40,4) ,CNAME(8s,4) ,ISTUDY,ANAME(3s,8) ,R DATE,
<br />+OWN(ll) ,OWNACC(7) ,GAIN(s) ,STORE(7) ,OWNMAX(6) ,OWNPRE(ll) ,PRE ACC(7)
<br />DIMENSION DDATA(28),PATHTR(12)
<br />DATA MM/31. ,30. ,31. ,31. ,28. ,31. ,30. ,31. ,30. ,31. ,31. ,30./
<br />C
<br />c(13): operates natural flow accounting in the period of May-Sept.
<br />LJ - 28
<br />JJ - 77
<br />DO 30 K - 1,28
<br />DDATA(K) - CDATA(K+s3)*MM(J)*.00198347l
<br />30 ONTINUE
<br />50 IF(J.LT.8) GO TO 1000
<br />PATHTR(J) ~ 0.0
<br />C NATURAL FLOW ABOVE ALCOVA
<br />c(14): copmpute total natural flow above Alcova
<br />FLOW(77+l,J) - HDATA(l,J)
<br />FLOW(77+2,J) - HDATA(4,J) + HDATA(6,J)
<br />FLOW(77+3,J) - FLOW(77+l,J)+FLOW(77+2,J)
<br />C STORAGE RELEASE FROM GRAY REEF
<br />c(ls): compute storage release from Gray Reef Reservoir and its
<br />c evaporation loss
<br />FLOW(77+4,J) - HDATA(9,J)
<br />FLOW(77+S,J) - AMINl(FLOW(77+4,J),FLOW(77+3,J))
<br />FLOW(77+s,J) = FLOW(77+4,J) - FLOW(77+s,J)
<br />C CALCULATE STORAGE EVAPORATION
<br />FLOW(77+6,J) = (FLOW(77+S,J)/FLOW(77+4,J))*ADATA(LJ+l,J)
<br />C GLENDO INFLOW
<br />c(16): compute natural flow into Glendo and Glenda storage release
<br />FLOW(77+7,J) - FLOW(77+4,J) + HDATA(lO,J) - FLOW(224,J)
<br />FLOW(77+8,J) - FLOW(77+7,J) - (FLOW(77+s,J)-FLOW(77+6,J))
<br />FLOW(77+8,J) - AMAXl(FLOW(77+8,J) ,0.0)
<br />C STORAGE RELEASE FROM GLENDO
<br />FLOW(77+9.J) = HDATA(ll,J)
<br />FLOW(77+l0,J) = AMINl(FLOW(77+8,J),HDATA(11,J))
<br />FLOW(77+l0,J) = HDATA(ll,J) - FLOW(77+l0,J)
<br />C ADJUST INDICATED STORAGE RELEASE FROM GLENDO BY AMOUNT STORED
<br />C ABOVE ALCOVA
<br />IF (FLOW(77+3,J).GT.FLOW(77+4,J)) THEN
<br />PATHTR(J) - FLOW(77+3,J) - FLOW(77+4,J)
<br />RESDL = AMINl(PATHTR(J),FLOW(77+l0,J))
<br />FLOW(77+l0,J) = FLOW(77+IO,J) - RESDL
<br />PATHTR(J) ~ PATHTR(J) - RESDL
<br />
<br />25
<br />
|