# stackbar_groups = [r"AR 1.1\,$R_\odot$" + "\n" + r"$\log n_e = 8.1$",
# r"AR 1.3\,$R_\odot$" + "\n" + r"$\log n_e = 7.7$",
# r"*AR 1.5\,$R_\odot$" + "\n" + r"$\log n_e = 7.3$",
# r"ST 1.1\,$R_\odot$" + "\n" + r"$\log n_e = 7.85$",
# r"ST 1.25\,$R_\odot$" + "\n" + r"$\log n_e = 7.5$",
# r"*ST 1.5\,$R_\odot$" + "\n" + r"$\log n_e = 7.0$"]
stackbar_groups = [r"AR 1.1\,$R_\odot$" + "\n" + r"$\log n_e = 8.1$",
r"AR 1.3\,$R_\odot$" + "\n" + r"$7.7$",
r"*AR 1.5\,$R_\odot$" + "\n" + r"$7.3$",
r"ST 1.1\,$R_\odot$" + "\n" + r"$7.85$",
r"ST 1.25\,$R_\odot$" + "\n" + r"$7.5$",
r"*ST 1.5\,$R_\odot$" + "\n" + r"$7.0$"]
stackbar_colors = ["#D0E6A5","#FFDD94","#FA897B","#CCABDB"]
fig = plt.figure(figsize=(9,9),constrained_layout=True)
subfigs = fig.subfigures(2,1,height_ratios=[9,7])
ax1,ax2,ax3 = subfigs[0].subplots(1,3,gridspec_kw={"width_ratios":[1,1,2]})
ax4 = subfigs[1].subplots()
im1 = ax1.pcolormesh(comp_xcoord, comp_ycoord, enintensity_1074, cmap=aia_193_cm,
norm=ImageNormalize(enintensity_1074),rasterized=True)
ax1.set_title(r"\textbf{Fe \textsc{xiii} 1074.7\,nm}" + "\n" + r"\textbf{Enhanced Intensity}",fontsize=14,
y=1.02)
ax1.set_ylabel(r"\textbf{SOLAR-Y [arcsec]}",fontsize=14)
ax1.set_xlabel(r"\textbf{SOLAR-X [arcsec]}",fontsize=14)
ax1.plot(AR_cut_solarx_1, AR_cut_solary_1, color="#77AADD",lw=4)
ax1.plot(AR_cut_solarx_2, AR_cut_solary_2, color="#44BB99",lw=4)
ax1.plot(ST_cut_solarx_1, ST_cut_solary_1, color="#EE6677",lw=4)
clb1, clb_ax1 = plot_colorbar(im1,ax1,width="10%",bbox_to_anchor=(0.63, -0.48, 1, 0.9))
clb_ax1.yaxis.set_major_locator(MaxNLocator(2))
ax1.text(0.05,0.98,r"\textbf{(a)}",fontsize=14,transform=ax1.transAxes,ha="left",va="top")
im2 = ax2.pcolormesh(comp_xcoord[:310],comp_ycoord,FeXIII_dens_1074_1079[:,:310],vmin=7.5,vmax=8.5,
cmap=cmcm.batlowK,rasterized=True)
ax2.set_xlabel(r"\textbf{SOLAR-X [arcsec]}",fontsize=14)
ax2.set_title(r"\textbf{Fe \textsc{xiii} 1074/1079}" "\n" + r"$\boldsymbol{\log n_e/\mathrm{[cm^{-3}]}}$",fontsize=14,
y=1.02)
ax2.tick_params(labelleft=False)
clb2, clb_ax2 = plot_colorbar(im2,ax2,width="10%",bbox_to_anchor=(0.63, -0.48, 1, 0.9))
clb_ax2.yaxis.set_major_locator(MaxNLocator(2))
ax2.text(0.05,0.98,r"\textbf{(b)}",fontsize=14,transform=ax2.transAxes,ha="left",va="top")
ax3.scatter(np.sqrt(AR_cut_solarx_1**2 + AR_cut_solary_1**2)/header_1074["RSUN"],AR_density_1,label=r"\textbf{AR Outskirts}",
color="#77AADD",marker="1",s=80)
ax3.scatter(np.sqrt(AR_cut_solarx_2**2 + AR_cut_solary_2**2)/header_1074["RSUN"],AR_density_2,label=r"\textbf{AR Core}",
color="#44BB99",marker="2",s=80)
ax3.scatter(np.sqrt(ST_cut_solarx_1**2 + ST_cut_solary_1**2)/header_1074["RSUN"],ST_density_1,label=r"\textbf{Streamer}",
color="#EE6677",marker="x")
ax3.yaxis.tick_right()
ax3.yaxis.set_label_position("right")
ax3.set_ylabel(r"$\boldsymbol{\log n_e/\mathrm{[cm^{-3}]}}$",fontsize=14,)
ax3.set_xlabel(r"\textbf{Heliocentric Distance} $\boldsymbol{r\,[R_\odot]}$",fontsize=14)
ax3.set_title(r"\textbf{Density Drop along Cuts}",fontsize=14,y=1.04)
ax3.legend(fontsize=14,frameon=False,loc="upper right",handletextpad=0)
ax3.text(0.03,0.02,r"\textbf{(c)}",fontsize=14,transform=ax3.transAxes,ha="left",va="bottom")
ax4.bar(stackbar_groups, pop_processes_df["ph"].values, color=stackbar_colors[0], label=r"\textbf{PE}",zorder=2,
)
ax4.bar(stackbar_groups, pop_processes_df["rd"].values, bottom=pop_processes_df["ph"].values,
color=stackbar_colors[1], label=r"\textbf{RD}",zorder=2)
ax4.bar(stackbar_groups, pop_processes_df["e"].values,
bottom=pop_processes_df["rd"].values+pop_processes_df["ph"].values,
color=stackbar_colors[2], label=r"$\boldsymbol{e^-}$",zorder=2)
ax4.bar(stackbar_groups, pop_processes_df["p"].values,
bottom=pop_processes_df["e"].values+pop_processes_df["rd"].values+pop_processes_df["ph"].values,
color=stackbar_colors[3], label=r"$\boldsymbol{p^+}$",zorder=2)
ax4.legend(bbox_to_anchor=(1, 1.05), fontsize=14,frameon=False)
ax4.set_title(r"\textbf{Population entering Fe \textsc{xiv}} $\boldsymbol{3s^2 3p\ \mathrm{^2 P_{3/2}}}$",fontsize=14,
y=1.02)
ax4.set_ylabel(r"\textbf{Contribution [\%]}",fontsize=14)
ax4.grid(axis="y",linestyle="dotted",zorder=-1)
ax4.text(0.005,0.02,r"\textbf{(d)}",fontsize=14,transform=ax4.transAxes,ha="left",va="bottom")
for ax_ in (ax1,ax2):
ax_.set_aspect("equal")
ax_.set_xlim(right=-500)
ax_.set_ylim(bottom=-1000,top=1000)
for ax_ in (ax1,ax2,ax3,ax4):
ax_.tick_params(labelsize=14)
ax_.tick_params(which="both",direction="in",top=True,right=True,left=True,bottom=True)
ax4.tick_params(top=False,bottom=False)
plt.savefig(fname="../../figs/ms/FeXIII_dens_FeXIV_pop.pdf",format="pdf",bbox_inches="tight",dpi=200)
