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

⟦81faaddff⟧ TextFile

    Length: 6659 (0x1a03)
    Types: TextFile
    Names: »S55198.f«

Derivation

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

TextFile

      SUBROUTINE S55198(H,PP,T,RHO,DVDHP,DVDPH,DTDH,DTDP)
      IMPLICITDOUBLEPRECISION(A-H,O-Z)
      DATA A00,A01,A02,A03,A04,A05,A06,A07,A08,A09/
     1     6170.78663883D+00, 489.886451913D+00,-1.71011901617D+00,
     2     .00310665805845D+00,-34.6610591433D-07 , 24.1062845537D-10 ,
     3    -10.411445494D-13  ,  27.0989223969D-17,-38.9010136647D-21 ,
     4     23.654728597D-25/
      DATA A10,A11,A12,A13,A14,A15,A16,A17,A18,A19/
     1    -27.2234538756D+00,-1.55945023546D+00, .0055688166592D+00,
     2    -10.1907487805D-06 , 11.3988461203D-09 ,-79.2905387618D-13 ,
     3     34.2019458067D-16 ,-88.8161818151D-20 , 12.7093552905D-23 ,
     4    -76.9751114847D-28/
      DATA A20,A21,A22,A23,A24,A25,A26,A27,A28,A29/
     1     .0400334696045D+00, .0019832294859D+00,-72.2377802606D-07 ,
     2     13.3067830456D-09 ,-14.917384566D-12  , 10.3760380082D-15 ,
     3    -44.6915456261D-19 , 11.5763883026D-22 ,-16.5086089907D-26 ,
     4     99.5533322201D-31/
      DATA A30,A31,A32,A33,A34,A35,A36,A37,A38,A39/
     1    -26.7334101475D-06 ,-12.5868723956D-07 , 46.6552793801D-10 ,
     2    -86.4549289781D-13 , 97.1047694917D-16 ,-67.5258399441D-19 ,
     3     29.036844039D-22  ,-75.0092619426D-26 , 10.6572651116D-29 ,
     4    -63.9675459557D-34/
      DATA A40,A41,A42,A43,A44,A45,A46,A47,A48,A49/
     1     87.7523620691D-10 , 39.8481772183D-11 ,-15.0016095749D-13 ,
     2     27.947539761D-16  ,-31.4412831D-19    , 21.8544044613D-22 ,
     3    -93.8062091533D-26 , 24.1619703535D-29 ,-34.1939584316D-33 ,
     4     20.421351707D-37/
      DATA A50,A51,A52,A53,A54,A55,A56,A57,A58,A59/
     1    -11.4250142505D-13 ,-50.3259455874D-15 , 19.2105899709D-17 ,
     2    -35.9600016574D-20 , 40.5103597804D-23 ,-28.1410621514D-26 ,
     3     12.0553990633D-29 ,-30.9558667049D-33 , 43.6260674D-37    ,
     4    -25.9155193335D-38/
      DATA B00,B01,B02,B03,B04,B05,B06,B07/
     1     445.665390762D+00,-2.22263915916D+00, .00472884372253D+00,
     2    -55.6382273478D-07 , 39.0978422275D-10 ,-16.4098861643D-13 ,
     3     38.0911781785D-17 ,-37.7244080484D-21/
      DATA B10,B11,B12,B13,B14,B15,B16,B17/
     1    -1093412.10859D+00, 5392.48720906D+00,-11.3504501908D+00,
     2     .0132165571029D+00,-91.9374308338D-07 , 38.2114755965D-10 ,
     3    -87.8718361625D-14 , 86.2568295293D-18/
      DATA B20,B21,B22,B23,B24,B25,B26,B27/
     1     473986630.09D+00 ,-2303848.6116D+00 , 4783.56767793D+00,
     2    -5.50018860399D+00, .00378235279254D+00,-15.556306052D-07  ,
     3     35.4309477645D-11 ,-34.4728839248D-15/
      DATA B30,B31,B32,B33,B34,B35,B36,B37/
     1    -59101208761.4D+00, 285687397.896D+00,-590083.283098D+00,
     2     675.104420712D+00,-.46204835791D+00 , .000189172522631D+00,
     3    -42.8993080892D-09 , 41.5668966327D-13/
      DATA B40,B41,B42,B43,B44,B45,B46,B47/
     1     3147752810980.D+00,-15146899832.4D+00, 31147138.3352D+00,
     2    -35480.7250254D+00, 24.1808570239D+00,-.00985946157921D+00,
     3     22.2693241999D-07 ,-21.4941089448D-11/
      DATA B50,B51,B52,B53,B54,B55,B56,B57/
     1    -75347916606600.D+00, 361189325093.D+00,-739923492.971D+00,
     2     839728.397005D+00,-570.189938069D+00, .231648962465D+00,
     3    -52.1369298193D-06 , 50.1481876596D-10/
      DATA B60,B61,B62,B63,B64,B65,B66,B67/
     1     663617137215000.D+00,-3171123872380.D+00, 6475932094.22D+00,
     2    -7326561.7806D+00 , 4959.5175126D+00 ,-2.00874787721D+00,
     3    .000450750907579D+00,-43.2282398141D-09/
      P      = PP
      IF( P.LT.25. ) P= 25.D+00
      A0   = A00 +P*(A01+P*(A02+P*(A03+P*(A04+P*(A05+P*(A06+P*(A07+P*
     1       (A08+P*(A09)))))))))
      A1   = A10+P*(A11+P*(A12+P*(A13+P*(A14+P*(A15+P*(A16+P*(A17+P*
     1       (A18+P*(A19)))))))))
      A2   = A20+P*(A21+P*(A22+P*(A23+P*(A24+P*(A25+P*(A26+P*(A27+P*
     1       (A28+P*(A29)))))))))
      A3   = A30+P*(A31+P*(A32+P*(A33+P*(A34+P*(A35+P*(A36+P*(A37+P*
     1       (A38+P*(A39)))))))))
      A4   = A40+P*(A41+P*(A42+P*(A43+P*(A44+P*(A45+P*(A46+P*(A47+P*
     1       (A48+P*(A49)))))))))
      A5   = A50+P*(A51+P*(A52+P*(A53+P*(A54+P*(A55+P*(A56+P*(A57+P*
     1       (A58+P*(A59)))))))))
      DA0=A01+P*(2.*A02+P*(3.*A03+P*(4.*A04+P*(5.*A05+P*(6.*A06+
     1 P*(7.*A07+P*(8.*A08+9.*A09*P)))))))
      DA1=A11+P*(2.*A12+P*(3.*A13+P*(4.*A14+P*(5.*A15+P*(6.*A16+
     1 P*(7.*A17+P*(8.*A18+9.*A19*P)))))))
      DA2=A21+P*(2.*A22+P*(3.*A23+P*(4.*A24+P*(5.*A25+P*(6.*A26+
     1 P*(7.*A27+P*(8.*A28+9.*A29*P)))))))
      DA3=A31+P*(2.*A32+P*(3.*A33+P*(4.*A34+P*(5.*A35+P*(6.*A36+
     1 P*(7.*A37+P*(8.*A38+9.*A39*P)))))))
      DA4=A41+P*(2.*A42+P*(3.*A43+P*(4.*A44+P*(5.*A45+P*(6.*A46+
     1 P*(7.*A47+P*(8.*A48+9.*A49*P)))))))
      DA5=A51+P*(2.*A52+P*(3.*A53+P*(4.*A54+P*(5.*A55+P*(6.*A56+
     1 P*(7.*A57+P*(8.*A58+9.*A59*P)))))))
      T    = A0+H*(A1+H*(A2+H*(A3+H*(A4+H*(A5)))))
      DTDH=A1+H*(2.*A2+H*(3.*A3+H*(4.*A4+5.*A5*H)))
      DTDP=DA0+H*(DA1+H*(DA2+H*(DA3+H*(DA4+H*DA5))))
      B0 = B00+H*(B01+H*(B02+H*(B03+H*(B04+H*(B05+H*(B06+H*(B07)))))))
      B1 = B10+H*(B11+H*(B12+H*(B13+H*(B14+H*(B15+H*(B16+H*(B17)))))))
      B2 = B20+H*(B21+H*(B22+H*(B23+H*(B24+H*(B25+H*(B26+H*(B27)))))))
      B3 = B30+H*(B31+H*(B32+H*(B33+H*(B34+H*(B35+H*(B36+H*(B37)))))))
      B4 = B40+H*(B41+H*(B42+H*(B43+H*(B44+H*(B45+H*(B46+H*(B47)))))))
      B5 = B50+H*(B51+H*(B52+H*(B53+H*(B54+H*(B55+H*(B56+H*(B57)))))))
      B6 = B60+H*(B61+H*(B62+H*(B63+H*(B64+H*(B65+H*(B66+H*(B67)))))))
      DB0   = B01+H*(2.*B02+H*(3.*B03+H*(4.*B04+H*(5.*B05+H*(6.*B06
     1           +H*(7.*B07))))))
      DB1   = B11+H*(2.*B12+H*(3.*B13+H*(4.*B14+H*(5.*B15+H*(6.*B16
     1           +H*(7.*B17))))))
      DB2   = B21+H*(2.*B22+H*(3.*B23+H*(4.*B24+H*(5.*B25+H*(6.*B26
     1           +H*(7.*B27))))))
      DB3   = B31+H*(2.*B32+H*(3.*B33+H*(4.*B34+H*(5.*B35+H*(6.*B36
     1           +H*(7.*B37))))))
      DB4   = B41+H*(2.*B42+H*(3.*B43+H*(4.*B44+H*(5.*B45+H*(6.*B46
     1           +H*(7.*B47))))))
      DB5   = B51+H*(2.*B52+H*(3.*B53+H*(4.*B54+H*(5.*B55+H*(6.*B56
     1           +H*(7.*B57))))))
      DB6   = B61+H*(2.*B62+H*(3.*B63+H*(4.*B64+H*(5.*B65+H*(6.*B66
     1           +H*(7.*B67))))))
      R   = 1./P
      V  = B0+R*(B1+R*(B2+R*(B3+R*(B4+R*(B5+R*(B6))))))
      DVDPH = -R**2*(B1+R*(2.*B2+R*(3.*B3+R*(4.*B4+R*(5.*B5+
     1                  R*(6.*B6))))))
      DVDHP = DB0+R*(DB1+R*(DB2+R*(DB3+R*(DB4+R*(DB5+R*(DB6))))))
      IF( PP.LT.25.) GO TO 200
      RHO  = 1./V
      RETURN
200   CONTINUE
      P      = PP
      IF( P.LT.0.1 ) P= 0.1
      V      = V*25./P
      DVDHP = DVDHP*25./P
      DVDPH = DVDPH*(25./P)**2
      RHO=1./V
      RETURN
      END