|
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