DataMuseum.dk

top - metrics - download

    Length: 3104 (0xc20)
    Types: TextFile
    Names: »S00061.f«

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

TextFile

      SUBROUTINE S00061(I,HV,HVT,HVME,HVMS,YNU,RE,TFL,DTFL,IOPT)
      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/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)
      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/AAA22/DH(22,2),DM(22,2),DPDT,DUU(22,2),DV(22,2),DNU(22,2)
      COMMON/AAA18/DDEB(21),DDDEB(21),DXAT(21),DXAG(21),DXBT(21),
     1             DXBG(21),DNUG(21),DNUF(21),DHE(22,2),DDEBO(21)
      COMMON/AAA66/TS,HFS,HGS,HFGS,DHFPS,DHGPS,VVFS,VVGS,DVVFPS,DVVGPS
     1             ,DHFGS,VV1,VV2,DVV1P,DVV2P,T1,T2
      DOUBLEPRECISION X(21),Y(21)
      DATA A,B,C,D,E/.76082D-15,-.74064D-11,.23241D-7,.14794D-4,
     1               .10467D-1/
      DATA R,S/1.7399D-2,5.4997D-5/
      CALL X21Y21 (X,Y)
      HVT=0.D+00
      HVME=0.D+00
      HVMS=0.D+00
      TFLA=TFL
      IF(TFL.GE.1600.)TFLA=1600.D+00
      TF=TFLA
      IE=I-1
      DEBE=DEBAV(IE)
      DEBS=DEBAV(I)
      DEVE=XA(IE)*DEBE+XB(IE)
      DEVS=XA(I)*DEBS+XB(I)
      SEC=SCAPA(I)
      DEBM=(DEVE+DEVS)/2.D+00
      EPS=1.D+00
      IF(DEBM.LT.0.)EPS=-1.D+00
      VIS=(R+S*TF)/3600.D+00
      RE=EPS*DEBM*DIAHY/(VIS*SEC)
      REO=RE
      IF(RE.GE.1000.)GOTO 100
      I1=IDINT(RE/200.)+1
      GOTO 200
  100 IF(RE.GE.4000.)GOTO 110
      I1=IDINT(RE/500.)+4
      GOTO 200
  110 IF(RE.GE.10000.)GOTO 120
      I1=IDINT(RE/2000.)+10
      GOTO 200
  120 IF(RE.GE.30000.)GOTO 130
      I1=IDINT(RE/5000.)+13
      GOTO 200
  130 IF(RE.GE.49999.)RE=49999.D+00
      I1=IDINT(RE/10000.)+16
  200 I2=I1+1
      Y1=Y(I1)
      Y2=Y(I2)
      X1=X(I1)
      X2=X(I2)
      RAP=(Y2-Y1)/(X2-X1)
      TMRE=Y1+RAP*(RE-X1)
      YNU=.023*TMRE
      YNUO=YNU
      IF(YNU.LT.3.5)YNU=3.5
      TF=TFLA
      CON=E+TF*(D+TF*(C+TF*(B+A*TF)))
      HV=YNU*CON/(3600.*DIAHY)
      IF(IOPT.LE.0 )RETURN
      DTFLA=DTFL
      IF(TFL.GE.1600.)DTFLA=0.D+00
      DNUT=0.D+00
      IF(YNUO.LT.3.5)GOTO 1000
      IF(REO.GT.49999.)GOTO 1000
      DQE=(DXAG(IE)*DEBE+XA(IE)+DXBG(IE))*DDEB(IE)+DXAT(IE)*DEBE+
     1 DXBT(IE)
      DQS=(DXAG(I)*DEBS+XA(I)+DXBG(I))*DDEB(I)+DXAT(I)*DEBS+DXBT(I)
      DQT=.5*(DQE+DQS)
      DMUT=S*DTFLA/3600.D+00
      DRET=EPS*DIAHY*(DQT-DEBM*DMUT/VIS)/(VIS*SEC)
      DNUT=.023*RAP*DRET
      YNUMAX=5.*YNUO
      ABNUT=DABS(DNUT)
      IF(ABNUT.GE.YNUMAX)DNUT=DNUT*YNUMAX/ABNUT
 1000 CONTINUE
      DLAT=(D+TF*(2.*C+TF*(3.*B+4.*A*TF)))*DTFLA
      HVT=(YNU*DLAT+DNUT*CON)/(3600.*DIAHY)
      RETURN
      END