FORTRAN IV G LEVEL 18 



MAIN 



DATE = 71201 



22/08/17 



0001 



0002 

 0003 

 0004 

 0005 



0006 

 0007 



0008 



0009 



0010 



0011 

 0012 

 0013 

 0014 

 0015 

 0016 

 0017 

 0018 

 0019 

 0C20 

 0021 

 0022 

 0023 

 002A 

 0C2 5 

 0026 



0C27 

 0028 

 0029 

 0030 

 0031 

 0C32 

 0033 

 0034 

 0035 

 0036 

 0037 

 0038 

 0039 

 0040 

 0041 

 0042 

 0043 

 0044 

 0045 

 0046 

 0047 

 0048 

 0049 

 0050 

 0051 

 0052 

 0C53 



TO 8001 

 GO TO 100 

 J = J^1 



I •/15X, 'STCPM NUfBEP' , 14) 



DIMENSION S( 9000) ,HCuR<9000) ,0FLL(300) ,TS ( 300 » ,T P ( 300 > , TE ( 300 ) , 

 2RS ( 300) ,RE( 300) ,0R1 (3 00I ,DR2 (300 ) , DT 1( 300) .0T2(300 ), RP (300 ), 

 3UP ( 300) ,0M( 300) , PERI M( 300) ,A( 300) ,PERFL (300 ) ,A2 (300) , ARISE (300 ) , 

 4AGLDN( 300) ,BASE ( 300) 



DIMENSION OFLPK( 3C0) ,QFPPK(3C0) 



DIMENSION TOTTM(300) ,0INCH(300) 



DIMENSION SUBFF ( 3 00) , RF L L ( 300 ) , S LB IN (300 ) ,RI NCH( 300) 

 D I MENS I ON PIFLO{ 300) , PRFL0(300) ,P0BPK(300 ) ,PCAPK (300 ) 

 C THF PROGRAM ACCOMODATES 300 STORMS. IF A LARGER NUMBER OF STORMS IS 

 C ANTICIPATED, THE DIMFNSION STATEMENT ABOVE SHOULD BE CHANGED, OR 

 C THE NUMBER OF DATA BEING ANALYZED AT ONE TIME SHOULD BE REDUCED. 

 DO 7 J=l,9000 



7 S(J)=-1.0 

 ENDS=- 1 . 



C IN THIS PROGRAM, S=STAGE. THE CONVERSION FROM STAGE TO CFS TAKES PLACE 

 C BELPW, UNDER '0= AND •RS(K)= '.CONVERSION CONSTANTS SHOULD BE CHANGED 

 C TO FIT THE WEIR BEING USFO. 

 READ (5,l,ENr=e)S 



C THE FOLLOWING FORMAT STATEMENT MLST FIT THE FORfAT IN WHICH THE INPUT 

 C DATA FOR STAGE ARE P'JNCHED. 



1 FORMAT ( lOX, 12F 4. 3 ,2?X) 



C IF HYDROGRAPH READINGS ARE TAKEN AT INTERVALS OTHER THAN ONE HOUR, 

 C THEN A SEPARATE 'READ' STATEMENT SHOULD BE MADE. 



8 DO 2 J=1,900C 



2 HOUR(J)=J 

 J= 1 



K= 1 



4 START=S(J) 

 IF ( S( J ) .LE .ENDS) GO 



3 IF( S(J*1).GT. START) 

 IF ( S( J-t-1 ) .LF .START) 

 GO TO 4 



100 WRITF(6,5) K 



5 FORMATI ' 1 ' ,50X, • TABLE 

 TS(K)=HOUR(J) 



0= 1. 522«( B. 0**0. 026) 



RS(K) = S( J)**0*'32.G 



WR I TE (6,6) TS( K ) ,RS ( K ) 



6 F0RMAT(2X, 'STOPM BEGINS AT • , F 1 . 2 , 3 X , • I N I T I A L RUN OF F RATE IS',F12.6 

 2/2X,'TIME BASFFLORATE RUNOFRATE BASEFLO TCTALFLO STORMFLO? 



3 STORMFLO INTERFLO RAPIDFIO I NTE RF LORAT E • ) 

 U=0 

 RC=0 



UCUM1=RS(K) 

 DFL AF=0 



QFLO=0 ' 

 RP ( K )=RS( K ) 

 TP ( K ) =HOUR( J ) 

 U2 = 

 SUBFL=0 

 RFLO=0 

 UC UM 3 = R S( K ) 

 R 1 = S( J )**0»3?.0 

 R2 = S( J+1 )**0*32.0 

 200 PIFP=ABS(P2-UCUM1) 



niFT = HOUR( J*n-HCUR( J) 

 T=HnuR ( J + 1 ) 



H=S0''T(nirR**2*DIFT**2) 

 IF ( H.LT.niFP )H-DIFP 

 IF ( H.LT.DIFT )H»OIF T 

 Rr-(«2*R 1)/2.0«P1FT + RC 

 AL''HA-ATAN(PIFR/DrFT) 

 RETA«(H«ALPHA-DIFR|/H 

 IF(BrTA.LF.O) BETA-C 

 C SA» l.O+SIN( ALPHA) 



|J»FXP(-C SA )*SOR T(DIFR*niFT*TAN ( BET A) ) 

 Ur UM2xU*-UC UMl 

 IF(R2.GT.PP( K) )TP( K) =T 



16 



