|
<br />L
<br />
<br />C
<br />
<br />C
<br />
<br /> 53
<br /> 55
<br />C
<br /> 60
<br />C
<br /> 70
<br /> 90
<br />C
<br />C
<br />C 92
<br />C 93
<br /> 92
<br /> 95
<br /> 100
<br />C
<br />
<br />C
<br />
<br />.. .
<br />
<br />: OOlHi71
<br />
<br />oo20J=1,1::
<br />20 HOLD(K,J)=O.O
<br />I=ISTART-IFRST+1
<br />IE=IEND-IFRST+1
<br />IFLAG=O
<br />30 CONTINUE
<br />ZERO FLOW ARRAY
<br />DO 40 K=l,NL
<br />DO 40 J=1,13
<br />40 FLOW(K,J)=O.O
<br />IF(INPT.GT.O}CALL READONE
<br />CALL COMPUTE
<br />SUM ANNUAL DATA
<br />DO 50 K=l,NL
<br />DO 50 J=1,12
<br />50 FLOW(K,13)=FLOW(K,13)+FLOW(K,J}
<br />SET END OF YEAR RESERVOIR CONTENT TO ZERO HERE
<br />IF(NRES.EQ.O)GO TO 55
<br />DO 53 I1=1,NRES
<br />FLOW(IRES(I1),13)=O.0
<br />CONTINUE
<br />CONTINUE
<br />SUM DATA INTO SUMMARY ARRAY
<br />DO 60 K=l,NL
<br />DO 60 J=1,13
<br />HOLD(K,J)=HOLD(K,J)+FLOW(K,J)
<br />SUM ANNUAL DATA IN TABLES
<br />DO 70 K=l,NT
<br />DO 70 J=1,12
<br />TABLE(K,I ,13}=TABLE(K,I,13)+TABLE(K,I,J)
<br />CALL WRITE 1 (IFLAG,ANS,GALF)
<br />IYEAR = I + IFRST - 1
<br />DO 93 LQ = 1, NH
<br />WRITE(9,92)IYEAR, (HDATA(LQ,J),J=1,12)
<br />FORMAT(I4,12F7.2)
<br />CONTINUE
<br />WRITE(9,92) IYEAR, (FLOW(75,J) ,J=1,12)
<br />FORMAT(I4,12F7.2)
<br />IF(I.EQ.IE)GO TO 100
<br />1=1+1
<br />GO TO 30
<br />CONTINUE
<br />AVERAGE SUMMARY DATA
<br />NYEAR=IEND-ISTART+1
<br />AYEAR=NYEAR
<br />DO 110 K=l, NL
<br />DO 110 J=l, 13
<br />FLOW(K,J)=HOLD(K,J)jAYEAR
<br />IFLAG=l
<br />CALL WRITE1(IFLAG,ANS,GALF)
<br />CONTINUE
<br />CALL WRITE2
<br />IF(NOUT.GT.O)CALL SAVTAB(NOUT,IFILE,ISAVE)
<br />CLOSE(6,STATUS='KEEP'}
<br />CLOSE(9,STATUS='KEEP')
<br />IF(IPLT.NE.O)CALL PLOTTER(IPLT,IPLOT)
<br />STOP
<br />END
<br />*******************
<br />
<br />110
<br />
<br />105
<br />107
<br />125
<br />
<br />130
<br />
<br />SUBROUTINE READS(INPT,NOUT,IFILE,ISAVE,IPLT,IPLOT,IRES)
<br />
|