[cig-commits] [commit] add_gibbs_energy: Corrected sign of S excess (43d6979)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/burnman

On branch  : add_gibbs_energy
Link       : https://github.com/geodynamics/burnman/compare/0000000000000000000000000000000000000000...2148b324d3e8aa7b527f831eb397590942563008

>---------------------------------------------------------------

commit 43d69799b72eaa7177d83bbd24fa74c844afe527
Author: Bob Myhill <myhill.bob at gmail.com>
Date:   Sat Aug 30 21:41:34 2014 +0200

    Corrected sign of S excess


>---------------------------------------------------------------

43d69799b72eaa7177d83bbd24fa74c844afe527
 burnman/solidsolution.py | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/burnman/solidsolution.py b/burnman/solidsolution.py
index 44198d6..18a3ee6 100644
--- a/burnman/solidsolution.py
+++ b/burnman/solidsolution.py
@@ -135,9 +135,8 @@ class SolidSolution(Mineral):
         phi=np.divide(phi, np.sum(phi))
 
         self.H_excess=np.dot(self.alpha.T,molar_fraction)*np.dot(phi.T,np.dot(self.Wh,phi))
-        self.S_deficit=np.dot(self.alpha.T,molar_fraction)*np.dot(phi.T,np.dot(self.Ws,phi))
+        self.S_excess=np.dot(self.alpha.T,molar_fraction)*np.dot(phi.T,np.dot(self.Ws,phi))
         self.V_excess=np.dot(self.alpha.T,molar_fraction)*np.dot(phi.T,np.dot(self.Wv,phi))
-        self.S_excess=0.0-S_deficit
 
     def set_state(self, pressure, temperature, molar_fraction):
         # Set the state of all the endmembers