Changeset 5587 for anuga_work/development/anuga_1d/shallow_water_domain.py
- Timestamp:
- Jul 30, 2008, 5:03:47 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/anuga_1d/shallow_water_domain.py
r5563 r5587 59 59 tagged_elements) 60 60 61 from config import minimum_allowed_height, g 61 from config import minimum_allowed_height, g, h0 62 62 self.minimum_allowed_height = minimum_allowed_height 63 63 self.g = g 64 self.h0 = h0 64 65 65 66 #forcing terms not included in 1d domain ?WHy? … … 308 309 """ 309 310 310 from config import g, epsilon 311 from config import g, epsilon, h0 311 312 from math import sqrt 312 313 from Numeric import array … … 333 334 h_left = 0.0 334 335 else: 335 u_left = uh_left/h_left 336 u_left = uh_left/(h_left + h0/h_left) 337 338 339 uh_left = u_left*h_left 336 340 337 341 … … 339 343 h_right = w_right-z 340 344 uh_right = q_right[1] 341 342 345 343 346 if h_right < epsilon: … … 345 348 h_right = 0.0 346 349 else: 347 u_right = uh_right/h_right 350 u_right = uh_right/(h_right + h0/h_right) 351 352 uh_right = u_right*h_right 348 353 349 354 #vh_left = q_left[2]
Note: See TracChangeset
for help on using the changeset viewer.