|
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: 6522 (0x197a)
Types: TextFile
Names: »S55199.f«
└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S55199.f«
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S55199.f«
SUBROUTINE S55199(H,P,T,RHO,DVDHP,DVDPH,DTDH,DTDP)
IMPLICITDOUBLEPRECISION(A-H,O-Z)
DATA A00,A01,A02,A03,A04,A05,A06,A07,A08,A09,A010,A011,A012 /
* 1.0421485616135D-02 ,
* 3.7989819258871D-04 ,
* -1.0809520514321D-05 ,
* 1.7199071099232D-07 ,
* -1.7003453661702D-09 ,
* 1.1093079504746D-11 ,
* -4.9305639341012D-14 ,
* 1.5141687819488D-16 ,
* -3.2078176823915D-19 ,
* 4.5968793700556D-22 ,
* -4.2500133969565D-25 ,
* 2.2858451991157D-28 ,
* -5.4289561756668D-32
1 /
DATA A10,A11,A12,A13,A14,A15,A16,A17,A18,A19,A110,A111,A112 /
* 7.1558907160166D-06 ,
* -4.8989819387822D-07 ,
* 1.3966895742219D-08 ,
* -2.2153628233107D-10 ,
* 2.1877549314087D-12 ,
* -1.4259859281335D-14 ,
* 6.3323408405408D-17 ,
* -1.9428241441369D-19 ,
* 4.1119914445076D-22 ,
* -5.8869455392068D-25 ,
* 5.4376738094709D-28 ,
* -2.9220477630923D-31 ,
* 6.9343597264571D-35
1 /
DATA A20,A21,A22,A23,A24,A25,A26,A27,A28,A29,A210,A211,A212 /
* -3.0284448323971D-09 ,
* 2.0572908444785D-10 ,
* -5.8631773583883D-12 ,
* 9.2947600855704D-14 ,
* -9.1733812676273D-16 ,
* 5.9751462062171D-18 ,
* -2.6513621151927D-20 ,
* 8.1280449373158D-23 ,
* -1.7188520355861D-25 ,
* 2.4587005585713D-28 ,
* -2.2691496114378D-31 ,
* 1.2183868386189D-34 ,
* -2.8891726411793D-38
1 /
DATA A30,A31,A32,A33,A34,A35,A36,A37,A38,A39,A310,A311,A312 /
* 4.1603162475165D-13 ,
* -2.8253313464467D-14 ,
* 8.0498643149541D-16 ,
* -1.2756275155692D-17 ,
* 1.2583652213261D-19 ,
* -8.1918809930522D-22 ,
* 3.6327419331754D-24 ,
* -1.1129150871768D-26 ,
* 2.3518675463763D-29 ,
* -3.3618179351310D-32 ,
* 3.1004787014482D-35 ,
* -1.6636314001942D-38 ,
* 3.9424525577961D-38
1 /
DATA B00,B01,B02,B03,B04,B05,B06,B07,B08,B09,B010,B011,B012 /
* 3.4799516636947D+01 ,
* 8.0207756895533D-01 ,
* 5.7134938990324D-03 ,
* -9.0888455214870D-05 ,
* 9.0343489356475D-07 ,
* -5.9586500435765D-09 ,
* 2.6874936213899D-11 ,
* -8.4045091912731D-14 ,
* 1.8181932356958D-16 ,
* -2.6657692257362D-19 ,
* 2.5243994056927D-22 ,
* -1.3912049726064D-25 ,
* 3.3848430431244D-29
1 /
DATA B10,B11,B12,B13,B14,B15,B16,B17,B18,B19,B110,B111,B112 /
* -3.3870909274776D-03 ,
* 3.3749976806899D-05 ,
* -8.3841532461841D-07 ,
* 1.3908406425265D-08 ,
* -1.5041874666294D-10 ,
* 1.1068510439562D-12 ,
* -5.6404562404400D-15 ,
* 1.9963935983267D-17 ,
* -4.8613722666270D-20 ,
* 7.9478277453504D-23 ,
* -8.3010018206858D-26 ,
* 4.9898102992625D-29 ,
* -1.3103660877630D-32
1 /
DATA B20,B21,B22,B23,B24,B25,B26,B27,B28,B29,B210,B211,B212 /
* -2.2063755868780D-07 ,
* 1.4268968388826D-08 ,
* -3.9416351705981D-10 ,
* 5.9606582726984D-12 ,
* -5.4871340648895D-14 ,
* 3.2304069082462D-16 ,
* -1.2421935589492D-18 ,
* 3.1127728952968D-21 ,
* -4.9182825188505D-24 ,
* 4.4404482566015D-27 ,
* -1.5665047149903D-30 ,
* -5.6075256366923D-34 ,
* 4.6890466767990D-37
1 /
C0= A00 + H*(A01+H*(A02+H*(A03+H*(A04+H*(A05+H*(A06+H*
1 (A07+H*(A08+H*(A09+H*(A010+H*(A011+H*A012)))))))))))
C1= A10 + H*(A11+H*(A12+H*(A13+H*(A14+H*(A15+H*(A16+H*
1 (A17+H*(A18+H*(A19+H*(A110+H*(A111+H*A112)))))))))))
C2= A20 + H*(A21+H*(A22+H*(A23+H*(A24+H*(A25+H*(A26+H*
1 (A27+H*(A28+H*(A29+H*(A210+H*(A211+H*A212)))))))))))
C3= A30 + H*(A31+H*(A32+H*(A33+H*(A34+H*(A35+H*(A36+H*
1 (A37+H*(A38+H*(A39+H*(A310+H*(A311+H*A312)))))))))))
D0= B00 + H*(B01+H*(B02+H*(B03+H*(B04+H*(B05+H*(B06+H*
1 (B07+H*(B08+H*(B09+H*(B010+H*(B011+H*B012)))))))))))
D1= B10 + H*(B11+H*(B12+H*(B13+H*(B14+H*(B15+H*(B16+H*
1 (B17+H*(B18+H*(B19+H*(B110+H*(B111+H*B112)))))))))))
D2= B20 + H*(B21+H*(B22+H*(B23+H*(B24+H*(B25+H*(B26+H*
1 (B27+H*(B28+H*(B29+H*(B210+H*(B211+H*B212)))))))))))
DD0=B01+H*(2.*B02+H*(3.*B03+H*(4.*B04+H*(5.*B05+H*(6.*B06+
1 H*(7.*B07+H*(8.*B08+H*(9.*B09+H*(10.*B010+H*(11.*B011+
2 12.*B012*H))))))))))
DD1=B11+H*(2.*B12+H*(3.*B13+H*(4.*B14+H*(5.*B15+H*(6.*B16+
1 H*(7.*B17+H*(8.*B18+H*(9.*B19+H*(10.*B110+H*(11.*B111+
2 12.*B112*H))))))))))
DD2=B21+H*(2.*B22+H*(3.*B23+H*(4.*B24+H*(5.*B25+H*(6.*B26+
1 H*(7.*B27+H*(8.*B28+H*(9.*B29+H*(10.*B210+H*(11.*B211+
2 12.*B212*H))))))))))
DVDH1 = A01 + H*(2.*A02+H*(3.*A03+H*(4.*A04+H*(5.*A05+H*(
1 6.*A06+H*(7.*A07+H*(8.*A08+H*(9.*A09+H*(10.*A010+H*(
2 11.*A011+H*12.*A012))))))))))
DVDH2 = A11 + H*(2.*A12+H*(3.*A13+H*(4.*A14+H*(5.*A15+H*(
1 6.*A16+H*(7.*A17+H*(8.*A18+H*(9.*A19+H*(10.*A110+H*(
2 11.*A111+H*12.*A112))))))))))
DVDH3 = A21 + H*(2.*A22+H*(3.*A23+H*(4.*A24+H*(5.*A25+H*(
1 6.*A26+H*(7.*A27+H*(8.*A28+H*(9.*A29+H*(10.*A210+H*(
2 11.*A211+H*12.*A212))))))))))
DVDH4 = A31 + H*(2.*A32+H*(3.*A33+H*(4.*A34+H*(5.*A35+H*(
1 6.*A36+H*(7.*A37+H*(8.*A38+H*(9.*A39+H*(10.*A310+H*(
2 11.*A311+H*12.*A312))))))))))
DVDHP = DVDH1 + P*(DVDH2+P*(DVDH3+P*DVDH4))
DVDPH = C1 + 2. * P * C2 + 3. * P**2 * C3
V = C0 + P*(C1+P*(C2+P*C3))
T = D0 + P*(D1+P*D2)
DTDH=DD0+P*(DD1+P*DD2)
DTDP=D1+2.*D2*P
RHO = 1. / V
RETURN
END