Changeset 4720 for anuga_work/development/demos/island.py
 Timestamp:
 Sep 10, 2007, 6:32:46 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/demos/island.py
r4631 r4720 3 3 Island surrounded by water. 4 4 This example investigates onshore 'creep' 5 6 Creep is defined as water moving uphill over many timesteps. 7 The initial adjustment at time = 0 is not 'creep' 5 8 6 9 """ … … 33 36 'top': [2], 34 37 'left': [3]}, 35 maximum_triangle_area = 5,38 maximum_triangle_area = 2, 36 39 filename = 'island.msh' , 37 40 interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 10.0)] … … 60 63 #domain.alpha_balance = 10.0 61 64 #domain.minimum_allowed_height = 1.0e4 62 #domain.maximum_allowed_speed = 100.065 #domain.maximum_allowed_speed = 0 63 66 #domain.minimum_storable_height = 1.0e4 64 67 65 domain.beta_h = 0.0 66 domain.tight_slope_limiters = 1 68 # 69 # Combinations with creep 70 # 71 72 #domain.maximum_allowed_speed = 0 73 #domain.beta_h = 0.0 74 #domain.tight_slope_limiters = 0 75 76 #domain.maximum_allowed_speed = 0 77 #domain.beta_h = 0.0 78 #domain.H0 = 0.01 79 #domain.tight_slope_limiters = 1 80 81 # 82 # Combinations with less creep ?? 83 # 84 85 domain.maximum_allowed_speed = 1 86 domain.beta_h = 0.0 87 domain.tight_slope_limiters = 0 88 89 #domain.maximum_allowed_speed = 10 90 #domain.beta_h = 0.0 91 #domain.H0 = 0.01 92 #domain.tight_slope_limiters = 1 93 94 95 96 67 97 68 98 # … … 86 116 return z 87 117 88 stage_value = 1 5.0118 stage_value = 10.0 89 119 #domain.set_quantity('friction', 0.1) #Honky dory 90 120 domain.set_quantity('friction', 0.01) #Creep … … 110 140 111 141 import time 112 for t in domain.evolve(yieldstep = 1, finaltime = 5):142 for t in domain.evolve(yieldstep = 1, finaltime = 25): 113 143 domain.write_time() 114 144 #if allclose(t, 100):
Note: See TracChangeset
for help on using the changeset viewer.