plot_fgout.py.html
Source file:   plot_fgout.py
Directory:   /Users/rjl/clawpack_src/clawpack_master/geoclaw/examples/tsunami/chile2010_fgmax-fgout
Converted:   Mon Feb 19 2024 at 14:27:24   using clawcode2html
This documentation file will not reflect any later changes in the source file.

 

from pylab import *
import os
from clawpack.visclaw import plottools, geoplot
from clawpack.visclaw import animation_tools
from matplotlib import animation, colors
from clawpack.geoclaw import fgout_tools
from datetime import timedelta 

fgno = 1
outdir = '_output'
output_format = 'binary'

# Instantiate object for reading fgout frames:
fgout_grid = fgout_tools.FGoutGrid(fgno, outdir, output_format)

# Plot one frame of fgout data 
fgframe = 20
fgout = fgout_grid.read_frame(2)

fgout = fgout_grid.read_frame(fgframe)

figure(1, figsize=(8,8))
imshow(flipud(fgout.B.T), extent=fgout.extent_edges,
       cmap=geoplot.land_colors)

clim(0,100)

eta_water = where(fgout.h > 0, fgout.eta, nan)
imshow(flipud(eta_water.T), extent=fgout.extent_edges,
       cmap=geoplot.tsunami_colormap)

clim(-0.2, 0.2)

title('Surface at time %s' % timedelta(seconds=fgout.t))

fname = 'fgout_frame%s.png' % str(fgframe).zfill(4)
savefig(fname)
print('Created ',fname)