Changeset 6695
 Timestamp:
 Apr 1, 2009, 8:55:37 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_work/development/anuga_1d/dam_h_sudi.py
r5587 r6695 2 2 from math import sqrt 3 3 #from shallow_water_h import * 4 from shallow_water_domain import *4 from shallow_water_domain_suggestion2 import * 5 5 from Numeric import zeros, Float 6 6 from analytic_dam_sudi import AnalyticDam … … 32 32 33 33 L=2000.0 34 N= 40034 N=200 35 35 36 36 cell_len=L/N … … 42 42 domain=Domain(points) 43 43 44 domain. default_order = 245 domain. default_time_order = 144 domain.order = 2 45 domain.set_timestepping_method('rk2') 46 46 domain.cfl = 1.0 47 47 domain.limiter = "vanleer" … … 60 60 return y 61 61 62 domain.set_quantity('stage',height) #('height', height) 63 domain.order=domain.default_order 64 print "domain order", domain.order 62 domain.set_quantity('stage',height) 65 63 66 64 domain.set_boundary({'exterior':Reflective_boundary(domain)}) … … 74 72 import time 75 73 t0=time.time() 76 yieldstep= 30.074 yieldstep=20.0 77 75 finaltime=20.0 78 76 print "integral", domain.quantities['stage'].get_integral() #['height'].get_integral()
Note: See TracChangeset
for help on using the changeset viewer.