compare_gauges.py.html
Source file:   compare_gauges.py
Directory:   /var/www/html/clawpack/links/awr11/radial-ocean-island
Converted:   Sat Mar 5 2011 at 20:51:33   using clawcode2html
This documentation file will not reflect any later changes in the source file.

 

import pylab
from pyclaw.plotters.data import ClawPlotData

outdir1 = '_output_260'
outdir2 = '_output_220'
gaugenos = [1,2]

plotdata = ClawPlotData()
    
pylab.figure(301)
pylab.clf()

plotdata.outdir = outdir1
for gaugeno in gaugenos:
    gs = plotdata.getgauge(gaugeno)
    pylab.plot(gs.t, gs.q[:,3], 'b',linewidth=2)


plotdata.outdir = outdir2
for gaugeno in gaugenos:
    gs = plotdata.getgauge(gaugeno)
    pylab.plot(gs.t, gs.q[:,3], 'r--',linewidth=2)

pylab.xlim([7500,14000])
pylab.ylim([-2.5,4])
pylab.xticks(fontsize=15)
pylab.yticks(fontsize=15)

pylab.annotate('Gauge 1',[8200,1.9],[7700,2.5],arrowprops={'width':1,'color':'k'})
pylab.annotate('Gauge 2',[9300,2.7] ,[9700,3.1],arrowprops={'width':1,'color':'k'})

#pylab.savefig('gauges1-2.tif')
pylab.savefig('gauges1-2.png')

#print "Created gauges1-2.png and gauges1-2.tif"

#===============================

gaugenos = [3,4]
pylab.figure(302)
pylab.clf()

plotdata.outdir = outdir1
for gaugeno in gaugenos:
    gs = plotdata.getgauge(gaugeno)
    pylab.plot(gs.t, gs.q[:,3], 'b',linewidth=2)


plotdata.outdir = outdir2
for gaugeno in gaugenos:
    gs = plotdata.getgauge(gaugeno)
    pylab.plot(gs.t, gs.q[:,3], 'r--',linewidth=2)

pylab.xlim([7500,14000])
pylab.ylim([-2.5,4])
pylab.xticks(fontsize=15)
pylab.yticks(fontsize=15)

pylab.annotate('Gauge 3',[12500,2.6],[13000,3.1],arrowprops={'width':1,'color':'k'})
pylab.annotate('Gauge 4',[11550,1.0] ,[11200,2.0],arrowprops={'width':1,'color':'k'})

#pylab.savefig('gauges3-4.tif')
pylab.savefig('gauges3-4.png')

#print "Created gauges3-4.png and gauges3-4.tif"