def test_eta(self): f = AgDegBW(friction = "Manning", Qf = 400, If = 0.1, S = 0.0002, initdepth = 15, B = 60, L = 50000, dx = 833, Gtf = 187059, D90 = 0.6, porosity = 0.4, kc = 75, phi = 1, alpha_r = 8.1, alpha_t = 8., tau_c = 0.047, n = 1.5, Cz = 22, alpha_u = 1, dt = 0.1, totalT = 120, print_dt = 120, g = 9.81, rho_w = 1000, rho_s = 2650, verbose = False, save_output = False) f.run() self.assertEqual(f.bed_elevation[0], 16.150678907133734) self.assertEqual(f.bed_elevation[-2], 0.23969483388856572)
def test_shear_stress(self): f = AgDegBW(friction = "Manning", Qf = 400, If = 0.1, S = 0.0002, initdepth = 15, B = 60, L = 50000, dx = 833, Gtf = 187059, D90 = 0.6, porosity = 0.4, kc = 75, phi = 1, alpha_r = 8.1, alpha_t = 8., tau_c = 0.047, n = 1.5, Cz = 22, alpha_u = 1, dt = 0.1, totalT = 120, print_dt = 120, g = 9.81, rho_w = 1000, rho_s = 2650, verbose = False, save_output = False) f.run() self.assertEqual(f.bed_shear_stress[0], 0.89439215587050513) self.assertEqual(f.bed_shear_stress[-2], 0.055014859454631616)
def test_water_surface(self): f = AgDegBW(friction = "Manning", Qf = 400, If = 0.1, S = 0.0002, initdepth = 15, B = 60, L = 50000, dx = 833, Gtf = 187059, D90 = 0.6, porosity = 0.4, kc = 75, phi = 1, alpha_r = 8.1, alpha_t = 8., tau_c = 0.047, n = 1.5, Cz = 22, alpha_u = 1, dt = 0.1, totalT = 120, print_dt = 120, g = 9.81, rho_w = 1000, rho_s = 2650, verbose = False, save_output = False) f.run() self.assertEqual(f.water_surface_elevation[0], 20.619136144306204) self.assertEqual(f.water_surface_elevation[-2], 15.002769439775431)