<s> r"""<STR_LIT>""" <EOL> import scipy as sp <EOL> def water ( phase , <EOL> pore_T = '<STR_LIT>' , <EOL> pore_salinity = '<STR_LIT>' , <EOL> ** kwargs ) : <EOL> r"""<STR_LIT>""" <EOL> T = phase [ '<STR_LIT>' ] <EOL> try : <EOL> S = phase [ '<STR_LIT>' ] <EOL> except : <EOL> S = <NUM_LIT:0> <EOL> sigma_w = <NUM_LIT> * ( ( <NUM_LIT:1> - ( T / <NUM_LIT> ) ) ** <NUM_LIT> ) * ( <NUM_LIT:1> - <NUM_LIT> * ( <NUM_LIT:1> - ( T / <NUM_LIT> ) ) ) <EOL> a1 = <NUM_LIT> <EOL> a2 = <NUM_LIT> <EOL> a3 = <NUM_LIT> <EOL> TC = T - <NUM_LIT> <EOL> sigma_sw = sigma_w * ( <NUM_LIT:1> + ( a1 * TC + a2 ) * sp . log ( <NUM_LIT:1> + a3 * S ) ) <EOL> value = sigma_sw <EOL> return value <EOL> def eotvos ( phase , k , <EOL> pore_T = '<STR_LIT>' , <EOL> pore_Tc = '<STR_LIT>' , <EOL> pore_molar_density = '<STR_LIT>' , ** kwargs ) : <EOL> r"""<STR_LIT>""" <EOL> Tc = phase [ pore_Tc ] <EOL> T = phase [ pore_T ] <EOL> Vm = <NUM_LIT:1> / phase [ pore_molar_density ] <EOL> value = k * ( Tc - T ) / ( Vm ** ( <NUM_LIT:2> / <NUM_LIT:3> ) ) <EOL> return value <EOL> def guggenheim_katayama ( phase , K2 , n , <EOL> pore_T = '<STR_LIT>' , <EOL> pore_Tc = '<STR_LIT>' , <EOL> pore_Pc = '<STR_LIT>' , <EOL> ** kwargs ) : <EOL> r"""<STR_LIT>""" <EOL> T = phase [ pore_T ] <EOL> Pc = phase [ pore_Pc ] <EOL> Tc = phase [ pore_Tc ] <EOL> sigma_o = K2 * Tc ** ( <NUM_LIT:1> / <NUM_LIT:3> ) * Pc ** ( <NUM_LIT:2> / <NUM_LIT:3> ) <EOL> value = sigma_o * ( <NUM_LIT:1> - T / Tc ) ** n <EOL> return value <EOL> def brock_bird_scaling ( phase , sigma_o , To , <EOL> pore_T = '<STR_LIT>' , <EOL> pore_Tc = '<STR_LIT>' , <EOL> ** params ) : <EOL> r"""<STR_LIT>""" <EOL> Tc = phase [ pore_Tc ] <EOL> Ti = phase [ pore_T ] <EOL> Tro = To / Tc <EOL> Tri = Ti / Tc <EOL> value = sigma_o * ( <NUM_LIT:1> - Tri ) ** ( <NUM_LIT:11> / <NUM_LIT:9> ) / ( <NUM_LIT:1> - Tro ) ** ( <NUM_LIT:11> / <NUM_LIT:9> ) <EOL> return value </s>
