def make_fgmax(FG):
print "---------------------------------------------- "
x1,x2 = FG.x1, FG.x2
y1,y2 = FG.y1, FG.y2
point_style = FG.point_style
if point_style not in [1,2]:
raise NotImplementedError("make_fgmax not implemented for point_style %i" \
% point_style)
if point_style == 2:
# 2d grid of points
if FG.nx is None:
dx = FG.dx
nx = int(round((x2-x1)/dx)) + 1
if abs((nx-1)*dx + x1 - x2) > 1e-6:
print "Warning: abs((nx-1)*dx + x1 - x2) = ", \
abs((nx-1)*dx + x1 - x2)
print " old x2: %22.16e" % x2
x2 = x1 + dx*(nx-1)
print " resetting x2 to %22.16e" % x2
else:
nx = FG.nx
dx = (x2-x1)/(nx+1.)
if FG.dx is not None:
print "*** Warning: dx specified over-ridden by: ",dx
if FG.ny is None:
dy = FG.dy
if dy is None:
dy = dx
ny = int(round((y2-y1)/dy)) + 1
if abs((ny-1)*dy + y1 - y2) > 1e-6:
print "Warning: abs((ny-1)*dy + y1 - y2) = ", \
abs((ny-1)*dy + y1 - y2)
print " old y2: %22.16e" % y2
y2 = y1 + dy*(ny-1)
print " resetting y2 to %22.16e" % y2
else:
ny = FG.ny
dy = (y2-y1)/(ny+1.)
if FG.dy is not None:
print "*** Warning: dy specified over-ridden by: ",dy
npts = nx*ny
fid = open(FG.fname,'w')
fid.write("%16.10e # tstart_max\n" % FG.tstart_max)
fid.write("%16.10e # tend_max\n" % FG.tend_max)
fid.write("%16.10e # dt_check\n" % FG.dt_check)
fid.write("%i %s # min_level_check\n" \
% (FG.min_level_check,16*" "))
fid.write("%16.10e # arrival_tol\n" % FG.arrival_tol)
fid.write("%i %s # point_style\n" \
% (FG.point_style,16*" "))
fid.write("%i %i %s # nx,ny\n" \
% (nx,ny,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1,y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2,y2))
fid.close()
print "Created file ", FG.fname
print " specifying fixed grid with shape %i by %i, with %i points" \
% (nx,ny,npts)
print " lower left = (%15.10f,%15.10f)" % (x1,y1)
print " upper right = (%15.10f,%15.10f)" % (x2,y2)
print " dx = %15.10e, dy = %15.10e" % (dx,dy)
xy = [x1,x2,y1,y2]
fname_root = os.path.splitext(FG.fname)[0]
kml_file = fname_root + '.kml'
kmltools.box2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style==1:
# 1d transect of points
if FG.npts is None:
dx = FG.dx
npts = int(round(sqrt((x2-x1)**2 + (y2-y1)**2)/dx)) + 1
if abs((npts-1)*dx + x1 - x2) > 1e-6:
print "Warning: abs((npts-1)*dx + x1 - x2) = ", \
abs((npts-1)*dx + x1 - x2)
x2 = x1 + dx*(npts-1)
y2 = y1 + dx*(npts-1)
print " resetting x2 to %g" % x2
print " resetting y2 to %g" % y2
else:
npts = FG.npts
dx = sqrt((x2-x1)**2 + (y2-y1)**2)/(npts+1.)
if FG.dx is not None:
print "*** Warning: dx specified over-ridden by: ",dx
print "Creating 1d fixed grid with %s points" % npts
print " dx = %g" % dx
fid = open(FG.fname,'w')
fid.write("%g # tstart_max\n" % FG.tstart_max)
fid.write("%g # tend_max\n" % FG.tend_max)
fid.write("%g # dt_check\n" % FG.dt_check)
fid.write("%i # min_level_check\n" % FG.min_level_check)
fid.write("%g # arrival_tol\n" % FG.arrival_tol)
fid.write("%g # point_style\n" % FG.point_style)
fid.write("%i # npts\n" % (npts))
fid.write("%g %g # x1, y1\n" % (x1,y1))
fid.write("%g %g # x2, y2\n" % (x2,y2))
fid.close()
print "Created file ", FG.fname
print " specifying fixed grid with %i points equally spaced from " \
% npts
print " (%g,%g) to (%g,%g)" % (x1,y1,x2,y2)
"""
Code to make the file CCregion.kml. You can only run this if you have Clawpack installed.
The CCregion.kml file can be opened in Google Earth and should look like the image in
CCregion.png.
Note that longitude must be shifted by 360 degrees to work on GE.
"""
from clawpack.geoclaw import kmltools
kmltools.box2kml((235.79781-360, 235.82087-360, 41.739671,41.762726), fname='CCregion.kml',
name='Crescent City subset', color='FFFFFF')
"""
Code to make the file CCregion.kml. You can only run this if you have Clawpack installed.
The CCregion.kml file can be opened in Google Earth and should look like the image in
CCregion.png.
Note that longitude must be shifted by 360 degrees to work on GE.
"""
from clawpack.geoclaw import kmltools
kmltools.box2kml((235.79781 - 360, 235.82087 - 360, 41.739671, 41.762726),
fname='CCregion.kml',
name='Crescent City subset',
color='FFFFFF')
def make_fgmax(FG):
"""Deprecated -- use FGmaxGrid.write_input_data"""
print "---------------------------------------------- "
x1, x2 = FG.x1, FG.x2
y1, y2 = FG.y1, FG.y2
point_style = FG.point_style
if point_style not in [1, 2, 3]:
raise NotImplementedError("make_fgmax not implemented for point_style %i" \
% point_style)
if point_style == 2:
# 2d grid of points
if FG.nx is None:
dx = FG.dx
nx = int(round((x2 - x1) / dx)) + 1
if abs((nx - 1) * dx + x1 - x2) > 1e-6:
print "Warning: abs((nx-1)*dx + x1 - x2) = ", \
abs((nx-1)*dx + x1 - x2)
print " old x2: %22.16e" % x2
x2 = x1 + dx * (nx - 1)
print " resetting x2 to %22.16e" % x2
else:
nx = FG.nx
dx = (x2 - x1) / (nx + 1.)
if FG.dx is not None:
print "*** Warning: dx specified over-ridden by: ", dx
if FG.ny is None:
dy = FG.dy
if dy is None:
dy = dx
ny = int(round((y2 - y1) / dy)) + 1
if abs((ny - 1) * dy + y1 - y2) > 1e-6:
print "Warning: abs((ny-1)*dy + y1 - y2) = ", \
abs((ny-1)*dy + y1 - y2)
print " old y2: %22.16e" % y2
y2 = y1 + dy * (ny - 1)
print " resetting y2 to %22.16e" % y2
else:
ny = FG.ny
dy = (y2 - y1) / (ny + 1.)
if FG.dy is not None:
print "*** Warning: dy specified over-ridden by: ", dy
npts = nx * ny
fid = open(FG.fname, 'w')
fid.write("%16.10e # tstart_max\n" % FG.tstart_max)
fid.write("%16.10e # tend_max\n" % FG.tend_max)
fid.write("%16.10e # dt_check\n" % FG.dt_check)
fid.write("%i %s # min_level_check\n" \
% (FG.min_level_check,16*" "))
fid.write("%16.10e # arrival_tol\n" % FG.arrival_tol)
fid.write("%i %s # point_style\n" \
% (FG.point_style,16*" "))
fid.write("%i %i %s # nx,ny\n" \
% (nx,ny,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1, y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2, y2))
fid.close()
print "Created file ", FG.fname
print " specifying fixed grid with shape %i by %i, with %i points" \
% (nx,ny,npts)
print " lower left = (%15.10f,%15.10f)" % (x1, y1)
print " upper right = (%15.10f,%15.10f)" % (x2, y2)
print " dx = %15.10e, dy = %15.10e" % (dx, dy)
xy = [x1, x2, y1, y2]
fname_root = os.path.splitext(FG.fname)[0]
kml_file = fname_root + '.kml'
kmltools.box2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style == 1:
# 1d transect of points
if FG.npts is None:
dx = FG.dx
npts = int(round(sqrt((x2 - x1)**2 + (y2 - y1)**2) / dx)) + 1
if abs((npts - 1) * dx + x1 - x2) > 1e-6:
print "Warning: abs((npts-1)*dx + x1 - x2) = ", \
abs((npts-1)*dx + x1 - x2)
x2 = x1 + dx * (npts - 1)
y2 = y1 + dx * (npts - 1)
print " resetting x2 to %g" % x2
print " resetting y2 to %g" % y2
else:
npts = FG.npts
dx = sqrt((x2 - x1)**2 + (y2 - y1)**2) / (npts + 1.)
if FG.dx is not None:
print "*** Warning: dx specified over-ridden by: ", dx
print "Creating 1d fixed grid with %s points" % npts
print " dx = %g" % dx
fid = open(FG.fname, 'w')
fid.write("%g # tstart_max\n" % FG.tstart_max)
fid.write("%g # tend_max\n" % FG.tend_max)
fid.write("%g # dt_check\n" % FG.dt_check)
fid.write("%i # min_level_check\n" %
FG.min_level_check)
fid.write("%g # arrival_tol\n" % FG.arrival_tol)
fid.write("%g # point_style\n" % FG.point_style)
fid.write("%i # npts\n" % (npts))
fid.write("%g %g # x1, y1\n" % (x1, y1))
fid.write("%g %g # x2, y2\n" % (x2, y2))
fid.close()
print "Created file ", FG.fname
print " specifying fixed grid with %i points equally spaced from " \
% npts
print " (%g,%g) to (%g,%g)" % (x1, y1, x2, y2)
# not yet implemented:
#fname_root = os.path.splitext(FG.fname)[0]
#kml_file = fname_root + '.kml'
#kmltools.line2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style == 3:
# arbitrary quadrilateral
x3, x4 = FG.x3, FG.x4
y3, y4 = FG.y3, FG.y4
if FG.n12 is None:
raise NotImplementedError("Need to set n12 and n23")
else:
npts = FG.n12 * FG.n23
fid = open(FG.fname, 'w')
fid.write("%g # tstart_max\n" % FG.tstart_max)
fid.write("%g # tend_max\n" % FG.tend_max)
fid.write("%g # dt_check\n" % FG.dt_check)
fid.write("%i # min_level_check\n" %
FG.min_level_check)
fid.write("%g # arrival_tol\n" % FG.arrival_tol)
fid.write("%g # point_style\n" % FG.point_style)
fid.write("%i %i %s # FG.n12,FG.n23\n" \
% (FG.n12,FG.n23,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1, y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2, y2))
fid.write("%16.10e %20.10e # x3, y3\n" % (x3, y3))
fid.write("%16.10e %20.10e # x4, y4\n" % (x4, y4))
fid.close()
print "Created file ", FG.fname
print " specifying fixed grid as a quadrilateral"
print " %i by %i, with %i points" \
% (FG.n12,FG.n23,npts)
print " corner 1 = (%15.10f,%15.10f)" % (x1, y1)
print " corner 2 = (%15.10f,%15.10f)" % (x2, y2)
print " corner 3 = (%15.10f,%15.10f)" % (x3, y3)
print " corner 4 = (%15.10f,%15.10f)" % (x4, y4)
xy = [x1, y1, x2, y2, x3, y3, x4, y4]
fname_root = os.path.splitext(FG.fname)[0]
kml_file = fname_root + '.kml'
kmltools.quad2kml(xy, kml_file, fname_root, color='8888FF')
def write_input_data(self, input_file_name=None):
r"""
Write input data from a file like *fgmax.txt* that
will be read in by GeoClaw as input data.
"""
if input_file_name is not None:
self.input_file_name = input_file_name
print "---------------------------------------------- "
point_style = self.point_style
if point_style not in [0, 1, 2, 3]:
raise NotImplementedError("make_fgmax not implemented for point_style %i" \
% point_style)
# write header, independent of point_style:
fid = open(self.input_file_name, 'w')
fid.write("%16.10e # tstart_max\n" % self.tstart_max)
fid.write("%16.10e # tend_max\n" % self.tend_max)
fid.write("%16.10e # dt_check\n" % self.dt_check)
fid.write("%i %s # min_level_check\n" \
% (self.min_level_check,16*" "))
fid.write("%16.10e # arrival_tol\n" % self.arrival_tol)
fid.write("%i %s # point_style\n" \
% (self.point_style,16*" "))
if point_style == 0:
# list of points
npts = self.npts
print "Creating unstructured grid of %s points" % npts
fid.write("%i # npts\n" % (npts))
for k in range(npts):
fid.write("%g %g \n" % (self.X[k], self.Y[k]))
fid.close()
print "Created file ", self.input_file_name
elif point_style == 1:
# 1d transect of points
x1, x2 = self.x1, self.x2
y1, y2 = self.y1, self.y2
if self.npts is None:
dx = self.dx
npts = int(round(sqrt((x2 - x1)**2 + (y2 - y1)**2) / dx)) + 1
if abs((npts - 1) * dx + x1 - x2) > 1e-6:
print "Warning: abs((npts-1)*dx + x1 - x2) = ", \
abs((npts-1)*dx + x1 - x2)
x2 = x1 + dx * (npts - 1)
y2 = y1 + dx * (npts - 1)
print " resetting x2 to %g" % x2
print " resetting y2 to %g" % y2
else:
npts = self.npts
dx = sqrt((x2 - x1)**2 + (y2 - y1)**2) / (npts + 1.)
if self.dx is not None:
print "*** Warning: dx specified over-ridden by: ", dx
print "Creating 1d fixed grid with %s points" % npts
print " dx = %g" % dx
fid.write("%i # npts\n" % (npts))
fid.write("%g %g # x1, y1\n" % (x1, y1))
fid.write("%g %g # x2, y2\n" % (x2, y2))
fid.close()
print "Created file ", self.input_file_name
print " specifying fixed grid with %i points equally spaced from " \
% npts
print " (%g,%g) to (%g,%g)" % (x1, y1, x2, y2)
# not yet implemented:
#fname_root = os.path.splitext(self.input_file_name)[0]
#kml_file = fname_root + '.kml'
#kmltools.line2kml(xy, kml_file, fname_root, color='8888FF')
if point_style == 2:
# 2d grid of points
x1, x2 = self.x1, self.x2
y1, y2 = self.y1, self.y2
if self.nx is None:
dx = self.dx
nx = int(round((x2 - x1) / dx)) + 1
if abs((nx - 1) * dx + x1 - x2) > 1e-6:
print "Warning: abs((nx-1)*dx + x1 - x2) = ", \
abs((nx-1)*dx + x1 - x2)
print " old x2: %22.16e" % x2
x2 = x1 + dx * (nx - 1)
print " resetting x2 to %22.16e" % x2
else:
nx = self.nx
dx = (x2 - x1) / (nx + 1.)
if self.dx is not None:
print "*** Warning: dx specified over-ridden by: ", dx
if self.ny is None:
dy = self.dy
if dy is None:
dy = dx
ny = int(round((y2 - y1) / dy)) + 1
if abs((ny - 1) * dy + y1 - y2) > 1e-6:
print "Warning: abs((ny-1)*dy + y1 - y2) = ", \
abs((ny-1)*dy + y1 - y2)
print " old y2: %22.16e" % y2
y2 = y1 + dy * (ny - 1)
print " resetting y2 to %22.16e" % y2
else:
ny = self.ny
dy = (y2 - y1) / (ny + 1.)
if self.dy is not None:
print "*** Warning: dy specified over-ridden by: ", dy
npts = nx * ny
fid.write("%i %i %s # nx,ny\n" \
% (nx,ny,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1, y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2, y2))
fid.close()
print "Created file ", self.input_file_name
print " specifying fixed grid with shape %i by %i, with %i points" \
% (nx,ny,npts)
print " lower left = (%15.10f,%15.10f)" % (x1, y1)
print " upper right = (%15.10f,%15.10f)" % (x2, y2)
print " dx = %15.10e, dy = %15.10e" % (dx, dy)
xy = [x1, x2, y1, y2]
fname_root = os.path.splitext(self.input_file_name)[0]
kml_file = fname_root + '.kml'
kmltools.box2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style == 3:
# arbitrary quadrilateral
x1, x2 = self.x1, self.x2
y1, y2 = self.y1, self.y2
x3, x4 = self.x3, self.x4
y3, y4 = self.y3, self.y4
if self.n12 is None:
raise NotImplementedError("Need to set n12 and n23")
else:
npts = self.n12 * self.n23
fid.write("%i %i %s # self.n12,self.n23\n" \
% (self.n12,self.n23,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1, y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2, y2))
fid.write("%16.10e %20.10e # x3, y3\n" % (x3, y3))
fid.write("%16.10e %20.10e # x4, y4\n" % (x4, y4))
fid.close()
print "Created file ", self.input_file_name
print " specifying fixed grid as a quadrilateral"
print " %i by %i, with %i points" \
% (self.n12,self.n23,npts)
print " corner 1 = (%15.10f,%15.10f)" % (x1, y1)
print " corner 2 = (%15.10f,%15.10f)" % (x2, y2)
print " corner 3 = (%15.10f,%15.10f)" % (x3, y3)
print " corner 4 = (%15.10f,%15.10f)" % (x4, y4)
xy = [x1, y1, x2, y2, x3, y3, x4, y4]
fname_root = os.path.splitext(self.input_file_name)[0]
kml_file = fname_root + '.kml'
kmltools.quad2kml(xy, kml_file, fname_root, color='8888FF')
def make_fgmax(FG):
"""Deprecated -- use FGmaxGrid.write_input_data"""
print("---------------------------------------------- ")
x1,x2 = FG.x1, FG.x2
y1,y2 = FG.y1, FG.y2
point_style = FG.point_style
if point_style not in [1,2,3]:
raise NotImplementedError("make_fgmax not implemented for point_style %i" \
% point_style)
if point_style == 2:
# 2d grid of points
if FG.nx is None:
dx = FG.dx
nx = int(round((x2-x1)/dx)) + 1
if abs((nx-1)*dx + x1 - x2) > 1e-6:
print("Warning: abs((nx-1)*dx + x1 - x2) = ", \
abs((nx-1)*dx + x1 - x2))
print(" old x2: %22.16e" % x2)
x2 = x1 + dx*(nx-1)
print(" resetting x2 to %22.16e" % x2)
else:
nx = FG.nx
dx = (x2-x1)/(nx+1.)
if FG.dx is not None:
print("*** Warning: dx specified over-ridden by: ",dx)
if FG.ny is None:
dy = FG.dy
if dy is None:
dy = dx
ny = int(round((y2-y1)/dy)) + 1
if abs((ny-1)*dy + y1 - y2) > 1e-6:
print("Warning: abs((ny-1)*dy + y1 - y2) = ", \
abs((ny-1)*dy + y1 - y2))
print(" old y2: %22.16e" % y2)
y2 = y1 + dy*(ny-1)
print(" resetting y2 to %22.16e" % y2)
else:
ny = FG.ny
dy = (y2-y1)/(ny+1.)
if FG.dy is not None:
print("*** Warning: dy specified over-ridden by: ",dy)
npts = nx*ny
fid = open(FG.fname,'w')
fid.write("%16.10e # tstart_max\n" % FG.tstart_max)
fid.write("%16.10e # tend_max\n" % FG.tend_max)
fid.write("%16.10e # dt_check\n" % FG.dt_check)
fid.write("%i %s # min_level_check\n" \
% (FG.min_level_check,16*" "))
fid.write("%16.10e # arrival_tol\n" % FG.arrival_tol)
fid.write("%i %s # point_style\n" \
% (FG.point_style,16*" "))
fid.write("%i %i %s # nx,ny\n" \
% (nx,ny,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1,y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2,y2))
fid.close()
print("Created file ", FG.fname)
print(" specifying fixed grid with shape %i by %i, with %i points" \
% (nx,ny,npts))
print(" lower left = (%15.10f,%15.10f)" % (x1,y1))
print(" upper right = (%15.10f,%15.10f)" % (x2,y2))
print(" dx = %15.10e, dy = %15.10e" % (dx,dy))
xy = [x1,x2,y1,y2]
fname_root = os.path.splitext(FG.fname)[0]
kml_file = fname_root + '.kml'
kmltools.box2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style==1:
# 1d transect of points
if FG.npts is None:
dx = FG.dx
npts = int(round(sqrt((x2-x1)**2 + (y2-y1)**2)/dx)) + 1
if abs((npts-1)*dx + x1 - x2) > 1e-6:
print("Warning: abs((npts-1)*dx + x1 - x2) = ", \
abs((npts-1)*dx + x1 - x2))
x2 = x1 + dx*(npts-1)
y2 = y1 + dx*(npts-1)
print(" resetting x2 to %g" % x2)
print(" resetting y2 to %g" % y2)
else:
npts = FG.npts
dx = sqrt((x2-x1)**2 + (y2-y1)**2)/(npts+1.)
if FG.dx is not None:
print("*** Warning: dx specified over-ridden by: ",dx)
print("Creating 1d fixed grid with %s points" % npts)
print(" dx = %g" % dx)
fid = open(FG.fname,'w')
fid.write("%g # tstart_max\n" % FG.tstart_max)
fid.write("%g # tend_max\n" % FG.tend_max)
fid.write("%g # dt_check\n" % FG.dt_check)
fid.write("%i # min_level_check\n" % FG.min_level_check)
fid.write("%g # arrival_tol\n" % FG.arrival_tol)
fid.write("%g # point_style\n" % FG.point_style)
fid.write("%i # npts\n" % (npts))
fid.write("%g %g # x1, y1\n" % (x1,y1))
fid.write("%g %g # x2, y2\n" % (x2,y2))
fid.close()
print("Created file ", FG.fname)
print(" specifying fixed grid with %i points equally spaced from " \
% npts)
print(" (%g,%g) to (%g,%g)" % (x1,y1,x2,y2))
# not yet implemented:
#fname_root = os.path.splitext(FG.fname)[0]
#kml_file = fname_root + '.kml'
#kmltools.line2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style==3:
# arbitrary quadrilateral
x3,x4 = FG.x3, FG.x4
y3,y4 = FG.y3, FG.y4
if FG.n12 is None:
raise NotImplementedError("Need to set n12 and n23")
else:
npts = FG.n12 * FG.n23
fid = open(FG.fname,'w')
fid.write("%g # tstart_max\n" % FG.tstart_max)
fid.write("%g # tend_max\n" % FG.tend_max)
fid.write("%g # dt_check\n" % FG.dt_check)
fid.write("%i # min_level_check\n" % FG.min_level_check)
fid.write("%g # arrival_tol\n" % FG.arrival_tol)
fid.write("%g # point_style\n" % FG.point_style)
fid.write("%i %i %s # FG.n12,FG.n23\n" \
% (FG.n12,FG.n23,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1,y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2,y2))
fid.write("%16.10e %20.10e # x3, y3\n" % (x3,y3))
fid.write("%16.10e %20.10e # x4, y4\n" % (x4,y4))
fid.close()
print("Created file ", FG.fname)
print(" specifying fixed grid as a quadrilateral")
print(" %i by %i, with %i points" \
% (FG.n12,FG.n23,npts))
print(" corner 1 = (%15.10f,%15.10f)" % (x1,y1))
print(" corner 2 = (%15.10f,%15.10f)" % (x2,y2))
print(" corner 3 = (%15.10f,%15.10f)" % (x3,y3))
print(" corner 4 = (%15.10f,%15.10f)" % (x4,y4))
xy = [x1,y1,x2,y2,x3,y3,x4,y4]
fname_root = os.path.splitext(FG.fname)[0]
kml_file = fname_root + '.kml'
kmltools.quad2kml(xy, kml_file, fname_root, color='8888FF')
def write_input_data(self, input_file_name=None):
r"""
Write input data from a file like *fgmax.txt* that
will be read in by GeoClaw as input data.
"""
if input_file_name is not None:
self.input_file_name = input_file_name
print("---------------------------------------------- ")
point_style = self.point_style
if point_style not in [0,1,2,3]:
raise NotImplementedError("make_fgmax not implemented for point_style %i" \
% point_style)
# write header, independent of point_style:
fid = open(self.input_file_name,'w')
fid.write("%16.10e # tstart_max\n" % self.tstart_max)
fid.write("%16.10e # tend_max\n" % self.tend_max)
fid.write("%16.10e # dt_check\n" % self.dt_check)
fid.write("%i %s # min_level_check\n" \
% (self.min_level_check,16*" "))
fid.write("%16.10e # arrival_tol\n" % self.arrival_tol)
fid.write("%i %s # point_style\n" \
% (self.point_style,16*" "))
if point_style == 0:
# list of points
npts = self.npts
print("Creating unstructured grid of %s points" % npts)
fid.write("%i # npts\n" % (npts))
for k in range(npts):
fid.write("%g %g \n" % (self.X[k],self.Y[k]))
fid.close()
print("Created file ", self.input_file_name)
elif point_style==1:
# 1d transect of points
x1,x2 = self.x1, self.x2
y1,y2 = self.y1, self.y2
if self.npts is None:
dx = self.dx
npts = int(round(sqrt((x2-x1)**2 + (y2-y1)**2)/dx)) + 1
if abs((npts-1)*dx + x1 - x2) > 1e-6:
print("Warning: abs((npts-1)*dx + x1 - x2) = ", \
abs((npts-1)*dx + x1 - x2))
x2 = x1 + dx*(npts-1)
y2 = y1 + dx*(npts-1)
print(" resetting x2 to %g" % x2)
print(" resetting y2 to %g" % y2)
else:
npts = self.npts
dx = sqrt((x2-x1)**2 + (y2-y1)**2)/(npts+1.)
if self.dx is not None:
print("*** Warning: dx specified over-ridden by: ",dx)
print("Creating 1d fixed grid with %s points" % npts)
print(" dx = %g" % dx)
fid.write("%i # npts\n" % (npts))
fid.write("%g %g # x1, y1\n" % (x1,y1))
fid.write("%g %g # x2, y2\n" % (x2,y2))
fid.close()
print("Created file ", self.input_file_name)
print(" specifying fixed grid with %i points equally spaced from " \
% npts)
print(" (%g,%g) to (%g,%g)" % (x1,y1,x2,y2))
# not yet implemented:
#fname_root = os.path.splitext(self.input_file_name)[0]
#kml_file = fname_root + '.kml'
#kmltools.line2kml(xy, kml_file, fname_root, color='8888FF')
if point_style == 2:
# 2d grid of points
x1,x2 = self.x1, self.x2
y1,y2 = self.y1, self.y2
if self.nx is None:
dx = self.dx
nx = int(round((x2-x1)/dx)) + 1
if abs((nx-1)*dx + x1 - x2) > 1e-6:
print("Warning: abs((nx-1)*dx + x1 - x2) = ", \
abs((nx-1)*dx + x1 - x2))
print(" old x2: %22.16e" % x2)
x2 = x1 + dx*(nx-1)
print(" resetting x2 to %22.16e" % x2)
else:
nx = self.nx
dx = (x2-x1)/(nx+1.)
if self.dx is not None:
print("*** Warning: dx specified over-ridden by: ",dx)
if self.ny is None:
dy = self.dy
if dy is None:
dy = dx
ny = int(round((y2-y1)/dy)) + 1
if abs((ny-1)*dy + y1 - y2) > 1e-6:
print("Warning: abs((ny-1)*dy + y1 - y2) = ", \
abs((ny-1)*dy + y1 - y2))
print(" old y2: %22.16e" % y2)
y2 = y1 + dy*(ny-1)
print(" resetting y2 to %22.16e" % y2)
else:
ny = self.ny
dy = (y2-y1)/(ny+1.)
if self.dy is not None:
print("*** Warning: dy specified over-ridden by: ",dy)
npts = nx*ny
fid.write("%i %i %s # nx,ny\n" \
% (nx,ny,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1,y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2,y2))
fid.close()
print("Created file ", self.input_file_name)
print(" specifying fixed grid with shape %i by %i, with %i points" \
% (nx,ny,npts))
print(" lower left = (%15.10f,%15.10f)" % (x1,y1))
print(" upper right = (%15.10f,%15.10f)" % (x2,y2))
print(" dx = %15.10e, dy = %15.10e" % (dx,dy))
xy = [x1,x2,y1,y2]
fname_root = os.path.splitext(self.input_file_name)[0]
kml_file = fname_root + '.kml'
kmltools.box2kml(xy, kml_file, fname_root, color='8888FF')
elif point_style==3:
# arbitrary quadrilateral
x1,x2 = self.x1, self.x2
y1,y2 = self.y1, self.y2
x3,x4 = self.x3, self.x4
y3,y4 = self.y3, self.y4
if self.n12 is None:
raise NotImplementedError("Need to set n12 and n23")
else:
npts = self.n12 * self.n23
fid.write("%i %i %s # self.n12,self.n23\n" \
% (self.n12,self.n23,10*" "))
fid.write("%16.10e %20.10e # x1, y1\n" % (x1,y1))
fid.write("%16.10e %20.10e # x2, y2\n" % (x2,y2))
fid.write("%16.10e %20.10e # x3, y3\n" % (x3,y3))
fid.write("%16.10e %20.10e # x4, y4\n" % (x4,y4))
fid.close()
print("Created file ", self.input_file_name)
print(" specifying fixed grid as a quadrilateral")
print(" %i by %i, with %i points" \
% (self.n12,self.n23,npts))
print(" corner 1 = (%15.10f,%15.10f)" % (x1,y1))
print(" corner 2 = (%15.10f,%15.10f)" % (x2,y2))
print(" corner 3 = (%15.10f,%15.10f)" % (x3,y3))
print(" corner 4 = (%15.10f,%15.10f)" % (x4,y4))
xy = [x1,y1,x2,y2,x3,y3,x4,y4]
fname_root = os.path.splitext(self.input_file_name)[0]
kml_file = fname_root + '.kml'
kmltools.quad2kml(xy, kml_file, fname_root, color='8888FF')
