|
141 | 141 | "metadata": {}, |
142 | 142 | "outputs": [], |
143 | 143 | "source": [ |
144 | | - "from gadopt import *" |
| 144 | + "from gadopt import *\n", |
| 145 | + "from firedrake.pyplot import tripcolor" |
145 | 146 | ] |
146 | 147 | }, |
147 | 148 | { |
|
290 | 291 | "Ra, delta_t = Constant(1e4), Constant(1e-3)\n", |
291 | 292 | "\n", |
292 | 293 | "approximation = BoussinesqApproximation(Ra)\n", |
293 | | - "t_adapt = TimestepAdaptor(delta_t, V, maximum_timestep=0.1, increase_tolerance=1.5)" |
| 294 | + "t_adapt = TimestepAdaptor(delta_t, u, V, maximum_timestep=0.1, increase_tolerance=1.5)" |
294 | 295 | ] |
295 | 296 | }, |
296 | 297 | { |
|
377 | 378 | "gd = GeodynamicalDiagnostics(u, p, T, bottom, top)\n", |
378 | 379 | "\n", |
379 | 380 | "for timestep in range(0, no_timesteps+1):\n", |
380 | | - " dt = t_adapt.update_timestep(u)\n", |
| 381 | + " dt = t_adapt.update_timestep()\n", |
381 | 382 | " time += dt\n", |
382 | 383 | "\n", |
383 | 384 | " stokes_solver.solve()\n", |
|
521 | 522 | "gd = GeodynamicalDiagnostics(u, p, T, bottom, top)\n", |
522 | 523 | "\n", |
523 | 524 | "for timestep in range(0, no_timesteps+1):\n", |
524 | | - " dt = t_adapt.update_timestep(u)\n", |
| 525 | + " dt = t_adapt.update_timestep()\n", |
525 | 526 | " time += dt\n", |
526 | 527 | "\n", |
527 | 528 | " stokes_solver.solve()\n", |
|
0 commit comments