Changeset 5180


Ignore:
Timestamp:
Mar 28, 2008, 2:13:13 PM (16 years ago)
Author:
steve
Message:

Added an example which has soliton type solutions (tilted.py)

Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/test_eq.py

    r5179 r5180  
    2323        test = -2205.904774487
    2424
    25         print 'depth', depth
     25        #print 'depth', depth
    2626        z = Okada_func(length=length, width=width, dip=dip, \
    2727                       x0=x0, y0=y0, strike=strike, depth=depth, \
     
    3838        assert allclose(z.rake, rake)
    3939
    40         print 'in test', z.test
     40        #print 'in test', z.test
    4141        assert allclose(z.test, -2205.904774487)
    4242       
    43         print 'hello finished okada'
     43        #print 'hello finished okada'
    4444
    4545    def test_earthquake_tsunami(self):
  • anuga_validation/convergence_study/wave.py

    r5175 r5180  
    4242# Setup domain
    4343#------------------------------------------------------------------------------
    44 dx = 200.
     44dx = 100.
    4545dy = dx
    4646L = 100000.
     
    7676domain.beta_h    = 0.0
    7777
     78interactive_visualisation = False
     79
    7880
    7981#------------------------------------------------------------------------------
     
    107109
    108110#===============================================================================
    109 from anuga.visualiser import RealtimeVisualiser
    110 vis = RealtimeVisualiser(domain)
    111 vis.render_quantity_height("stage", zScale =10000, dynamic=True)
    112 vis.colour_height_quantity('stage', (1.0, 0.5, 0.5))
    113 vis.start()
     111if interactive_visualisation:
     112    from anuga.visualiser import RealtimeVisualiser
     113    vis = RealtimeVisualiser(domain)
     114    vis.render_quantity_height("stage", zScale =10000, dynamic=True)
     115    vis.colour_height_quantity('stage', (1.0, 0.5, 0.5))
     116    vis.start()
    114117#===============================================================================
    115118
     
    121124for t in domain.evolve(yieldstep = 50.0, finaltime = 60*60.):
    122125    domain.write_time()
    123     vis.update()
    124    
    125 vis.evolveFinished()
     126    if interactive_visualisation:
     127        vis.update()
    126128
     129if interactive_visualisation:
     130    vis.evolveFinished()
     131
Note: See TracChangeset for help on using the changeset viewer.