wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-08-11 11:26:01 +08:00
Code Issues Packages Projects Releases Wiki Activity

Make working in log conductivity a bit more explicit.

Browse Source 
Should be fixed in the actual implementation, and use SimPEG maps
instead of forcing the user to work in weird units.
This commit is contained in:
rowanc1
2014-07-03 09:56:05 -07:00
parent 23141c682f
commit f747d724c9
2 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
docs/examples/richards_comparisonToCeiliaEtAl1990.py
+1
View File
@@ -9,6 +9,7 @@ import matplotlib.pyplot as plt
M = Mesh.TensorMesh([np.ones(40)])
M.setCellGradBC('dirichlet')
params = Richards.Empirical.HaverkampParams().celia1990
params['Ks'] = np.log(params['Ks'])
E = Richards.Empirical.Haverkamp(M, **params)
bc = np.array([-61.5,-20.7])
simpegFLOW/Richards/Empirical.py
+3 -3
View File
@@ -153,7 +153,7 @@ class HaverkampParams(object):
"""
return {'alpha':1.611e+06, 'beta':3.96,
'theta_r':0.075, 'theta_s':0.287,
'Ks':np.log(9.44e-03), 'A':1.175e+06,
'Ks':9.44e-03, 'A':1.175e+06,
'gamma':4.74}
@@ -243,7 +243,7 @@ class Haverkamp(RichardsMap):
theta_r = _ModelProperty('theta_r', ['thetaModel'], default=0.075)
theta_s = _ModelProperty('theta_s', ['thetaModel'], default=0.287)
Ks = _ModelProperty('Ks', ['kModel'], default=np.log(24.96))
Ks = _ModelProperty('Ks', ['kModel'], default=24.96)
A = _ModelProperty('A', ['kModel'], default=1.1750e+06)
gamma = _ModelProperty('gamma', ['kModel'], default=4.74)
@@ -357,7 +357,7 @@ class VanGenuchten(RichardsMap):
alpha = _ModelProperty('alpha', ['thetaModel', 'kModel'], default=0.036)
n = _ModelProperty('n', ['thetaModel', 'kModel'], default=1.560)
Ks = _ModelProperty('Ks', ['kModel'], default=np.log(24.96))
Ks = _ModelProperty('Ks', ['kModel'], default=24.96)
I = _ModelProperty('I', ['kModel'], default=0.500)
def __init__(self, mesh, **kwargs):