 Timestamp:
 Jun 22, 2010, 5:30:32 PM (13 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/2010projects/anuga_1d/sww/run_dry_dam.py
r7860 r7868 12 12 13 13 h1 = 10.0 14 h0 = 0. 014 h0 = 0.1 15 15 16 16 def analytical_sol(C,t): … … 81 81 k = 0 82 82 83 N = 80083 N = 3200 84 84 print "Evaluating domain with %d cells" %N 85 domain = Domain(*uniform_mesh(N ))85 domain = Domain(*uniform_mesh(N,x_1=L)) 86 86 87 87 domain.set_quantity('stage', stage) … … 90 90 91 91 domain.set_boundary({'left': Br, 'right' : Br}) 92 domain. order = 293 domain.set_timestepping_method(' euler')92 domain.set_spatial_order(2) 93 domain.set_timestepping_method('rk2') 94 94 domain.set_CFL(1.0) 95 domain.set_limiter("minmod ")95 domain.set_limiter("minmod_kurganov") 96 96 #domain.h0=0.0001 97 97 … … 136 136 137 137 pylab.plot(X.flat,VelocityV.flat) 138 plot2.set_ylim([ 20,10])138 plot2.set_ylim([15,15]) 139 139 140 140 pylab.xlabel('Position')
Note: See TracChangeset
for help on using the changeset viewer.