Beispiel #1
0
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')
Beispiel #2
0
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')
Beispiel #3
0
    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')
Beispiel #4
0
    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')