def write_data_objects(self):
super(PlaneWaveTest, self).write_data_objects()
# Write out bathy file
mx = self.rundata.clawdata.num_cells[0]
my = self.rundata.clawdata.num_cells[1]
xlower = self.rundata.clawdata.lower[0]
xupper = self.rundata.clawdata.upper[0]
ylower = self.rundata.clawdata.lower[1]
yupper = self.rundata.clawdata.upper[1]
dx = (xupper - xlower) / mx
dy = (yupper - ylower) / my
d = min(dx,dy)
mx = int((xupper - xlower) / d) + 8
my = int((yupper - ylower) / d) + 8
xlower = xlower - d*4.0
ylower = ylower - d*4.0
xupper = xupper + d*4.0
yupper = yupper + d*4.0
step = lambda x,y: setrun.bathy_step(x, y, location=self.bathy_location,
angle=self.bathy_angle)
tt.topo2writer('./topo.tt2', step, xlower, xupper, ylower, yupper,
mx, my, nodata_value=-99999)
# Write out simple bathy geometry file for communication to the plotting
with open("./bathy_geometry.data", 'w') as bathy_geometry_file:
bathy_geometry_file.write("%s\n%s" % (self.bathy_location,
self.bathy_angle) )
def write_data_objects(self):
super(PlaneWaveTest, self).write_data_objects()
# Write out bathy file
mx = self.rundata.clawdata.num_cells[0]
my = self.rundata.clawdata.num_cells[1]
xlower = self.rundata.clawdata.lower[0]
xupper = self.rundata.clawdata.upper[0]
ylower = self.rundata.clawdata.lower[1]
yupper = self.rundata.clawdata.upper[1]
dx = (xupper - xlower) / mx
dy = (yupper - ylower) / my
d = min(dx,dy)
mx = int((xupper - xlower) / d) + 8
my = int((yupper - ylower) / d) + 8
xlower = xlower - d*4.0
ylower = ylower - d*4.0
xupper = xupper + d*4.0
yupper = yupper + d*4.0
step = lambda x,y: setrun.bathy_step(x, y, location=self.bathy_location,
angle=self.bathy_angle)
tt.topo2writer('./topo.tt2', step, xlower, xupper, ylower, yupper,
mx, my, nodata_value=-99999)
# Write out simple bathy geometry file for communication to the plotting
with open("./bathy_geometry.data", 'w') as bathy_geometry_file:
bathy_geometry_file.write("%s\n%s" % (self.bathy_location,
self.bathy_angle) )
def write_topo_file(run_data, out_file, **kwargs):
# Write out bathy file
mx = run_data.clawdata.num_cells[0]
my = run_data.clawdata.num_cells[1]
xlower = run_data.clawdata.lower[0]
xupper = run_data.clawdata.upper[0]
ylower = run_data.clawdata.lower[1]
yupper = run_data.clawdata.upper[1]
dx = (xupper - xlower) / mx
dy = (yupper - ylower) / my
d = min(dx,dy)
mx = int((xupper - xlower) / d) + 8
my = int((yupper - ylower) / d) + 8
xlower = xlower - d*4.0
ylower = ylower - d*4.0
xupper = xupper + d*4.0
yupper = yupper + d*4.0
step = lambda x,y: bathy_step(x, y, **kwargs)
tt.topo2writer(out_file, step, xlower, xupper, ylower, yupper, mx, my, nodata_value=-99999)
# Write out simple bathy geometry file for communication to the plotting
with open("./bathy_geometry.data", 'w') as bathy_geometry_file:
if kwargs.has_key("location"):
location = kwargs['location']
else:
location = 0.15
if kwargs.has_key("angle"):
angle = kwargs['angle']
else:
angle = 0.0
bathy_geometry_file.write("%s\n%s" % (location, angle) )
def maketopo1b(path):
"""
Output topography file for the entire domain
"""
nxpoints = 101
nypoints = 71
xlower = 0.0
xupper = 1.0
ylower = 0.5
yupper = 1.2
topotools.topo2writer(path, topo1, xlower, xupper, ylower, yupper, nxpoints, nypoints)
def maketopo_flat():
"""
Output topography file for the entire domain
"""
nxpoints = 201
nypoints = 301
xlower = 204.812
xupper = 205.012
ylower = 19.7
yupper = 20.0
outfile= "flat.tt2"
topotools.topo2writer(outfile,topo_flat,xlower,xupper,ylower,yupper,nxpoints,nypoints)
def maketopo1a(path):
"""
Output topography file for the entire domain
"""
nxpoints = 101
nypoints = 76
xlower = -1.5
xupper = 2.5
ylower = -1.0
yupper = 2.0
topotools.topo2writer(path, topo1, xlower, xupper, ylower, yupper, nxpoints, nypoints)
def maketopo1b(path):
"""
Output topography file for the entire domain
"""
nxpoints = 101
nypoints = 71
xlower = 0.0
xupper = 1.0
ylower = 0.5
yupper = 1.2
topotools.topo2writer(path,topo1,xlower,xupper,ylower,yupper,\
nxpoints,nypoints)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints = 1001
nypoints = 21
xlower = -50.e0
xupper = 50.e0
yupper = 1.e0
ylower = -1.e0
outfile= "tank.tt2"
topotools.topo2writer(outfile,topo,xlower,xupper,ylower,yupper,nxpoints,nypoints)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints = 2001
nypoints = 2001
xlower = 0
xupper = 200000.e0
yupper = 200000.e0
ylower = 0
outfile= "bowl.topotype2"
topotools.topo2writer(outfile,topo,xlower,xupper,ylower,yupper,nxpoints,nypoints)
def maketopo_flat():
"""
Output topography file for the entire domain
"""
nxpoints = 201
nypoints = 301
xlower = 204.812
xupper = 205.012
ylower = 19.7
yupper = 20.0
outfile = "flat.tt2"
topotools.topo2writer(outfile, topo_flat, xlower, xupper, ylower, yupper,
nxpoints, nypoints)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints = 101
nypoints = 101
xlower = -100.e0
xupper = 100.e0
yupper = 100.e0
ylower = -100.e0
outfile = "bowl.topotype2"
topotools.topo2writer(outfile, topo, xlower, xupper, ylower, yupper,
nxpoints, nypoints)
def maketopo1a(path):
"""
Output topography file for the entire domain
"""
nxpoints = 101
nypoints = 76
xlower = -1.5
xupper = 2.5
ylower = -1.
yupper = 2.
topotools.topo2writer(path,topo1,xlower,xupper,ylower,yupper,\
nxpoints,nypoints)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints = 1001
nypoints = 21
xlower = -50.e0
xupper = 50.e0
yupper = 1.e0
ylower = -1.e0
outfile = "tank.tt2"
topotools.topo2writer(outfile, topo, xlower, xupper, ylower, yupper,
nxpoints, nypoints)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints=201
nypoints=201
xlower = -10e0
xupper= 10e0
ylower = -10e0
yupper= 10e0
outfile= "topo1.topotype2"
topo2writer(outfile,topo,xlower,xupper,ylower,yupper,nxpoints,nypoints)
dx = (xupper-xlower)/(nxpoints-1)
dy = (yupper-ylower)/(nypoints-1)
print "==> topo1 has dx = %g, dy = %g" % (dx,dy)
def maketopo2():
"""
Another topo file for testing
"""
nxpoints=21
nypoints=51
xlower = -0.5
xupper= -0.3
ylower = -0.1
yupper= 0.4
outfile= "topo2.topotype2"
topo2writer(outfile,topo2,xlower,xupper,ylower,yupper,nxpoints,nypoints)
dx = (xupper-xlower)/(nxpoints-1)
dy = (yupper-ylower)/(nypoints-1)
print "==> topo2 has dx = %g, dy = %g" % (dx,dy)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints=200
nypoints=200
xupper=2.e0
yupper=2.e0
xlower = -2.e0
ylower = -2.e0
outfile= "bowl.topotype2"
topo2writer(outfile,topo,xlower,xupper,ylower,yupper,nxpoints,nypoints)
# convert to netcdf form:
topo4 = topotools.Topography()
topo4.read(outfile,2)
topo4.write("bowl.nc",topo_type=4)
def maketopo():
"""
Output topography file for the entire domain
"""
nxpoints=200
nypoints=200
xupper=2.e0
yupper=2.e0
xlower = -2.e0
ylower = -2.e0
outfile= "bowl.topotype2"
topo2writer(outfile,topo,xlower,xupper,ylower,yupper,nxpoints,nypoints)
# convert to netcdf form:
topo4 = topotools.Topography()
topo4.read(outfile,2)
topo4.write("bowl.nc",topo_type=4)
def write_topo_file(run_data, out_file, **kwargs):
# Write out bathy file
mx = run_data.clawdata.num_cells[0]
my = run_data.clawdata.num_cells[1]
xlower = run_data.clawdata.lower[0]
xupper = run_data.clawdata.upper[0]
ylower = run_data.clawdata.lower[1]
yupper = run_data.clawdata.upper[1]
dx = (xupper - xlower) / mx
dy = (yupper - ylower) / my
d = min(dx, dy)
mx = int((xupper - xlower) / d) + 8
my = int((yupper - ylower) / d) + 8
xlower = xlower - d * 4.0
ylower = ylower - d * 4.0
xupper = xupper + d * 4.0
yupper = yupper + d * 4.0
step = lambda x, y: bathy_step(x, y, **kwargs)
tt.topo2writer(out_file,
step,
xlower,
xupper,
ylower,
yupper,
mx,
my,
nodata_value=-99999)
# Write out simple bathy geometry file for communication to the plotting
with open("./bathy_geometry.data", 'w') as bathy_geometry_file:
if kwargs.has_key("location"):
location = kwargs['location']
else:
location = 0.15
if kwargs.has_key("angle"):
angle = kwargs['angle']
else:
angle = 0.0
bathy_geometry_file.write("%s\n%s" % (location, angle))
def maketopo():
"""
Output topography file for the entire domain
and near-shore in one location.
"""
nxpoints = 400
nypoints = 400
xlower = -30.e0
xupper = 30.e0
ylower = 20.e0
yupper = 60.e0
outfile = "ocean.topotype2"
topotools.topo2writer(outfile, topo, xlower, xupper, ylower, yupper,
nxpoints, nypoints)
nxpoints = 200
nypoints = 200
xlower = xisland - 1.
xupper = xisland + 1.
ylower = yisland - 1.
yupper = yisland + 1.
outfile = "island.topotype2"
topotools.topo2writer(outfile, topo, xlower, xupper, ylower, yupper,
nxpoints, nypoints)
