Commit 48eb988b authored by Henning Francke's avatar Henning Francke
Browse files

clean up

parent 0f4b751f
......@@ -12,14 +12,6 @@ protected
//print("x_H2_0: "+String(x_H2_0");
Real[2] a = {0.018519, -0.30185103};
algorithm
/* if AssertLevel>0 then
assert(ignoreTlimit or ignoreLimitH2_T or ((if molalities[iNaCl] > 0 then 273.15 else 323.15)<T and T<373.15), "Temperature out of validity range: T=" + String(T - 273.15) + ".\nTo ignore set ignoreLimitN2_T=true",aLevel);
assert(ignoreLimitH2_p or ((if molalities[iNaCl] > 0 then 10e5 else 1e5)<p and p<203e5),"Pressure out of validity rangep=" + String(p/1e5) + " bar.\nTo ignore set ignoreLimitN2_p=true",aLevel);
assert(ignoreLimitH2_b or molalities[iNaCl]<5,"Molality out of validity range: mola[NaCl]=" + String(molalities[iNaCl]) + " mol/kg.\nTo ignore set ignoreLimitH2_b=true",aLevel);
end if;
*/
//Debug.Print "ln(xH2/xH2°)"; a1 * b_NaCl ^ 2 + a2 * b_NaCl 'equ. 12
y_H2 :=y_H2_0*exp(a[1]*b_NaCl^2 + a[2]*b_NaCl); //x_H2 (equ. 12)
// print("y_H2: " + String( y_H2)); //equ. 13
end solubility_H2_pTb_Chabab2020_y;
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment