def test_mass_balance(self): f = DeltaBW(friction = "Chezy", flood_discharge = 6., flood_intermittency = 1., bedload_discharge_input = 0.001, grain_size = 0.5, coeff_Chezy = 15., exponent_load_relation = 2.5, critical_Shields_stress = 0., foreset_elev__top__init = 3., foreset_elev__base__init = 0., bed_slope__init = 0.00025, basin_slope = 0., domain_length = 10000., domain_length_max = 500000., foreset_slope = 0.2, bed_porosity = 0.4, coeff_Manning_ks = 0.3, coeff_bedload = 7.2, coeff_Manning = 8.1, dx = 500., dt__days = 0.182625, num_iterations = 10000, print_step = 5000, density_sediment = 2650., density_water = 1000., basin_water_surface_elevation = 8.5, verbose = False, save_output = False) f.run() f.finalize() self.assertEqual(f.sed_volume__init, 26338.5) self.assertEqual(f.sed_volume__final, 190680.45389578832) self.assertEqual(f.sed_volume__feed, 157788.00000000003) self.assertEqual(f.mass_balance_error, 0.034371398650908352)
def test_stats(self): f = DeltaBW(friction = "Chezy", flood_discharge = 6., flood_intermittency = 1., bedload_discharge_input = 0.001, grain_size = 0.5, coeff_Chezy = 15., exponent_load_relation = 2.5, critical_Shields_stress = 0., foreset_elev__top__init = 3., foreset_elev__base__init = 0., bed_slope__init = 0.00025, basin_slope = 0., domain_length = 10000., domain_length_max = 500000., foreset_slope = 0.2, bed_porosity = 0.4, coeff_Manning_ks = 0.3, coeff_bedload = 7.2, coeff_Manning = 8.1, dx = 500., dt__days = 0.182625, num_iterations = 10000, print_step = 5000, density_sediment = 2650., density_water = 1000., basin_water_surface_elevation = 8.5, verbose = False, save_output = False) f.run() self.assertEqual(f.domain_length, 34481.071182531079) self.assertEqual(f.bed_length, 34500.977813788544) self.assertEqual(f.eta_upstream, 15.154555352651583) self.assertEqual(f.eta[-2], 3.9813262514556356) self.assertEqual(f.eta[-1], 0.0)