Ignore:
Timestamp:
Sep 5, 2007, 5:17:30 PM (16 years ago)
Author:
ole
Message:

Some documentation of ticket:192 plus storage of restricting polygon and time_interval for extrema computations.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/test_data_manager.py

    r4704 r4705  
    274274                                               'ymomentum',
    275275                                               'stage-elevation'])
     276
     277        assert domain.monitor_polygon is None
     278        assert domain.monitor_time_interval is None
     279
     280
     281        domain.set_quantities_to_be_monitored(['xmomentum',
     282                                               'ymomentum',
     283                                               'stage-elevation'],
     284                                              polygon=domain.get_boundary_polygon(),
     285                                              time_interval=[0,1])
     286       
    276287       
    277288        assert len(domain.quantities_to_be_monitored) == 3
     
    279290        assert domain.quantities_to_be_monitored.has_key('xmomentum')               
    280291        assert domain.quantities_to_be_monitored.has_key('ymomentum')       
    281         assert domain.monitor_polygon is None
    282         assert domain.monitor_time_interval is None       
     292
    283293       
    284294        sww = get_dataobject(domain)
     
    287297            pass
    288298            #print domain.timestepping_statistics()
    289             #print domain.quantity_statistics(precision = '%.8f')
     299            domain.quantity_statistics(precision = '%.8f') # Silent
    290300
    291301           
     
    312322        assert allclose(extrema,[0.00, 0.06241221])       
    313323
     324        time_interval = fid.variables['extrema:time_interval'][:]
     325        assert allclose(time_interval, [0,1])
     326       
     327        polygon = fid.variables['extrema:polygon'][:]       
     328        assert allclose(polygon, domain.get_boundary_polygon())
    314329       
    315330        fid.close()
Note: See TracChangeset for help on using the changeset viewer.