Ignore:
Timestamp:
Jun 22, 2010, 5:30:32 PM (13 years ago)
Author:
steve
Message:

Added in some more c based limiters

File:
1 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/development/2010-projects/anuga_1d/sww/run_dry_dam.py

    r7860 r7868  
    1212
    1313h1 = 10.0
    14 h0 = 0.0
     14h0 = 0.1
    1515
    1616def analytical_sol(C,t):
    … …  
    8181k = 0
    8282
    83 N = 800
     83N = 3200
    8484print "Evaluating domain with %d cells" %N
    85 domain = Domain(*uniform_mesh(N))
     85domain = Domain(*uniform_mesh(N,x_1=L))
    8686   
    8787domain.set_quantity('stage', stage)
    … …  
    9090
    9191domain.set_boundary({'left': Br, 'right' : Br})
    92 domain.order = 2
    93 domain.set_timestepping_method('euler')
     92domain.set_spatial_order(2)
     93domain.set_timestepping_method('rk2')
    9494domain.set_CFL(1.0)
    95 domain.set_limiter("minmod")
     95domain.set_limiter("minmod_kurganov")
    9696#domain.h0=0.0001
    9797
    … …  
    136136
    137137pylab.plot(X.flat,VelocityV.flat)
    138 plot2.set_ylim([-20,10])
     138plot2.set_ylim([-15,15])
    139139
    140140pylab.xlabel('Position')
Note: See TracChangeset for help on using the changeset viewer.