Changeset 4992


Ignore:
Timestamp:
Feb 6, 2008, 9:15:00 AM (17 years ago)
Author:
steve
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed

  • anuga_validation/convergence_study/convergence_structured.py

    r4838 r4992  
    2323# Setup computational domain
    2424#------------------------------------------------------------------------------
    25 dx = 100.
     25dx = 1000.
    2626dy = dx
    2727L = 100000.
    … …  
    2929
    3030# structured mesh
    31 points, vertices, boundary = rectangular_cross(int(L/dx), int(W/dy),
    32                                                L, W, (0.0, -W/2)) # Basic mesh
    33 #points, vertices, boundary = rectangular_cross(666, 3, 100000, 3000, (0.0, -0.0)) # Basic mesh
    34 #points, vertices, boundary = rectangular_cross(530, 10, 5300, 100, (-5000.0, -50.0)) # Basic mesh
    35 #points, vertices, boundary = rectangular_cross(1000, 100, 20, 3) # Basic mesh
     31points, vertices, boundary = rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
     32
    3633domain = Domain(points, vertices, boundary)
    37 
    38 
    39 ## # unstructured mesh
    40 ## poly_domain = [[0,-W],[0,W],[L,W],[L,-W]]
    41 ## meshname = 'test.msh'
    42 ## from anuga.pmesh.mesh_interface import create_mesh_from_regions
    43 ## # Create mesh
    44 ## create_mesh_from_regions(poly_domain,
    45 ##                          boundary_tags={'left': [0], 'top': [1],
    46 ##                                         'right': [2], 'bottom': [3]},
    47 ##                          maximum_triangle_area = 1000,
    48 ##                          filename=meshname)
    49 
    50 ## # Create domain
    51 ## domain = Domain(meshname, use_cache=True, verbose = True)
    5234
    5335domain.set_timestepping_method('euler')
    … …  
    5638domain.set_datadir('.')                     # Use current directory for output
    5739
    58 domain.beta_w      = 100.0
     40domain.beta_w      = 1.0
    5941domain.beta_w_dry  = 0.2
    60 domain.beta_uh     = 100.0
     42domain.beta_uh     = 1.0
    6143domain.beta_uh_dry = 0.2
    62 domain.beta_vh     = 100.0
     44domain.beta_vh     = 1.0
    6345domain.beta_vh_dry = 0.2
    64 domain.beta_h      = 100.0
     46domain.beta_h      = 1.0
    6547
    6648#------------------------------------------------------------------------------
    … …  
    9476domain.set_boundary({'left': Bw, 'right': Bt, 'top': Br, 'bottom': Br})
    9577
    96 ## from anuga.visualiser import RealtimeVisualiser
    97 ## vis = RealtimeVisualiser(domain)
    98 ## vis.render_quantity_height("elevation", dynamic=False)
    99 ## vis.render_quantity_height("stage", dynamic=True, zScale=100.0)
    100 ## #vis.colour_height_quantity('stage', (0.0, 0.0, 0.8))
    101 ## vis.colour_height_quantity('stage', (lambda q:q['stage'], -1, 1))
    102 ## vis.start()
     78
     79#===============================================================================
     80from anuga.visualiser import RealtimeVisualiser
     81vis = RealtimeVisualiser(domain)
     82#vis.render_quantity_height("elevation", zScale=1, offset = 5.0, dynamic=False)
     83vis.render_quantity_height("stage", zScale =10000, dynamic=True)
     84#vis.colour_height_quantity('stage', (lambda q:q['stage'], -1.0, 1.0))
     85vis.colour_height_quantity('stage', (1.0, 0.5, 0.5))
     86vis.start()
     87#===============================================================================
     88
    10389
    10490#------------------------------------------------------------------------------
    … …  
    10692#------------------------------------------------------------------------------
    10793
    108 for t in domain.evolve(yieldstep = 20.0, finaltime = 40*60.):
     94for t in domain.evolve(yieldstep = 20.0, finaltime = 10*40*60.):
    10995    domain.write_time()
    110     #vis.update()
     96    vis.update()
    11197   
    112 #vis.evolveFinished()
    113 #vis.join()
     98vis.evolveFinished()
     99
Note: See TracChangeset for help on using the changeset viewer.