Ignore:
Timestamp:
Mar 16, 2009, 11:06:22 AM (14 years ago)
Author:
rwilson
Message:

Hand-merged recent changes in main trunk. Still work to be done in shallow_water.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/numpy/anuga/geospatial_data/geospatial_data.py

    r6441 r6517  
    18061806    points = G_small.get_data_points()
    18071807
     1808    #FIXME: Remove points outside boundary polygon
     1809#    print 'new point',len(points)
     1810#
     1811#    new_points=[]
     1812#    new_points=array([],typecode=Float)
     1813#    new_points=resize(new_points,(len(points),2))
     1814#    print "BOUNDARY", boundary_poly
     1815#    for i,point in enumerate(points):
     1816#        if is_inside_polygon(point,boundary_poly, verbose=True):
     1817#            new_points[i] = point
     1818#            print"WOW",i,new_points[i]
     1819#    points = new_points
     1820
    18081821    if verbose: print "Number of points in sample to compare: ", len(points)
    18091822
Note: See TracChangeset for help on using the changeset viewer.