Changeset 1744


Ignore:
Timestamp:
Aug 23, 2005, 3:04:27 PM (19 years ago)
Author:
ole
Message:

lwru2.py and future_directions.txt

Location:
inundation
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • inundation/validation/LWRU2/lwru2.py

    r1742 r1744  
    134134#                                              verbose=True)
    135135
    136 #vertex_attributes = cache(least_squares.fit_to_mesh,
    137 #                          (domain.coordinates, domain.triangles,
    138 #                           points, elevation) ,
    139 #                          {'alpha': 0.0001, 'verbose': True})
     136vertex_attributes = cache(least_squares.fit_to_mesh,
     137                          (domain.coordinates, domain.triangles,
     138                           points, elevation) ,
     139                          {'alpha': 0.0001, 'verbose': True})
    140140
    141141print 'Initial values'
    142 #domain.set_quantity('elevation', vertex_attributes)
    143 domain.set_quantity('elevation', points=points, attributes=attributes['elevation'], alpha=0.0001)
     142domain.set_quantity('elevation', vertex_attributes)
     143#domain.set_quantity('elevation', points=points, attributes=attributes['elevation'], alpha=0.0001)
    144144domain.set_quantity('friction', 0.0)
    145145domain.set_quantity('stage', 0.0)
  • inundation/wiki/future_directions.txt

    r1739 r1744  
    4646     
    4747Georeferencing
    48   Code works in UTM coordinates. Zone boundaries cannot be crossed
     48  Code works in UTM coordinates. Zone boundaries cannot be crossed: Introduce nonstandard contral meridians if necessary
     49 
    4950  Could code be cast in terms of lat/lons?
    5051  make sure that georeferencing is used consistently (domain, pmesh, sww files, ...) (Duncan) 
     
    5657        Give feedback when
    5758                points outside mesh
    58                 mesh outside points boundary
     59                mesh outside points boundary  !
    5960                if mesh is corrupted (see lwru2.py and create_mesh.py in debug dir)
    6061                if segments are repeated (maybe in mesh.check_integrity)   
Note: See TracChangeset for help on using the changeset viewer.