[cig-commits] [commit] update_examples: Reverted changes to solidsolution.py (315226f)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Wed Dec 31 15:38:35 PST 2014


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

On branch  : update_examples
Link       : https://github.com/geodynamics/burnman/compare/bf39df300a99fe98a6add654cf40cc69ed4206af...315226f591a5794325c8c13b2db9571ce90ac130

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

commit 315226f591a5794325c8c13b2db9571ce90ac130
Author: Bob Myhill <myhill.bob at gmail.com>
Date:   Wed Dec 31 23:38:12 2014 +0000

    Reverted changes to solidsolution.py


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

315226f591a5794325c8c13b2db9571ce90ac130
 burnman/solidsolution.py | 11 +++--------
 1 file changed, 3 insertions(+), 8 deletions(-)

diff --git a/burnman/solidsolution.py b/burnman/solidsolution.py
index ff90993..1b785c7 100644
--- a/burnman/solidsolution.py
+++ b/burnman/solidsolution.py
@@ -81,14 +81,9 @@ class SolidSolution(Mineral):
         self.partial_gibbs = np.array([self.base_material[i][0].gibbs for i in range(self.n_endmembers)]) + self.excess_partial_gibbs
         self.gibbs= sum([ self.base_material[i][0].gibbs * self.molar_fraction[i] for i in range(self.n_endmembers) ]) + self.excess_gibbs
 
-        if self.solution_model.__class__.__name__ == 'IdealSolution':
-            self.excess_enthalpy = 0.
-            self.excess_entropy = 0.
-            self.excess_volume = 0.
-        else:
-            self.excess_enthalpy = self.solution_model.excess_enthalpy( pressure, temperature, self.molar_fraction)
-            self.excess_entropy = self.solution_model.excess_entropy( pressure, temperature, self.molar_fraction)
-            self.excess_volume = self.solution_model.excess_volume( pressure, temperature, self.molar_fraction)
+        self.excess_enthalpy = self.solution_model.excess_enthalpy( pressure, temperature, self.molar_fraction)
+        self.excess_entropy = self.solution_model.excess_entropy( pressure, temperature, self.molar_fraction)
+        self.excess_volume = self.solution_model.excess_volume( pressure, temperature, self.molar_fraction)
 
         self.H = sum([ self.base_material[i][0].H * self.molar_fraction[i] for i in range(self.n_endmembers) ]) + self.excess_enthalpy
         self.S = sum([ self.base_material[i][0].S * self.molar_fraction[i] for i in range(self.n_endmembers) ]) + self.excess_entropy



More information about the CIG-COMMITS mailing list