Changeset 7924
 Timestamp:
 Aug 5, 2010, 8:32:33 PM (12 years ago)
 Location:
 trunk/anuga_work/development/sudi/sw_1d/periodic_waves
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

trunk/anuga_work/development/sudi/sw_1d/periodic_waves/cg/runprogram.py
r7922 r7924 297 297 298 298 299 from pylab import clf,plot,title,xlabel,ylabel,legend,s avefig,show,hold,subplot,ion299 from pylab import clf,plot,title,xlabel,ylabel,legend,show,hold,subplot,ion#,savefig 300 300 hold(False) 301 301 clf() 302 302 303 303 plot1 = subplot(311) 304 plot(Ver /1e+4,AnalitW_V,'b', Ver/1e+4,StageV,'g', Ver/1e+4,ElevationV,'k')304 plot(Ver.flat/1e+4,AnalitW_V.flat,'b', Ver.flat/1e+4,StageV.flat,'g', Ver.flat/1e+4,ElevationV.flat,'k') 305 305 #plot(Ver,StageV, Ver,ElevationV) 306 306 #plot(Ver/L,StageV/h_0, Ver/L,ElevationV/h_0) … … 313 313 314 314 plot2 = subplot(312) 315 plot(Ver /1e+4,AnalitP_V,'b', Ver/1e+4,MomV,'g')315 plot(Ver.flat/1e+4,AnalitP_V.flat,'b', Ver.flat/1e+4,MomV.flat,'g') 316 316 #plot(Ver/L, VelV/sqrt(g*h_0)) 317 317 #xlabel('Position') … … 322 322 323 323 plot3 = subplot(313) 324 plot(Ver /1e+4,AnalitU_V,'b', Ver/1e+4,VelV,'g')324 plot(Ver.flat/1e+4,AnalitU_V.flat,'b', Ver.flat/1e+4,VelV.flat,'g') 325 325 #plot(Ver/L, VelV/sqrt(g*h_0)) 326 326 xlabel('Position / 10,000') … … 332 332 #filename = "%s%04i%s%f%s" %("numerical_cg_", counter,"_", domain.time, ".png") 333 333 #savefig(filename) 334 #show()334 show() 335 335 #raw_input("Press ENTER to continue") 336 336 
trunk/anuga_work/development/sudi/sw_1d/periodic_waves/johns/runprogram.py
r7922 r7924 299 299 300 300 #put this 301 from pylab import clf,plot,title,xlabel,ylabel,legend,s avefig,show,hold,subplot,ion301 from pylab import clf,plot,title,xlabel,ylabel,legend,show,hold,subplot,ion#,savefig 302 302 hold(False) 303 303 clf() 304 304 305 305 plot1 = subplot(311) 306 plot(Ver /1e+4,AnalitW_V,'b', Ver/1e+4,StageV,'g', Ver/1e+4,ElevationV,'k')306 plot(Ver.flat/1e+4,AnalitW_V.flat,'b', Ver.flat/1e+4,StageV.flat,'g', Ver.flat/1e+4,ElevationV.flat,'k') 307 307 #plot(Ver,StageV, Ver,ElevationV) 308 308 #plot(Ver/L,StageV/h_0, Ver/L,ElevationV/h_0) … … 315 315 316 316 plot2 = subplot(312) 317 plot(Ver /1e+4,AnalitP_V,'b', Ver/1e+4,MomV,'g')317 plot(Ver.flat/1e+4,AnalitP_V.flat,'b', Ver.flat/1e+4,MomV.flat,'g') 318 318 #plot(Ver/L, VelV/sqrt(g*h_0)) 319 319 #xlabel('Position') … … 324 324 325 325 plot3 = subplot(313) 326 plot(Ver /1e+4,AnalitU_V,'b', Ver/1e+4,VelV,'g')326 plot(Ver.flat/1e+4,AnalitU_V.flat,'b', Ver.flat/1e+4,VelV.flat,'g') 327 327 #plot(Ver/L, VelV/sqrt(g*h_0)) 328 328 xlabel('Position / 10,000') … … 332 332 'upper center', shadow=False) 333 333 334 filename = "%s%04i%s%f%s" %("numerical_johns_", counter,"_", domain.time, ".png")335 savefig(filename)336 #show()334 #filename = "%s%04i%s%f%s" %("numerical_johns_", counter,"_", domain.time, ".png") 335 #savefig(filename) 336 show() 337 337 #raw_input("Press ENTER to continue") 338 #put this339 338 counter = counter+1 340 339 finaltime = finaltime + waktu
Note: See TracChangeset
for help on using the changeset viewer.