Python topo2writerの例

コード例 #1

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)

コード例 #2

    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)

コード例 #3

ファイル: maketopo.py プロジェクト: dlgeorge/flume

        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'

コード例 #4

ファイル: maketopo.py プロジェクト: dlgeorge/flume

    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)

コード例 #5

ファイル: maketopo.py プロジェクト: dlgeorge/flume

    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)

コード例 #6

ファイル: maketopo.py プロジェクト: dlgeorge/flume

        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'

コード例 #7

    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)

コード例 #8

ファイル: maketopo.py プロジェクト: dlgeorge/flume

    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)

コード例 #9

        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'

コード例 #10

ファイル: makeaux.py プロジェクト: dlgeorge/flume

    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)

コード例 #11

    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)

コード例 #12

ファイル: maketopo.py プロジェクト: dlgeorge/flume

    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)

コード例 #13

ファイル: maketopo.py プロジェクト: dlgeorge/flume

    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)