From 8d295d2ba1b9d1eb4e7d8743d0fef6cbb8bfad55 Mon Sep 17 00:00:00 2001
From: Jean-Baptiste Keck <Jean-Baptiste.Keck@imag.fr>
Date: Fri, 14 Dec 2018 22:43:19 +0100
Subject: [PATCH] update levelset example
---
examples/sediment_deposit/sediment_deposit_levelset.py | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/examples/sediment_deposit/sediment_deposit_levelset.py b/examples/sediment_deposit/sediment_deposit_levelset.py
index 53e73d40a..2aadf45fa 100644
--- a/examples/sediment_deposit/sediment_deposit_levelset.py
+++ b/examples/sediment_deposit/sediment_deposit_levelset.py
@@ -8,7 +8,7 @@ TANK_RATIO = 1
FILL_PCT = 0.25
SEDIMENT_RADIUS = 5e-3
SEDIMENT_COUNT = int(1.15*(FILL_PCT*TANK_RATIO / (np.pi*SEDIMENT_RADIUS**2)))
-DISCRETIZATION = 2048
+DISCRETIZATION = 256
BLOB_INIT = False
NB = 5
@@ -322,7 +322,7 @@ def compute(args):
**extra_op_kwds)
#> Operator to compute and save mean fields
- axes = list(range(0, dim-1))
+ axes = list(range(1, dim))
view = [slice(None,None,None),]*dim
view = tuple(view)
io_params = IOParams(filename='horizontally_averaged_profiles', frequency=0)
@@ -334,7 +334,7 @@ def compute(args):
### Adaptive timestep operator
adapt_dt = AdaptiveTimeStep(dt, equivalent_CFL=True,
name='merge_dt', pretty_name='dt',
- max_dt=1)
+ max_dt=1e-1)
dt_cfl = adapt_dt.push_cfl_criteria(cfl=args.cfl,
Finf=min_max_U.Finf,
equivalent_CFL=True,
@@ -430,8 +430,8 @@ if __name__=='__main__':
npts=(TANK_RATIO*DISCRETIZATION+1,DISCRETIZATION+1),
box_origin=(0.0,), box_length=(1.0,),
tstart=0.0, tend=100.1,
- dt=1e-6, cfl=8.00, lcfl=0.50,
- dump_times=tuple(float(5*x) for x in range(20)),
+ dt=1e-6, cfl=0.50, lcfl=0.50,
+ dump_times=tuple(float(1*x) for x in range(100)),
dump_freq=0)
parser.run(compute)
--
GitLab