예제 #1
0
    def test_shear_stress(self):
    
        f = AgDegNormalSub(friction = 'Manning',
							 Qf = 70,
							 If = 0.03,
							 L = 10000,
							 dx = 100,
							 S = 0.008,
							 Gtf = 700000,
							 D = 30,
							 porosity = 0.35,
							 kc = 75,
							 phi = 1,
							 alpha_r = 8.1,
							 alpha_t = 8.,
							 tau_c = 0.047,
							 nt = 1.5,
							 Cz = 22,
							 alpha_u = 0.5,
							 dt = 0.01,
							 totalT = 120,
							 print_dt = 120,
							 g = 9.81,
							 rho_w = 1000,
							 rho_s = 2650,
							 subsidence = 5,
							 channel_width = 25,
							 depositional_width = 25*60,
							 washload_fraction = 1.5,
							 sinuosity = 1.5,
							 verbose = False,
							 save_output = False)
        f.run()
        self.assertEqual(f.bed_shear_stress[0], 0.21046460651748258)
        self.assertEqual(f.bed_shear_stress[-2], 0.14610015954984537)
예제 #2
0
    def test_water_surface(self):
    
        f = AgDegNormalSub(friction = 'Manning',
							 Qf = 70,
							 If = 0.03,
							 L = 10000,
							 dx = 100,
							 S = 0.008,
							 Gtf = 700000,
							 D = 30,
							 porosity = 0.35,
							 kc = 75,
							 phi = 1,
							 alpha_r = 8.1,
							 alpha_t = 8.,
							 tau_c = 0.047,
							 nt = 1.5,
							 Cz = 22,
							 alpha_u = 0.5,
							 dt = 0.01,
							 totalT = 120,
							 print_dt = 120,
							 g = 9.81,
							 rho_w = 1000,
							 rho_s = 2650,
							 subsidence = 5,
							 channel_width = 25,
							 depositional_width = 25*60,
							 washload_fraction = 1.5,
							 sinuosity = 1.5,
							 verbose = False,
							 save_output = False)
        f.run()
        self.assertEqual(f.water_surface_elevation[0], 103.00768207498517)
        self.assertEqual(f.water_surface_elevation[-2], 1.100884370946358)
예제 #3
0
    def test_eta(self):
    
        f = AgDegNormalSub(friction = 'Manning',
							 Qf = 70,
							 If = 0.03,
							 L = 10000,
							 dx = 100,
							 S = 0.008,
							 Gtf = 700000,
							 D = 30,
							 porosity = 0.35,
							 kc = 75,
							 phi = 1,
							 alpha_r = 8.1,
							 alpha_t = 8.,
							 tau_c = 0.047,
							 nt = 1.5,
							 Cz = 22,
							 alpha_u = 0.5,
							 dt = 0.01,
							 totalT = 120,
							 print_dt = 120,
							 g = 9.81,
							 rho_w = 1000,
							 rho_s = 2650,
							 subsidence = 5,
							 channel_width = 25,
							 depositional_width = 25*60,
							 washload_fraction = 1.5,
							 sinuosity = 1.5,
							 verbose = False,
							 save_output = False)
                 
        f.run()
        self.assertEqual(f.bed_elevation[0], 102.26911259423767)
        self.assertEqual(f.bed_elevation[-2], 0.23724567562372906)