|
setplot_kml.py.html |
 |
|
Source file: setplot_kml.py
|
|
Directory: /Users/rjl/clawpack_src/clawpack_master/geoclaw/examples/tsunami/chile2010
|
|
Converted: Mon Feb 19 2024 at 14:25:29
using clawcode2html
|
|
This documentation file will
not reflect any later changes in the source file.
|
"""
Set up the plot figures, axes, and items to be done for each frame.
This module is imported by the plotting routines and then the
function setplot is called to set the plot parameters.
"""
import numpy as np
import matplotlib.pyplot as plt
from clawpack.geoclaw import topotools, kmltools
try:
TG32412 = np.loadtxt('32412_notide.txt')
except:
print("*** Could not load DART data file")
#--------------------------
def setplot(plotdata):
#--------------------------
"""
Specify what is to be plotted at each frame.
Input: plotdata, an instance of pyclaw.plotters.data.ClawPlotData.
Output: a modified version of plotdata.
"""
from clawpack.visclaw import colormaps, geoplot
from numpy import linspace
plotdata.clearfigures() # clear any old figures,axes,items data
plotdata.verbose = False
# To plot gauge locations on pcolor or contour plot, use this as
# an afteraxis function:
def addgauges(current_data):
from clawpack.visclaw import gaugetools
gaugetools.plot_gauge_locations(current_data.plotdata, \
gaugenos='all', format_string='ko', add_labels=True)
#-----------------------------------------
# Some global kml flags
#-----------------------------------------
plotdata.kml_name = "Chile 2010"
plotdata.kml_starttime = [2010,2,27,6,34,0] # Time of event in UTC [None]
plotdata.kml_tz_offset = 3 # Time zone offset (in hours) of event. [None]
plotdata.kml_index_fname = "Chile_2010" # name for .kmz and .kml files ["_GoogleEarth"]
# Set to a URL where KMZ file will be published.
# plotdata.kml_publish = None
# Colormap range
cmin = -0.2
cmax = 0.2
cmap = geoplot.googleearth_transparent
#-----------------------------------------------------------
# Figure for KML files
#----------------------------------------------------------
plotfigure = plotdata.new_plotfigure(name='Sea Surface',figno=1)
plotfigure.show = True
plotfigure.use_for_kml = True
plotfigure.kml_use_for_initial_view = True
plotfigure.kml_show_figure = True
# These override any axes limits set below in plotaxes
plotfigure.kml_xlimits = [-120,-60]
plotfigure.kml_ylimits = [-60, 0.0];
# Resolution (should be consistent with data)
# Refinement levels : [2,6]; max level = 3; num_cells = [30,30]
# rcl : resolution of the coarsest level in this figure
rcl = 5 # choose so rcl*figsize = integer multiple of num_cells
plotfigure.kml_figsize = [30.0,30.0]
plotfigure.kml_dpi = rcl*2*6 # Resolve all three levels
plotfigure.kml_tile_images = False # Tile images for faster loading. Requires GDAL
# Water
plotaxes = plotfigure.new_plotaxes('kml')
plotitem = plotaxes.new_plotitem(plot_type='2d_pcolor')
plotitem.plot_var = geoplot.surface_or_depth
plotitem.pcolor_cmap = cmap
plotitem.pcolor_cmin = cmin
plotitem.pcolor_cmax = cmax
def kml_colorbar(filename):
kmltools.kml_build_colorbar(filename,cmap=cmap,cmin=cmin,cmax=cmax,
label='meters', title='surface',
extend='both')
plotfigure.kml_colorbar = kml_colorbar
#-----------------------------------------------------------
# Figure for KML files (zoom)
#----------------------------------------------------------
plotfigure = plotdata.new_plotfigure(name='Sea Surface (zoom)',figno=2)
#plotfigure.show = False
plotfigure.use_for_kml = True
plotfigure.kml_use_for_initial_view = False # Use large plot for view
plotfigure.kml_show_figure = True
# Set Google Earth bounding box and figure size
plotfigure.kml_xlimits = [-84,-74]
plotfigure.kml_ylimits = [-18,-4]
plotfigure.kml_figsize = [10,14] # inches.
# Resolution
rcl = 10 # Over-resolve the coarsest level
plotfigure.kml_dpi = rcl*2*6 # Resolve all three levels
plotfigure.kml_tile_images = False # Tile images for faster loading.
plotaxes = plotfigure.new_plotaxes('kml')
plotitem = plotaxes.new_plotitem(plot_type='2d_pcolor')
plotitem.plot_var = geoplot.surface_or_depth
plotitem.pcolor_cmap = cmap
plotitem.pcolor_cmin = cmin
plotitem.pcolor_cmax = cmax
#-----------------------------------------
# Figures for gauges
#-----------------------------------------
plotfigure = plotdata.new_plotfigure(name='Surface at gauges', figno=300, \
type='each_gauge')
plotfigure.clf_each_gauge = True
# Set up for axes in this figure:
plotaxes = plotfigure.new_plotaxes()
plotaxes.xlimits = 'auto'
plotaxes.ylimits = 'auto'
# plotaxes.title = 'Surface'
# Plot surface as blue curve:
plotitem = plotaxes.new_plotitem(plot_type='1d_plot')
plotitem.plot_var = 3
plotitem.plotstyle = 'b-'
# Plot topo as green curve:
plotitem = plotaxes.new_plotitem(plot_type='1d_plot')
plotitem.show = True
def gaugetopo(current_data):
q = current_data.q
h = q[0,:]
eta = q[3,:]
topo = eta - h
return topo
plotitem.plot_var = gaugetopo
plotitem.plotstyle = 'g-'
def add_zeroline(current_data):
from pylab import plot, legend, xticks, floor, axis, xlabel
t = current_data.t
gaugeno = current_data.gaugeno
if gaugeno == 32412:
try:
plot(TG32412[:,0], TG32412[:,1], 'r')
legend(['GeoClaw','Obs'],loc='lower right')
except: pass
axis((0,t.max(),-0.3,0.3))
plot(t, 0*t, 'k')
n = int(floor(t.max()/3600.) + 2)
xticks([3600*i for i in range(n)], ['%i' % i for i in range(n)])
xlabel('time (hours)')
plotaxes.afteraxes = add_zeroline
#-----------------------------------------
# Parameters used only when creating html and/or latex hardcopy
# e.g., via pyclaw.plotters.frametools.printframes:
plotdata.printfigs = True # print figures
plotdata.print_format = 'png' # file format
plotdata.print_framenos = 'all' # list of frames to print
plotdata.print_gaugenos = 'all' # list of gauges to print
plotdata.print_fignos = 'all' # list of figures to print
plotdata.html = False # create html files of plots?
plotdata.html_movie = None # create html files of plots?
plotdata.html_homelink = '../README.html' # pointer for top of index
plotdata.latex = False # create latex file of plots?
plotdata.latex_figsperline = 2 # layout of plots
plotdata.latex_framesperline = 1 # layout of plots
plotdata.latex_makepdf = False # also run pdflatex?
plotdata.kml = True
return plotdata