def test_eta(self):
    
        f = AgDegNormalFault(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,
							 fault_position = 0.5,
							 fault_height = 15,
							 fault_time = 1,
							 verbose = False,
							 save_output = False)
                 
        f.run()
        self.assertEqual(f.bed_elevation[0], 126.64142774887191)
        self.assertEqual(f.bed_elevation[-2], -13.585970469445467)
    def test_shear_stress(self):
    
        f = AgDegNormalFault(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,
							 fault_position = 0.5,
							 fault_height = 15,
							 fault_time = 1,
							 verbose = False,
							 save_output = False)
        f.run()
        self.assertEqual(f.bed_shear_stress[0], 0.21151184144832058)
        self.assertEqual(f.bed_shear_stress[-2], 0.21082629624709234)
    def test_water_surface(self):
    
        f = AgDegNormalFault(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,
							 fault_position = 0.5,
							 fault_height = 15,
							 fault_time = 1,
							 verbose = False,
							 save_output = False)
        f.run()
        self.assertEqual(f.water_surface_elevation[0], 127.37842780835773)
        self.assertEqual(f.water_surface_elevation[-2], -12.847944287798242)