[cig-commits] [commit] solid_solution_molar_mass: Add arithmetic and harmonic averages for G in solidsolution (thanks to Mark Panning for pointing out they were missing) (30107f8)

[cig_noreply at geodynamics.org]

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

On branch  : solid_solution_molar_mass
Link       : https://github.com/geodynamics/burnman/compare/0000000000000000000000000000000000000000...fe211b4b277a753e04e153250e4ed8b270dee9b7

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

commit 30107f8bfcb75e1414ee388fa9390cee3e55b617
Author: ian-r-rose <ian.r.rose at gmail.com>
Date:   Tue Dec 16 20:58:46 2014 -0800

    Add arithmetic and harmonic averages for G in solidsolution (thanks to Mark Panning for pointing out they were missing)


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

30107f8bfcb75e1414ee388fa9390cee3e55b617
 burnman/solidsolution.py | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/burnman/solidsolution.py b/burnman/solidsolution.py
index 6b4a905..1287d23 100644
--- a/burnman/solidsolution.py
+++ b/burnman/solidsolution.py
@@ -98,6 +98,8 @@ class SolidSolution(Mineral):
         self.C_p = sum([ self.base_material[i][0].C_p * self.molar_fraction[i] for i in range(self.n_endmembers) ])
         self.alpha = (1./self.V) * sum([ self.base_material[i][0].alpha * self.base_material[i][0].V * self.molar_fraction[i] for i in range(self.n_endmembers) ])
         self.K_T = self.V * 1./(sum([ self.base_material[i][0].V / (self.base_material[i][0].K_T)  * self.molar_fraction[i] for i in range(self.n_endmembers) ]))
+#        self.G = self.V * 1./(sum([ self.base_material[i][0].V / (self.base_material[i][0].G)  * self.molar_fraction[i] for i in range(self.n_endmembers) ]))
+#        self.G = (1./self.V) * sum([ self.base_material[i][0].G * self.base_material[i][0].V * self.molar_fraction[i] for i in range(self.n_endmembers) ])
 
         # Derived properties
         self.C_v = self.C_p - self.V*temperature*self.alpha*self.alpha*self.K_T