In [ ]:
colorbar_width = "5%"
def create_rec_eqs():
return patches.Rectangle((eis_eqs_xstart, eis_eqs_ystart),
eis_eqs_fovx, eis_eqs_fovy,linewidth=2,edgecolor="red",
facecolor="none",alpha=0.6)
def create_limb_circle(rsun):
return patches.Circle((0,0),rsun,linewidth=2,edgecolor="grey",
facecolor="none",alpha=0.6,ls="--")
def create_lunar_mask(rmoon):
return patches.Circle((0,0),rmoon,linewidth=0,edgecolor="none",
facecolor="white",alpha=1)
fig, ax = plt.subplots(figsize=(8,8),constrained_layout=True,gridspec_kw={'left':0.4,'top':0.8,'bottom':0.2})
green_where_disk_ext = np.where((x_2d_grid_green_arcsec_rot_ext + slit_center_x_green)**2 + \
(y_2d_grid_green_arcsec_rot_ext + slit_center_y_green)**2 < 940**2)
# green_line_int_masked = np.copy(green_fit_matrix_ext[1,:,:])
# green_line_int_masked[green_where_disk_ext] = np.nan
# norm_green_line_int = ImageNormalize(green_line_int_masked,interval=ManualInterval(0,850),
# stretch=SqrtStretch())
im = ax.pcolormesh(x_2d_grid_green_arcsec_rot_ext + slit_center_x_green,
y_2d_grid_green_arcsec_rot_ext + slit_center_y_green,
green_line_cont_ratio_masked,cmap=cmr.jungle_r,rasterized=True,vmin=0,vmax=0.16)
# ax.contour(img_xarcsec_array,img_yarcsec_array,FeXIV_line_cont_image_rot_scipy,levels=[0.4,0.65,0.9],alpha=0.6,
# colors=['grey'])
# plot_colorbar(im, ax,width=colorbar_width,fontsize=16)
ax.add_patch(create_lunar_mask(rsun_arcsec*np.float64(FeXI_line_cont_frame.header["MOONR"])/np.float64(FeXI_line_cont_frame.header["SUNR"])))
norm_aia_211 = ImageNormalize(aia_211_image,interval=AsymmetricPercentileInterval(0.5,99.97),stretch=LogStretch())
ax.pcolormesh(aia_211_xcoord, aia_211_ycoord, aia_211_image,cmap=aia_211_cm,norm=norm_aia_211,
rasterized=True,shading="auto")
ax.set_xlabel(r"\textbf{Solar-X [arcsec]}",fontsize=18,x=0.4,y=-0.3)
ax.set_ylabel(r"\textbf{Solar-Y [arcsec]}",fontsize=18,labelpad=-0.05)
xlim_zoomin = [-1500,-400]
ylim_zoomin = [-1000,900]
ax.set_aspect(1)
# ax.add_patch(create_rec_eqs())
ax.add_patch(create_limb_circle(rsun_arcsec))
ax.set_xlim(xlim_zoomin)
ax.set_ylim(ylim_zoomin)
ax.tick_params(labelsize=16)
feed_params = (x_2d_grid_green_arcsec_bin_rot_ext + slit_center_x_green,
y_2d_grid_green_arcsec_bin_rot_ext + slit_center_y_green,
green_fit_filename_index, green_path,
totality_green_df_ext,FeXIV_63rd_xslice,flatfield_1d_FeXIV_63rd,63,
green_fit_matrix_bin_ext, green_fit_matrix_bin_ext_err,
green_vlos_masked, green_vlos_masked_err, green_veff_masked, green_veff_masked_err)
spectral_line_zoom_in(-1020,42,fig,ax,[-1.3,0,0.95,0.45],(1,0.5),
*feed_params)
spectral_line_zoom_in(-1424,154,fig,ax,[-1.3,0.55,0.95,0.45],(1,0.5),
*feed_params)
spectral_line_zoom_in(-1103,-740,fig,ax,[-1.1,-0.60,0.95,0.45],(1,1),
*feed_params,label_axes=True,label_fontsize=18)
spectral_line_zoom_in(-976,-694,fig,ax,[0.1,-0.60,0.95,0.45],(1,0.5),
*feed_params)
spectral_line_zoom_in(-885,443,fig,ax,[-1.3,1.1,0.95,0.45],(1,0.1),
*feed_params)
spectral_line_zoom_in(-735,855,fig,ax,[0,1.1,0.95,0.45],(0.5,0),
*feed_params)
ax.text(-0.45,1.58,r"\textbf{Fe \textsc{xiv} 530.3 nm}",fontsize=30,transform=ax.transAxes)
fig.savefig(fname="../../figs/ms/FeXIV_zoomin.pdf",format="pdf",dpi=300,bbox_inches="tight")
# fig.savefig(fname="../../slides/reg_meet/figs/FeXIV_zoomin.png",format="png",dpi=300,bbox_inches="tight")
