- Timestamp:
- Aug 5, 2010, 7:16:33 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_work/development/sudi/sw_1d/avalanche/B_momentum/shallow_water_domain_avalanche.py
r7915 r7921 673 673 import sys 674 674 from Numeric import Float 675 from numpy import zeros675 from numpy import array, zeros 676 676 from config import epsilon, h0 677 677 … … 691 691 h_C = Height.centroid_values 692 692 u_C = Velocity.centroid_values 693 694 #for i in range(N):695 # h_C[i] = w_C[i] - z_C[i]696 # if h_C[i] <= epsilon :697 # uh_C[i] = 0.0698 # u_C[i] = 0.0699 # w_C[i] = z_C[i]700 #u_C[:] = uh_C/(h_C + h0/h_C)701 693 702 694 for i in range(N): … … 727 719 728 720 for i in range(N): 729 if h_V[i]<= 0.0:730 h_V[i] = 0.0721 if min(h_V[i]) <= 0.0: 722 h_V[i] = array([0.0, 0.0]) 731 723 stage_V[i] = bed_V[i] 732 xmom_V[i] = 0.0724 xmom_V[i] = array([0.0, 0.0]) 733 725 734 726 u_V[:] = xmom_V/(h_V + h0/h_V)
Note: See TracChangeset
for help on using the changeset viewer.