Ask Question

Name:
Title:
Your Question:

Answer Question

Name:
Your Answer:
User Submitted Source Code!


Description:
  WATERFLOW
Language: FORTRAN
Code:
C
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C  WATERFLOW.FOR
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C
C  WATER FLOW INTO OV FOR GIVEN LEAK DIAMETERS
C
C  BEGUN ON FEB 5, 2013
C  LAST EDITED ON 2/5/13
C
C  R.F. BOURQUE
C
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
C
     OPEN(UNIT=21, FILE='DATAOUT.TXT', STATUS='UNKNOWN')
C     OPEN(UNIT=24, FILE='DATAIN.TXT', STATUS='OLD')
C
     PI = 355./113.
     C = 39.37
C
C  INPUT:
C
     DOVIN = 110.
     DIVIN = 60.
     AREAIV = 28.
     QMAX = 25000.
C
C  OV BOTTOM SIMPLIFIED TO FLAT
C  IV BOTTOM HEIGHT ABOVE OF BOTTOM:
C
     AIN = 22.
C
C  LEAK DIAMETER:
C
     DLEAKIN = 3.
     PAIRPSIA = 20.
     HEADFT = 5.
C
C  CONVERT TO METRIC:
C
     ROV = DOVIN/(2.*C)
     AOV = PI*ROV**2
     RIV =DIVIN/(2.*C)
     ASMOOTH = 4*PI*RIV**2
     FMULT = AREAIV/ASMOOTH
     A = AIN/C
     DLEAK = DLEAKIN/C
     ALEAK = (PI/4.)*DLEAK**2
     PAIR = PAIRPSIA*6895.
     PHEAD = HEADFT*(62.4/144.)*6895.
     PTOT = PAIR+PHEAD
     VELH2O = SQRT(2.*PTOT/1000.)
     VOLRATE = ALEAK*VELH2O
C
C  Y IS VERTICAL HEIGHT FROM BOTTOM OF OV
C
     Y = 0.
     TIME = 0.
     VH201 = 0.
     VH202 = 0.
     VH20TOT = 0.
     ASURF = 0.
     QTOT = 0.
     YMAX = A+2.*RIV
     DY = .01
C
   1  CONTINUE
C
     WRITE(21,*) TIME,Y,QKW
C
     Y = Y+DY
C
     IF(Y .LE. A) VH2O1 = AOV*Y
     IF(Y .LE. A) TIME = VH2O1/VOLRATE
C  TIME IT TAKES FOR WATER TO CONTACT IV:
     IF(Y .LE.A) TTOUCH = TIME
     IF(Y .LE. A) GO TO 1
C
     H = Y-A
     X = RIV-H
C  VOLUME OF IV IN WATER:
     VIV = (2.*PI/3.)*RIV**3-(PI/3.)*X**2*(3.*RIV-X)
     VH2O2 = AOV*H-VIV
     VH2OTOT = VH2O1+VH2O2
     TIME = VH2OTOT/VOLRATE
     AIVSMUTH = 2.*PI*RIV*H
     ASURF = FMULT*AIVSMUTH
     QTOT = QMAX*ASURF
     QKW = QTOT/1000.
C
     IF(Y .GT. YMAX) GO TO 2
     GO TO 1
C
   2  CONTINUE
C
     END
C

          
Comments: