|
|
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: 2148 (0x864)
Types: TextFile
Names: »S33019.f«
└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S33019.f«
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S33019.f«
SUBROUTINE S33019(I,SEC)
IMPLICITDOUBLEPRECISION(A-H,O-Z)
COMMON/AAA26/DEB(21),DEBV(21),DEBL(21),ALFAD(21),ALFB(21),
2ALFT(21),HA(21),HB(21),HE(21,2),HEM(21),VE(21,2),
3 XA(21),XB(21),XD(21),XEM(21),SDSC
4 ,TJONC(21),HJONB(21,2),HJONT(21,2),DEBAV(21)
5 ,GLINF(21),GLSUP(21),PP1,PP2,ALNIV
6 ,VGJO(21),GLINFP(21),V55198(21),ICAS(21),ICZW(21)
COMMON/AAA55/KPRESS,KDEB,KCON,KTGAI,KT876,ICAT(3),KIMP,KGLISS
COMMON/AAA44/SC(21),SCAPA(22),SHTC(22),VC(22),XL(22)
1,VFIX,ZCOT(21),I6FT,J6FT,I8FT,J8FT,NC,NC1,NC2,IBCH,IHCH
COMMON/AAA13/TEM,P,DT,HNIV(22),VITESS
1 ,ZNIV,RI(2),ISS,NITERA,NIV,ISLNI,IBTR,IHTR,ITYP(22)
COMMON/AAA66/TS,HFS,HGS,HFGS,DHFPS,DHGPS,VVFS,VVGS,DVVFPS,DVVGPS
1 ,DHFGS,VV1,VV2,DVV1P,DVV2P,T1,T2
COMMON/AAA10/V(22,2),VV(22,2),H(22,2),HF(22,2),HG(22,2),
1 DNUHP(22,2),U(22,2),XM(22,2),XT(22,2),ALFA(22,2),
2 T(22,2),DVVH(22,2),DELH(22,2),DVVP(22,2),DELP(22,2)
3 ,AA(22,2),BB(22,2),DTDH(22,2),DTDP(22,2),DTSAT,
4 DNUPP(22,2),A(22),B(22),C(22),TVID(22,12),M(22,2),IST(22,2)
DOUBLEPRECISION M
VF=VE(I,1)
VG=VE(I,2)
G=DEB(I)/SEC
TJON=TJONC(I)
ICHOIX=KGLISS
GOTO (10,20,30,40),ICHOIX
10 ALF=ALFB(I)
IC=1
IF(G.GT.0.)GOTO 15
ALF=ALFT(I)
IC=3
15 XAA=ALF*VF/(ALF*VF+(1.-ALF)*VG)
XBB=0.
ALFAD(I)=ALF
GLINF(I)=0.
GLSUP(I)=0.
GOTO 100
20 CONTINUE
RETURN
30 GOTO 90
40 CONTINUE
ALB=ALFB(I)
ALT=ALFT(I)
42 CALL F88345(G,ALB,ALT,VG,VF,TJON,X1,X2,X3,ALFAD(I),X4,PE,
1 CST,IC,GL1,GL2,ICZW(I))
90 GLINF(I)=GL2
GLSUP(I)=GL1
VGJO(I)=X4
XAA=PE
XBB=CST*SEC
100 ICAS(I)=IC
XA(I)=XAA
XB(I)=XBB
DELL=(1.-XAA)*DEB(I)-XBB
DEVV=XAA*DEB(I)+XBB
HEL=HJONT(I,1)
HEV=HJONB(I,2)
IF(DELL.GT.0.)HEL=HJONB(I,1)
IF(DEVV.LT.0.)HEV=HJONT(I,2)
HE(I,1)=HEL
HE(I,2)=HEV
HA(I)=XAA*HEV+(1.-XAA)*HEL
HB(I)=XBB*(HEV-HEL)
RETURN
END