fig, (ax1,ax2,ax3) = plt.subplots(3,1,figsize=(6,12),constrained_layout=True,gridspec_kw={"height_ratios":[3,3,1]})
im1 = ax1.pcolormesh(img_xarcsec_array,img_yarcsec_array,FeXIV_line_cont_image_rot_scipy,vmin=0.1,vmax=1.3,
cmap=cmr.jungle_r,shading="auto",rasterized=True)
im2 = ax2.pcolormesh(img_xarcsec_array,img_yarcsec_array,FeXI_line_cont_image_rot_scipy,vmin=1,vmax=4,
cmap=cmcm.lajolla,shading="auto",rasterized=True)
ax1_xlim = ax1.get_xlim()
ax1_ylim = ax1.get_ylim()
ax2_xlim = ax2.get_xlim()
ax2_ylim = ax2.get_ylim()
starttime = datetime(2017,8,21,17,45,36)
for ii, green_row_ in totality_green_df_ext.iterrows():
date_obs = green_row_["date-obs"]
exptime = np.float64(green_row_["exptime"])
time_difference = date_obs - datetime(2017,8,21,17,45,36)
startindex = time_difference.total_seconds()
endindex = startindex + exptime
if startindex < 5:
slit_xshift_start = sun_x_fitpoly(5) - slit_pos
else:
slit_xshift_start = sun_x_fitpoly(startindex) - slit_pos
if endindex < 5:
slit_xshift_end = sun_x_fitpoly(5) - slit_pos
else:
slit_xshift_end = sun_x_fitpoly(endindex) - slit_pos
if startindex == endindex:
pass
else:
slit_start_x = - slit_xshift_start/pixel_ratio_to_arcsec*np.cos(np.deg2rad(np.abs(rotate_angle_context)))
slit_start_y = slit_xshift_start/pixel_ratio_to_arcsec*np.sin(np.deg2rad(np.abs(rotate_angle_context)))
slit_end_x = - slit_xshift_end/pixel_ratio_to_arcsec*np.cos(np.deg2rad(np.abs(rotate_angle_context)))
slit_end_y = slit_xshift_end/pixel_ratio_to_arcsec*np.sin(np.deg2rad(np.abs(rotate_angle_context)))
slit_center_x = (slit_start_x + slit_end_x)/2.
slit_center_y = (slit_start_y + slit_end_y)/2.
slit_plot_interval = (slit_end_y - slit_center_y)/(np.sin(np.deg2rad(np.abs(rotate_angle_context)))**2)
slit_plot_top_y = slit_center_y + slit_plot_interval
slit_plot_bottom_y = slit_center_y - slit_plot_interval
fill_between_x = np.linspace(1000,-2500)
fill_between_y_top = straight_line_slope(fill_between_x, slit_center_x, slit_plot_top_y,
1./np.tan(np.deg2rad(np.abs(rotate_angle_context))))
fill_between_y_bottom = straight_line_slope(fill_between_x, slit_center_x, slit_plot_bottom_y,
1./np.tan(np.deg2rad(np.abs(rotate_angle_context))))
ax1.fill_between(fill_between_x, fill_between_y_bottom, fill_between_y_top, alpha=0.4,color="grey",
linewidth=0,label=r"\textbf{3PAMIS}",zorder=3)
ax1.plot(fill_between_x, fill_between_y_bottom, alpha=1, color="#373C38",linewidth=0.06)
ax1.plot(fill_between_x, fill_between_y_top, alpha=1, color="#373C38",linewidth=0.06)
# plot_eis_fov(ax1,eis_atlas_30_fexiv_intmap_shifted,color="red",alpha=0.5,
# label=r"\textbf{EIS}")
# plot_eis_fov(ax1,eis_ch_fexiv_intmap_shifted,color="red",alpha=0.5,
# label=r"\textbf{EIS}")
# plot_comp_fov(ax1,zorder=1,alpha=0.4,color="#B481BB",label=r"\textbf{CoMP}")
# plot_comp_fov(ax2,zorder=1,color="#B481BB",alpha=0.4)
for ii, red_row_ in totality_red_df_ext.iterrows():
date_obs = red_row_["date-obs"]
exptime = np.float64(red_row_["exptime"])
time_difference = date_obs - datetime(2017,8,21,17,45,36+2)
startindex = time_difference.total_seconds()
endindex = startindex + exptime
if startindex < 5:
slit_xshift_start = sun_x_fitpoly(5) - slit_pos
else:
slit_xshift_start = sun_x_fitpoly(startindex) - slit_pos
if endindex < 5:
slit_xshift_end = sun_x_fitpoly(5) - slit_pos
else:
slit_xshift_end = sun_x_fitpoly(endindex) - slit_pos
if startindex == endindex:
pass
else:
slit_start_x = - slit_xshift_start/pixel_ratio_to_arcsec*np.cos(np.deg2rad(np.abs(rotate_angle_context)))
slit_start_y = slit_xshift_start/pixel_ratio_to_arcsec*np.sin(np.deg2rad(np.abs(rotate_angle_context)))
slit_end_x = - slit_xshift_end/pixel_ratio_to_arcsec*np.cos(np.deg2rad(np.abs(rotate_angle_context)))
slit_end_y = slit_xshift_end/pixel_ratio_to_arcsec*np.sin(np.deg2rad(np.abs(rotate_angle_context)))
slit_center_x = (slit_start_x + slit_end_x)/2.
slit_center_y = (slit_start_y + slit_end_y)/2.
slit_plot_interval = (slit_end_y - slit_center_y)/(np.sin(np.deg2rad(np.abs(rotate_angle_context)))**2)
slit_plot_top_y = slit_center_y + slit_plot_interval
slit_plot_bottom_y = slit_center_y - slit_plot_interval
fill_between_x = np.linspace(1000,-2500)
fill_between_y_top = straight_line_slope(fill_between_x, slit_center_x, slit_plot_top_y,
1./np.tan(np.deg2rad(np.abs(rotate_angle_context))))
fill_between_y_bottom = straight_line_slope(fill_between_x, slit_center_x, slit_plot_bottom_y,
1./np.tan(np.deg2rad(np.abs(rotate_angle_context))))
ax2.fill_between(fill_between_x, fill_between_y_bottom, fill_between_y_top, alpha=0.4,facecolor="grey",
linewidth=0,label=r"\textbf{3PAMIS}",zorder=3)
ax2.plot(fill_between_x, fill_between_y_bottom, alpha=1, color="#373C38",linewidth=0.06)
ax2.plot(fill_between_x, fill_between_y_top, alpha=1, color="#373C38",linewidth=0.06)
# plot_eis_fov(ax2,eis_atlas_30_fex_intmap_shifted,color="red",alpha=0.5,
# label=r"\textbf{EIS}")
# plot_eis_fov(ax2,eis_ch_fex_intmap_shifted,color="red",alpha=0.5,
# label=r"\textbf{EIS}")
xlim_zoomin = (-1600,1600)
ylim_zoomin = (-1300,1300)
handles, labels = ax1.get_legend_handles_labels()
temp_dict = {k:v for k,v in zip(labels, handles)}
ax1.legend(temp_dict.values(), temp_dict.keys(),fontsize=18,frameon=False,
bbox_to_anchor=(0,1.05,1,0.1),mode="expand",ncol=3)
ax3.step(green_frame_wavelength/63./10., green_frame_example.data[170,:],
where="mid",color="#787D7B")
ax1.text(0.02,0.97,r"\textbf{(a)}",fontsize=18,ha="left",va="top",
transform=ax1.transAxes)
ax2.text(0.02,0.97,r"\textbf{(b)}",fontsize=18,ha="left",va="top",
transform=ax2.transAxes)
ax3.text(0.02,0.93,r"\textbf{(c)}",fontsize=18,ha="left",va="top",
transform=ax3.transAxes)
for ax_ in (ax1,ax2):
ax_.set_aspect(1)
ax_.tick_params(labelsize=18,direction="in",top=True,right=True)
ax_.tick_params(which="major",length=6)
ax_.set_ylabel(r"\textbf{SOLAR-Y [arcsec]}",fontsize=18)
ax1.set_xlim(xlim_zoomin)
ax1.set_ylim(ylim_zoomin)
ax2.set_xlim(xlim_zoomin)
ax2.set_ylim(ylim_zoomin)
ax2.set_xlabel(r"\textbf{SOLAR-X [arcsec]}",fontsize=18)
ax3.set_ylim(bottom=2e3)
ax3_62 = ax3.secondary_xaxis('bottom', functions=(order63to62, order62to63))
ax3_62.spines['bottom'].set_position(('outward', 25))
ax3_61 = ax3.secondary_xaxis('bottom', functions=(order63to61, order61to63))
ax3_61.spines['bottom'].set_position(('outward', 50))
ax3.text(530.3*61/63,1.3e4,r"\textbf{61}",fontsize=18,color="#C1272D",va="center",ha="center")
ax3.text(530.3*62/63,2e4,r"\textbf{62}",fontsize=18,color="#008176",va="center",ha="center")
ax3.text(530.3,1.2e4,r"\textbf{63}",fontsize=18,va="center",ha="center")
ax3_61.set_xlabel(r"\textbf{Wavelength [nm]}",fontsize=18)
for ax_ in (ax3,ax3_62,ax3_61):
ax_.tick_params(labelsize=18,direction="in",which="major",length=6)
ax3_62.tick_params(colors="#008176")
ax3_62.spines['bottom'].set_color("#008176")
ax3_61.tick_params(colors="#C1272D")
ax3_61.spines['bottom'].set_color("#C1272D")
ax3.set_ylabel(r"\textbf{Counts}",fontsize=18)
plt.savefig("../../figs/ms/obs_summary_PaMIS.pdf",dpi=300,format="pdf",
bbox_inches="tight")
