def create_sediment_model(self, mesh, bathymetry): self.P1DG = FunctionSpace(mesh, "DG", 1) self.P1_vec_dg = VectorFunctionSpace(mesh, "DG", 1) #if uv_init is None: self.uv_d = Function(self.P1_vec_dg).project(self.uv_init) #if elev_init is None: self.eta_d = Function(self.P1DG).project(self.elev_init) self.sediment_model = SedimentModel( ModelOptions2d, suspendedload=self.suspended, convectivevel=self.convective_vel_flag, bedload=self.bedload, angle_correction=self.angle_correction, slope_eff=self.slope_eff, seccurrent=False, sediment_slide=self.sediment_slide, mesh2d=mesh, bathymetry_2d=bathymetry, dt=self.dt, uv_init=self.uv_d, elev_init=self.eta_d, ks=self.ks, average_size=self.average_size, porosity=self.porosity, max_angle=self.max_angle, meshgrid_size=self.meshgrid_size, cons_tracer=self.use_tracer_conservative_form, morfac=self.morphological_acceleration_factor, wetting_and_drying=self.wetting_and_drying, wetting_alpha=self.wetting_and_drying_alpha)
def create_sediment_model(self, mesh, bathymetry): self.P1DG = FunctionSpace(mesh, "DG", 1) self.P1_vec_dg = VectorFunctionSpace(mesh, "DG", 1) # TODO: Do these really need to be attributes? Would prefer to avoid stashing func spaces self.uv_d = Function(self.P1_vec_dg).project(self.uv_init) self.eta_d = Function(self.P1DG).project(self.elev_init) self.sediment_model = SedimentModel( # TODO: Use kwargs dict format ModelOptions2d, suspendedload=self.suspended, convectivevel=self.convective_vel_flag, bedload=self.bedload, angle_correction=self.angle_correction, slope_eff=self.slope_eff, seccurrent=False, mesh2d=mesh, bathymetry_2d=bathymetry, uv_init=self.uv_d, elev_init=self.eta_d, ks=self.ks, average_size=self.average_size, cons_tracer=self.use_tracer_conservative_form, wetting_and_drying=self.wetting_and_drying, wetting_alpha=self.wetting_and_drying_alpha)
def create_sediment_model(self, mesh, bathymetry): self.P1DG = FunctionSpace(mesh, "DG", 1) self.P1_vec_dg = VectorFunctionSpace(mesh, "DG", 1) self.uv_d = Function(self.P1_vec_dg).project(self.uv_init) self.eta_d = Function(self.P1DG).project(self.elev_init) self.sediment_model = SedimentModel(ModelOptions2d, suspendedload=self.suspended, convectivevel=self.convective_vel_flag, bedload=self.bedload, angle_correction=self.angle_correction, slope_eff=self.slope_eff, seccurrent=False, mesh2d=mesh, bathymetry_2d=bathymetry, uv_init=self.uv_d, elev_init=self.eta_d, ks=self.ks, average_size=self.average_size, cons_tracer=self.use_tracer_conservative_form, wetting_and_drying=self.wetting_and_drying, wetting_alpha=self.wetting_and_drying_alpha)