Ignore:
Timestamp:
Mar 30, 2006, 10:13:03 PM (18 years ago)
Author:
steve
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/examples/island.py

    r2636 r2638  
    22
    33Island surrounded by water.
    4 This example investigates onshore 'creep' 
     4This example investigates onshore 'creep'
    55
    66"""
     
    3737                          , interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 3)]
    3838                          )
    39                          
    40                                                    
    41                          
     39
     40
     41
    4242#Create shallow water domain
    4343domain = Domain(mesh_filename = 'island.msh')
     
    6666        z[i] = 8*exp( -((x[i]-50)**2 + (y[i]-50)**2)/100 )
    6767
    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
    7070    return z
    7171
     
    7373    z = 0*x
    7474    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
    7777    return z
    7878
    79 #domain.set_quantity('friction', 0.1)  #Honky dory
    80 domain.set_quantity('friction', 2)     #Creep
     79domain.set_quantity('friction', 0.1)  #Honky dory
     80#domain.set_quantity('friction', 2)     #Creep
    8181domain.set_quantity('elevation', island)
    8282domain.set_quantity('stage', 1)
Note: See TracChangeset for help on using the changeset viewer.