### 2 panels plot ###

import matplotlib.pyplot as plt
import matplotlib.image as mpimg

## Load your two saved plots
#img1 = mpimg.imread("/home/taneja/Documents/kanika/work/Paper_2026/supplementary_figures/hyytiala_lidar_updraft_lognormal_fit_JulAug2022.png")
#img2 = mpimg.imread("/home/taneja/Documents/kanika/work/Paper_2026/supplementary_figures/sigma_cb_lognormal_fit_IQR_UPDRAFTS_SGP.png")

img1 = mpimg.imread("/media/taneja/Expansion/Paper_2026/figures/hyytiälä/EGU_lnN100_vs_T_dailyMean_density_summerSquares_oneLegend_hyytiala_1.0_2012_2023_adiab.png")
img2 = mpimg.imread("/home/taneja/Documents/kanika/work/NEW_2025/SGP_new_local_filter/results_LWP_TM/1.0deg_adiab/TM_revised_lnN100_vs_T_dailyMean_density_summerSquares_oneLegend_SGP_1.0_2017_2023_adiab.png")


fig, axes = plt.subplots(1, 2, figsize=(12, 5))

for ax, img, label in zip(axes, [img1, img2], ["(a)", "(b)"]):
    ax.imshow(img)
    ax.axis("off")
    ax.text(
        0.15, 0.90, label,
        transform=ax.transAxes,
        fontsize=16,
##        fontweight="bold",
        va="top",
        ha="right",
        bbox=dict(facecolor="white", edgecolor="none", alpha=0.7, pad=2)
    )

plt.tight_layout()
#plt.savefig("/home/taneja/Documents/kanika/work/Paper_2026/figures/combined_hyttiälä_sgp_updrafts_lognormal_fit.pdf", bbox_inches="tight")
#plt.savefig("/home/taneja/Documents/kanika/work/Paper_2026/figures/combined_hyttiälä_sgp_updrafts_lognormal_fit.png", dpi=300, bbox_inches="tight")
plt.savefig("/media/taneja/Expansion/Paper_2026/figures/combined_plots/combined_hyttiälä_sgp_lnN100_vs_T_dailyMean_density_summerSquares_revised_with_multivariate_fit.png", dpi=300, bbox_inches="tight")
#plt.show()


##### 4 plots in one figure ####

import matplotlib.pyplot as plt
import matplotlib.image as mpimg

# ============================================================
# Load your four saved plots
# ============================================================

img1 = mpimg.imread(
    "/home/taneja/Documents/kanika/work/Paper_2026/supplementary_figures/"
    #"cpm_vs_modis_cdnc_n100_sat_0.5deg_localFilter_NOadiab_2012_2023_hyytiala_surfTemp.png"
    "cpm_vs_modis_cdnc_n100_sat_0.5deg_localFilter_NOadiab_2017_2023_SGP_LWP_ARMtemp_halfadiab.png"
)

img2 = mpimg.imread(
    "/home/taneja/Documents/kanika/work/Paper_2026/supplementary_figures/"
    #"cpm_vs_modis_cdnc_n100_sat_0.5deg_localFilter_adiab_2012_2023_hyytiala_surfTemp.png"
    "cpm_vs_modis_cdnc_n100_sat_0.5deg_localFilter_adiab_2017_2023_SGP_LWP_ARMtemp_halfadiab.png"
)

img3 = mpimg.imread(
    "/home/taneja/Documents/kanika/work/Paper_2026/supplementary_figures/"
    #"cpm_vs_modis_cdnc_n100_sat_2.0deg_localFilter_NOadiab_2012_2023_hyytiala_surfTemp.png"
    "cpm_vs_modis_cdnc_n100_sat_2.0deg_localFilter_NOadiab_2017_2023_SGP_LWP_ARMtemp_halfadiab.png"
)

img4 = mpimg.imread(
    "/home/taneja/Documents/kanika/work/Paper_2026/supplementary_figures/"
    #"cpm_vs_modis_cdnc_n100_sat_2.0deg_localFilter_adiab_2012_2023_hyytiala_surfTemp.png"
    "cpm_vs_modis_cdnc_n100_sat_2.0deg_localFilter_adiab_2017_2023_SGP_LWP_ARMtemp_halfadiab.png"
)

# ============================================================
# Arrange images and labels
# ============================================================

images = [img1, img2, img3, img4]
labels = ["(a)", "(b)", "(c)", "(d)"]

# ============================================================
# Create 2 x 2 panel figure
# ============================================================

fig, axes = plt.subplots(2, 2, figsize=(12, 10))

# Flatten axes for easy looping
axes = axes.flatten()

for ax, img, label in zip(axes, images, labels):
    ax.imshow(img)
    ax.axis("off")

    # Panel label at bottom right
    ax.text(
        0.75, 0.12, label,
        transform=ax.transAxes,
        fontsize=16,
        va="bottom",
        ha="right",
        bbox=dict(
            facecolor="white",
            edgecolor="none",
            alpha=0.7,
            pad=2
        )
    )

# ============================================================
# Adjust spacing
# ============================================================

plt.tight_layout()

# If the panels are still too far apart, use this instead:
# plt.subplots_adjust(wspace=0.02, hspace=0.02)

# ============================================================
# Save combined figure
# ============================================================

#plt.savefig(
#    "/home/taneja/Documents/kanika/work/Paper_2026/figures/"
#    "combined_cpm_vs_modis_hyytiala_sgp_4panels_other_domains_sgp.pdf",
#    bbox_inches="tight"
#)

#plt.savefig(
#    "/home/taneja/Documents/kanika/work/Paper_2026/figures/"
#    "combined_cpm_vs_modis_hyytiala_sgp_4panels_other_domains_sgp.png",
#    dpi=300,
#    bbox_inches="tight"
#)

#plt.savefig(
#    "/home/taneja/Documents/kanika/work/Paper_2026/figures/"
#    "combined_cpm_vs_modis_4panels_other_domains_hyytiala_TM.pdf",
#    bbox_inches="tight"
#)

#plt.savefig(
#   "/home/taneja/Documents/kanika/work/Paper_2026/figures/combined_cpm_vs_modis_4panels_other_domains_sgp_TM.pdf",
#   bbox_inches="tight"
#


#plt.show()