- Timestamp:
- Mar 31, 2008, 4:39:49 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/fit_interpolate/test_interpolate.py
r4820 r5181 1582 1582 return -x 1583 1583 1584 # create mesh1584 # Create mesh 1585 1585 mesh_file = tempfile.mktemp(".tsh") 1586 1586 points = [[0.0,0.0],[6.0,0.0],[6.0,6.0],[0.0,6.0]] … … 1591 1591 m.export_mesh_file(mesh_file) 1592 1592 1593 # Create shallow water domain1593 # Create shallow water domain 1594 1594 domain = Domain(mesh_file) 1595 1595 os.remove(mesh_file) 1596 1596 1597 domain.default_order =21597 domain.default_order = 2 1598 1598 domain.beta_h = 0 1599 1599 1600 #Set some field values 1600 # This test was made before tight_slope_limiters were introduced 1601 # Since were are testing interpolation values this is OK 1602 domain.tight_slope_limiters = 0 1603 1604 # Set some field values 1601 1605 domain.set_quantity('elevation', elevation_function) 1602 1606 domain.set_quantity('friction', 0.03) … … 1622 1626 sww.store_connectivity() 1623 1627 sww.store_timestep(['stage', 'xmomentum', 'ymomentum']) 1624 domain.set_quantity('stage', 10.0) # This is automatically lim mited1625 # so it will not be less than the elevation1628 domain.set_quantity('stage', 10.0) # This is automatically limited 1629 # So it will not be less than the elevation 1626 1630 domain.time = 2. 1627 1631 sww.store_timestep(['stage', 'xmomentum', 'ymomentum']) 1628 1632 1629 # test the function1633 # Test the function 1630 1634 points = [[5.0,1.],[0.5,2.]] 1631 1635 depth_file = tempfile.mktemp(".csv")
Note: See TracChangeset
for help on using the changeset viewer.