small update fig 4

Author: pdnooteboom

figures/figure4.py

@@ -270,7 +270,8 @@ def check_lens(dic):
 
 cbar_ax.set_visible(False)
 #    cbar.ax.xaxis.set_label_position('right')
-cbar.ax.set_ylabel('$^{\circ}$C', fontsize=fs, rotation=-90)
+cbar.ax.set_ylabel('$^{\circ}$C', fontsize=fs, rotation='horizontal')
+cbar.ax.yaxis.set_label_coords(-0.1, -0.3)
 cbar.ax.tick_params(labelsize=fs)
 
 #legend

figures/merTf.py

@@ -132,7 +132,7 @@ def nearest_int(T, latm, lonm, lond, latd, res='HR'):
     Td = griddata(points, values, (lond, latd), method='nearest')
     return Td.flatten()
 
-def model_uncertainty(T, latm, lonm, lond, latd, deg = 2, res='HR'):
+def model_uncertainty(T, latm, lonm, lond, latd, deg = 3, res='HR'):
     # deg determines the degree box that is used to calculate the uncertainty
     # in the model due to the paleolocation
     # first nearest interpolate land values
@@ -157,9 +157,12 @@ def model_uncertainty(T, latm, lonm, lond, latd, deg = 2, res='HR'):
     for i in range(len(resl)):
         distances = (lond[i]-lonmf)**2+(latd[i]-latmf)**2
         arg = np.argmin(distances)
-        idx = np.where(np.logical_and(np.logical_and(np.logical_and(lonmf>=lonmf[arg]-deg,lonmf<=lonmf[arg]+deg),
-                                                     latmf>=latmf[arg]-deg),
-                latmf<=latmf[arg]+deg))
+ #       idx = np.where(np.logical_and(np.logical_and(np.logical_and(lonmf>=lonmf[arg]-deg,lonmf<=lonmf[arg]+deg),
+#                                                     latmf>=latmf[arg]-deg),
+#                latmf<=latmf[arg]+deg))
+        dist = np.sqrt((lonmf-lonmf[arg])**2+(latmf-latmf[arg])**2)
+        
+        idx = np.where(dist<=deg)
         resl[i] = np.min(T[idx])
         resh[i] = np.max(T[idx])