50    if (temperature > H2O::criticalTemperature())
 
   51        temperature = H2O::criticalTemperature();
 
   53    Scalar Tr = temperature/H2O::criticalTemperature();
 
   56    static const Scalar c[6] = {
 
   57        1.99274064, 1.09965342, -0.510839303,
 
   58        -1.75493479,-45.5170352, -6.7469445e5
 
   60    static const Scalar d[6] = {
 
   62        16/3.0, 43/3.0, 110/3.0
 
   64    static const Scalar q = -0.023767;
 
   68    for (
int i = 0; i < 6; ++i) {
 
   69        f += c[i]*pow(tau, d[i]);
 
   78        exp((H2O::tripleTemperature() - temperature)/100);
 
   81    return exp(exponent)*H2O::vaporPressure(temperature);
 
 
Scalar henryIAPWS(Scalar E, Scalar F, Scalar G, Scalar H, Scalar temperature)
The Henry constants in liquid water using the IAPWS 2004 formulation.
Definition henryiapws.hh:40