|
<br />~.,' ,; ,
<br />
<br />OOOG88
<br />
<br />+IFRST ,'IGROUP (40) , HOLD(224, 13) , ADATA (35,12) ,ISTART, lEND, J, I,
<br />+NT,NCL,IE,NG,NL,NC,NA,NH,NYI,NRES,TITLE(2,8),ATHEAD(l0,8),
<br />+ALHEAD(224,4) ,AGHEAD(40,4),CNAME(85,4),ISTUDY,ANAME(35,8),RDATE,
<br />+OWN(ll) ,OWNACC(7) ,GAIN(5),STORE(7) ,OWNMAX(6),OWNPRE(ll),PREACC(7)
<br />DIMENSION DDATA(28) ,PATHTR(12)
<br />DATA M11/31., 30. , 31 . , 3 1 . , 28 . , 31 . , 30. , 31. , 30 . , 31 . , 31. , 30 . /
<br />
<br />C
<br />
<br />LJ = 28
<br />JJ = 77
<br />DO 30 K = 1,28
<br />DDATA(K) = CDATA(K+53)*MM(J)*.001983471
<br />30 CONTINUE
<br />SO IF(J.LT.8) GO TO 1000
<br />PATHTR(J) = 0.0
<br />C NATURAL FLOW ABOVE ALCOVA
<br />FLOW(77+1,J) = HDATA(l,J)
<br />FLOW(77+2,J) = HDATA(4,J) + HDATA(6,J)
<br />FLOW(77+3,J) = FLOW(77+1,J)+FLOW(77+2,J)
<br />C STORAGE RELEASE FROM GRAY REEF
<br />FLOW(77+4,J) = HDATA(9,J)
<br />FLOW(77+S,J) = AMIN1(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+1,J)
<br />C GLENDO INFLOW
<br />FLOW(77+7,J) = FLOW(77+4,J) + HDATA(10,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) = N1AX1(FLOW(77+8,J),0.0)
<br />C STORAGE RELEASE FROM GLENDO
<br />FLOW(77+9,J) = HDATA(ll,J)
<br />FLOW(77+10,J) = AMIN1(FLOW(77+8,J),HDATA(ll,J))
<br />FLOW(77+10,J) = HDATA(ll,J) - FLOW(77+10,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 = AMIN1(PATHTR(J),FLOW(77+10,J))
<br />FLOW(77+10,J) = FLOW(77+10,J) - RESDL
<br />PATHTR(J) = PATHTR(J) - RESDL
<br />ENDIF
<br />C NATURAL FLOW AT GUERNSEY
<br />FLOW(77+11,J) = HDATA(ll,J) +,HDATA(12,J)
<br />FLOW(77+12,J) = FLOW(77+11,J) - FLOW(77+10,J)
<br />FLOW(77+12,J) = AMAX1(FLOW(77+12,J),0.0)
<br />C STORAGE RELEASE FROM GUERNSEY
<br />FLOW(77+13,J) = HDATA(13,J)
<br />FLOW(77+14,J) = AMIN1(FLOW(77+12,J),HDATA(13,J))
<br />FLOW(77+14,J) = HDATA(13,J) - FLOW(77+14,J)
<br />IF(PATHTR(J).GT.FLOW(77+14,J)) THEN
<br />RESDL = AMIN1(PATHTR(J),FLOW(77+14,J))
<br />FLOW(77+14,J) = FLOW(77+14,J) - RESDL
<br />PATHTR(J) = PATHTR(J) - RESDL
<br />ENDIF
<br />C PROPORTION GUERNSEY-WHALEN LOSS
<br />FLOW(77+15,J) = AMAX1(0.0,HDATA(13,J) + HDATA(14,J))
<br />IF (HDATA(14,J) .LT.O.O) THEN
<br />FLOW(77+17,J)=FLOW(77+14,J)+((FLOW(77+14,J)/HDATA(13,J)) "
<br />1*HDATA(14,J))
<br />ELSE
<br />FLOW(77+17,J) = FLOW(77+14,J)
<br />
|