# display_name: Python 3 # language: python # name: python3 # --- # %% from useful_scit.plot.fig_manip import subp_insert_abc # %% from sectional_v2.constants import get_plotpath from sectional_v2.util.practical_functions import make_folders from IPython.display import clear_output plot_path = get_plotpath('comparison') +'prof_map_abs' print(plot_path) make_folders(plot_path) fn_base = plot_path + '/abs_' print(fn_base) make_folders(fn_base) # %% from sectional_v2.util.plot.plot_maps import plot_map_diff, fix_axis4map_plot, plot_map_abs_abs_diff, plot_map,plot_map_diff_only from useful_scit.imps import (np, xr, plt, pd) from sectional_v2.util.imports import get_averaged_fields from sectional_v2.util.plot.plot_profiles import plot_profile, set_legend_area_profs, set_scalar_formatter # load and autoreload from useful_scit.plot import get_cmap_dic from IPython import get_ipython
# %% [markdown] # ## Cases # %% cases_sec = ['SECTv21_ctrl'] cases_orig = ['noSECTv21_default_dd', 'noSECTv21_ox_ricc_dd'] # cases_orig =['noSECTv21_ox_ricc'] cases = cases_orig + cases_sec # %% # %% version = 'v21dd' plot_path = get_plotpath('levlat') filen_base = plot_path + '/_%s' % version # print(plot_path) make_folders(plot_path) # %% # %% def load_and_plot_diffs(varl, case_ctrl, case_other, start_time, end_time, pressure_coords=True, relative=False, cbar_orient='vertical', asp_ratio=2, subfig_size=3, ncol=None, ylim=None, yscale='log',
startyear = '2008-01' endyear = '2010-12' p_level = 1013. pmin = 850. # minimum pressure level avg_over_lev = True # True#True#False#True pressure_adjust = True # Can only be false if avg_over_lev false. Plots particular hybrid sigma lev if avg_over_lev: pressure_adjust = True p_levels = [1013., 900., 800., 700., 600.] # used if not avg # %% from sectional_v2.constants import get_plotpath from sectional_v2.util.practical_functions import make_folders version = 'v21dd_both' plot_path = get_plotpath('maps') filen_base = plot_path + '/_%s' % version #print(plot_path) make_folders(plot_path) # %% [markdown] # ## Cases # %% #cases_sec = ['SECTv21_ctrl'] #cases_orig =['noSECTv21_default'] #cases_orig =['noSECTv21_ox_ricc'] to_case = 'SECTv21_ctrl_koagD' from_cases = ['noSECTv21_default_dd', 'noSECTv21_ox_ricc_dd'] cases = [to_case] + from_cases
# noinspection PyBroadException try: _ipython = get_ipython() _magic = _ipython.magic _magic('load_ext autoreload') _magic('autoreload 2') except: pass # %% from sectional_v2.constants import get_outdata_path path_in = get_outdata_path('eusaar') version ='_noresmv21_dd' file_in = path_in + 'Nd_cat_sources_timeseries%s.csv'%version plot_path = get_plotpath('eusaar') version ='_noresmv21dd_both' # %% # %% # case_ns = 'noSECTv11_ctrl_fbvoc' # case_sec = 'SECTv11_ctrl_fbvoc' case_sec='SECTv21_ctrl_koagD' #'SECTv11_ctrl_fbvoc']#'SECTv11_ctrl']#,'SECTv11_ctrl_fbvoc']#'SECTv11_ctrl'] case_ns ='noSECTv21_ox_ricc' #'noSECTv11_ctrl_fbvoc'] #/no SECTv11_ctrl cases_ns = ['noSECTv21_default_dd','noSECTv21_ox_ricc_dd'] cases_s = [case_sec] cases = cases_ns + cases_s eusaar='eusaar' # %% [markdown] #
# format_version: '1.3' # jupytext_version: 1.3.2 # kernelspec: # display_name: Python 3 # language: python # name: python3 # --- # %% from sectional_v2.constants import get_plotpath from sectional_v2.util.practical_functions import make_folders from IPython.display import clear_output plot_path = get_plotpath('comparison') print(plot_path) make_folders(plot_path) fn_base = plot_path + '/prof_map_' print(fn_base) # %% from sectional_v2.util.plot.plot_maps import plot_map_diff, fix_axis4map_plot, plot_map_abs_abs_diff, plot_map,plot_map_diff_only from useful_scit.imps import (np, xr, plt, pd) from sectional_v2.util.imports import get_averaged_fields from sectional_v2.util.plot.plot_profiles import plot_profile, set_legend_area_profs, set_scalar_formatter # load and autoreload from useful_scit.plot import get_cmap_dic from IPython import get_ipython
try: _ipython = get_ipython() _magic = _ipython.magic _magic('load_ext autoreload') _magic('autoreload 2') except: pass # %% model = 'NorESM' # %% from sectional_v2.constants import get_plotpath from sectional_v2.util.practical_functions import make_folders plot_path = get_plotpath('comparison') + '/scatter/' print(plot_path) make_folders(plot_path) # %% from_time = '2008-01' to_time = '2008-12' pmin = 850. # minimum pressure level avg_over_lev = True # True#True#False#True pressure_adjust = False # Can only be false if avg_over_lev false. Plots particular hybrid sigma lev #if avg_over_lev: # pressure_adjust = True p_levels = [1013., 900., 800., 700., 600.] # used if not avg lev_lim = 0. # %%