Changeset 5742 for anuga_work/development/anuga_1d/dam_h_elevation.py
- Timestamp:
- Sep 6, 2008, 4:30:16 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/anuga_1d/dam_h_elevation.py
r5741 r5742 1 1 import os 2 2 from math import sqrt, sin, cos, pi, exp 3 from shallow_water_domain _newimport *3 from shallow_water_domain import * 4 4 from Numeric import zeros, Float 5 5 #from analytic_dam_sudi import AnalyticDam … … 31 31 32 32 L=2000.0 33 N= 40033 N=200 34 34 35 35 cell_len=L/N … … 61 61 y=zeros(len(x), Float) 62 62 for i in range(len(x)): 63 if 1100.0 <= x[i] <=1200.0:63 if 800.0 <= x[i] <=1200.0: 64 64 y[i]=10.0 65 65 else: 66 y[i]= 10.066 y[i]=5.0 67 67 return y 68 68 … … 108 108 109 109 domain.set_quantity('stage',stage_perturb) 110 domain.set_quantity('elevation', elevation_ parabol)110 domain.set_quantity('elevation', elevation_box) 111 111 #domain.set_quantity('xmomentum', xmom_sincos) 112 112 domain.order=domain.default_order … … 152 152 153 153 plot(X,ElevationQ,X,HeightQ) 154 plot1.set_ylim([9.99,10.01])154 #plot1.set_ylim([9.99,10.01]) 155 155 xlabel('Position') 156 156 ylabel('Stage') … … 180 180 181 181 raw_input("Press any key") 182
Note: See TracChangeset
for help on using the changeset viewer.