diff --git a/PE_CBL.py b/PE_CBL.py
index 5e346ac8fe5789cd44e7feaa12a751ad9d4b5b91..941c268f332cb4e18d7882d5cbf9bb655b21e8e6 100644
--- a/PE_CBL.py
+++ b/PE_CBL.py
@@ -668,17 +668,17 @@ if __name__ == '__main__':
         fig.set_size_inches(6,6)
         #
         [ax1,ax0],c0,c1 = plot_CBL_identifiability(exp_A,da_settings_A['obs_error_sdev_assimilate'][0],ax=[ax1,ax0])
-        ax0.set_title(r'a) Exp. A, $\delta\overline{y}$')
+        ax0.set_title(r'a) Exp. A, $\delta\overline{y}$ (K)')
         ax0.set_xlabel('Time (h)')
         ax0.set_ylabel('Height (m)')
-        ax1.set_title(r'b) Exp. A, $\sigma_{p\prime\prime}/\sigma_{y^b}$')
+        ax1.set_title(r'b) Exp. A, $\sigma_{p\prime\prime}/\sigma_{y^b}$ (K$^{-1}$)')
         ax1.set_xlabel('Time (h)')
         ax1.set_ylabel('Height (m)')
         [ax3,ax2],c2,c3 = plot_CBL_identifiability(exp_D,da_settings_D['obs_error_sdev_assimilate'][0],ax=[ax3,ax2])
-        ax2.set_title(r'c) Exp. D, $\delta\overline{y}$')
+        ax2.set_title(r'c) Exp. D, $\delta\overline{y}$ (K)')
         ax2.set_xlabel('Time (h)')
         ax2.set_ylabel('Height (m)')
-        ax3.set_title(r'd) Exp. D, $\sigma_{p\prime\prime}/\sigma_{y^b}$')
+        ax3.set_title(r'd) Exp. D, $\sigma_{p\prime\prime}/\sigma_{y^b}$ (K$^{-1}$)')
         ax3.set_xlabel('Time (h)')
         ax3.set_ylabel('Height (m)')
         p.colorbar(c0,orientation='horizontal')