[cig-commits] [commit] master: add moduli to seismic models (e0a47a1)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/burnman/compare/d4a00047c3fa7097111958d0b97cb77140ee94e3...4f4886fe01e0d6065e0a8ade5a079caacfe81cc9

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

commit e0a47a1b69ad6d38d4132908e3f2daf92dfa9366
Author: sannecottaar <sanne.cottaar at gmail.com>
Date:   Tue Dec 9 09:32:50 2014 -0800

    add moduli to seismic models


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

e0a47a1b69ad6d38d4132908e3f2daf92dfa9366
 burnman/seismic.py | 20 ++++++++++++++++++++
 1 file changed, 20 insertions(+)

diff --git a/burnman/seismic.py b/burnman/seismic.py
index 5da0382..0fe0a44 100644
--- a/burnman/seismic.py
+++ b/burnman/seismic.py
@@ -138,6 +138,26 @@ class Seismic1DModel:
         raise ValueError, "not implemented"
         return 0
     
+    
+    def G(self, depth):
+        """
+        Parameters
+        ----------
+        depth : float or array of floats
+        Shear modulus at given for depth(s) in [Pa].
+        """
+        return np.power(self.v_s(depth),2.) * self.density(depth)
+    
+    def K(self, depth):
+        """
+        Parameters
+        ----------
+        depth : float or array of floats
+        Bulk modulus at given for depth(s) in [Pa]
+        """
+        return np.power(self.v_phi(depth),2.) * self.density(depth)
+
+    
     def depth(self, pressure):
         """
         Parameters