Changeset 4572


Ignore:
Timestamp:
Jul 2, 2007, 11:10:51 AM (18 years ago)
Author:
duncan
Message:

small fixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/fit_interpolate/fit.py

    r4569 r4572  
    471471def _fit_to_mesh(vertex_coordinates,
    472472                 triangles,
    473                  point_coordinates, # this can also be a .csv/.txt file name
     473                 point_coordinates, # this can also be a points file name
    474474                 point_attributes=None,
    475475                 alpha=DEFAULT_ALPHA,
     
    497497
    498498          point_coordinates: List of coordinate pairs [x, y] of data points
    499           (or an nx2 Numeric array). This can also be a .csv/.txt file name
    500           FIXME (Ole): Can't it be a pts file too?
     499          (or an nx2 Numeric array). This can also be a .csv/.txt/.pts
     500          file name.
    501501
    502502          alpha: Smoothing parameter.
     
    601601    if verbose: print 'tsh file %s loaded' %mesh_file
    602602
    603     # load in the .pts file
     603    # load in the points file
    604604    try:
    605         ###point_dict = import_points_file(point_file, verbose=verbose)
    606605        geo = Geospatial_data(point_file, verbose=verbose)
    607606    except IOError,e:
     
    611610
    612611    point_coordinates = geo.get_data_points(absolute=True)
    613     #point_dict['pointlist']
    614     #
    615     # return list of attribute titles, array of attributes
    616     # title_list,point_attributes = concatinate_attributelist(point_dict['attributelist'])
    617     title_list,point_attributes = concatinate_attributelist(geo.get_all_attributes())
    618 
    619     if mesh_dict.has_key('geo_reference') and not mesh_dict['geo_reference'] is None:
     612    title_list,point_attributes = concatinate_attributelist( \
     613        geo.get_all_attributes())
     614
     615    if mesh_dict.has_key('geo_reference') and \
     616           not mesh_dict['geo_reference'] is None:
    620617        mesh_origin = mesh_dict['geo_reference'].get_origin()
    621618    else:
    622         mesh_origin = None #(56, 0, 0) #FIXME(DSG-DSG)
     619        mesh_origin = None
    623620
    624621    if verbose: print "points file loaded"
     
    649646        mesh_dict['vertex_attribute_titles'] = title_list
    650647
    651     #FIXME (Ole): Remember to output mesh_origin as well
    652648    if verbose: print "exporting to file ", mesh_output_file
    653649
Note: See TracChangeset for help on using the changeset viewer.