# 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.
# %%
