Python plot_distributed_thickness Examples

Example #1

File: test_graphics.py Project: oberrauch/oggm

def test_thick_elev_bands():
    fig, ax = plt.subplots()
    gdir = init_columbia_eb(dir_name='test_thick_eb')
    workflow.inversion_tasks(utils.tolist(gdir))
    inversion.distribute_thickness_per_altitude(gdir)
    graphics.plot_distributed_thickness(gdir, ax=ax)
    fig.tight_layout()
    return fig

Example #2

File: test_graphics.py Project: zxl951128/oggm

def test_thick_interp():
    fig, ax = plt.subplots()
    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, ax=ax, varname_suffix='_interp')
    fig.tight_layout()
    return fig

Example #3

File: run_itmix_synth.py Project: JohannesUIBK/oggm

if do_itmix:

    done = False

    # V1
    distrib = partial(distribute_thickness, how='per_altitude',
                      add_slope=True,
                      smooth=True)
    execute_entity_task(distrib, gdirs)
    pdir = os.path.join(PLOTS_DIR, 'out_dis') + '/'
    if not os.path.exists(pdir):
        os.mkdir(pdir)
    for gd in gdirs:
        itmix.write_itmix_ascii(gd, 1)
        graphics.plot_distributed_thickness(gd)
        plt.savefig(pdir + gd.name + '_' + gd.rgi_id + '_d1.png')
        plt.close()

    # V2
    distrib = partial(distribute_thickness, how='per_altitude',
                      add_slope=False,
                      smooth=True)
    execute_entity_task(distrib, gdirs)
    for gd in gdirs:
        itmix.write_itmix_ascii(gd, 2)
        graphics.plot_distributed_thickness(gd)
        plt.savefig(pdir + gd.name + '_' + gd.rgi_id + '_d2.png')
        plt.close()

Example #4

def test_thick_interp():
    fig, ax = plt.subplots()
    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, ax=ax, how='per_interpolation')
    fig.tight_layout()
    return fig

Example #5

def test_thick_alt():
    fig, ax = plt.subplots()
    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, ax=ax, how='per_altitude')
    fig.tight_layout()
    return fig

Example #6

File: run_itmix.py Project: anton-ub/oggm

        _addt = addt if 'Columbia' in gd.name else ''
        graphics.plot_inversion(gd, add_title_comment=_addt)
        plt.savefig(pdir + gd.name + '_' + gd.rgi_id + '_inv.png')
        plt.close()

    # V1
    distrib = partial(distribute_thickness, how='per_altitude',
                      add_slope=True,
                      smooth=True)
    execute_entity_task(distrib, gdirs)
    pdir = os.path.join(PLOTS_DIR, 'out_dis') + '/'
    if not os.path.exists(pdir):
        os.mkdir(pdir)
    for gd in gdirs:
        itmix.write_itmix_ascii(gd, 1)
        graphics.plot_distributed_thickness(gd)
        plt.savefig(pdir + gd.name + '_' + gd.rgi_id + '_d1.png')
        plt.close()

    # V2
    distrib = partial(distribute_thickness, how='per_altitude',
                      add_slope=False,
                      smooth=True)
    execute_entity_task(distrib, gdirs)
    for gd in gdirs:
        itmix.write_itmix_ascii(gd, 2)
        graphics.plot_distributed_thickness(gd)
        plt.savefig(pdir + gd.name + '_' + gd.rgi_id + '_d2.png')
        plt.close()

    # Columbia

Example #7

File: test_graphics.py Project: fmaussion/oggm

def test_plot_distrib():

    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, how='per_altitude')
    graphics.plot_distributed_thickness(gdir, how='per_interpolation')

Example #8

File: test_graphics.py Project: MachineAi/oggm

def test_plot_distrib():

    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, how='per_altitude')
    graphics.plot_distributed_thickness(gdir, how='per_interpolation')

Example #9

File: test_graphics.py Project: JohannesUIBK/oggm

def test_thick_interp():
    fig, ax = plt.subplots()
    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, ax=ax, how='per_interpolation')
    fig.tight_layout()
    return fig

Example #10

File: test_graphics.py Project: JohannesUIBK/oggm

def test_thick_alt():
    fig, ax = plt.subplots()
    gdir = init_hef()
    graphics.plot_distributed_thickness(gdir, ax=ax, how='per_altitude')
    fig.tight_layout()
    return fig