import numpy as np import pylab import calc_prism_class as cpc p_dict = {} p_dict['g'] = 9.81 p_dict['b'] = 6.1 p_dict['Q'] = 17.0 p_dict['n'] = 0.014 p_dict['s'] = 0.001 p_dict['kn'] = 1.0 p_dict['Y'] = np.arange(0.001,10.0,0.001) rect = cpc.rect(p_dict) #--normal depth y0 = rect.y0() #--critical depth yc = rect.yc() print 'y0,yc',y0,yc #--M3 curve from gate to HJ start = 0.47 end = rect.yc()*1.001 Y_m3 = np.linspace(start,end,200) m3_results = rect.direct_step(Y_m3)
import math import numpy as np import pylab import calc_prism_class as cpc p_dict = {} p_dict['g'] = 9.81 p_dict['b'] = 6.1 p_dict['Q'] = 17.0 p_dict['n'] = 0.014 p_dict['s'] = 0.001 p_dict['kn'] = 1.0 p_dict['Y'] = np.arange(0.001, 10.0, 0.001) rect = cpc.rect(p_dict) #--normal depth y0 = rect.y0() #--critical depth yc = rect.yc() print 'y0,yc', y0, yc #--M3 curve from gate to HJ start = 0.47 end = rect.yc() * 1.001 Y_m3 = np.linspace(start, end, 200) m3_results = rect.direct_step(Y_m3) #for x,y in zip(m3_results[0],Y_m3):
import copy import numpy as np import pylab import calc_prism_class as cpc p_dict = {} p_dict['g'] = 9.81 p_dict['b'] = 10.0 p_dict['Q'] = 101.0 p_dict['n'] = 0.03 p_dict['s'] = 0.005 p_dict['kn'] = 1.0 p_dict['Y'] = np.arange(0.001, 10.0, 0.001) rect = cpc.rect(p_dict) #--normal depth y0 = rect.y0() #--critical depth yc = rect.yc() print 'y0,yc', y0, yc #--steep s_dict = copy.deepcopy(p_dict) s_dict['s'] = 0.02 rect_s = cpc.rect(s_dict) #--normal depth y0_s = rect_s.y0()
import numpy as np import pylab import calc_prism_class as cpc p_dict = {} p_dict['g'] = 9.81 p_dict['b'] = 10.0 p_dict['Q'] = 101.0 p_dict['n'] = 0.03 p_dict['s'] = 0.005 p_dict['kn'] = 1.0 p_dict['Y'] = np.arange(0.001,10.0,0.001) rect = cpc.rect(p_dict) #--normal depth y0 = rect.y0() #--critical depth yc = rect.yc() print 'y0,yc',y0,yc #--steep s_dict = copy.deepcopy(p_dict) s_dict['s'] = 0.02 rect_s = cpc.rect(s_dict) #--normal depth y0_s = rect_s.y0()