Changeset 2207


Ignore:
Timestamp:
Jan 13, 2006, 4:53:30 PM (18 years ago)
Author:
steve
Message:
 
Location:
inundation/pyvolution
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/pyvolution/netherlands.py

    r2105 r2207  
    1515     Transmissive_boundary, Constant_height, Constant_stage
    1616
    17 from mesh_factory import rectangular
     17from mesh_factory import rectangular_cross
    1818from Numeric import array
    1919
     
    8383N = 130 #size = 33800
    8484N = 600 #Size = 720000
    85 N = 50
     85N = 20
    8686
    8787
     
    9090print 'Creating domain'
    9191#Create basic mesh
    92 points, elements, boundary = rectangular(N, N)
     92points, elements, boundary = rectangular_cross(N, N)
    9393
    9494#Create shallow water domain
     
    168168#V.update_quantity('elevation')
    169169
    170 for t in domain.evolve(yieldstep = 0.005, finaltime = 15.0):
     170
     171
     172
     173for t in domain.evolve(yieldstep = 0.05, finaltime = 1.0):
    171174    domain.write_time()
    172175    #domain.write_boundary()
     
    174177    print domain.quantities['stage'].get_values(location='centroids',
    175178                                                indices=[0])
     179    #domain.visualiser.update_quantity('elevation')
     180    #time.sleep(0.1)
     181    #raw_input('pause>')
    176182    #V.update_quantity('stage')
    177183    #rpdb.set_active()
Note: See TracChangeset for help on using the changeset viewer.