CE=7.E-5
CH=3.4E-5
AEDA=0.9
TSD=120.
TRD=60.
CALL WTRTBL(TRD,PRD,RHORD,RHD,RMUD)
CALL WTRTBL(TSD,PSD,RHOSD,SHD,SMUD)
RHOD=(RHOSD+RHORD)/2.
TS=120.
TA=20.
A12=A*PVFE*PL*(EPS**B)*(1.273240**(2+C))
QF=(0.425*COS(6.283*T/8760.))+0.575
CALL CLMTD(TS,TA,QF,FRF,TR,TRG,FRFG)
FT=FRF
C CALCULATE QUANITIES WHICH MAY BE A FUNCTIONS OF TIME.
CALL WTRTBL(TS,PS,RHOS,SH,SMU)
CALL WTRTBL(TR,PR,RHOR,RH,RMU)
A7=(((1./(SMU**C))/(RHOS*RHOS))+((1./(RMU**C))/(RHOR*RHOR)))/2.
RHOA=(RHOS+RHOR)/2.
FUN=A7*((((CE*RHOA)/(AEDA*RHOD))*(FT* *(2+C)))-((PVFH/PVFE)
&*CH*(FT**(3+C))))
FUN=FUN*A12*(FRD**(3+C))
PRINT 102, T,FT,FUN
102 FORMAT(5X,3E 15.4)
RETURN
END
Program I2
PROGRAM I2
EXTERNAL FUN
CALL ROMBRG(FUN,0.,4380.,0.,1.E-7,PI2)
PRINT 101,PI2+2.
101 FORMAT(l0X,F15.5)
STOP
END
Program I2-C-GMT
PROGRAM I2-C-GMT
SUBROUTINE ROMBRG(FUN,A,B,C,ERR,RES)
C THIS SUBROUTINE COMPUTES INTEGRALS OF A USER SUPPLIED
FUNCTION
C USING ROMBERG'S METHOD. THIS SUBROUTINE IS FROM
"NUMERICAL
C METHODS FOR ENGINEERING APPLICATION",J.H. FERZIGER, JOHN
WILEY
C AND SONS, 1981. THE ARGUMENTS ARE:
C FUN = THE FUNCTION TO BE INTEGRATED
C A = LOWER LIMIT
88
Previous Page
