tmpa1 = - 171.9065 - (0.077993*tk) + (2839.319/tk)
& + (71.595*log10(tk))
tmpa2 = +(-0.77712 + (0.0028426*tk) + (178.34/tk) )*sqrts
tmpa1 = - 171.945 - (0.077993*tk) + (2903.293/tk)
tmpa2 = +(-0.068393 + (0.0017276*tk) + (88.135/tk) )*sqrts
_RL tk
tk = 273.15 _d 0 + t
tk100 = tk/100. _d 0
invtk=1.0 _d 0/tk
dlogtk=log(tk)
RT = Rgas*tk
delta = (57.7 _d 0 - 0.118 _d 0*tk)
B1 = -1636.75 _d 0 + 12.0408 _d 0*tk - 0.0327957 _d 0*tk*tk
B = B1 + 3.16528 _d 0*tk*tk*tk*(0.00001 _d 0)
& dlogtk + 0.053105 _d 0*sqrts*tk)
tk = 273.15 + t
tk100 = tk/100.0
invtk=1.0/tk
& exp( (24.2-0.085*t)*(pressc-1.0)/(83.143*tk) )
& exp( (16.4-0.040*t)*(pressc-1.0)/(83.143*tk) )
& dlogtk + 0.053105*sqrts*tk)
pfactor = - (dv/(bigR*tk))*pressc
& + (0.5*dk/(bigR*tk))*pressc*pressc