|
|
DataMuseum.dkPresents historical artifacts from the history of: DKUUG/EUUG Conference tapes |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about DKUUG/EUUG Conference tapes Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - downloadIndex: S T
Length: 1757 (0x6dd)
Types: TextFile
Names: »S22202.f«
└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S22202.f«
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S22202.f«
DOUBLEPRECISIONFUNCTION S22202(AL,DAL,TW,DTW,TF,DTF,P,DPDT,
* DIAHY,IOPT)
IMPLICITDOUBLEPRECISION(A-H,O-Z)
DOUBLEPRECISION R(15),F(15)
DATA R /.15,.25,.35,.45,.65,.85,1.05,1.25,1.45,1.95,2.45,
1 2.95,3.45,3.95,4.45/
DATA F /.4258,.5359,.6235,.6981,.8238,.9295,1.0222,
1 1.1056,1.1820,1.3506,1.4967,1.6271,1.7459,1.8555,1.9578/
TWRA=TW/1000.+.45967
TFRA=TF/1000.+.45967
TFC=TFRA*TFRA
TWC=TWRA*TWRA
TS=TWRA+TFRA
TSC=TS*TS
HRT=TS*(TWC+TFC)
HRFT=TSC+2.*TFC
HRWT=TSC+2.*TWC
IF(IOPT.EQ.3)GOTO 2000
RR=TFRA/TWRA
IF(RR.GE..45)GOTO 200
I1= IDINT(RR/.1-1.5)+1
GOTO 500
200 IF(RR.GE.1.45)GOTO 300
I1=IDINT(RR/.2-2.25)+4
GOTO 500
300 I1=IDINT(RR/.5-2.9)+9
500 I2=I1+1
R1=R(I1)
R2=R(I2)
F1=F(I1)
F2=F(I2)
PENT=(F2-F1)/(R2-R1)
FF=F1+(RR-R1)*PENT
DTWA=DTW/1000.D+00
DTFA=DTF/1000.D+00
XT=1.D-4
PX=P+XT*DPDT
TWX=TWRA+DTWA*XT
TFX=TFRA+DTFA*XT
PWLX=(PX/14.696)*DIAHY*TWX/TFX
PWL=(P/14.696)*DIAHY*TWRA/TFRA
EWALL=.6
TT=TWRA*1000.D+00
TTX=TWX*1000.D+00
EH2O=S00093(PWL,TT)
EH2OX=S00093(PWLX,TTX)
DEH2O=(EH2OX-EH2O)/XT
AH2O=2.*EH2O*FF
DRR=(DTFA-RR*DTWA)/TWRA
DF=PENT*DRR
DAH2O=2.*(EH2O*DF+DEH2O*FF)
DEN=1./EWALL+1./AH2O-1.D+00
DCC=1.73*DAH2O/(DEN*DEN*AH2O*AH2O)
CC=1.73/DEN
S22202=(DCC*HRT+CC*(DTWA*HRWT+DTFA*HRFT))/3600.D+00
RETURN
2000 CC=1.73/3600.D+00
DEN=1.25-.15*AL
DHDTF=CC*(1.-.6*AL)*HRFT/DEN
DHDTW=CC*(1.-.6*AL)*HRWT/DEN
DHDAL=CC*HRT*(-.6)/(DEN*DEN)
S22202=DHDTF*DTF+DHDTW*DTW+DHDAL*DAL
RETURN
END