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()
for gd in gdirs:
_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()
for gd in gdirs:
_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()
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()