def make_data(): #phi file outfile= 'FlumePhi.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 140.0 ylower = -15.0 yupper = 15.0 nxpoints = int((xupper-xlower)/1.0) + 1 nypoints = int((yupper-ylower)/1.0) + 1 gt.topo2writer(outfile,flume_phi,xlower,xupper,ylower,yupper,nxpoints,nypoints) #phi file outfile= 'FlumeThetaRamp.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 140.0 ylower = -15.0 yupper = 15.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,flume_theta_ramp,xlower,xupper,ylower,yupper,nxpoints,nypoints) #phi file outfile= 'FlumeTheta.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 140.0 ylower = -15.0 yupper = 15.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,flume_theta,xlower,xupper,ylower,yupper,nxpoints,nypoints)
Z[np.ix_(yind,x3ind)] = depth Z[np.ix_(yind,xposind)] = (xpos -X[np.ix_(yind,xposind)])*np.tan(alpha3) return Z #flat topo outfile= 'ZeroTopo.tt2' outfile = os.path.join('topo',outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #wall topo outfile= 'Wall1Topo.tt2' outfile = os.path.join('topo',outfile) xlower = -15.0 xupper = 90.0 ylower = -0.5 yupper = 0.0 nxpoints = int((xupper-xlower)/0.05) + 1 nypoints = int((yupper-ylower)/0.05) + 1 gt.topo2writer(outfile,wallzero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #wall topo outfile= 'Wall2Topo.tt2' outfile = os.path.join('topo',outfile)
yind, x0ind)]) / flume_eta_res(X[np.ix_(yind, x0ind)], Y[np.ix_( yind, x0ind)]) return Z #flat topo outfile = 'ZeroTopo.tt2' outfile = os.path.join('topo', outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper - xlower) / 0.1) + 1 nypoints = int((yupper - ylower) / 0.1) + 1 gt.topo2writer(outfile, zero, xlower, xupper, ylower, yupper, nxpoints, nypoints) #flat topo outfile = 'ZeroTopoGate.tt2' outfile = os.path.join('topo', outfile) xlower = -8.0 xupper = 10.0 ylower = 0.0 yupper = 2.0 nxpoints = int((xupper - xlower) / 0.005) + 1 nypoints = int((yupper - ylower) / 0.005) + 1 gt.topo2writer(outfile, zero, xlower, xupper, ylower, yupper, nxpoints, nypoints) #wall topo outfile = 'Wall1Topo.tt2'
Z[np.ix_(yind,x1ind)] = m0 Z[np.ix_(yind,x0ind)] = m0*flume_eta(X[np.ix_(yind,x0ind)],Y[np.ix_(yind,x0ind)])/flume_eta_res(X[np.ix_(yind,x0ind)],Y[np.ix_(yind,x0ind)]) return Z #flat topo outfile= 'ZeroTopo.tt2' outfile = os.path.join('topo',outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #flat topo outfile= 'ZeroTopoGate.tt2' outfile = os.path.join('topo',outfile) xlower = -8.0 xupper = 10.0 ylower = 0.0 yupper = 2.0 nxpoints = int((xupper-xlower)/0.005) + 1 nypoints = int((yupper-ylower)/0.005) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #wall topo outfile= 'Wall1Topo.tt2' outfile = os.path.join('topo',outfile)
Z[np.ix_(yind,x1ind)] = m0 Z[np.ix_(yind,x0ind)] = m0*flume_eta(X[np.ix_(yind,x0ind)],Y[np.ix_(yind,x0ind)])/flume_eta_res(X[np.ix_(yind,x0ind)],Y[np.ix_(yind,x0ind)]) return Z #flat topo outfile= 'ZeroTopo.tt2' outfile = os.path.join('topo',outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #flat topo outfile= 'ZeroTopoGate.tt2' outfile = os.path.join('topo',outfile) xlower = -8.0 xupper = 10.0 ylower = 0.0 yupper = 2.0 nxpoints = int((xupper-xlower)/0.005) + 1 nypoints = int((yupper-ylower)/0.005) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #wall topo outfile= 'Wall1Topo.tt2' outfile = os.path.join('topo',outfile)
yind, x0ind)]) / flume_eta_res(X[np.ix_(yind, x0ind)], Y[np.ix_( yind, x0ind)]) return Z #flat topo outfile = 'ZeroTopo.tt2' outfile = os.path.join('topo', outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper - xlower) / 0.1) + 1 nypoints = int((yupper - ylower) / 0.1) + 1 gt.topo2writer(outfile, zero, xlower, xupper, ylower, yupper, nxpoints, nypoints) #flat topo outfile = 'ZeroTopoGate.tt2' outfile = os.path.join('topo', outfile) xlower = -8.0 xupper = 10.0 ylower = 0.0 yupper = 2.0 nxpoints = int((xupper - xlower) / 0.005) + 1 nypoints = int((yupper - ylower) / 0.005) + 1 gt.topo2writer(outfile, zero, xlower, xupper, ylower, yupper, nxpoints, nypoints) #wall topo outfile = 'Wall1Topo.tt2'
Z[np.ix_(yind,x2ind)] = etax1 - (x1-X[np.ix_(yind,x2ind)])*np.sin(backangle) return Z #pad outfile = 'FlumeRunoutPad_10cm.tt2' outfile = os.path.join('topo',outfile) xlower = 65.0 xupper = 180.00 ylower = -10.0 yupper = 12.0 nxpoints = int((xupper-xlower)/0.05) + 1 nypoints = int((yupper-ylower)/0.05) + 1 gt.topo2writer(outfile,bed_curve,xlower,xupper,ylower,yupper,nxpoints,nypoints) #hill outfile = 'FlumeHillside_1m.tt2' outfile = os.path.join('topo',outfile) xlower = -10.0 xupper = 72.0 ylower = -10.0 yupper = 12.0 nxpoints = int((xupper-xlower)/1.0) + 1 nypoints = int((yupper-ylower)/1.0) + 1 gt.topo2writer(outfile,hill,xlower,xupper,ylower,yupper,nxpoints,nypoints) #bed outfile = 'FlumeHopper_-10.0m_4.4m_1cm.tt2' outfile = os.path.join('topo',outfile)
Z[np.ix_(yind, x5ind)] = 0.0 pdb.set_trace() return Z outfile = 'Rotating_Wall_10.tt2' outfile = os.path.join('topo', outfile) xlower = 90.0 xupper = 120.0 ylower = -4.0 yupper = 6.0 nxpoints = int((xupper - xlower) / 0.05) + 1 nypoints = int((yupper - ylower) / 0.05) + 1 gt.topo2writer(outfile, adverse_ramp_prog, xlower, xupper, ylower, yupper, nxpoints, nypoints) #flat topo outfile = 'ZeroTopo.tt2' outfile = os.path.join('topo', outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper - xlower) / 0.1) + 1 nypoints = int((yupper - ylower) / 0.1) + 1 #gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #flat topo outfile = 'ZeroTopoGate.tt2' outfile = os.path.join('topo', outfile)
yind, x0ind)]) / flume_eta_res(X[np.ix_(yind, x0ind)], Y[np.ix_( yind, x0ind)]) return Z #flat topo outfile = 'ZeroTopo.tt2' outfile = os.path.join('topo', outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper - xlower) / 0.1) + 1 nypoints = int((yupper - ylower) / 0.1) + 1 gt.topo2writer(outfile, zero, xlower, xupper, ylower, yupper, nxpoints, nypoints) #flat topo outfile = 'ZeroTopoGate.tt2' outfile = os.path.join('topo', outfile) xlower = -8.0 xupper = 10.0 ylower = 0.0 yupper = 2.0 nxpoints = int((xupper - xlower) / 0.005) + 1 nypoints = int((yupper - ylower) / 0.005) + 1 gt.topo2writer(outfile, zero, xlower, xupper, ylower, yupper, nxpoints, nypoints) #wall topo outfile = 'Wall1Topo.tt2'
Z = deg2rad*Z return Z #phi file outfile= 'FlumePhi.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 180.0 ylower = -15.0 yupper = 15.0 nxpoints = int((xupper-xlower)/1.0) + 1 nypoints = int((yupper-ylower)/1.0) + 1 gt.topo2writer(outfile,phi,xlower,xupper,ylower,yupper,nxpoints,nypoints) #phi file outfile= 'FlumeTheta.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 180.0 ylower = -15.0 yupper = 15.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,flume_theta,xlower,xupper,ylower,yupper,nxpoints,nypoints)
Z = deg2rad*Z return Z #phi file outfile= 'FlumePhi.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 160.0 ylower = -4.0 yupper = 6.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,phi,xlower,xupper,ylower,yupper,nxpoints,nypoints) #theta file outfile= 'FlumeTheta.tt2' outfile = os.path.join('aux',outfile) xlower = -15.0 xupper = 160.0 ylower = -4.0 yupper = 6.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,flume_theta,xlower,xupper,ylower,yupper,nxpoints,nypoints)
Z[np.ix_(yind,x1ind)] = m0 Z[np.ix_(yind,x0ind)] = m0*flume_eta(X[np.ix_(yind,x0ind)],Y[np.ix_(yind,x0ind)])/flume_eta_res(X[np.ix_(yind,x0ind)],Y[np.ix_(yind,x0ind)]) return Z #flat topo outfile= 'ZeroTopo.tt2' outfile = os.path.join('topo',outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #flat topo outfile= 'ZeroTopoGate.tt2' outfile = os.path.join('topo',outfile) xlower = -8.0 xupper = 10.0 ylower = 0.0 yupper = 2.0 nxpoints = int((xupper-xlower)/0.005) + 1 nypoints = int((yupper-ylower)/0.005) + 1 gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #wall topo outfile= 'Wall1Topo.tt2' outfile = os.path.join('topo',outfile)
Z[np.ix_(yind,x4ind)] = z5 + np.tan(thetaramp*deg2rad)*(xend-x5) Z[np.ix_(yind,x5ind)] = 0.0 pdb.set_trace() return Z outfile= 'Rotating_Wall_10.tt2' outfile = os.path.join('topo',outfile) xlower = 90.0 xupper = 120.0 ylower = -4.0 yupper = 6.0 nxpoints = int((xupper-xlower)/0.05) + 1 nypoints = int((yupper-ylower)/0.05) + 1 gt.topo2writer(outfile,adverse_ramp_prog,xlower,xupper,ylower,yupper,nxpoints,nypoints) #flat topo outfile= 'ZeroTopo.tt2' outfile = os.path.join('topo',outfile) xlower = -10.0 xupper = 180 ylower = -10.0 yupper = 10.0 nxpoints = int((xupper-xlower)/0.1) + 1 nypoints = int((yupper-ylower)/0.1) + 1 #gt.topo2writer(outfile,zero,xlower,xupper,ylower,yupper,nxpoints,nypoints) #flat topo outfile= 'ZeroTopoGate.tt2' outfile = os.path.join('topo',outfile)