[cig-commits] [commit] inversion, master, validate_MT_params: Added comment for origin of formula in benchmark (c9e4070)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Dec 12 18:26:15 PST 2014
Repository : https://github.com/geodynamics/burnman
On branches: inversion,master,validate_MT_params
Link : https://github.com/geodynamics/burnman/compare/80c2a295c42dfdb38f83f6c1334bf7d8f97a8463...409647ff05dfad6a686198cac1481bd46b5e2e62
>---------------------------------------------------------------
commit c9e4070c2c8c56b6284c03f52859aeb32578bb14
Author: Bob Myhill <myhill.bob at gmail.com>
Date: Thu Sep 4 00:17:11 2014 +0200
Added comment for origin of formula in benchmark
>---------------------------------------------------------------
c9e4070c2c8c56b6284c03f52859aeb32578bb14
solidsolution_benchmarks.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/solidsolution_benchmarks.py b/solidsolution_benchmarks.py
index 3b9deb2..f254a85 100644
--- a/solidsolution_benchmarks.py
+++ b/solidsolution_benchmarks.py
@@ -56,11 +56,11 @@ for i,c in enumerate(comp):
sp.set_composition( np.array(molar_fractions) )
sp.set_state( 1e5, 298.15 )
sp_entropies[i] = sp.solution_model.configurational_entropy( molar_fractions )
- sp_S[i] = -8.3145*(c*np.log(c) + (1.-c)*np.log(1.-c) + c*np.log(c/2.) + (2.-c)*np.log(1.-c/2.))
+ sp_entropies_NK1967[i] = -8.3145*(c*np.log(c) + (1.-c)*np.log(1.-c) + c*np.log(c/2.) + (2.-c)*np.log(1.-c/2.)) # eq. 7 in Navrotsky and Kleppa, 1967.
#fig1 = mpimg.imread('configurational_entropy.png') # Uncomment these two lines if you want to overlay the plot on a screengrab from SLB2011
#plt.imshow(fig1, extent=[0.0, 1.0,0.,17.0], aspect='auto')
-plt.plot( comp, sp_S, 'b-', linewidth=3.)
+plt.plot( comp, sp_entropies_NK1967, 'b-', linewidth=3.)
plt.plot( comp, sp_entropies, 'r--', linewidth=3.)
plt.xlim(0.0,1.0)
plt.ylim(0.,17.0)
More information about the CIG-COMMITS
mailing list