TSD=120.
TRD=55.
TA=20.
TS=120.
TGMTD=59.1608
RN1=1.3
A13=1.0
PL=1000.
DMF=4.
TF=(0.425*COS(6.283 *T/8760.))+0.575
TOP=(TSD-TRD)*TF
BOTTOM=TS-TA-((1/(TS-TA))*(TGMTD*TGMTD)*((TF/A13)**(2/RN1)))
A=TOP/BOTTOM
C "A" is the normalized mass flow rate, m/md.
FUN=(2.043E-8*DMF*((DMF*A)**1.9432))-(8.924E-9*((A*DMF)**2.9432))
FUN=FUN*PL/1000
PRINT 102, T,A,FUN
102 FORMAT(5X,3E15.4)
RETURN
END
Program I2
PROGRAM I2
EXTERNAL FUN
CALL ROMBRG(FUN,0.,4380.,0.,1.E-9,PI2)
PRINT 101,PI2*2.
101 FORMAT(l0X,E15.5)
STOP
END
Program I1EQ3-26
PROGRAM I1EQ3-26
FUNCTION FUN1(T,Z)
C THIS FUNCTION CONTAINS THE INTEGRAND OF THE I1 PARAMETER.
C It uses equation 3.26 which was derived using the GMTD appx.
C Modified from "IlFT.for" on 1/12/94.
RN1=1.3
A14=0.575
A15=0.425
A13=1.0
FUN1=(((A15*COS(6.283*T/8760.))+A14)/A13)**(2/RN1)
PRINT 102, T,FUN1
102 FORMAT(5X,2E15.4)
RETURN
END
90
Previous Page
