DataMuseum.dk

Presents historical artifacts from the history of:

DKUUG/EUUG Conference tapes

This is an automatic "excavation" of a thematic subset of
artifacts from Datamuseum.dk's BitArchive.

See our Wiki for more about DKUUG/EUUG Conference tapes

Excavated with: AutoArchaeologist - Free & Open Source Software. 

top - metrics - download
Index: S T

⟦46a74c21d⟧ TextFile

    Length: 3686 (0xe66)
    Types: TextFile
    Names: »S00009.f«

Derivation

└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
    └─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S00009.f« 
    └─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S00009.f«

TextFile

      SUBROUTINE S00009
      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/AAA22/DH(22,2),DM(22,2),DPDT,DUU(22,2),DV(22,2),DNU(22,2)
      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/AAA88/FLUP(22,2),FLUV(22,2),FLUL(22,2),FLUIV(22,2),
     2DEBIL(22,2),FLUIL(22,2),DEBI(22),DEBHI(22),DEBIV(22,2)
      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/AAA66/TS,HFS,HGS,HFGS,DHFPS,DHGPS,VVFS,VVGS,DVVFPS,DVVGPS
     1             ,DHFGS,VV1,VV2,DVV1P,DVV2P,T1,T2
      COMMON/AAA55/KPRESS,KDEB,KCON,KTGAI,KT876,ICAT(3),KIMP,KGLISS
      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/AAA16/DBENT(20),TENT(20),HENT(20),DEBE,ENTE,
     2        DBSOR(20),TSOR(20),DEBS,DPDTT(20),TDPDT(20)
     3  ,HSOR(20),THSOR(20),ENTS,NDBENT,NDPDT,NDBSOR,NHSOR
      DOUBLEPRECISION M,CC(22,2)
      COEF=144.D+00/778.D+00
      DO 100 I=2,NC1
      J=I-1
      ITY=ITYP(I)
      J1=2
      IF(ITY.EQ.1)J1=1
      GOTO(101,102,103),ITY
 101  CONTINUE
      DM(I,1)=DEB(J)-DEB(I)
      DUU(I,1)=DEBL(J)*HE(J,1)+DEBV(J)*HE(J,2)-DEBL(I)*HE(I,1)-DEBV(I)*
     2HE(I,2)+FLUP(I,1)
      GO TO 1000
 102  CONTINUE
       DM(I,1)=DEBL(J)-DEBL(I)-DEBI(I)
      DM(I,2)=DEBV(J)-DEBV(I)+DEBI(I)
      DUU(I,1)=DEBL(J)*HE(J,1)-DEBL(I)*HE(I,1)
     1    +FLUL(I,1)+FLUIL(I,1)
     1   -DEBIL(I,1)*HF(I,1)+DEBIV(I,2)*HFS
     3    -DEBIV(I,1)*HGS+DEBIL(I,2)*HFS
      DUU(I,2)=DEBV(J)*HE(J,2)-DEBV(I)*HE(I,2)
     1    +FLUV(I,2)+FLUIV(I,2)
     2    +DEBIL(I,1)*HGS-DEBIV(I,2)*HG(I,2)
     3    +DEBIV(I,1)*HGS-DEBIL(I,2)*HFS
      GO TO 1000
 103  CONTINUE
      DM(I,1)=DEB(J)-DEBI(I)
      DM(I,2)=-DEB(I)+DEBI(I)
      DUU(I,1)=DEBL(J)*HE(J,1)+DEBV(J)*HE(J,2)-DEBHI(I)+FLUP(I,1)
      DUU(I,2)=-DEBL(I)*HE(I,1)-DEBV(I)*HE(I,2)+DEBHI(I)+FLUP(I,2)
 1000 CONTINUE
      DO 500 K=1,J1
      CC(I,K)=DUU(I,K)-H(I,K)*DM(I,K)
      AA(I,K)=VV(I,K)*DM(I,K)+DVVH(I,K)*CC(I,K)
      DV(I,K)=AA(I,K)+BB(I,K)*DPDT
      IF(DABS(DV(I,K)).LT.1.D-10)DV(I,K)=0.D+00
      EPSM=DELM*VC(I )/VVFS
      YM=M(I,K)
      IF(YM.LT.EPSM)YM=EPSM
      DH(I,K)=(CC(I,K)+COEF*V(I,K)*DPDT)/YM
      DNU(I,K)=DVVH(I,K)*DH(I,K)+DVVP(I,K)*DPDT
  500 CONTINUE
 100  CONTINUE
      TAFT=TEM+1.D-4
      CALL S00934(THSOR,HSOR,HSOR,TAFT,HAFT,X,2)
      DH(NC2,1)=(HAFT-H(NC2,1))/1.D-4
      DNU(NC2,1)=DVVH(NC2,1)*DH(NC2,1)+DVVP(NC2,1)*DPDT
      CALL S00934(TENT,HENT,HENT,TAFT,HAFT,X,20)
      DH(1,1)=(HAFT-H(1,1))/1.D-4
      DNU(1,1)=DVVH(1,1)*DH(1,1)+DVVP(1,1)*DPDT
      IF(ICAT(2).LE.2)RETURN
      ISS=1
      VITESS=DV(NIV,1)
      IF(VITESS.LT.0.)ISS=2
      RETURN
      END