Changeset 6713 for anuga_work/development/anuga_1d/discontinuous_z.py
- Timestamp:
- Apr 3, 2009, 12:34:47 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/anuga_1d/discontinuous_z.py
r6694 r6713 19 19 domain.set_timestepping_method('rk2') 20 20 domain.cfl = 1.0 21 domain.limiter = " minmod"21 domain.limiter = "vanleer" 22 22 23 23 def stage_flat(x): … … 85 85 domain.set_quantity('elevation',elevation_arbitrary) 86 86 domain.set_boundary({'exterior':Reflective_boundary(domain)}) 87 X=domain.vertices 88 C=domain.centroids 87 X=domain.vertices.flat 88 C=domain.centroids.flat 89 89 90 90 t0=time.time() … … 95 95 print "t=",t 96 96 print "integral", domain.quantities['height'].get_integral() 97 StageQ=domain.quantities['stage'].vertex_values 98 MomentumQ=domain.quantities['xmomentum'].vertex_values 99 ElevationQ=domain.quantities['elevation'].vertex_values 100 VelocityQ=domain.quantities['velocity'].vertex_values 97 StageQ=domain.quantities['stage'].vertex_values.flat 98 MomentumQ=domain.quantities['xmomentum'].vertex_values.flat 99 ElevationQ=domain.quantities['elevation'].vertex_values.flat 100 VelocityQ=domain.quantities['velocity'].vertex_values.flat 101 101 102 102 hold(False)
Note: See TracChangeset
for help on using the changeset viewer.