def obsolete_plot_hazard_map(indir, datafile, out_dir, output_file, save_file=None,
                    title=None, np_posn=None, s_posn=None, cb_label=None,
                    cb_steps=None, invert=False, colourmap=None,
                    annotate=None):
    """Plot a hazard map from gridded data.

    indir       input directory
    data        name of the data file to plot (in indir)
    out_dir     general output directory
    output_file name of map output file to create in 'out_dir' directory
    save_file   name of map data output file to create in 'out_dir' directory
    title       title to put on the graph
    np_posn     place to put a north pointer symbol at
    s_posn      place to put a length scale at
    cb_label    label text to put on the colour bar
                (if not defined, no colour bar)
    cb_steps    sequence of required values in colourbar
    colourmap   the base colourmap to use
    annotate    user annotations required

    """

    if cb_steps is None:
        cb_steps = [0.1,0.15,0.20,0.30,0.40,0.50]

    data = cld.calc_load_data(datafile, invert=invert)
    pgxc.plot_gmt_xyz_contour(data, output_file,
                              title=title, np_posn=np_posn,
                              s_posn=s_posn, cb_label=cb_label,
                              cb_steps=cb_steps, annotate=annotate)
def do_example(datafile, outfile, invert=False):
#    ignore_value = 0.025
#
    annotations = []
#    if ignore_value:
#        annotations.append(('text', (143.6,-37.10),
#                            'Values <= %.2f are ignored' % ignore_value))

    data = cld.calc_load_data(datafile, invert=invert)
    # TODO : add call to calc_ignore
    pgxc.plot_gmt_xyz_contour(data, outfile,
                              title='Bedrock Hazard, RP=?, RSA=?',
                              np_posn='NE', s_posn='SE',
                              cb_label='Acceleration (g)',
                              cb_steps=[0.1,0.15,0.20,0.30,0.40,0.50],
                              annotate=annotations, linewidth=0.0)
def do_example(datafile, outfile, invert=False):
    #    ignore_value = 0.025
    #
    annotations = []
    #    if ignore_value:
    #        annotations.append(('text', (143.6,-37.10),
    #                            'Values <= %.2f are ignored' % ignore_value))

    data = cld.calc_load_data(datafile, invert=invert)
    # TODO : add call to calc_ignore
    pgxc.plot_gmt_xyz_contour(data,
                              outfile,
                              title='Bedrock Hazard, RP=?, RSA=?',
                              np_posn='NE',
                              s_posn='SE',
                              cb_label='Acceleration (g)',
                              cb_steps=[0.1, 0.15, 0.20, 0.30, 0.40, 0.50],
                              annotate=annotations,
                              linewidth=0.0)