Changeset 1014
- Timestamp:
- Mar 4, 2005, 7:57:38 PM (20 years ago)
- Location:
- inundation/ga/storm_surge/pyvolution
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/pyvolution/test_shallow_water.py
r1013 r1014 1539 1539 domain = Domain(points, vertices, boundary) 1540 1540 domain.smooth = False 1541 domain.default_order= 11541 domain.default_order=2 1542 1542 1543 1543 #IC … … 1592 1592 domain = Domain(points, vertices, boundary) 1593 1593 domain.smooth = False 1594 domain.default_order= 11594 domain.default_order=2 1595 1595 1596 1596 #IC … … 1643 1643 domain = Domain(points, vertices, boundary) 1644 1644 domain.smooth = False 1645 domain.default_order=1 1645 domain.default_order=1 #FIXME: 2 will fail 1646 domain.beta_h = 0 1646 1647 1647 1648 #IC … … 1655 1656 if 0.5 <= x[i] < 0.7: 1656 1657 z[i] = 0.4 1658 #z[i] = 0.7 1657 1659 if 0.7 <= x[i]: 1658 1660 z[i] = x[i]/3 … … 1689 1691 1690 1692 #Evolution 1691 for t in domain.evolve(yieldstep = 0.05, finaltime = 5.0):1693 for t in domain.evolve(yieldstep = 0.05, finaltime = 10.0): 1692 1694 volume = domain.quantities['stage'].get_integral() 1693 1695 #print t, volume, initial_volume -
inundation/ga/storm_surge/pyvolution/twolevels.py
r911 r1014 36 36 # beta_h == 0.0 reveals the problem 37 37 # beta_h > 0.2 alleviates it 38 domain.beta_h = 0. 038 domain.beta_h = 0.2 39 39 40 40 #IC
Note: See TracChangeset
for help on using the changeset viewer.