update example to work with new problem formulation in simpeg

2026-06-30 22:03:56 +08:00 · 2014-07-03 07:50:25 -07:00
parent 0718c2c50a
commit 23141c682f
1 changed files with 2 additions and 1 deletions
@@ -16,7 +16,8 @@ h = np.zeros(M.nC) + bc[0]


 def getFields(timeStep,method):
-    prob = Richards.RichardsProblem(M, mapping=E, timeStep=timeStep, timeEnd=360,
+    timeSteps = np.ones(360/timeStep)*timeStep
+    prob = Richards.RichardsProblem(M, mapping=E, timeSteps=timeSteps,
                                    boundaryConditions=bc, initialConditions=h,
                                    doNewton=False, method=method)
    return prob.fields(params['Ks'])