Ignore:
Timestamp:
Mar 29, 2006, 4:33:30 PM (18 years ago)
Author:
ole
Message:

Work on creep demo

File:
1 edited

Legend:

Unmodified
Added
Removed

  • inundation/examples/island.py

    r2621 r2632  
    77
    88
    9 #-------------------------------------------------------------------------------
     9#------------------------------------------------------------------------------
    1010# Import necessary modules
    11 #-------------------------------------------------------------------------------
     11#------------------------------------------------------------------------------
    1212
    1313# Standard modules
    … …  
    2121
    2222
    23 #-------------------------------------------------------------------------------   
     23#------------------------------------------------------------------------------
    2424# Setup computational domain
    25 #-------------------------------------------------------------------------------
     25#------------------------------------------------------------------------------
    2626
    2727#Create basic mesh
    … …  
    3737                          maximum_triangle_area = 25,
    3838                          filename = 'island.msh',
    39                           interior_regions=[ ([[55,25], [75,25],
    40                                                [75,75], [55,75]], 3)]
     39                          interior_regions=[ ([[50,25], [70,25],
     40                                               [70,75], [50,75]], 3)]
    4141                          )
    4242                         
    … …  
    6363
    6464
    65 #-------------------------------------------------------------------------------
     65#------------------------------------------------------------------------------
    6666# Setup initial conditions
    67 #-------------------------------------------------------------------------------
     67#------------------------------------------------------------------------------
    6868
    6969def island(x, y):
    … …  
    7373    return z
    7474
    75 domain.set_quantity('friction', 0.0)
     75domain.set_quantity('friction', 2.0)
    7676domain.set_quantity('elevation', island)
    7777domain.set_quantity('stage', 1)
    7878
    7979
    80 #-------------------------------------------------------------------------------
     80#------------------------------------------------------------------------------
    8181# Setup boundary conditions (all reflective)
    82 #-------------------------------------------------------------------------------
     82#------------------------------------------------------------------------------
    8383
    8484Br = Reflective_boundary(domain)
    … …  
    8888
    8989
    90 #-------------------------------------------------------------------------------
     90#------------------------------------------------------------------------------
    9191# Evolve system through time
    92 #-------------------------------------------------------------------------------
     92#------------------------------------------------------------------------------
    9393
    9494import time
    95 for t in domain.evolve(yieldstep = 0.5, finaltime = 100):
     95for t in domain.evolve(yieldstep = 1, finaltime = 100):
    9696    domain.write_time()
     97#for t in domain.evolve(yieldstep = 0.05, finaltime = 5):
     98#    domain.write_time()
     99#    #print '    Volume: ', domain.get_quantity('stage').get_integral()
    97100
    98101print 'Done'
Note: See TracChangeset for help on using the changeset viewer.