old_contents: import pylab
  import numpy
  from util import runparams
  
  def dovis(myData, n):
  
      pylab.clf()
  
      dens = myData.getVarPtr("density")
  
      myg = myData.grid
  
      pylab.imshow(numpy.transpose(dens[myg.ilo:myg.ihi+1,myg.jlo:myg.jhi+1]), 
                   interpolation="nearest", origin="lower",
                   extent=[myg.xmin, myg.xmax, myg.ymin, myg.ymax])
  
      pylab.xlabel("x")
      pylab.ylabel("y")
      pylab.title("density")
  
      pylab.colorbar()
  
      pylab.draw()
  
  
new_contents: import pylab
  import numpy
  from util import runparams
  
  def dovis(myData, n):
  
      pylab.clf()
  
      dens = myData.getVarPtr("density")
  
      myg = myData.grid
  
      pylab.imshow(numpy.transpose(dens[myg.ilo:myg.ihi+1,myg.jlo:myg.jhi+1]), 
                   interpolation="nearest", origin="lower",
                   extent=[myg.xmin, myg.xmax, myg.ymin, myg.ymax])
  
      pylab.xlabel("x")
      pylab.ylabel("y")
      pylab.title("density")
  
      pylab.colorbar()
  
      pylab.figtext(0.05,0.0125, "t = %10.5f" % myData.t)
  
      pylab.draw()
  
  
subject: add time