plot_dart.py.html
Source file:   plot_dart.py
Directory:   /var/www/html/clawpack/links/awr11/chile2010awr
Converted:   Sat Feb 5 2011 at 15:39:50   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_4levels"
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')