Commit 87e30ca2 authored by David Vincent's avatar David Vincent

fix some minor problem for sediment 2D

parent 0d2ab455
Pipeline #1881 passed with stage
in 42 minutes and 57 seconds
......@@ -817,8 +817,6 @@ class ShallowWaterTracer2d:
- float [in seconds]
- date string in format "year-month-day hours:minutes:seconds" (e.g. "2000-06-20 22:42:17")
* is_sediment
use shallow water tracer law to predict the sediment transport (default: False)
"""
fmt = '%Y-%m-%d %H:%M:%S'
if initial_time is None:
......@@ -995,9 +993,12 @@ class ShallowWaterTracer2d:
* windV
netcdf or .msh file containing the surface wind velocity along the y-axis in the local basis [in m*s^-1].
"""
self._is_sediment = True
self._sediment_bottom = dgpy.dgDofContainer(self._domain._groups, 1)
slim_private._load(self._sediment_bottom, initial_bottom_concentration)
self._wind = slim_private._check_vector(windU, windV, None, self._domain)
_windU=slim_private._load_function(windU, self._domain._groups)
_windV=slim_private._load_function(windV, self._domain._groups)
self._wind = slim_private._check_vector(_windU, _windV, None, self._domain)
def compute_sediment(self,dt):
"""Compute the sediment flux
......
......@@ -36,7 +36,7 @@ class _sum:
bath = cm.get(self._bath)
nPt = sol.shape
val = np.zeros([nPt[0], 1])
val[:,0] = sol[:,0] + bath[:]
val[:,0] = sol[:,0] + bath[:,0]
return val
class _get_x:
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment