def MeshPlot(self, Lon, Lat, Elev, Cmap=[], Clim=[], Mesh=True):
    from matplotlib.colors import BoundaryNorm
    from matplotlib.ticker import MaxNLocator
    (x, y) = self._map(Lon, Lat)
    z = Elev
    if (not Cmap):
        Cmap = cm.jet
    if (not Clim):
        Clim = [z.min(), z.max()]
    levels = MaxNLocator(nbins=16).tick_values(Clim[0], Clim[1])
    norm = BoundaryNorm(levels, ncolors=Cmap.N, clip=True)
    if (not Mesh):
        self._zo += 1
        h = plt.scatter(x, y, c=z, s=20, marker='s', cmap=Cmap, vmin=Clim[0], vmax=Clim[1], lw=0, alpha=1.0, zorder=self._zo)
    else:
        self._zo += 1
        z = z[:(- 1), :(- 1)]
        h = plt.pcolormesh(x, y, z, cmap=Cmap, norm=norm, vmin=Clim[0], vmax=Clim[1], lw=0, alpha=1.0, zorder=self._zo)
    clb = plt.gcf().colorbar(h, orientation='vertical')
    clb.outline.set_linewidth(1)
    clb.ax.tick_params(labelsize=14)
    clb.set_label('Spectral Acceleration ($g$)', size=12)
def MeshPlot(self, Lon, Lat, Elev, Cmap=[], Clim=[], Mesh=True):
    from matplotlib.colors import BoundaryNorm
    from matplotlib.ticker import MaxNLocator
    (x, y) = self._map(Lon, Lat)
    z = Elev
    if (not Cmap):
        Cmap = cm.jet
    if (not Clim):
        Clim = [z.min(), z.max()]
    levels = MaxNLocator(nbins=16).tick_values(Clim[0], Clim[1])
    norm = BoundaryNorm(levels, ncolors=Cmap.N, clip=True)
    if (not Mesh):
        self._zo += 1
        h = plt.scatter(x, y, c=z, s=20, marker='s', cmap=Cmap, vmin=Clim[0], vmax=Clim[1], lw=0, alpha=1.0, zorder=self._zo)
    else:
        self._zo += 1
        z = z[:(- 1), :(- 1)]
        h = plt.pcolormesh(x, y, z, cmap=Cmap, norm=norm, vmin=Clim[0], vmax=Clim[1], lw=0, alpha=1.0, zorder=self._zo)
    clb = plt.gcf().colorbar(h, orientation='vertical')
    clb.outline.set_linewidth(1)
    clb.ax.tick_params(labelsize=14)
    clb.set_label('Spectral Acceleration ($g$)', size=12)