Changeset 5737 for anuga_work/development/anuga_1d/dam_h_elevation.py
- Timestamp:
- Sep 5, 2008, 1:59:15 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/anuga_1d/dam_h_elevation.py
r5733 r5737 61 61 for i in range(len(x)): 62 62 if 1100.0 <= x[i] <=1200.0: 63 y[i]=10.0 63 y[i]=10.001 64 64 else: 65 65 y[i]=10.0 … … 106 106 return y 107 107 108 domain.set_quantity('stage',stage )109 domain.set_quantity('elevation', elevation_ box)108 domain.set_quantity('stage',stage_perturb) 109 domain.set_quantity('elevation', elevation_parabol) 110 110 #domain.set_quantity('xmomentum', xmom_sincos) 111 111 domain.order=domain.default_order … … 122 122 import time 123 123 t0=time.time() 124 yieldstep= 0.5#30.0124 yieldstep=10.0 #30.0 125 125 finaltime=10.0 #20.0 126 126 print "integral", domain.quantities['stage'].get_integral()
Note: See TracChangeset
for help on using the changeset viewer.