Changeset 5643
- Timestamp:
- Aug 13, 2008, 8:14:27 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/util.py
r5482 r5643 2262 2262 vel = 0.0 2263 2263 else: 2264 momentum = sqrt(point_quantities[2]**2 +\ 2265 point_quantities[3]**2) 2266 # vel = momentum/depth 2267 vel = momentum/(point_quantities[0] - point_quantities[1]) 2268 # vel = momentum/(depth + 1.e-6/depth) 2269 2264 if point_quantities[2] < 1.0e6: 2265 momentum = sqrt(point_quantities[2]**2 +\ 2266 point_quantities[3]**2) 2267 # vel = momentum/depth 2268 vel = momentum/(point_quantities[0] - point_quantities[1]) 2269 # vel = momentum/(depth + 1.e-6/depth) 2270 else: 2271 momentum = 0 2272 vel = 0 2273 2270 2274 points_list.append(vel) 2271 2275
Note: See TracChangeset
for help on using the changeset viewer.