plot_dart.py.html CLAWPACK  
 Source file:   plot_dart.py
 Directory:   /Users/rjl/git/rjleveque/clawpack-4.6.3/apps/tsunami/chile2010
 Converted:   Mon Jan 21 2013 at 20:15:27   using clawcode2html
 This documentation file will not reflect any later changes in the source file. 

 
"""
Compare DART buoy results to NOAA results.
"""

import pylab 
from pyclaw.plotters.data import ClawPlotData
from matplotlib import image

dartpng = image.imread('dart32412_comp-2.png')
pylab.figure(400)
pylab.clf()
pylab.imshow(dartpng)
pylab.hold(True)

# origin of plot in pixels on image:
t0 = 195.
y0 = 221.
tscale = 120.0
yscale = -446.66666666666669
plotdata = ClawPlotData()
plotdata.outdir = "_output"
gaugedata = plotdata.getgauge(32412)
t = gaugedata.t
eta = gaugedata.q[:,3]
t = t0 + tscale * (t / 3600.)
y = y0 + yscale*eta

pylab.plot(t,y,'b')
pylab.plot([160,225],[550,550],'b')
pylab.text(240,560,'GeoClaw (added to NOAA original)',fontsize=8)
pylab.text(100,600,'NOAA Original from http://nctr.pmel.noaa.gov/chile20100227/dart32412_comp-2.pdf',fontsize=8)


pylab.xlim([50,1300])
pylab.ylim([640,0])
pylab.axis('off')

pylab.savefig('dart.png')