Laserfiche WebLink
<br />OD(J52$ <br /> <br />SUMDEM = HDATA(15.J)+HDATA(16,J)+HDATA(17,J)+HDATA(lB,J) <br />SUMDEM - SUMDEM+HDATA(20.J)+HDATA(22.J)+HDATA(26.J)+HDATA(27.J) <br />SUMDEM = SUMDEM+HDATA(2B,J)+HDATA(29.J)+HDATA(30.J)+HDATA(3l.J) <br />SUMDEM - SUMDEM+HDATA(39,J)+HDATA(40,J) <br />KTDEM(B) - HDATA(19,J)+HDATA(2l.J)+HDATA(32,J)+HDATA(34.J) <br />KTDEM(B) - KTDEM(B) + AMAXl(O. .SUMDEM-SUMlNF) <br />REi1FL - O. <br />C ************************************************************************ <br />C MEET REACH DEMANDS ,FILL THE STORAGES BY INFLOW, DECREMENT <br />C BY EVAPORATION AND COMPUTE RESIDUAL <br />c(5): Meet reach demand, compute residual demand and remaining flow, <br />c record reach flow, compute evaporation loss, fill reservoir by <br />c remaining inflow and make minimum release, update reservoir <br />c storage. And compute the below Guernsey reservoir reach's <br />c flow and residual demand. <br />DO 210 K-l, 7 <br />RDEM(K) = KTDEM(K) <br />TDEM = TDEM + ROEM(K) <br />REMFL - REMFL + INF(K) <br />C *********TAKE AS MUCH OF THE REACH DEMAND FROM THE SPILL AND INF AS POSSIBLE <br />DEL = AMINl(RDEM(K),REMFL) <br />RDEM(K) - RDEM(K) - DEL <br />REMFL - REMFL - DEL <br />C *********INITIALIZE THE REACH FLOW <br />RFLO(K)=REMFL <br />C*********FILL THE RESERVOIR AND MAKE MINIMUM RELEASE <br />C EVAPORATION AT BEGINING OF MONTH <br />TEVAP = O. <br />STl - RST(K) <br />IF (STl .LE. 0.) GO TO 60 <br />AR-O <br />AR - COEFFS(1,K)*STI**COEFFS(2,K) <br />TEVAP- AR*AVEV(J,K) <br />60 CONTINUE <br />REVAP(K)-TEVAP <br />RST(K)-RST(K)-REVAP(K) <br />DEL - AMAXl( -RST(K),AMINl(REMFL-YRE(K,J),XST(K)-RST(K))) <br />RST(K)-RST(K)+DEL <br />REMFL - REMFL - DEL <br />C ********INITIALIZE RUNNING RELEASE <br />RRE(K) - REMFL <br />c Adding test code here by NRC <br />if(((xst(K)-rst(k)).lt.(remfl-yre(k,j))).and. <br />1 ((remfl.gt.xre(k)))) then <br />print *, 'There is a spill !' <br />endif <br />210 CONTINUE <br />C ********************************************************************* <br />C ADD <br />C <br />C <br />C <br />C <br />C <br /> <br />FLOW BELOW GUERNSEY, SATISFY REACH <br />IF((SUMDEM-SUMINF).LT.O.O) THEN <br />SUMINF - SUMINF - SUMDEM <br />ELSE <br />SUMINF = 0.0 <br />ENDIF <br />REMFL = REMFL + INF(B) <br /> <br />6 DEMANDS FROM SPILL <br /> <br />15 <br />