|
|
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: 3291 (0xcdb)
Types: TextFile
Names: »S00022.f«
└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S00022.f«
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S00022.f«
SUBROUTINE S00022
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/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/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)
COMMON/AAA13/TEM,P,DT,HNIV(22),VITESS
1 ,ZNIV,RI(2),ISS,NITERA,NIV,ISLNI,IBTR,IHTR,ITYP(22)
COMMON/AAA77/QINIT,TPUI(20),QPUI(20),RPUI,PLMOY,
1 XHTC(22,12),PL(22,12),QREPA(22,12),
2 XMCGA(22,12),XMCUO(22,12),QCED(22,12),QCON(22,12)
3 ,QUGA(22,12),QFOU(22,12),QCEI(2),QCONI(2),QUGII(2)
4 ,TMI(2),DTMI(2),QTOTAL,QTUG,QTGF,TERC(2)
4 ,TGAI(22,12),TGAII(22,12),T876(22,12),TUMAX(22,12)
5 ,TUMIN(22,12),DTGAI(22,12),DT876(22,12)
6 ,TERCO(22,12),TMU(2),DTMU(2),NPUI,NASL(22),NSL(22)
COMMON/AAA11/V0011N(20),V0012X(20),V0011,V0012,
1TLIM,ZLIMIN,ZLIMAX,DTVOI,TIMIMP(20),DIFREF(20),DIFMIC(20),
3V00001,V22202,TEMIMP,TEMMIC,DELM,DELV,DELVI,DELVS,
6V22201,V22203,DTPRE,DTPHY,DTTEMP,DTCHAN,DTDIS,DTBEF,DTORG
7,NIMIC,NMULTR,IZONE,N00011(20),ISORT,NECRIT,NSTAR,IPLOT,NZONE
+,IFRE(20),ICO,IIMP,NPAS1(20),NPAS2(20)
COMMON/AAA55/KPRESS,KDEB,KCON,KTGAI,KT876,ICAT(3),KIMP,KGLISS
COMMON/AAA66/TS,HFS,HGS,HFGS,DHFPS,DHGPS,VVFS,VVGS,DVVFPS,DVVGPS
1 ,DHFGS,VV1,VV2,DVV1P,DVV2P,T1,T2
COMMON/AAA99/HLCF,HVCFO,HVAP,HCON,HVIN,HLIN,R876,RGAI
1 ,XM876,XMGAI,XL876,XLGAI,CP876,CPGAI,XMC876,XMCGAI
2 ,DIAHY,HGAPP(22),HCG(22),HTC(22,12),HVCFF(22)
DOUBLEPRECISION M
IF(ICAT(2).LE.2)GOTO 900
Z1=XL(NIV)/2.
Z2=HNIV(NIV)/2.
IF(NIV.GT.2)Z1=XL(NIV-1)/2.
HJ=F00111(NIV-1)
X=(HJ-HFS)/(HGS-HFS)
AL3=X/(X+(1.-X)*VVFS/VVGS)
IF(AL3.LT..05)Z1=Z1*AL3/.05
PP1=(ALFA(NIV,1)-AL3)/(Z1+Z2)
IF(PP1.LT.0.)PP1=0.
PP2=PP1
900 CONTINUE
CALL S00012
DO 910 I=1,NC1
DEBAV(I)=0.
XA(I)=0.
910 XB(I)=0.
CALL S00013
DO 2000 I=1,NC1
IDODUC=ITYP(I)
GOTO(1010,1020,1030),IDODUC
1010 HA(I)=H(I,1)
XA(I)=(HA(I)-HF(I,1))/(HG(I,1)-HF(I,1))
ALFAD(I)=ALFA(I,1)
HE(I,1)=HF(I,1)
HE(I,2)=HG(I,1)
GOTO 1000
1020 HA(I)=F00111(I)
XA(I)=(HA(I)-H(I,1))/(H(I,2)-H(I,1))
HE(I,1)=H(I,1)
HE(I,2)=H(I,2)
ALFAD(I)=(ALFA(I,1)*V(I,1)+ALFA(I,2)*V(I,2))/VC(I)
GOTO 1000
1030 HA(I)=H(I,2)
XA(I)=(HA(I)-HF(I,2))/(HG(I,2)-HF(I,2))
ALFAD(I)=ALFA(I,2)
HE(I,1)=HF(I,2)
HE(I,2)=HG(I,2)
1000 HB(I)=0.
XB(I)=0.
2000 CONTINUE
CALL S00014
RETURN
END