Changeset 2638 for inundation/examples/island.py
- Timestamp:
- Mar 30, 2006, 10:13:03 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/examples/island.py
r2636 r2638 2 2 3 3 Island surrounded by water. 4 This example investigates onshore 'creep' 4 This example investigates onshore 'creep' 5 5 6 6 """ … … 37 37 , interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 3)] 38 38 ) 39 40 41 39 40 41 42 42 #Create shallow water domain 43 43 domain = Domain(mesh_filename = 'island.msh') … … 66 66 z[i] = 8*exp( -((x[i]-50)**2 + (y[i]-50)**2)/100 ) 67 67 68 #z[i] += 0.5*exp( -((x[i]-10)**2 + (y[i]-10)**2)/50 ) 69 68 #z[i] += 0.5*exp( -((x[i]-10)**2 + (y[i]-10)**2)/50 ) 69 70 70 return z 71 71 … … 73 73 z = 0*x 74 74 for i in range(len(x)): 75 z[i] -= 0.7*exp( -((x[i]-10)**2 + (y[i]-10)**2)/200 ) 76 75 z[i] -= 0.7*exp( -((x[i]-10)**2 + (y[i]-10)**2)/200 ) 76 77 77 return z 78 78 79 #domain.set_quantity('friction', 0.1) #Honky dory 80 domain.set_quantity('friction', 2) #Creep79 domain.set_quantity('friction', 0.1) #Honky dory 80 #domain.set_quantity('friction', 2) #Creep 81 81 domain.set_quantity('elevation', island) 82 82 domain.set_quantity('stage', 1)
Note: See TracChangeset
for help on using the changeset viewer.