Changeset 5738 for anuga_work/development/anuga_1d/dam_h_elevation.py
- Timestamp:
- Sep 5, 2008, 3:06:01 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/development/anuga_1d/dam_h_elevation.py
r5737 r5738 42 42 43 43 domain.default_order = 2 44 domain.default_time_order = 244 domain.default_time_order = 1 45 45 domain.cfl = 1.0 46 46 domain.limiter = "vanleer" … … 61 61 for i in range(len(x)): 62 62 if 1100.0 <= x[i] <=1200.0: 63 y[i]=10.0 0163 y[i]=10.0 64 64 else: 65 65 y[i]=10.0 … … 140 140 from pylab import clf,plot,title,xlabel,ylabel,legend,savefig,show,hold,subplot,ion 141 141 #print 'Test1' 142 ion() 142 143 hold(False) 143 144 #clf() … … 150 151 151 152 plot(X,ElevationQ,X,HeightQ) 152 plot1.set_ylim([ -1,12])153 plot1.set_ylim([9.99,10.01]) 153 154 xlabel('Position') 154 155 ylabel('Stage') … … 161 162 plot2 = subplot(212) 162 163 plot(X,MomentumQ) 163 plot2.set_ylim([-1,1])164 #plot2.set_ylim([-1,1]) 164 165 165 166 #legend( ('Numerical Solution', 'for momentum'), … … 173 174 file += ".eps" 174 175 #savefig(file) 175 show()176 #show() 176 177 177 178 print 'That took %.2f seconds'%(time.time()-t0) 179 180 raw_input("Press any key")
Note: See TracChangeset
for help on using the changeset viewer.