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()