compare_gauges.py.html
Source file:   compare_gauges.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.

 

import pylab
from pyclaw.plotters.data import ClawPlotData

outdir1 = '_output_noamr'
outdir2 = '_output'
outdir3 = '_output_4levels'
gaugenos = [32412]

plotdata = ClawPlotData()
    
pylab.figure(301,figsize=(10,5))
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=3)

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

n = pylab.floor(gs.t.max()/3600.) + 2
pylab.xticks([3600*i for i in range(n)], [str(i) for i in range(n)],fontsize=15)
d = pylab.loadtxt('32412_dart.txt')
pylab.plot(d[:,0],d[:,1],'k-')
pylab.legend(['Uniform grid','3 levels of AMR','4 levels of AMR','DART data'])

pylab.xlim([9000,18800])
pylab.ylim([-0.2,0.3])
pylab.yticks(fontsize=15)
pylab.xlabel("Hours after earthquake",fontsize=15)
pylab.ylabel("Meters",fontsize=15)
pylab.title("Surface elevation at DART 32412",fontsize=15)

pylab.savefig('dart_gauges.eps')
pylab.savefig('dart_gauges.png')

print "Created plots dart_gauges"